Assign palette (e.g. rainbow) to a series of points on 1 plot
One of the reasons we ask for a *reproducible* example, is that
it allows us to test our ideas and make sure that all the details
are taken care of. Here's a reproducible example that may help
solve your problem:
lx100=c(1,1,1,.8,.5,.4,.2,0)
day100=c(0,1,2,3,4,5,6,7)
lx90=c(1,1,1,1,.9,.8,.6,.4,.2,.1,0)
day90=c(0,1,2,3,4,5,6,7,8,9,10)
lx0=c(1,1,1,1,1,1,.9,.9,.8,.8,.6,.5,.4,.3,.2,.1,.1,0)
day0=c(0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17)
lx = list(lx0=lx0,lx90=lx90,lx100=lx100)
day = list(day0=day0,day90=day90,day100=day100)
plot(range(unlist(day)),range(unlist(lx)), type="n", xlab="Day of adult life",ylab='lx',lwd=2.2)
mapply(function(x,y,col)points(x,y,type='b',col=col),day,lx,rainbow(length(lx)))
legend('topright',names(lx),pch=1,lty=1,col=rainbow(length(lx)))
title('Survival vs. Time')
- Phil Spector
Statistical Computing Facility
Department of Statistics
UC Berkeley
spector at stat.berkeley.edu
On Sun, 29 Nov 2009, Frostygoat wrote:
I have 11 vectors representing insect survival probabilities in response to different levels of toxins at 10 concentrations lx100=c(1,1,1,.8,.5,.4,.2,0) day100=c(0,1,2,3,4,5,6,7,8) lx90=c(1,1,1,1,.9,.8,.6,.4,.2,.1,0) day90=c(0,1,2,3,4,5,6,7,8,9,10) #...and so on10% and a zero (control) series lx0=c(1,1,1,1,1,1,.9,.9,.8,.8,.6,.5,.4,.3,.2,.1,.1,0) day0=c(0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17) I want to plot them on one plot with a palette color scheme such as rainbow or topo.colors, and I want one color per concentration on both the point and line. I have found a number of ways to plot them: 1. Plot a blank frame big enough to accommodate the control data and then add the x,y coords using points. plot(x=day0,y=lx0, type="n", xlab="Day of adult life",ylab=lx,lwd=2.2,ylim=c(0.0,1),col=cols) points(x=day100,y=lx100,type="b",col="#FF8B00") points(x=day90,y=lx90,type="b",col= "#E8FF00") ... This is harder than it should be, I tracked down the color names generated with rainbow(11) and individually name each points command. 2. Bind the respective x and y coords into 2 respective matrices and plot. lxs=matrix(c(lx100,lx90,lx80,lx70,lx60...)) days=matrix(c(day100,day90,day80,day70,day60...)) plot(x=days,y=lxs, col=rainbow(11), type="b") The points rainbow (not as series) and the lines are red. I tried various methods of binding data and plotting data frames without success. I would appreciate it if someone would kindly put me on the trail. Thank you for your time.
______________________________________________ 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.