R Dataframe новый столбец с скользящей корреляцией двух других столбцов - PullRequest
Новая
       10

R Dataframe новый столбец с скользящей корреляцией двух других столбцов

1 голос
/ 30 апреля 2019

Я застрял около часа с чем-то простым. У меня есть фрейм данных, и мне нужно будет добавить новый столбец, содержащий соотношение двух других столбцов. Важно отметить, что мне нужно, чтобы вывод был новым столбцом данных, содержащим все начальные строки, в том числе без корреляции. Вот мой код: 

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 строки показывают NA в столбце «out», но я ожидаю, что значение NA будет в первых семи строках столбца «out». 

> 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')

это устранило проблему наличия NA в первых 7 значениях, но все равно дало кучу строк NA в хвостовой части кадра данных. Любое предложение, пожалуйста? Спасибо

...