[R-br] Desdobrar interação dupla e tripla com glm

Maurício Lordêlo mslordelo em gmail.com
Quarta Outubro 29 10:46:43 BRST 2014


Saudações a todos!
No mês passado, solicitei um auxílio sobre como analisar os dados de
um experimento
 realizado com explantes de bananeiras considerando três fatores. Na
ocasião descrevi o problema mais ou menos da seguinte forma (ocorreram
algumas mudanças):

"*Experimento  realizado com explantes de bananeiras considerando três
fatores:*
*1. Genótipo, com dois níveis;*
*2. Meio,  com sete níveis;*
*3. Tipo de corte, com oito níveis.*
*O objetivo é verificar qual tipo de corte induziu melhor o explante na
formação do embrião considerando os setes meios e os dois genótipos.*
*Para isso foram preparadas duas placas (caracterizando o número de
repetições) com uma quantidade de explantes idênticos de acordo com o corte
e o genótipo (esta quantidade de explantes em cada placa variou de 6 a 11).
Cada placa com estes explantes foi introduzida nos sete diferentes meios.
As respostas foram obtidas registrando-se, em relação ao total de cada
placa, quantos explantes oxidaram, quantos formaram calos, quantos não
desenvolveram, etc.*
*Os dados com uma variável resposta está neste link:*

*https://www.dropbox.com/s/ci4u8c4vn4cmk4l/dados.txt?dl=0
<https://www.dropbox.com/s/ci4u8c4vn4cmk4l/dados.txt?dl=0>  *"

Obtive retorno do Walmes que contribui com valiosas sugestões.
*O que preciso agora é realizar o desdobramento de interações duplas e
triplas utilizando "glm"*.
No CRM abaixo, apenas a interação dupla foi significativa. Porém, tenho
outras variáveis neste experimento cuja interação tripla foi significativa
e será necessário também desdobrar.
Na busca que fiz, encontrei apenas comandos de desdobramentos para modelos
ajustados com "aov".
Antecipadamente, agradeço quem possa contribuir.

Maurício


dados = read.table("dados.txt",header=TRUE, sep="\t")
str(dados)
attach(dados)

# Modelo fatorial completo, glm (quasi)binomial para variável resposta
"n_ex_ox"
model1 = glm(cbind(good=n_ex_ox, bad=n_explant-n_ex_ox)~genotipo*meio*tipo,
data=dados, family=quasibinomial)
anova(model1,test="F")
summary(model1)
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20141029/afe61395/attachment.html>


Mais detalhes sobre a lista de discussão R-br