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

 

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

_______________________________________________
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.