我的数据看起来像
pos year A B
1 2012 1 1
1 2013 1 NA
2 2012 NA NA
2 2013 NA 1
只有当 A 和 B 都不是 NA 时,我才尝试计算每个 pos 的 A 和 B 之间的百分比差异。
ddply(x, .(pos), summarize, diff = ifelse(is.na(A)==FALSE & is.na(B)==FALSE,
(rowsum(A, pos, na.rm=TRUE)-rowsum(B, pos, na.rm=TRUE))/rowsum(A, pos, na.rm=TRUE),""))
所以结果应该是
pos diff
1 0.5
2 NA
我不知道哪一部分是错的,但我的代码会产生
pos diff
1 0.5
1
2
2
任何想法,将不胜感激。谢谢!