[R-br] ADAPTAR UMA MATRIZ NO R

[David Feitosa]

Tendi foi nada...



Atenciosamente,

David Feitosa

(\_(\
(=°;°)
(("")("")

Em 15 de outubro de 2014 19:29, Tiago Souza Marçal <
tiagosouzamarcal em hotmail.com> escreveu:

> 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 em 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.
>
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20141015/82112d18/attachment.html>