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