Bernardo Powaga <bspowaga <at> gmail.com> writes: > I would like to fit a generalized linear model for the binomial family > with some non standard link functions. For instance, this is the > Aranda-Ordaz link: > > η = ln( ( (1 - π)^-α - 1 )/α) > >>> snip <<<< > Is there any way to tell glm() to add this > parameter in the estimation or do I have to write my own estimator with > optim()?
If the parameter cannot be made into a coefficient of the linear predictor, then I'm afraid that you will have to roll your own. > Thanks, > > BP -- Ken Knoblauch Inserm U846 Stem-cell and Brain Research Institute Department of Integrative Neurosciences 18 avenue du Doyen Lépine 69500 Bron France tel: +33 (0)4 72 91 34 77 fax: +33 (0)4 72 91 34 61 portable: +33 (0)6 84 10 64 10 http://www.sbri.fr/members/kenneth-knoblauch.html ______________________________________________ R-help@r-project.org 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.