我有一个像
df<-data.frame(date=c(rep("1/27/2010",times=30)),
loc1=c(rep(9:13,each=6)),
loc2=c(rep(c("N","E","W"),each=2)),
loc3=c(rep(c(1,2))),
tr1=c(rep(c(0,1),each=15)),
tr2=c(0,0,0,0,0,0,0,1,1,1,1,1,1,1,1,0,0,0,0,0,0,0,0,1,1,1,1,1,1,1),
tr3=c(1,1,1,1,1,1,1,2,2,2,2,2,2,2,2,3,3,3,3,3,3,3,4,4,4,4,4,4,4,4),
Birth=c(sample(c("early","late"),30,replace=TRUE,prob=c(0.5,0.5))),
Species=c(rep(c("A","B"),times=15)),
Status=c(sample(c(0,1),30,replace=TRUE,prob=c(0.7,0.3))))
df<-rbind(df,df)
我想为 loc3 的每个值创建单独的列,行由 loc1、loc2、tr1、tr2、tr3、Birth 和 Species 定义。我想“计算”所有共享这些值的观察的状态,并按 loc3 对计数进行分组。
我计划使用 reshape2 包中的 dcast 。
我写了一个函数来执行我想要的“计数”。我是 R 新手,虽然我确信有一个函数可以做到这一点,但我无法立即找到它,而且鉴于任务的简单性,尝试自己编写脚本似乎是一个值得练习的练习。
d.count<-function(x){
j=0
for (i in 1:length(x))
if (is.na(x{i])){
j<-j+0
}else if(x[i]==0){
j<-j+1
} else if(x[i]==1){
j<-j+0
}
return(j)
}
0 应该增加计数,而 1 和 NA 不应该。
所以
df_1<-dcast(df,date+loc1+loc2+tr1+tr2+tr3+Birth+Species~loc3,value.var="Status",fun.aggregate=d.count)
我得到错误
Error in if (is.na(x[i])) { : argument is of length zero
这让我觉得我不明白 dcast 是如何对待 fun.aggregate...
谢谢您的帮助!-JJE