Olá, acho que é uma pergunta bem simple... Estou utilizando o pacote glm, e queria saber uma maneira fácil de se obter os IC´s para os parâmetros, dados que o summary passa o erro padrão....fiz na mão, calculando uma aproximação pra Normal, mas queria conhecer um método mais elegante. Em tempo, estou utilizando um modelo quasipoisson, e já queria as estimativas exponenciadas, para obter as taxas. Segue o exemplo do que estou usando. OBS: A classe svyglm deriva da glm glm2.coautL <- svyglm(coaut_lattes ~ biota + tematico, family = quasipoisson, design = design.ps ) summary(glm2.coautL) Obrigado -- Fernando A.B. Colugnati Pesquisador Associado Instituto de Pesquisas em Tecnologia e Inovação - IPTI Tel. 55 11 8704-9812 www.ipti.org.br fernando@ipti.org.br
ve se o confint() funciona 2011/9/20 Fernando Colugnati <fernando@ipti.org.br>
Olá, acho que é uma pergunta bem simple...
Estou utilizando o pacote glm, e queria saber uma maneira fácil de se obter os IC´s para os parâmetros, dados que o summary passa o erro padrão....fiz na mão, calculando uma aproximação pra Normal, mas queria conhecer um método mais elegante.
Em tempo, estou utilizando um modelo quasipoisson, e já queria as estimativas exponenciadas, para obter as taxas. Segue o exemplo do que estou usando. OBS: A classe svyglm deriva da glm
glm2.coautL <- svyglm(coaut_lattes ~ biota + tematico, family = quasipoisson, design = design.ps ) summary(glm2.coautL)
Obrigado
-- Fernando A.B. Colugnati Pesquisador Associado
Instituto de Pesquisas em Tecnologia e Inovação - IPTI Tel. 55 11 8704-9812 www.ipti.org.br fernando@ipti.org.br
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
Ótimo, funcionou! *Só um aviso para quem for utilizar o svyglm (e imagino que para qualquer outro método com svy)*...Por exemplo: mod1 <- svyglm(x ~ y, ...) confint(mod1) <> confint(svyglm(x~y, ...)) O termo da direita estima os ICs corretos, respeitando o desenho amostral. Abs Em 20 de setembro de 2011 12:38, Rodrigo Coster <rcoster@gmail.com>escreveu:
ve se o confint() funciona
2011/9/20 Fernando Colugnati <fernando@ipti.org.br>
Olá, acho que é uma pergunta bem simple...
Estou utilizando o pacote glm, e queria saber uma maneira fácil de se obter os IC´s para os parâmetros, dados que o summary passa o erro padrão....fiz na mão, calculando uma aproximação pra Normal, mas queria conhecer um método mais elegante.
Em tempo, estou utilizando um modelo quasipoisson, e já queria as estimativas exponenciadas, para obter as taxas. Segue o exemplo do que estou usando. OBS: A classe svyglm deriva da glm
glm2.coautL <- svyglm(coaut_lattes ~ biota + tematico, family = quasipoisson, design = design.ps ) summary(glm2.coautL)
Obrigado
-- Fernando A.B. Colugnati Pesquisador Associado
Instituto de Pesquisas em Tecnologia e Inovação - IPTI Tel. 55 11 8704-9812 www.ipti.org.br fernando@ipti.org.br
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
-- Fernando A.B. Colugnati Pesquisador Associado Instituto de Pesquisas em Tecnologia e Inovação - IPTI Tel. 55 11 8704-9812 www.ipti.org.br fernando@ipti.org.br
para objectos glm(), o metodo confint() deve bastar... para derivados dele, deveria ser automatico, mas depende do autor da funcao. library(MASS) confint(glm(y~x)) b
participantes (3)
-
Benilton Carvalho -
Fernando Colugnati -
Rodrigo Coster