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
____________________________________________________________...