<div dir="ltr"><div class="gmail_default" style="font-family:trebuchet ms,sans-serif">Bem, eu confesso que também não teria entendido nada se eu tivesse um tiquinho de conhecimento bem pequeno em experimento dialelo. Do pouco que sei e nas vezes que precisei, tive que gerar a matriz que você menciona e fiz isso a partir da informação dos genitores. O CMR abaixo mostra como obter uma matrix com 1 e -1 em cada linha para genitor I e genitor II (não deveria ser 0.5 e 0.5 haja visto que o genótipo resultante tem como esperado 0.5 do efeito aditivo de cada genitor? Enfim...).<br><br>## 5 genótipos. Dialelo completo.<br>gen <- gl(5,1)<br><br>## m <- diag(nlevels(gen))<br>## U <- upper.tri(m)<br>## cbind(col(m)[U], row(m)[U])<br>da <- as.data.frame(t(combn(5, 2))) ## Equivalente porém mais conveniente.<br>names(da) <- c("genitor1","genitor2")<br>da <- transform(da,<br>                genitor1=factor(genitor1, levels(gen)),<br>                genitor2=factor(genitor2, levels(gen)))<br>levels(da$genitor1)<br>levels(da$genitor2)<br><br>X1 <- model.matrix(~0+genitor1, da)<br>X2 <- model.matrix(~0+genitor2, da)<br><br>X <- X1-X2<br><br>da$y <- rnorm(nrow(X))<br><br>lm(da$y~0+X)<br></div>​<br>À<div class="gmail_default" style="font-family:trebuchet ms,sans-serif;display:inline">​ disposição.<br></div><div class="gmail_default" style="font-family:trebuchet ms,sans-serif;display:inline">Walmes.​</div><br><br></div>