Message-ID: <alpine.LFD.2.02.1106101944380.5166@gannet.stats.ox.ac.uk>
Date: 2011-06-10T18:45:27Z
From: Brian Ripley
Subject: Plotting NLS profiles
In-Reply-To: <alpine.OSX.1.00.1106101933540.41859@tystie.local>
On Fri, 10 Jun 2011, Prof Brian Ripley wrote:
> 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.
But as they hardcoded the colours to numbers, see ?palette .
>> 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
--
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