我有一些带有动物园对象的文件,看起来像这样(每个文件在不同的日期开始和结束):
code pp
1942-06-01 4016 0
1942-06-02 4016 NA
1942-06-03 4016 0
1942-06-04 4016 0
1942-06-05 NA 0
1942-06-06 NA 0
我想在 9 月、10 月和 11 月的所有文件的 pp 之间做一个相关矩阵(显示代码,以便我可以识别谁是谁)。由于代码列中的 NA,我无法使用 list.files 函数(由 Joran 在不同文件之间的相关矩阵中提供的代码)。所以我想出了以下代码:
files <- list.files(pattern=".csv")
xx<-read.zoo(files[1],sep=",", header=TRUE,index.column=1)
name<- as.name(xx$code[[1]])
colnames(xx) <- c("code", name)
x<-xx[months(time(xx), TRUE) %in% c("Sep", "Oct", "Nov")]
yy<-read.zoo(files[2],sep=",", header=TRUE,index.column=1)
name<- as.name(yy$code[1])
colnames(yy) <- c("code", name)
y<-yy[months(time(yy), TRUE) %in% c("Sep", "Oct", "Nov")]
CET<-merge(x, y, all = TRUE, fill = NA, check.names=FALSE)
for (i in 3:length(files))
{
z<-read.zoo(files[i],sep=",", header=TRUE,index.column=1)
name<- as.name(z$code[1])
colnames(z) <- c("code", name)
CET<-merge(CET, z, all = TRUE, fill = NA, check.names=FALSE)
}
a<-1:(dim(CET)[2])
even <- a[ a%%2 == 0 ]
# saves the precipitation column (even numbers) and discards the code ones
dat<-CET[,even]
c.mat<-cor(dat,use="pairwise.complete.obs" )
但是有些问题:在相关矩阵中,一些列/行名称有一个额外的“.z”或“.CET”,最重要的是相关系数不正确!我找不到问题所以任何帮助找到问题或提出更简单的代码来做到这一点将不胜感激!