<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 3.2//EN">
<HTML>
<HEAD>
<META HTTP-EQUIV="Content-Type" CONTENT="text/html; charset=iso-8859-1">
<META NAME="Generator" CONTENT="MS Exchange Server version 6.5.7653.38">
<TITLE>RE: [adegenet-forum] question colorplot</TITLE>
</HEAD>
<BODY>
<!-- Converted from text/plain format -->
<P><FONT SIZE=2>Michael C Fontaine wrote:<BR>
> Thanks a lot for your answer. So if I'm drawing just the first PC, how can I display the scale using colorplot?<BR>
><BR>
> Thanks again for the help!<BR>
> Michaël<BR>
><BR>
> On 13 janv. 09, at 18:08, Jombart, Thibaut J X wrote:<BR>
<BR>
In such case I would not recommand using a colorplot: s.value, or s.image (if you need interpolation) from the ade4 package would be more appropriate.<BR>
<BR>
But if you want to add a legend, you could use something like:<BR>
##<BR>
par(mar=c(3,3,3,5))<BR>
xy <- matrix(runif(120),ncol=2)<BR>
z <- matrix(rnorm(60), ncol=1)<BR>
colorplot(xy,z, cex=3, xlab="", ylab="")<BR>
<BR>
<BR>
<BR>
<BR>
addLeg <- function(coords, z){<BR>
dev.w <- abs(diff(par("usr")[1:2]))<BR>
dev.h <- abs(diff(par("usr")[3:4]))<BR>
x <- coords$x<BR>
y <- coords$y<BR>
<BR>
res.xl <- rep(x, 10)<BR>
res.xr <- rep(x + 0.05 * dev.w, 10)<BR>
temp <- seq(y, dev.h*0.7, le=11)<BR>
res.yb <- temp[1:10]<BR>
res.yt <- temp[2:11]<BR>
<BR>
<BR>
rect(res.xl, res.yb, res.xr, res.yt, col=col)<BR>
<BR>
txt.x <- res.xr + 0.05*dev.w<BR>
txt.y <- 0.5*(res.yb+res.yt)<BR>
txt.lab <- seq(min(z), max(z), le=10)<BR>
text(txt.x, txt.y, lab=round(txt.lab, 3))<BR>
<BR>
return()<BR>
}<BR>
<BR>
<BR>
par(xpd=TRUE)<BR>
addLeg(locator(1),z) # here clic for the bottom left corner of the legend<BR>
<BR>
##<BR>
<BR>
Cheers,<BR>
<BR>
Thibaut.<BR>
<BR>
<BR>
--<BR>
######################################<BR>
Dr Thibaut JOMBART<BR>
MRC Centre for Outbreak Analysis and Modelling<BR>
Department of Infectious Disease Epidemiology<BR>
Imperial College - Faculty of Medicine<BR>
St Mary’s Campus<BR>
Norfolk Place<BR>
London W2 1PG<BR>
United Kingdom<BR>
Tel. : 0044 (0)20 7594 3658<BR>
jombart@biomserv.univ-lyon1.fr<BR>
<A HREF="http://biomserv.univ-lyon1.fr/%7Ejombart/">http://biomserv.univ-lyon1.fr/%7Ejombart/</A><BR>
<A HREF="http://adegenet.r-forge.r-project.org/">http://adegenet.r-forge.r-project.org/</A><BR>
</FONT>
</P>
</BODY>
</HTML>