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:
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)