[R] plotting prcomp
Gavin Simpson
gavin.simpson at ucl.ac.uk
Sat Jun 28 13:25:51 CEST 2008
On Fri, 2008-06-27 at 21:43 +0200, Dan Bolser wrote:
> 2008/6/27 Prof Brian Ripley <ripley at stats.ox.ac.uk>:
<snip />
> >
> > I forgot to answer that: just ask R to show you stats:::biplot.prcomp
>
> Ah! I never knew the ::: trick!
In this case you also need to know that biplot.prcomp is in stats.
getAnywhere() is a more general solution:
getAnywhere(biplot.prcomp)
The following R News article by Uwe Ligges covers aspects of finding the
sources for R functions:
Uwe Ligges. R help desk: Accessing the sources. R News, 6(4):43-45,
October 2006.
HTH
G
>
> I got it working after looking at the code -
>
>
> my.df <-
> data.frame(A=(x <- rnorm(100,mean=100, sd=10)),
> B=(y <- x + rnorm(100,mean=10, sd=10)))
>
> plot(my.df)
>
> my.pc <-
> prcomp(my.df, center=TRUE, scale=TRUE)
>
> biplot(my.pc)
>
> attributes(my.pc)
>
> head(scores <- my.pc$x)
>
> lam <- my.pc$sdev
>
> n <- nrow(my.pc$x)
>
> lam <- lam * sqrt(n)
>
> my.x <- (my.pc$x)[,1] / lam[1]
> my.y <- (my.pc$x)[,2] / lam[2]
>
> plot(my.x, my.y, type="n")
> text(my.x, my.y, labels=1:100)
>
>
>
> >>> What am I missing?
> >>
> >> Reading the references on ?biplot.prcomp, or a good book on the subject
> >> (such as MASS or Gower & Hand).
> >>
> >> A biplot is a specific type of plot and the help page does not spell out
> >> the rather mathematical details (and there are several variants).
> >>
> >>>
> >>> Dan.
> >>>
> >>> [[alternative HTML version deleted]]
> >>>
> >>> ______________________________________________
> >>> 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.
> >>>
> >>
> >> --
> >> Brian D. Ripley, ripley at stats.ox.ac.uk
> >> Professor of Applied Statistics, http://www.stats.ox.ac.uk/~ripley/
> >> University of Oxford, Tel: +44 1865 272861 (self)
> >> 1 South Parks Road, +44 1865 272866 (PA)
> >> Oxford OX1 3TG, UK Fax: +44 1865 272595
> >>
> >
> > --
> > Brian D. Ripley, ripley at stats.ox.ac.uk
> > Professor of Applied Statistics, http://www.stats.ox.ac.uk/~ripley/
> > University of Oxford, Tel: +44 1865 272861 (self)
> > 1 South Parks Road, +44 1865 272866 (PA)
> > Oxford OX1 3TG, UK Fax: +44 1865 272595
> >
>
>
>
--
%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%
Dr. Gavin Simpson [t] +44 (0)20 7679 0522
ECRC, UCL Geography, [f] +44 (0)20 7679 0565
Pearson Building, [e] gavin.simpsonATNOSPAMucl.ac.uk
Gower Street, London [w] http://www.ucl.ac.uk/~ucfagls/
UK. WC1E 6BT. [w] http://www.freshwaters.org.uk
%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%
More information about the R-help
mailing list