[R] transfer function estimation
Liaw, Andy
andy_liaw at merck.com
Wed Jan 26 03:02:47 CET 2005
> From: Peter Dalgaard
>
> "Kemp S E (Comp)" <sekemp at glam.ac.uk> writes:
>
> > I have got as far as being able to compute the residual noise, a_t.
> > However, I am slightly confused about what to do next. Reading
> > Box-Jenkins, 1976 (pp. 391) they state the following
> >
> > "....However, it seems simplest to work with a standard nonlinear
> > least squares computer program in which the derivatives are
> > determined numerically and an option is available of 'constrained
> > iteration' to prevent instability. It is then only necessary to
> > program the computation of a_t itself......."
> >
> > I know that there is a 'nls' function in R but I really do not have
> > a clue about how to use it in this situation. Perhaps Box-Jenkins
> > are confusing me with their last sentance with regard to the a_t's -
> > is this really possible?
>
> I'm not sure I'd trust any computer recommendation from 1976, no
> matter how famous the authors are. However, the hint would lead me to
> consider the optim() function.
My copy of Box/Jenkins/Reinsel (1994, 3rd ed.) has exactly the same passage
on page 430 (except an reference to chapter 7)...
Andy
> --
> O__ ---- Peter Dalgaard Blegdamsvej 3
> c/ /'_ --- Dept. of Biostatistics 2200 Cph. N
> (*) \(*) -- University of Copenhagen Denmark Ph:
> (+45) 35327918
> ~~~~~~~~~~ - (p.dalgaard at biostat.ku.dk) FAX:
> (+45) 35327907
>
> ______________________________________________
> R-help at stat.math.ethz.ch mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide!
> http://www.R-project.org/posting-guide.html
>
>
More information about the R-help
mailing list