[Rd] nls function does not use subset argument (PR#9290)
Prof Brian Ripley
ripley at stats.ox.ac.uk
Wed Oct 11 11:42:31 CEST 2006
This is not the place to ask a question: do read the FAQ.
?nls shows that nls does have a 'subset' argument, and it does work.
Compare
> ?nls
> x <- 1:10
> y <- 2*x + 3 # perfect fit
> yeps <- y + rnorm(length(y), sd = 0.01) # added noise
> nls(yeps ~ a + b*x, start = list(a = 0.12345, b = 0.54321),
+ trace = TRUE)
1360.779 : 0.12345 0.54321
0.000778241 : 3.002946 1.999563
0.000778241 : 3.002946 1.999563
Nonlinear regression model
model: yeps ~ a + b * x
data: parent.frame()
a b
3.002946 1.999563
residual sum-of-squares: 0.000778241
> nls(yeps ~ a + b*x, start = list(a = 0.12345, b = 0.54321),
+ trace = TRUE, subset=1:9)
1056.760 : 0.12345 0.54321
0.0007055261 : 3.000838 2.000138
0.0007055261 : 3.000838 2.000138
Nonlinear regression model
model: yeps ~ a + b * x
data: parent.frame()
a b
3.000838 2.000138
residual sum-of-squares: 0.0007055261
I find your comments baffling: please provide a reproducible example of
your claims.
On Wed, 11 Oct 2006, tadakado at gmail.com wrote:
> Full_Name: Tadashi Kadowaki
> Version: 2.4.0
> OS: Redhat Linux 9
> Submission from: (NULL) (58.12.166.67)
>
>
> Doesn't nls function support subset? It seems not to work.
> And, there are no information in the online help.
> Has it sunk into oblivion?
>
> ______________________________________________
> R-devel at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-devel
>
--
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
More information about the R-devel
mailing list