Boa noite, R's Estou com uma dúvida, Como eu faço para saber o tamanho da diferença entre dois tratamentos (no caso "Alto" e "Baixo") em um modelo linear. 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. 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. Alguém poderia me ajudar, por favor? O script é o mesmo que eu estou usando para gerar diferentes padrões e tentar entender os resultados. Abraços Luciano #-------------------------------------------------------------------------------- SD=0.1 n=50 #----- aA=0.8#intercepto Alto bA=-0#inclinação Alta aB=0.3#intercepto Baixo bB=-0#inclinação Baixo xA<-runif(n) yA<-(aA+(bA*xA)) yA<-yA+rnorm(n,0,sd=SD) xA<-xA*100 xB<-runif(n) yB<-(aB+(bB*xB)) yB<-yB+rnorm(n,0,sd=SD) xB<-xB*100 x<-c(xB,xA) y<-c(yB,yA) v<-rep(c("BAIXO","ALTO"),c(length(xB),length(xA))) plot(y~x,pch="",xlab="Eixo X",ylab="Eixo Y") points(xA,yA,col="pink",pch=16) points(xB,yB,col="grey",,pch=16) abline(lm(yB~xB),lwd=4,col="black") abline(lm(yA~xA),lwd=4,col="red") legend("topright", c("ALTO","BAIXO"), col = c("red","black"),pch=16,merge = TRUE, bg = 'gray90',lty=1) model<-lm(y~x+factor(v)) summary(model) #-------------------------------------------------------------------------------- -- Luciano F. Sgarbi Mestrando em Ecologia e Evolução - UFG Laboratório de Ecologia de Insetos Cel. (62)8174-2262 Lab. (62)3521-1732