[R] Plotting NLS profiles

[Prof Brian Ripley]

On Fri, 10 Jun 2011, R Help wrote:

> Hello list,
>
> I'm trying to plot nls profiles, but the plot.profile.nls function in
> R doesn't seem to accept any plot.default variables.  Specifically,
> I'd like to be able to change the x-axis title and the colors to black
> and white.  Has anyone had any luck with this?
>
> If not, is there a way to override to plotting colors, perhaps in par()?

No.  The authors hardcoded all these.

Take a copy of the function and modify it to suit your purposes.

>
> Thanks,
> Sam
>
> fm1 <- nls(demand ~ SSasympOrig(Time, A, lrc), data = BOD)
> pr1 <- profile(fm1, alpha = 0.05)
> opar <- par(mfrow = c(2,2), oma = c(1.1, 0, 1.1, 0), las = 1)
> plot(pr1, conf = c(95, 90, 80, 50)/100) # works fine
> plot(pr1, conf = c(95, 90, 80, 50)/100,xlab=expression(theta),col=1) #
> doesn't change
>
> ______________________________________________
> 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.
>

-- 
Brian D. Ripley,                  ripley at stats.ox.ac.uk
Professor of Applied Statistics,  http://www.stats.ox.ac.uk/~ripley/
University of Oxford,             Tel:  +44 1865 272861 (self)
1 South Parks Road,                     +44 1865 272866 (PA)
Oxford OX1 3TG, UK                Fax:  +44 1865 272595