[R] Unable to retrieve residual sum of squares from nls output
Dennis Fisher
fisher at plessthan.com
Fri Aug 13 23:21:31 CEST 2010
Colleagues,
I am using "nls" successfully (2.11.1, OS X) but I am having difficulties retrieving part of the output - residual sum of squares. I have assigned the output to FIT:
> > FIT
> Nonlinear regression model
> model: NEWY ~ PMESOR + PAMPLITUDE * cos(2 * pi * (NEWX - POFFSET)/PERIOD)
> data: parent.frame()
> PMESOR PAMPLITUDE POFFSET
> 1153.02 -1183.09 24.58
> residual sum-of-squares: 1815056
>
> Number of iterations to convergence: 8
> Achieved convergence tolerance: 1.643e-08
As you can see, the object contains "residual sum-of-squares: 1815056". However, I can't figure out how to grab that value. The tried using "str":
> > str(FIT)
> List of 6
> $ m :List of 16
> ..$ resid :function ()
> ..$ fitted :function ()
> ..$ formula :function ()
> ..$ deviance :function ()
> ..$ lhs :function ()
> ..$ gradient :function ()
> ..$ conv :function ()
> ..$ incr :function ()
> ..$ setVarying:function (vary = rep(TRUE, length(useParams)))
> ..$ setPars :function (newPars)
> ..$ getPars :function ()
> ..$ getAllPars:function ()
> ..$ getEnv :function ()
> ..$ trace :function ()
> ..$ Rmat :function ()
> ..$ predict :function (newdata = list(), qr = FALSE)
> ..- attr(*, "class")= chr "nlsModel"
> $ convInfo :List of 5
> ..$ isConv : logi TRUE
> ..$ finIter : int 8
> ..$ finTol : num 1.64e-08
> ..$ stopCode : int 0
> ..$ stopMessage: chr "converged"
> $ data : language parent.frame()
> $ call : language nls(formula = NEWY ~ PMESOR + PAMPLITUDE * cos(2 * pi * (NEWX - POFFSET)/PERIOD), start = list(PMESOR = MESOR, PAMPLITUDE = AMPLITUDE, ...
> $ dataClasses: Named chr "numeric"
> ..- attr(*, "names")= chr "NEWX"
> $ control :List of 2
> ..$ maxiter : num 1000
> ..$ warnOnly: logi TRUE
> - attr(*, "class")= chr "nls"
However, none of these elements appears to contain the residual sum of squares.
Help! Thanks.
Dennis
Dennis Fisher MD
P < (The "P Less Than" Company)
Phone: 1-866-PLessThan (1-866-753-7784)
Fax: 1-866-PLessThan (1-866-753-7784)
www.PLessThan.com
More information about the R-help
mailing list