如何将循环中的数据子集,然后对r中生成的数据帧数据执行循环
您介意帮我运行从循环中的子集生成的数据帧吗。我想对每个子集进行Mann-Kendall测试(库(rkt)),数据如下:如何将循环中的数据子集,然后对r中生成的数据帧数据执行循环,r,loops,for-loop,subset,R,Loops,For Loop,Subset,您介意帮我运行从循环中的子集生成的数据帧吗。我想对每个子集进行Mann-Kendall测试(库(rkt)),数据如下: date Year Month Day sum-rainfall frac overall-slope Sens Slope Seq 1 1911-10-01 1911 10 1 216.3 1911.11 278.5459 a)1911-2019: -0.85 mm/yr 1 2 1912-10-01
date Year Month Day sum-rainfall frac overall-slope Sens Slope Seq
1 1911-10-01 1911 10 1 216.3 1911.11 278.5459 a)1911-2019: -0.85 mm/yr 1
2 1912-10-01 1912 10 1 275.3 1912.11 277.7006 a)1911-2019: -0.85 mm/yr 1
3 1913-10-01 1913 10 1 266.2 1913.11 276.8553 a)1911-2019: -0.85 mm/yr 1
.
.
78 1988-10-01 1988 10 1 225.6 1988.11 213.4582 a)1911-2019: -0.85 mm/yr 4
79 1989-10-01 1989 10 1 234.6 1989.11 212.6129 a)1911-2019: -0.85 mm/yr 4
.
.
108 2018-10-01 2018 10 1 206.8 2018.11 188.0994 a)1911-2019: -0.85 mm/yr 6
109 2019-10-01 2019 10 1 168.8 2019.11 187.2541 a)1911-2019: -0.85 mm/yr 6
for (i in rainfall_wet$Seq){
m<- assign(paste("segment_wet", i, sep = "_"),subset(rainfall_wet, Seq==i))
FD<-rkt(m$Year,m$sum_rainfall,m$Month,,TRUE)
print(FD)
FD$B
b0<- median(m$sum_rainfall)-(FD$B* median(m$Year))
m$overall_slope<- b0+(FD$B*m$Year)
length(m$Year)
a3<-as.character(m$Year[1])
b3<-as.character(m$Year[length(m$Year)])
z3<-as.character(lapply(FD$B, round,2))
m$`Sens Slpe`<- paste("b)",a3,"-",b3,":"," ",z3," ","mm/yr",sep='')
}
日期年月日总降雨量分形总坡度传感器坡度序列
1911-10-011911101216.31911.11278.5459 a)1911-2019:-0.85毫米/年
1912-10-01 1912 10 1 275.3 1912.11 277.7006 a)1911-2019:-0.85毫米/年1
3 1913-10-01 1913 10 1 266.2 1913.11 276.8553 a)1911-2019:-0.85毫米/年1
.
.
78 1988-10-01 1988 10 1 225.6 1988.11 213.4582 a)1911-2019:-0.85毫米/年4
1989年10月01日1989年10月234.6日1989年11月212.6129 a)1911年至2019年:-0.85毫米/年4
.
.
108 2018-10-01 2018 10 1 206.8 2018.11 188.0994 a)1911-2019:-0.85毫米/年6
109 2019-10-01 2019 10 1 168.8 2019.11 187.2541 a)1911-2019:-0.85毫米/年6
(i)在雨量中(湿$Seq){
像这样的
plyr::ddply(rainfall_wet, plyr::.(Seq), function(m){
FD<-rkt(m$Year,m$sum_rainfall,m$Month,,TRUE)
print(FD)
FD$B
b0<- median(m$sum_rainfall)-(FD$B* median(m$Year))
m$overall_slope<- b0+(FD$B*m$Year)
length(m$Year)
a3<-as.character(m$Year[1])
b3<-as.character(m$Year[length(m$Year)])
z3<-as.character(lapply(FD$B, round,2))
m$`Sens Slpe`<- paste("b)",a3,"-",b3,":"," ",z3," ","mm/yr",sep='')
return(m)
}
)
plyr::ddply(降雨湿,plyr::(Seq),函数(m){
FD