Thanks a lot, it seems to be correct, however...
load^-1 produces lots of infinite values.
summary(c((p$loadings)^-1))
Min. 1st Qu. Median Mean 3rd Qu. Max.
-3.094e+06 -3.747e+01 -6.409e+00 Inf 3.761e+01 Inf
Any idea why ?
Gottfried Helms wrote:
> Alexander Sirotkin schrieb:
>
>>hmm... isn't it just
>>
>>recoverfrompca <- function(pcaobject, pcascores) {
>> pcascores %*% t(pcaobject$loadings)
>>}
>>
>
> Hi -
>
> I don't know the syntax (and the functions) in S-Plus.
> But it is simply
> with the original data in zdat
> with the pca-loadings in load
> with the factor scores in fsc
>
> load*fsc = zdat // if zdat is unknown
> fsc = load^-1 * zdat // if fsc is unknown
>
> Or what else dir you look for?
>
> Gottfried Helms
>
> --
> Gottfried Helms
> Univ Kassel
.
.
=================================================================
Instructions for joining and leaving this list, remarks about the
problem of INAPPROPRIATE MESSAGES, and archives are available at:
. http://jse.stat.ncsu.edu/ .
=================================================================
