Displaying 1 result from an estimated 1 matches for "contrp".
Did you mean:
contr
2012 Nov 29
0
Simper analysis with Morisita-Horn
...2 * top/((da + db) * aN * bN)
#############
}
}
average <- colMeans(contr)
if (nperm > 0) {
if (trace)
cat("Permuting", paste(comp[i, 1], comp[i, 2],
sep = "_"), "\n")
contrp <- matrix(ncol = P, nrow = n.a * n.b)
for (p in 1:nperm) {
groupp <- group[perm[p, ]]
ga <- comm[groupp == comp[i, 1], ]
gb <- comm[groupp == comp[i, 2], ]
for (j in 1:n.b) {
for (k in 1:n.a) {...