[R-br] fatorial com efeito aleatório

[Maurício Lordêlo]

Caros "listeiros" do R,
Estou precisando analisar os dados de um experimento já realizado,
infelizmente sem a prévia consulta a um pesquisador da área de Estatística.
O experimento foi realizado com explantes de bananeiras de modo que foram
levados em consideração três fatores:
1. Genótipo, com dois níveis;
2. Doses ou meios,  com sete diferentes concentrações;
3. Tipos de cortes do explante, com dez diferentes tipos.
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
(doses). As respostas foram obtidas registrando-se, em relação ao total de
cada placa, quantos explantes oxidaram, quantos formaram calos, quantos não
desenvolveram, quantos apresentaram fungos, quantos apresentaram bactéria,
etc.
Os dados com apenas uma variável resposta estão neste link:
https://www.dropbox.com/s/23x1dr0qxoj6ghm/explantes.txt?dl=0

Comecei a esboçar um ajuste de modelo, porém surgiram dúvidas se estou
fazendo de forma correta.

explantes=read.table("explantes.txt", sep="", header=TRUE)
attach(explantes)
str(explantes)
meio=as.factor(meio1)
tipo=as.factor(tipo1)
prop=(n_oxidaram /n_explant)
tapply(prop,genotipo,mean)
tapply(prop,meio,mean)
tapply(prop,tipo,mean)

resp=cbind(n_oxidaram, n_explant - n_oxidaram)
ajuste1=glm(resp~(genotipo/meio):tipo,family=binomial())
summary(ajuste1)

Um fatorial 2 x 7 x 10 reconheço que não é adequado. Posso pensar em "tipo"
como efeito aleatório? Se sim, qual seria o comando para este caso e como
interpretaria os resultados? Existe alguma uma outra opção para a análise?

Agradeço qualquer tipo de ajuda.

Maurício
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20140917/0331b317/attachment.html>