R help - Jan 2005 - Colors and legend on a scatterplot?

Hello

I am new to R and I cannot overcome the 2 problems

1.Defining my own color scale for a value on a scatterplot
2.Adding legend showing values for colours

I have a set of points in area of my research with coordinates X and Y and
a value which I would like to show with colors.
 
#values for points are residuals from regression

stepmod = step(wmodzad, direction="both")

#I am adding residuals to original data to get coordinates for residuals

expzad = cbind(zadrz, residuals(stepmod))

#I want to have colors for classes defined as standard deviation from mean


std = sd(residuals(stepmod))
mn = mean(residuals(stepmod))
sdclas = seq(-3*std+mn, 3*std+mn, by=std)

#I would like to have colors from green to red, but I can't do it so I am
doing only

palette( rainbow(6) )

#I am dividing residuals into earlier defined classes by

residsd = cut(residuals(stepmod), sdclas, labels=FALSE) 

#or by (THEN POINTS DON'T APPEAR ON PLOT IN NEXT STEP)

residsd = cut(residuals(stepmod), sdclas)

# I am creating a plot

plot(expzad$X,zadrz$Y,col=residsd, pch=15 )

# I would like to add legend to plot showing values for colours from plot
but I can't do this. Simply help for 'legend' on my level of
knowledge on
R is not helpfull.

loc = locator(1)
legend(loc$x, loc$y, legend = levels(residsd) )

Thanks in advance for your help

Rafal Buczkowski

buczkowski at delta.sggw.waw.pl wrote:
> Hello
> 
> I am new to R and I cannot overcome the 2 problems
> 
> 1.Defining my own color scale for a value on a scatterplot
> 2.Adding legend showing values for colours
> 
> I have a set of points in area of my research with coordinates X and Y and
> a value which I would like to show with colors.
>  
> #values for points are residuals from regression
> 
> stepmod = step(wmodzad, direction="both")
> 
> #I am adding residuals to original data to get coordinates for residuals
> 
> expzad = cbind(zadrz, residuals(stepmod))
> 
> #I want to have colors for classes defined as standard deviation from mean
> 
> 
> std = sd(residuals(stepmod))
> mn = mean(residuals(stepmod))
> sdclas = seq(-3*std+mn, 3*std+mn, by=std)
> 
> #I would like to have colors from green to red, but I can't do it so I
am
> doing only
> 
> palette( rainbow(6) )
You could either use something like the RColorBrewer package, or specify 
the colors yourself using, e.g., rgb().


> #I am dividing residuals into earlier defined classes by
> 
> residsd = cut(residuals(stepmod), sdclas, labels=FALSE) 
> 
> #or by (THEN POINTS DON'T APPEAR ON PLOT IN NEXT STEP)
> 
> residsd = cut(residuals(stepmod), sdclas)
> 
> # I am creating a plot
> 
> plot(expzad$X,zadrz$Y,col=residsd, pch=15 )
> 
> # I would like to add legend to plot showing values for colours from plot
> but I can't do this. Simply help for 'legend' on my level of
knowledge on
> R is not helpfull.
> 
> loc = locator(1)
> legend(loc$x, loc$y, legend = levels(residsd) )
I guess you want something like

legend(loc$x, loc$y, col = residsd, legend=levels(residsd), pch=15)

Uwe Ligges
> Thanks in advance for your help
> 
> Rafal Buczkowski
> 
> ______________________________________________
> R-help at stat.math.ethz.ch mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide!
http://www.R-project.org/posting-guide.html