Dúvida em Modelo linear
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
Se x é não significativo você abandona o termo. A diferença de efeito entre níveis A e B saí do summary() do modelo, veja --- sessão carregada com objetos da mensagem anterior ---
model<-lm(y~factor(v)) # abandona x do modelo summary(model)
... Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) 0.80315 0.01388 57.85 <2e-16 *** factor(v)BAIXO -0.51275 0.01964 -26.11 <2e-16 *** ... Em vermelho a diferença estimada. À disposição. Walmes. ========================================================================== Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 VoIP: (3361 3600) 1053 1173 e-mail: walmes@ufpr.br skype: walmeszeviani twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================
Olá Walmes, E se eu tiver interesse no x (que em alguns casos isso acontece), a interpretação é a mesma? Qual seria a interpretação com o x? Abraços Luciano Em 26 de outubro de 2012 19:47, Walmes Zeviani <walmeszeviani@gmail.com>escreveu:
Se x é não significativo você abandona o termo. A diferença de efeito entre níveis A e B saí do summary() do modelo, veja
--- sessão carregada com objetos da mensagem anterior ---
model<-lm(y~factor(v)) # abandona x do modelo summary(model)
... Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) 0.80315 0.01388 57.85 <2e-16 *** factor(v)BAIXO -0.51275 0.01964 -26.11 <2e-16 *** ...
Em vermelho a diferença estimada.
À disposição. Walmes.
========================================================================== Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 VoIP: (3361 3600) 1053 1173 e-mail: walmes@ufpr.br skype: walmeszeviani twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================
_______________________________________________ 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.
-- Luciano F. Sgarbi Mestrando em Ecologia e Evolução - UFG Laboratório de Ecologia de Insetos Cel. (62)8174-2262 Lab. (62)3521-1732 https://groups.google.com/d/forum/r-ecologia
Recomendo a leitura destas fontes http://www.maths.bath.ac.uk/~jjf23/book/ http://www.leg.ufpr.br/~paulojus/embrapa/Rembrapa/ À disposição. Walmes. ========================================================================== Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 VoIP: (3361 3600) 1053 1173 e-mail: walmes@ufpr.br skype: walmeszeviani twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================
Olá Walmes, Muito obrigado pela força... Abraços, Luciano Em 27 de outubro de 2012 11:47, Walmes Zeviani <walmeszeviani@gmail.com>escreveu:
Recomendo a leitura destas fontes
http://www.maths.bath.ac.uk/~jjf23/book/ http://www.leg.ufpr.br/~paulojus/embrapa/Rembrapa/
À disposição. Walmes.
========================================================================== Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 VoIP: (3361 3600) 1053 1173 e-mail: walmes@ufpr.br skype: walmeszeviani twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================
_______________________________________________ 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.
-- Luciano F. Sgarbi Mestrando em Ecologia e Evolução - UFG Laboratório de Ecologia de Insetos Cel. (62)8174-2262 Lab. (62)3521-1732 https://groups.google.com/d/forum/r-ecologia
participantes (2)
-
Luciano F. Sgarbi -
Walmes Zeviani