[R] Modify the plot from countCDFxt

Michelle DePrenger-Levin DePrengM at botanicgardens.org
Fri Oct 16 22:53:23 CEST 2009


Hello,

I am running the plot from countCDFxt (popbio). I would like to report the
y-axis as a percent instead of the log scale (e^01...). I can add an axis
with axis(2, 0:1, line =2) but I'm having trouble understanding how to assign
the tic marks (with 'at ='). I'd like to tell it to make tics with the
percent value (or decimal if that's my only option) at the equivalent of
e^-1, e^-2, etc, but the values don't seem to line up correctly. I'd also
like to replace the provided axis with my axis instead of just placing my
axis out a line but I'm unsure how to do that. 

library(popbio)
logNE <- log(erbr$NEast[-1]/erbr$NEast[-6])
countCDFxt(mean(logNE), var(logNE), nt=5, Nc=1317, Ne=200)
axis(2, 0:1, labels = c(0.018, 0.050, 0.135, 0.368, 1), at = c(0.018, 0.050,
0.135, 0.368, 1), line =1)

Year        NEast NWest
2004        731   1732
2005        898   2004
2006        714   1130
2007        1748  1722
2008        1901  1661
2009        1317  1563

1. How do I tell axis() to report as a percent?
2. Why don't my values in label and at match the given axis?
3. How do I replace the axis label in countCDFxt?

Thanks for any help,

Michelle




More information about the R-help mailing list