Boa noite pessoal,

estou trabalhando com o desenvolvimento matricial dos efeitos dialélicos dos modelos de Griffing no R.

No entanto, estou tendo dificuldade para adaptar a matriz CR (mostrada abaixo) para a matriz CR1 (mostrada abaixo) que codifica

os efeitos recíprocos. Qualquer sugestão é bem vinda.

diall = structure(c(1L, 2L, 3L, 4L, 2L, 5L, 6L, 7L, 3L, 6L, 8L, 9L, 4L,

7L, 9L, 10L), .Label = c("1", "2", "3", "4", "5", "6", "7", "8",

"9", "10"), class = "factor")

mat <- model.matrix(~-1+diall)

mat[,colSums(mat)==1] = 0

CR = mat

CR1 <- structure(c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,

1, 0, 0, -1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0,

0, 0, 0, -1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0,

0, 0, 0, -1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,

0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, -1, 0, 0, 0, 0, 0, 0, 0,

0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, -1, 0, 0, 0, 0, 0, 0, 0,

0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,

0, 1, 0, 0, -1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,

0, 0), .Dim = c(16L, 10L), .Dimnames = list(c("1", "2", "3",

"4", "5", "6", "7", "8", "9", "10", "11", "12", "13", "14", "15",

"16"), c("dia1", "dia2", "dia3", "dia4", "dia5", "dia6", "dia7",

"dia8", "dia9", "dia10")), assign = c(1L, 1L, 1L, 1L, 1L, 1L,

1L, 1L, 1L, 1L), contrasts = structure(list(dia = "contr.treatment"), .Names = "dia"))

Att.

Tiago.

#################################################################

Tiago de Souza Marçal - Graduando em Agronomia pelo CCA-UFES

Bolsista de Iniciação Científica da área de Genética e Melhoramento de Plantas

#################################################################