[R-br] Quadrado médio de duas características QM(y1+y2)

[FHRB Toledo]

Assim, para o seu CMR, segue o código:

genotipo<-gl(10,3)
rep<-rep(c(1,2,3),10)
resp1<-
c(249.75,227.56,368.18,615.42,756.11,481.54,610.40,534.73,488.87,587.84,627.27,594.91,358.48,341.18,409.65,530.40,644.79,619.42,592.22,588.89,627.72,605.05
,620.72,
711.55,493.93,457.53,366.90,479.98,671.50,471.44)
resp2<-
 c(604.15,527.36,711.87,548.00,59.78,464.04,529.91,436.40,277.88,186.90,478.44,454.60,480.34,238.63,537.39,430.57,352.27,256.40,563.79,862.62,521.22,420.44,
519.95,375.29,425.79,450.88,264.96,611.18,372.21,284.37)

dados<-data.frame(genotipo ,rep,resp1,resp2)

ANOVA_resp1 <- anova(lm(resp1 ~ genotipo, data = dados)) # SQ da resposta 1
ANOVA_resp2 <- anova(lm(resp2 ~ genotipo, data = dados)) # SQ da resposta 2

ANOVA_soma <- anova(lm(I(resp1 + resp2) ~ genotipo, data = dados)) # SQ da
soma

Espero ter ajudado.

att,
FH

2013/2/27 Tiago Souza Marçal <tiagosouzamarcal em hotmail.com>

>  Fernando agradeço pela ajuda e como você pediu eu estou enviando um CMR
> para um delineamento em dic:
>
> genotipo<-gl(10,3)
> rep<-rep(c(1,2,3),10)
> resp1<-
> c(249.75,227.56,368.18,615.42,756.11,481.54,610.40,534.73,488.87,587.84,627.27,594.91,358.48,341.18,409.65,530.40,644.79,619.42,592.22,588.89,627.72,605.05
> ,620.72,
> 711.55,493.93,457.53,366.90,479.98,671.50,471.44)
> resp2<-
>  c(604.15,527.36,711.87,548.00,59.78,464.04,529.91,436.40,277.88,186.90,478.44,454.60,480.34,238.63,537.39,430.57,352.27,256.40,563.79,862.62,521.22,420.44,
> 519.95,375.29,425.79,450.88,264.96,611.18,372.21,284.37)
>
> dados<-data.frame(gen,rep,resp1,resp2)
>
> Att.
>
> Tiago.
>
> ------------------------------
> Date: Wed, 27 Feb 2013 20:47:01 -0300
> From: fernandohtoledo em gmail.com
> To: r-br em listas.c3sl.ufpr.br
> Subject: Re: [R-br] Quadrado médio de duas características QM(y1+y2)
>
>
> Tiago,
>
> Somar as duas variáveis é proibido?
>
> Se, por exemplo for um DIC e seus dados estão em um data frame chamado
> dados, com a variável resposta chamada Y, faça:
>
> lm(Y ~ a, data = dados) # soma de quadrados da variável a
> lm(Y ~ b, data = dados) # soma de quadrados da variável b
> lm(Y ~ I(a + b), data = dados) # soma de quadrados da soma
>
> Vale destacar que sem CMR isso é apenas um chute.
>
> att,
> FH
>
> 2013/2/27 Tiago Souza Marçal <tiagosouzamarcal em hotmail.com>
>
>  Boa noite pessoal,
>
> estou querendo estimar a correlação genética entre duas variáveis (y1 e
> y2). Entretanto, para isso eu preciso do quadrado médio de a mais b QM
> (a+b).
>
> Alguém sabe como estimar este quadrado médio para as duas variáveis no R.
>
> Att.
>
> Tiago.
>
> _______________________________________________
> 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.
>
>
>
> _______________________________________________ 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.
>
> _______________________________________________
> 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/20130228/5d9a9857/attachment-0001.html>