[R] What an element in loading is missing? (princomp)

Andris Jankevics andza at osi.lv
Mon Oct 19 16:51:55 CEST 2009


Hi,

?loadings

print (pca_result$loadings,cutoff=0)
unclass(pca_result$loadings)

-- 
Andris Jankevics
AIO
Groningen Bioinformatics Centre
Groningen Biomolecular Sciences and Biotechnology Institute
University of Groningen
Kerklaan 30, Haren, 9751 NN, The Netherlands

On Mon, Oct 19, 2009 at 4:34 PM, Peng Yu <pengyu.ut at gmail.com> wrote:
> On Mon, Oct 19, 2009 at 1:58 AM, Peter Ehlers <ehlers at ucalgary.ca> wrote:
>>
>> Peng Yu wrote:
>>>
>>> Please see below that [3,1] of loadings is not printed. I am wondering
>>> what the problem is?
>>
>> Not trying ?loadings, perhaps??
>
> 'loadings' gives me the same thing.
>
>> pca_result$loadings
>
> Loadings:
>     Comp.1 Comp.2 Comp.3 Comp.4
> [1,]  0.556  0.378 -0.467 -0.574
> [2,] -0.511  0.266 -0.761  0.299
> [3,]        -0.887 -0.425 -0.169
> [4,] -0.651         0.152 -0.744
>
>               Comp.1 Comp.2 Comp.3 Comp.4
> SS loadings      1.00   1.00   1.00   1.00
> Proportion Var   0.25   0.25   0.25   0.25
> Cumulative Var   0.25   0.50   0.75   1.00
>> loadings(pca_result)
>
> Loadings:
>     Comp.1 Comp.2 Comp.3 Comp.4
> [1,]  0.556  0.378 -0.467 -0.574
> [2,] -0.511  0.266 -0.761  0.299
> [3,]        -0.887 -0.425 -0.169
> [4,] -0.651         0.152 -0.744
>
>
>
>
>
>>>> set.seed(0)
>>>> m=10
>>>> n=4
>>>> X=replicate(n,rnorm(m))
>>>> pca_result=princomp(X)
>>>> svd_result=svd(apply(X,2,function(x){x-mean(x)}))
>>>> pca_result$loadings
>>>
>>> Loadings:
>>>     Comp.1 Comp.2 Comp.3 Comp.4
>>> [1,]  0.867  0.306 -0.281 -0.275
>>> [2,] -0.215  0.532 -0.614  0.541
>>> [3,]        -0.682 -0.713 -0.155
>>> [4,] -0.448  0.397 -0.186 -0.780
>>>
>>>               Comp.1 Comp.2 Comp.3 Comp.4
>>> SS loadings      1.00   1.00   1.00   1.00
>>> Proportion Var   0.25   0.25   0.25   0.25
>>> Cumulative Var   0.25   0.50   0.75   1.00
>>>>
>>>> svd_result$v
>>>
>>>            [,1]       [,2]      [,3]       [,4]
>>> [1,] -0.86704537 -0.3060826 0.2810007 -0.2749261
>>> [2,]  0.21540928 -0.5321956 0.6144209  0.5411595
>>> [3,]  0.03936291  0.6823365 0.7134127 -0.1546279
>>> [4,]  0.44752846 -0.3968605 0.1859243 -0.7795205
>>>
>>> ______________________________________________
>>> R-help at r-project.org mailing list
>>> https://stat.ethz.ch/mailman/listinfo/r-help
>>> PLEASE do read the posting guide
>>> http://www.R-project.org/posting-guide.html
>>> and provide commented, minimal, self-contained, reproducible code.
>>>
>>>
>>
>
> ______________________________________________
> R-help at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.
>




More information about the R-help mailing list