Displaying 1 result from an estimated 1 matches for "groupp".
Did you mean:
groups
2012 Nov 29
0
Simper analysis with Morisita-Horn
...s(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) {
########### Morisita-Horn
aNp <- sum(gpa[k, ])
bNp <- sum(gpb[j, ])
dap &...