文件名称：R语言 处理日期数据
文件大小：2KB
文件格式：R
更新时间：2020-09-19 12:57:25
R 日期数据
主要是用R语言，进行初期日期数据。 data1<-read.csv(file='C:/Users/li/Desktop/参数优化/a.csv',header = T) data2<-read.csv(file='C:/Users/li/Desktop/参数优化/b.csv',header = T) #提取日期数据 date1<-as.Date(data1$日期) dates<-paste(data2$年,data2$月,data2$日,sep='/') date2<-as.Date(dates,'%Y/%m/%d') xdate1<-date1[1:28]#tree 1 xdate2<-date1[29:54]#tree 2 xdate3<-date1[55:81]#tree 3 #构造09-01的日期形式 dates<-paste(data1$年-1,'09','01',sep='/') datestan<-as.Date(dates,'%Y/%m/%d') datestan1<-datestan[1:28]#tree 1 datestan2<-datestan[29:54]#tree 2 datestan3<-datestan[55:81]#tree 3 #提取气温数据 temp<-as.numeric(data2$平均气温) ##------------------------------------------------- #C t0 tb Tb Rc=function(c,Tb){ rc=0 ci<-datn1 while(rc!=c|ci<=datn1+140){ rc<-rc+sum(temp[ci]<=Tb) ci=ci+1 } print(ci) } # # Rf=function(f,Tb){ rf=0 fi<-tb while(rf<=f|fi<=tb+200){ if(temp[fi] <= Tb){ rf<-rf+0} else {rf<-rf+temp[fi]-Tb} fi=fi+1 } print(c(fi,rf)) } # #------------------------------ ypre<-datestan1 for(i in 1:28){ dat1<-datestan1[i] datn1<-which(date2==dat1) tb<-Rc(c=100,Tb=10) result<-Rf(f=100,Tb=10) y<-result[1] rf<-result[2] ypre[i]<-dat1+y-datn1 } r2<-mean(ypre-xdate1)