Boa noite, R's<div>Estou com uma dúvida,</div><div>Como eu faço para saber o tamanho da diferença entre dois tratamentos (no caso "Alto" e "Baixo") em um modelo linear.</div><div>No exemplo abaixo, não existe efeito do "x" sobre o "y", entretanto existe efeito de "t" (diferença entre os tratamentos "Alto" e "Baixo"). E é nisso que eu estou interessado.</div>

<div>Qual seria a forma mais adequada para encontrar essa diferença? Eu usei "lm(y~x+factor(t))" mas não consegui chegar a nada.</div><div>Alguém poderia me ajudar, por favor?</div><div><br></div><div>O script é o mesmo que eu estou usando para gerar diferentes padrões e tentar entender os resultados.</div>

<div>Abraços</div><div>Luciano</div><div><br></div><div><div>#--------------------------------------------------------------------------------</div><div>SD=0.1</div><div>n=50</div><div>#-----</div><div>aA=0.8#intercepto Alto</div>

<div>bA=-0#inclinação Alta</div><div>aB=0.3#intercepto Baixo</div><div>bB=-0#inclinação Baixo</div><div><br></div><div>xA<-runif(n)</div><div>yA<-(aA+(bA*xA))</div><div>yA<-yA+rnorm(n,0,sd=SD)</div><div>xA<-xA*100</div>

<div>xB<-runif(n)</div><div>yB<-(aB+(bB*xB))</div><div>yB<-yB+rnorm(n,0,sd=SD)</div><div>xB<-xB*100</div><div>x<-c(xB,xA)</div><div>y<-c(yB,yA)</div><div>v<-rep(c("BAIXO","ALTO"),c(length(xB),length(xA)))</div>

<div>plot(y~x,pch="",xlab="Eixo X",ylab="Eixo Y")</div><div>points(xA,yA,col="pink",pch=16)</div><div>points(xB,yB,col="grey",,pch=16)</div><div>abline(lm(yB~xB),lwd=4,col="black")</div>

<div>abline(lm(yA~xA),lwd=4,col="red")</div><div>legend("topright", c("ALTO","BAIXO"), col = c("red","black"),pch=16,merge = TRUE, bg = 'gray90',lty=1)</div>

<div><br></div><div><br></div><div>model<-lm(y~x+factor(v))</div><div>summary(model)</div></div><div><br></div><div>#--------------------------------------------------------------------------------</div><div><div><br>
</div>
-- <br><blockquote style="margin:0 0 0 40px;border:none;padding:0px"></blockquote><div>Luciano F. Sgarbi</div><div>Mestrando em Ecologia e Evolução - UFG</div><div>Laboratório de Ecologia de Insetos</div><div>Cel. (62)8174-2262   Lab. (62)3521-1732</div>

<div><br></div><div><br></div><br>
</div><div><br></div>