I don't know how to reparameterise the model. will you please give the detail for reparameterisation and delta method? Thanks
--- On Tue, 7/13/10, Rubén Roa <[email protected]> wrote: From: Rubén Roa <[email protected]> Subject: RE: [R] ed50 To: "Dipa Hari" <[email protected]>, [email protected] Date: Tuesday, July 13, 2010, 6:27 AM > -----Mensaje original----- > De: [email protected] > [mailto:[email protected]] En nombre de Dipa Hari > Enviado el: lunes, 12 de julio de 2010 22:19 > Para: [email protected] > Asunto: [R] ed50 > > I am using semiparametric Model > library(mgcv) > sm1=gam(y~x1+s(x2),family=binomial, f) > > How should I find out standard error for ed50 for the above > model ED50 =( -sm1$coef[1]-f(x2)) / sm1$coef [2] > > f(x2) is estimated value for non parametric term. > > Thanks Two ways, 1) Re-parameterise the model so that ed50 is an explicit parameter in the model, or 2) Taylor series (aka Delta method) using the standard errors of coef[1], coef[2] and their correlation. HTH ____________________________________________________________________________________ Dr. Rubén Roa-Ureta AZTI - Tecnalia / Marine Research Unit Txatxarramendi Ugartea z/g 48395 Sukarrieta (Bizkaia) SPAIN [[alternative HTML version deleted]]
______________________________________________ [email protected] mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
