[R-br] RES: Modelo Linear Generalizado

[Walmes Zeviani]

Podem haver outras abordagens, mas considere essa de início.

da <- with(subset(sub_data, tratamento=="A"),
           list(y=imoveis, n=moveis+imoveis))

library(bbmle)

fitA <- mle2(y~dbinom(prob=p, size=n),
             start=list(p=0.01),
             data=da,
             optimizer="nlminb",
             lower=c(p=10e-5))
confint(fitA)

prf <- profile(fitA, del=0.01)
plot(prf)

Eu tive que fazer uma restrição de borda à esquerda para que aconteça a
convergência (lower=10e-5). O valor estimado é justamente o valor que usei
na borda. Tentei valores menore mas ocorre falha (under/overflow numérico),
então mantive esse. Isso não é um problema, porque sabendo que todos os y
são zero, a estimativa amostral de p = y/n é zero (fora do espaço, claro).
De qualquer maneira, essa rotina é para calcular o IC e o confint retorna o
limite superior. Existe algumas preocupações sobre esses IC calculados para
estimativas de borda mas use em complemento ao fato da estimativa amostral
ter sido zero e que esse IC de verossimilhança é o melhor que se pode fazer
para reportar alguma incerteza sobre a estimativa.

​À
​ disposição.
Walmes.​
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20151006/1ee33cef/attachment.html>