Question

我希望在R中绘制残差的散点图矩阵，以便进行纵向数据分析。

我写了给定的代码：

head(LungCapData2)

数据如下：

     Age LungCap Height Gender Smoke
 1   9   3.124   57.0   female no
 2   8   3.172   67.5 female   no 
 3   7   3.160   54.5 female   no 
 4   9   2.674   53.0   male   no 
 5   9   3.685   57.0   male   no 
 6   8   5.008   61.0 female   no

我开始使用OLS获得残差：

lm1 <- lm(LungCap~Height, data=LungCapData)
lm1res<- resid(lm1)
roundht <- round(LungCapData2$Height)
Lung<-reshape(LungCapData2[c("Age","lm1res","roundht")],
          direction="wide",v.names="lm1res", timevar="roundht",idvar="Age")

问题在于将其重塑为矩阵的最后一段代码。我希望得到

的散点图

pairs(Lung)

但是我的肺数据如下：

head(Lung)
  Age
1   9
2   8
3   7
4   6
5   5
6   4
  lm1res.c(57, 68, 54, 53, 61, 58, 56, 60, 50, 59, 62, 49, 52, 48, 65, 66, 55, 51, 69, 63, 64, 47, 46, 70, 72, 71, 67, 74, 73)
1

不能弄清楚我哪里出了问题。

Answer 1

需要为lmframe和roundht分配给数据框列，然后再以重塑形式引用数据框：

    LungCapData$lm1res<- resid(lm1)
    LungCapData$roundht <- round(LungCapData$Height)
    Lung<-reshape(LungCapData[c("Age","lm1res","roundht")],
                  direction="wide", v.names="lm1res",timevar="roundht",idvar="Age")
    pairs(Lung)

获取残差散点图矩阵

1 个答案: