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

Tiago Souza Marçal tiagosouzamarcal em hotmail.com
Quinta Fevereiro 28 13:05:25 BRT 2013


Obrigado Fernando,
funcionou perfeitamente consegui estimar a correlação genética, fenotípica e ambiental.
Att
Tiago  

Date: Thu, 28 Feb 2013 10:59:20 -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)

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 alm(Y ~ b, data = dados) # soma de quadrados da variável blm(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.



_______________________________________________
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/f1fecd34/attachment.html>


Mais detalhes sobre a lista de discussão R-br