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 <tiagosouzamar...@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: fernandohtol...@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 <tiagosouzamar...@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.