Não foi possível enviar o arquivo. Será algum problema com as permissões?
Diferenças
Aqui você vê as diferenças entre duas revisões dessa página.
Ambos lados da revisão anterior Revisão anterior | Próxima revisão Ambos lados da revisão seguinte | ||
disciplinas:ce225-2012-02:historico [2013/02/15 21:16] paulojus [section 4] |
disciplinas:ce225-2012-02:historico [2013/02/15 21:21] paulojus |
||
---|---|---|---|
Linha 227: | Linha 227: | ||
M 375 134 | M 375 134 | ||
------------------------------ | ------------------------------ | ||
+ | </code> | ||
+ | |||
+ | <code R> | ||
+ | |||
+ | M <- cbind(c(435, 375), c(147, 134)) | ||
+ | dimnames(M) <- list(c("F", "M"), c("S", "N")) | ||
+ | M | ||
+ | |||
+ | chisq.test(M) | ||
+ | |||
+ | addmargins(M) | ||
+ | Mesp <- outer(rowSums(M), colSums(M))/sum(M) | ||
+ | Mesp | ||
+ | |||
+ | (Chi2 <- sum(((M - Mesp)^2)/Mesp)) | ||
+ | chisq.test(M, correct=F) | ||
+ | |||
+ | |||
+ | vam <- data.frame(Y=as.vector(M), Sexo = rownames(M), Acredita=rep(colnames(M), each=2)) | ||
+ | vam | ||
+ | |||
+ | ## Modelo 0 | ||
+ | |||
+ | ## Modelo 1: (independência) | ||
+ | ## E(Y) = mu = n * Sexo * Acredita | ||
+ | ## log(mu) = log(n) + log(Sexo) + log(Acredita) | ||
+ | mod1 <- glm(Y ~ Sexo + Acredita, family=poisson(link="log"), data=vam) | ||
+ | model.matrix(mod1) | ||
+ | mod1 | ||
+ | |||
+ | fitted(mod1) | ||
+ | Mesp | ||
+ | |||
+ | resid(mod1, type="pearson") | ||
+ | (M - Mesp)/sqrt(Mesp) | ||
+ | |||
+ | sum(resid(mod1, type="pearson")^2) | ||
+ | |||
+ | |||
+ | resid(mod1, type="deviance") | ||
+ | sum(resid(mod1, type="deviance")^2) | ||
+ | |||
+ | |||
+ | ## Modelo 2: não independência | ||
+ | mod2 <- glm(Y ~ Sexo * Acredita, family=poisson, data=vam) | ||
+ | model.matrix(mod2) | ||
+ | mod2 | ||
+ | |||
+ | fitted(mod2) | ||
+ | M | ||
+ | |||
+ | anova(mod1, mod2, test="Chisq") | ||
+ | |||
+ | ## extensível a várias dimensões - modelos log-lineares | ||
+ | |||
+ | # sob link canônico X'y = X'\hat{mu} | ||
+ | # - marginais preditas iguais a observadas em mod1 | ||
</code> | </code> | ||