R help - Nov 2003 - problem with nls()

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

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)

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.

Sincerely,
Siegfried Macho




=============================================Dr. Siegfried Macho
Department of Psychology
University of Fribourg
Rue de Faucigny 2
CH-1700 Fribourg

Tel.: +41-26-3007635
Fax.: +41-26-3009712
http://www.unifr.ch/psycho/general/macho.htm
=============================================	[[alternative HTML version
deleted]]