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@gmail.com
To: r-br@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@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@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@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.