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)