[R] plot legend from file

[Jim Lemon]

On 11/05/2009 09:10 PM, Kim Vanselow wrote:
> Dear r-helpers,
> I want to plot a legend (using legend())into an ordination plot (nmds produced by vegan/metaMDS). I used the output of a cluster analysis to create groups. I created a new .csv-file with 2 columns, one with the site ID used in the ordination, one with the group label derieved from the cluster dendrogram. I used this file to display the site-labels in differnt group-specified colours (using col=groupfile$grouplabel). This works. But now I want to create a legend with these groups and groupcolours from this file. I did not find anything in tutorials or help files about this. Maybe you have an idea?
> Thank you very much
> Kim
>    
Hi Kim,
Unless this plot is greatly different from others, you first need to 
know where you can put the legend. You can get the device limits in user 
coordinates from this function:

devUsr<-function() {
  usrcoord<-par("usr")
  plotin<-par("pin")
  devin<-par("din")
  plotmar<-par("mar")
  xcess<-(usrcoord[2]-usrcoord[1])*(devin[1]/plotin[1]-1)
  ycess<-(usrcoord[4]-usrcoord[3])*(devin[2]/plotin[2]-1)
  devusr<-usrcoord+c(-xcess*plotmar[2]/(plotmar[2]+plotmar[4]),
   xcess*plotmar[4]/(plotmar[2]+plotmar[4]),
   -ycess*plotmar[1]/(plotmar[1]+plotmar[3]),
   ycess*plotmar[3]/(plotmar[1]+plotmar[3]))
  return(devusr)
}

Your legend must be displayed within these coordinates. If you want it 
to extend beyond the plot, you must do this:

par(xpd=TRUE)

before the legend command. Then you should be able to use legend, 
probably using the "fill" argument to show the colors. Remember to:

par(xpd=FALSE)

when you have drawn the legend if it extends outside the plot.

Jim