Displaying 1 result from an estimated 1 matches for "netoboot1".
Did you mean:
netoboot
2011 Nov 05
1
testing significance of axis loadings from multivariate dudi.mix
...reve <- sign(diag(r[k,]))
sol <- pc$loadings[ ,k]
sol <- sweep(sol, 2, reve, "*")
out <- out + ifelse(res > 0, sol <= 0, sol >= 0)
}
out/permutations
}
I tried to aply it to the case of dudi.mix instead of princomp this way:
netoboot1<-function (x, permutations=1000,...) {
dudinull <- dudi.mix(x, scannf = FALSE, nf = 3)
res <- dudinull$c1
out <- matrix(0, nrow=nrow(res), ncol=ncol(res))
N <- nrow(x)
for (i in 1:permutations) {
dudi <- dudi.mix(x[sample(N, replace=TRUE), ], scannf = FA...