Thanks Andris, Michael and Petr for your prompt and kind feedbacks.
I will try generating my own biplot from low-level graph commands... I hope it
will work.
Best regards,
Marco
--
Marco Manca, MD
University of Maastricht
Faculty of Health, Medicine and Life Sciences (FHML)
Cardiovascular Research Institute (CARIM)
PO Box 616
6200 MD Maastricht
E-mail: m.manca at path.unimaas.nl
Office telephone: +31(0)433874633
Personal mobile: +31(0)626441205
Twitter: @markomanka
*********************************************************************************************************************
This email and any files transmitted with it are confidential and solely for the
use of the intended recipient.
It may contain material protected by privacy or attorney-client privilege. If
you are not the intended recipient or the person responsible for
delivering to the intended recipient, be advised that you have received this
email in error and that any use is STRICTLY PROHIBITED.
If you have received this email in error please notify us by telephone on
+31626441205 Dr Marco MANCA
*********************************************************************************************************************
________________________________________
Da: andris.jankevics at gmail.com [andris.jankevics at gmail.com] per conto di
Andris Jankevics [andza at osi.lv]
Inviato: mercoled? 2 settembre 2009 14.31
A: Manca Marco (PATH)
Cc: r-help at r-project.org
Oggetto: Re: [R] "biplot" graphical options?
Hi, You can make a biplot on Your own, it is not so hard. And in this
case You can change parameters for every low level function as You
wish.
PC <- prcomp (iris[,1:4])
lambda <- PC$sdev * sqrt(nrow(PC$x))
plot (t(t(PC$x)/lambda),pch=16,col=as.numeric(iris[,5]))
par (new=T)
Rot <- t(t(PC$rotation)*lambda)
XLIM <- c(-max(abs(Rot[,1])),max(abs(Rot[,1])))
XLIM <- XLIM+(XLIM*0.1)
plot(Rot,col=4,axes=FALSE,xlim=XLIM,ylim=XLIM,pch="")
arrows (rep(0,nrow(PC$rotation)),rep(0,nrow(PC$rotation)),Rot[,1],Rot[,2],col=4)
text (Rot[,1:2],rownames(Rot),col=6)
axis (3)
axis (4)
Best regards,
Andris
On Wed, Sep 2, 2009 at 1:02 PM, Manca Marco
(PATH)<m.manca at path.unimaas.nl> wrote:>
> Dear R-help fellows
>
> good afternoon.
>
> I am struggling in the attempt to impose some graphical conditions
(changing point symbols, colors, etc) to biplot function (I am using it to
visualize the results of princomp) but I can't apparently manage to change
anything but the axis... and I have been browsing manuals and vignettes without
finding any explicit suggestions on how to operate...
>
> Can anyone, please, point my attention to the relevant documentation?
>
> Thank you in advance and best regards,
> Marco
>
> --
> Marco Manca, MD
> University of Maastricht
> Faculty of Health, Medicine and Life Sciences (FHML)
> Cardiovascular Research Institute (CARIM)
> PO Box 616
> 6200 MD Maastricht
>
> E-mail: m.manca at path.unimaas.nl
> Office telephone: +31(0)433874633
> Personal mobile: +31(0)626441205
> Twitter: @markomanka
>
>
>
*********************************************************************************************************************
>
> This email and any files transmitted with it are confide...{{dropped:15}}
>
> ______________________________________________
> 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.
>
