我psych::principal
在另一个函数中使用,各种rotate
函数传递给principal
. (principal
提供许多旋转选项并将它们传递给不同的其他功能)。
我需要获得使用任何旋转过程的旋转矩阵,并实现。
所有下游轮换程序都提供此功能,但似乎没有return()
被principal
.
例如:
randomcor <- cor(matrix(data = rnorm(n = 100), nrow = 10))
library(psych)
principalres <- principal(r = randomcor, nfactors = 3, rotate = "none")
unrot.loa <- unclass(principalres$loadings)
principalrot <- principal(r = randomcor, nfactors = 3, rotate = "varimax") # there is no way to retrieve the rot.mat from principal
# but this CAN be done from the underlying varimax!
varimaxres <- varimax(x = unrot.loa)
varimaxres$rotmat # see, THIS is what I want!
我不愿意从principal
. (不要重复你自己或其他人,就像说的那样)。
有谁知道如何:
- 我可以优雅地,不知何故,神奇地,从检索 ,虽然它似乎没有返回它?
rotmat
principal()
- 或者,我可以估算任何
rotmat
必须“发生”的情况,因为我知道旋转和未旋转的载荷?