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