[BioC] screeplot in package pcaMethods fails

Yannick Wurm yannick.wurm at unil.ch
Fri Dec 12 09:37:24 CET 2008


Yes, I guess you're expected to use the following instead
  pcaMethods::plotR2

On Dec 12, 2008, at 0:41 , Wolfgang Huber wrote:

>
> Hi Yury,
>
> thank you. As far as I can tell, no "screeplot" method is defined  
> by the pcaMethods package, so you would be falling back on  
> stats:::screeplot.default, whose implementation indeed asks for x 
> $sdev and probably is what throws the error that you report  
> (remember next time to add the output of "traceback()" to your  
> error message).
>
> Perhaps the pcaMethods package documentation could be clarified  
> regarding the existence and/or usage of a "screeplot" method.
> (Typing "class ? pcaRes" currently tells you that such a method  
> exists.)
>
> Best wishes
>      Wolfgang
>
> ----------------------------------------------------
> Wolfgang Huber, EMBL-EBI, http://www.ebi.ac.uk/huber
>
>
> Yury Bukhman ha scritto:
>> Hi,  I have created a pcaRes object using pcaMethods.  slplot  
>> works OK, but screeplot fails with this error:
>> R> screeplot(kubo1.pca)
>> Error in x$sdev : $ operator not defined for this S4 class
>> More info below.
>> R> class(kubo1.pca)
>> [1] "pcaRes"
>> attr(,"package")
>> [1] "pcaMethods"
>> R> print(kubo1.pca)
>> svd calculated PCA
>> Importance of component(s):
>>                  PC1    PC2    PC3     PC4
>> R2            0.4585 0.2713 0.1381 0.04216
>> Cumulative R2 0.4585 0.7299 0.8680 0.91014
>> 11405 	Variables
>> 12 	Samples
>> 0 	NA's
>> 4 	Calculated component(s)
>> Data was mean centered before running PCA Scores structure:
>> [1] 12  4
>> Loadings structure:
>> [1] 11405     4
>> R> sessionInfo()
>> R version 2.8.0 (2008-10-20) x86_64-redhat-linux-gnu locale:
>> LC_CTYPE=en_US.UTF-8;LC_NUMERIC=C;LC_TIME=en_US.UTF-8;
>> LC_COLLATE=en_US.UTF-8;LC_MONETARY=en_US.UTF-8;
>> LC_MESSAGES=en_US.UTF-8;LC_PAPER=en_US.UTF-8;
>> LC_NAME=en_US.UTF-8;LC_ADDRESS=en_US.UTF-8;
>> LC_TELEPHONE=en_US.UTF-8;LC_MEASUREMENT=en_US.UTF-8;
>> LC_IDENTIFICATION=en_US.UTF-8
>> attached base packages:
>> [1] splines   tools     stats     graphics  grDevices utils      
>> datasets [8] methods   base     other attached packages:
>>  [1] simpleaffy_2.16.1    genefilter_1.20.1     
>> survival_2.34-1      [4] pcaMethods_1.18.0    pls_2.1-0             
>> MASS_7.2-44          [7] affyPLM_1.16.0       gcrma_2.12.1          
>> matchprobes_1.14.0  [10] affy_1.18.2          preprocessCore_1.2.1  
>> affyio_1.8.1        [13] Biobase_2.0.1        rkward_0.4.9         
>> loaded via a namespace (and not attached):
>> [1] annotate_1.18.0     AnnotationDbi_1.4.2 DBI_0.2-4          [4]  
>> grid_2.8.0          lattice_0.17-15     latticeExtra_0.5-4 [7]  
>> RSQLite_0.7-0      Thanks.
>> Yury
>
> _______________________________________________
> Bioconductor mailing list
> Bioconductor at stat.math.ethz.ch
> https://stat.ethz.ch/mailman/listinfo/bioconductor
> Search the archives: http://news.gmane.org/ 
> gmane.science.biology.informatics.conductor



More information about the Bioconductor mailing list