进行在线课程,我接到了以下任务:
load("skew.RData")
使用QQ图,比较矩阵每列的分布与法线。也就是说,在每列上使用qqnorm()。要快速完成此操作,您可以使用以下代码行为3x3 = 9绘图设置网格。 (" mfrow"意味着我们希望逐行填充多图形网格。另一种选择是mfcol。)
par(mfrow = c(3,3))
然后你可以使用for循环,遍历列,并一次显示一个qqnorm()图。您应该用您自己的代码替换**之间的文本。
for (i in 1:9) {
**put your qqnorm call here**
}
我有一个文件calld dat,有9列和大约1000行。
可以给我一个暗示QQ情节命令的样子吗?
答案 0 :(得分:2)
比较矩阵每列的分布与正常
在您的数据文件中读取后,假设您的矩阵为dat
load("skew.RData") ## read in matrix `dat`
dat <- scale(dat) ## standardization
par(mfrow = c(3, 3))
for (i in 1:9) {
qqnorm(dat[, i], main = paste0(i, "-th column"))
qqline(dat[, i])
}