迭代R中每列的圖

Question

我想為此數據集中的每一列（最后一列除外）創建一個qqplot：

data("pottery", package = "HSAUR2")

這是我到目前為止的內容，但僅列出了最后一列：

par(mfrow=c(3,3))
for (i in ncol(pottery[, 1:9])){
    qqnorm(pottery[, i], main = names(pottery[i]))
    qqline(x)
}

我不認為這應該很難，但是我不確定我缺少什么。

一如既往的感謝。

Answer 1

par(mfrow=c(3,3))
for (i in 1:ncol(pottery[,1: ncol(pottery) - 1 ])){  
qqnorm(pottery[, i], main = names(pottery[i]))
   qqline(pottery[, i])
}

我在您的腳本中更改了兩件事。 i in 1：ncol（pottery）表示您希望從開始到有多少列都跨越該范圍。 當您在ncol（pottery）中說“我”時，就是在最后一列中說“ i”。 然后我寫了pottery [，ncol（pottery）-1]，說您想忽略最后一列。

Answer 2

使用ggplot2和tidyr將數據tidyr轉換為長：

library(tidyverse)
gather(pottery, condition, measurement, Al2O3:BaO, factor_key = TRUE) %>%
  ggplot(aes(sample = measurement)) +
  facet_wrap(~ condition, scales = "free") +
  stat_qq() +
  stat_qq_line()