search for: pcasig

Displaying 2 results from an estimated 2 matches for "pcasig".

Did you mean: pasig
2007 Nov 26
0
writing statistical function for boot
...ne on m with a statistic function calculating the eigenvector bootstrap error ratio. If R = number of bootstrap replicates, then my function should look something like this where m.i is the ith bootstrap sample matrix (of course this is wrong, but I wrote it to give an idea of what I am after) pcasig <- function(m) { for (i in 1:R) { pca.i <- prcomp(m.i) eigen.i <- pca.i$rotation estim <- sum(eigen.1+….+eigen.R) / R a <- matrix(estim, p, n) se <- ((sum((eigen.i-estim)^2))/(R-1))^0.5 t.i <- a / se } } Thanks for your consideration, Monica ________________________________...
2007 Nov 27
0
Function to calculate eigenvector bootstrap error
...on m with a statistic function calculating the eigenvector bootstrap error ratio. If R = number of bootstrap replicates, then my function should look something like this where m.i is the ith bootstrap sample matrix (of course this is wrong, but I wrote it to give an idea what I actually want) pcasig <- function(m) { for (i in 1:R) { pca.i <- prcomp(m.i) eigen.i <- pca.i$rotation estim <- sum(eigen.1+….+eigen.R)/R se <- ((sum((eigen.i-estim)^2))/(R-1))^0.5 t.i <- estim/se } } Thanks for your consideration, Monica ____________________________________________________________...