Displaying 1 result from an estimated 1 matches for "svycor".
Did you mean:
svrcore
2016 Apr 30
0
Unexpected scores from weighted PCA with svyprcomp()
...at the principal
component returned by svyprcomp() was very different from from the
reconstituted indicator as well as from the indicator returned by
predict(). "Different" here means weak correlation and different
distributions.
I hope the following code illustrates what I mean:
=====
svycor <- function(formula, design) {
# https://stat.ethz.ch/pipermail/r-help/2003-July/036645.html
stopifnot(require(survey))
covariance.matrix <- svyvar(formula, design)
variables <- diag(covariance.matrix)
correlation.matrix <- covariance.matrix / sqrt(variables %*%
t(variables)...