【文件属性】:
文件名称: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)