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.Пожалуйста, научите меня, как это исправить.
Заранее большое спасибо!