[R] problem with nls()

Prof Brian Ripley ripley at stats.ox.ac.uk
Fri Nov 28 14:18:08 CET 2003


On Fri, 28 Nov 2003, Siegfried.Macho wrote:

> I wanted to use the nls() module to solve a Problem from Sen & Srivastava 
> (1990, p.209). Here is the (basic) code used to perform the estimation:
> 
> library(SenSrivastava)
> library(nls)
> data(E9.8)
> 
> # Use Linear Least Square for estimating start values
> lm.obj <- lm(R.1 ~ I.1 + S.1, data = E9.8)

Those are not the column names in E9.8!

> nls1.obj <- nls(R.1 ~ b.0 + b.1*(I.1^a.1-1)/a.1 + b.2*(S.1^a.2-1)/a.2,
>                          data = E9.8,
>                          start = list(b.0 = coef(lm.obj)[1], b.1 = 
> coef(lm.obj)[2], b.2 = coef(lm.obj)[3],                       a.1 = 1, a.2 = 1)
>                     )
> 
> unfortunately, I receive the following message from nls(): "Error in 
> object[[i]] : object is not subsettable"
> 
> Why do I get this message, what does it mean, and how can I avoid it ?
> Note 1: Enclosing expressions by using the function I() does not help either.
> Note 2: To avoid confusion I have changed the names I, R, S to I.1, R.1, S.1.

Perhaps that is your problem, as without changing the names this works for 
me in R 1.8.1.

If you still have a problem using the actual names and R 1.8.1, please
try using R's debugging facilities to find the problem on your computer.

-- 
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-help mailing list