[英]r - calculate using next non-na value in data frame column
我在數據框中有一些數據,我想計算month
值之間的百分比變化。 問題是我在某些條目中不NA
,並且引發了計算。
irm code price pct.change
1 201807 511130F075A04 4.6600 2.192982
2 201806 511130F075A04 4.5600 1.333333
3 201805 511130F075A04 4.5000 -13.461538
4 201804 511130F075A04 5.2000 NA
5 201803 511130F075A04 NA NA
6 201802 511130F075A04 4.9100 1.867220
7 201801 511130F075A04 4.8200 -5.304519
8 201712 511130F075A04 5.0900 2.414487
9 201711 511130F075A04 4.9700 -3.307393
10 201710 511130F075A04 5.1400 NA
11 201709 511130F075A04 NA NA
12 201708 511130F075A04 5.2900 2.918288
13 201707 511130F075A04 5.1400 66.553255
14 201706 511130F075A04 3.0861 -10.664351
15 201705 511130F075A04 3.4545 -7.241824
問題出在pct.change
列的第4行和第10行。 它們是NA
但我希望使用不是NA
的最新price
值來計算它們。 所需的輸出將是(請參閱第4和10行):
irm code price pct.change
1 201807 511130F075A04 4.6600 2.192982
2 201806 511130F075A04 4.5600 1.333333
3 201805 511130F075A04 4.5000 -13.461538
**4 201804 511130F075A04 5.2000 5.906314**
5 201803 511130F075A04 NA NA
6 201802 511130F075A04 4.9100 1.867220
7 201801 511130F075A04 4.8200 -5.304519
8 201712 511130F075A04 5.0900 2.414487
9 201711 511130F075A04 4.9700 -3.307393
**10 201710 511130F075A04 5.1400 -2.835539**
11 201709 511130F075A04 NA NA
12 201708 511130F075A04 5.2900 2.918288
13 201707 511130F075A04 5.1400 66.553255
14 201706 511130F075A04 3.0861 -10.664351
15 201705 511130F075A04 3.4545 -7.241824
我已經嘗試過標准(x/lead(x) - 1)*100
和使用(x/lag(which(!is.na(lead(x))
幾種變體,但我似乎缺少了一些東西。 base
或什dplyr
是dplyr
? 我不想替換NA,我想保留它們。
@LAP的評論可能是最好的方法。 使用data.table
的語法要好data.table
library(data.table)
setDT(df)
df[!is.na(price), pct.change := 100*(price/shift(price, type = 'lead') - 1)]
# irm code price pct.change
# 1: 201807 511130F075A04 4.6600 2.192982
# 2: 201806 511130F075A04 4.5600 1.333333
# 3: 201805 511130F075A04 4.5000 -13.461538
# 4: 201804 511130F075A04 5.2000 5.906314
# 5: 201803 511130F075A04 NA NA
# 6: 201802 511130F075A04 4.9100 1.867220
# 7: 201801 511130F075A04 4.8200 -5.304519
# 8: 201712 511130F075A04 5.0900 2.414487
# 9: 201711 511130F075A04 4.9700 -3.307393
# 10: 201710 511130F075A04 5.1400 -2.835539
# 11: 201709 511130F075A04 NA NA
# 12: 201708 511130F075A04 5.2900 2.918288
# 13: 201707 511130F075A04 5.1400 66.553255
# 14: 201706 511130F075A04 3.0861 -10.664351
# 15: 201705 511130F075A04 3.4545 NA
在Base R中,您可以決定替換:
a = which(is.na(df$price))-1
transform(df,pct.change=replace(pct.change,a,100*(price[a]/price[a+2]-1)))
irm code price pct.change
1 201807 511130F075A04 4.6600 2.192982
2 201806 511130F075A04 4.5600 1.333333
3 201805 511130F075A04 4.5000 -13.461538
4 201804 511130F075A04 5.2000 5.906314
5 201803 511130F075A04 NA NA
6 201802 511130F075A04 4.9100 1.867220
7 201801 511130F075A04 4.8200 -5.304519
8 201712 511130F075A04 5.0900 2.414487
9 201711 511130F075A04 4.9700 -3.307393
10 201710 511130F075A04 5.1400 -2.835539
11 201709 511130F075A04 NA NA
12 201708 511130F075A04 5.2900 2.918288
13 201707 511130F075A04 5.1400 66.553255
14 201706 511130F075A04 3.0861 -10.664351
15 201705 511130F075A04 3.4545 -7.241824
聲明:本站的技術帖子網頁,遵循CC BY-SA 4.0協議,如果您需要轉載,請注明本站網址或者原文地址。任何問題請咨詢:yoyou2525@163.com.