Warning: file_get_contents(/data/phpspider/zhask/data//catemap/4/r/65.json): failed to open stream: No such file or directory in /data/phpspider/zhask/libs/function.php on line 167

Warning: Invalid argument supplied for foreach() in /data/phpspider/zhask/libs/tag.function.php on line 1116

Notice: Undefined index: in /data/phpspider/zhask/libs/function.php on line 180

Warning: array_chunk() expects parameter 1 to be array, null given in /data/phpspider/zhask/libs/function.php on line 181
如何将循环中的数据子集,然后对r中生成的数据帧数据执行循环_R_Loops_For Loop_Subset - Fatal编程技术网

如何将循环中的数据子集,然后对r中生成的数据帧数据执行循环

如何将循环中的数据子集,然后对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

您介意帮我运行从循环中的子集生成的数据帧吗。我想对每个子集进行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  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