menhinick <- function(x) {
sum(x>0)/sqrt(sum(x))
}
ddply(df,~Newcode,function(x) {
data.frame(RICHNESS=menhinick(x[-1]))
})
Newcode RICHNESS
1 BL01 1.4508215
2 BL02 1.4050583
3 BL03 1.5321132
4 BL04 1.4142136
5 BL05 0.6479552
6 BL06 1.2880334
7 BL07 0.2820650
8 BL08 1.0684120
9 BL09 0.3252526
10 BL10 0.7363117
11 BL11 1.8349137
12 BL12 1.2174252
13 BL13 1.0563106
14 BL14 1.6669913
15 BL15 1.6175783
16 BL16 1.3056432
17 BL17 1.6294405
18 BL18 0.9885727
向量Newcode
已经在这里。但是我找不到RICHNESS
。
请教我如何解决。
非常感谢您!