[R] Odp: "biplot" graphical options?

Petr PIKAL petr.pikal at precheza.cz
Wed Sep 2 14:21:46 CEST 2009


Hi
I had similar problem some time ago. I was advised to use eqscplot  but it 
did not suite my purpose so I used a little twist of biplot. (I added a 
result)


fit<-princomp(some.data, cor=T)

biplot(fit, xlabs=rep("", no.of.poins)) #biplot without points

#some scaling of plot parameters

rrr<-apply(fit$scores[,1:2],2, range)
(abs(rrr)+.1)*sign(rrr)
par(usr=as.vector(rrr))

# plot the points according to some leves
points(fit$scores[,1:2], 
col=jet.colors(11)[(as.numeric(kvalita.ag$typ[cc]))], pch=20)

So my advice is to plot without points and then add points from fit and 
treat them in usual way (pch, col, cex, ...)

Or you could modify source code for biplot and create your own, which 
could by quite tricky.

Regards
Petr




r-help-bounces at r-project.org napsal dne 02.09.2009 13:02:27:

> 
> 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.


More information about the R-help mailing list