我被简单的事情困住了大约一个小时。 我有一个数据框,我需要添加一个新列,其中包含其他两个列的相关性。 重要的是要注意,我需要输出为包含所有初始行(包括不相关的行)的新数据框列。 这是我的代码:
data <- mtcars[1:50,]
data$out <- rollapply(data.frame(data$mpg, data$wt), 8 ,function(x) cor(x[,1],x[,2]), by.column=FALSE, fill=NA)
如果检查输出,出于某种奇怪的原因,尾巴上有很多NA行:
> tail(data)
mpg cyl disp hp drat wt qsec vs am gear carb new_col out
NA.12 NA NA NA NA NA NA NA NA NA NA NA <NA> NA
NA.13 NA NA NA NA NA NA NA NA NA NA NA <NA> NA
NA.14 NA NA NA NA NA NA NA NA NA NA NA <NA> NA
NA.15 NA NA NA NA NA NA NA NA NA NA NA <NA> NA
NA.16 NA NA NA NA NA NA NA NA NA NA NA <NA> NA
NA.17 NA NA NA NA NA NA NA NA NA NA NA <NA> NA
为什么会这样? 另外,如果我检查头部,则仅前3行在“ out”列中显示NA,但我希望在“ out”列的前7行中具有NA值。
> head(data,10)
mpg cyl disp hp drat wt qsec vs am gear carb new_col out
Mazda RX4 21.0 6 160.0 110 3.90 2.620 16.46 0 1 4 4 6 -- 2.62 NA
Mazda RX4 Wag 21.0 6 160.0 110 3.90 2.875 17.02 0 1 4 4 6 -- 2.875 NA
Datsun 710 22.8 4 108.0 93 3.85 2.320 18.61 1 1 4 1 4 -- 2.32 NA
Hornet 4 Drive 21.4 6 258.0 110 3.08 3.215 19.44 1 0 3 1 6 -- 3.215 -0.6294594
Hornet Sportabout 18.7 8 360.0 175 3.15 3.440 17.02 0 0 3 2 8 -- 3.44 -0.6231828
Valiant 18.1 6 225.0 105 2.76 3.460 20.22 1 0 3 1 6 -- 3.46 -0.6374638
Duster 360 14.3 8 360.0 245 3.21 3.570 15.84 0 0 3 4 8 -- 3.57 -0.9454915
Merc 240D 24.4 4 146.7 62 3.69 3.190 20.00 1 0 4 2 4 -- 3.19 -0.7318425
Merc 230 22.8 4 140.8 95 3.92 3.150 22.90 1 0 4 2 4 -- 3.15 -0.7428921
Merc 280 19.2 6 167.6 123 3.92 3.440 18.30 1 0 4 4 6 -- 3.44 -0.7784487
接下来,我意识到如果我使用
data$out <- rollapply(data.frame(data$mpg, data$wt), 8 ,function(x) cor(x[,1],x[,2]), by.column=FALSE, fill=NA, align='right')
这解决了在前7个值中包含NA的问题,但仍然在数据帧的尾部给出了一堆NA行。 有什么建议吗?谢谢