[R] User defined link function with extra parameters

[ken knoblauch]

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