[R-br] regressão linear pelo origem

Gilson Sanchez gilsonsch em gmail.com
Sexta Julho 8 15:54:24 BRT 2011


Pessoal mais uma vez procuroa juda de vc's, bom a verdade é pra um colega
ele é novo com o R, se alguem puder ajudar.

Estou precisando da ajuda de você, sou novo no grupo e estou tentando fazer
um grafico de regressao linear simples ajustando o intercepto da reta deve
passar pela origem dos eixos (0,0) mas não estou conseguindo.

O CMR que estou trabalhando é o seguinte.

# REGRESSÃO EI30 VS DS

require(XLConnect)
reg<-loadWorkbook("Estadistica Final.xls", create= FALSE)
reg1<-readWorksheet(reg, sheet="Plan1")
reg1
attach(reg1)
cor.test(EI30, DS, alternative = c("two.sided"),
         method = c("pearson"), conf.level = 0.95)
reglin1<-lm(DS~EI30, data=reg1)
anova(reglin1)
reglin1
summary(reglin1)

# GRAFICOS
par(mfrow = c(3,1), mar=c(4.5,5,2,2))
plot(EI30,DS, type="p", pch=15, lwd=1.5,
 ylim=c(0,100), xlim=c(0,50), cex.main=1.2, cex.lab=0.9, cex.axis=0.9,
 xlab="",
 ylab="",cex.lab=2,
 main="Descoberto", font.main=2)
abline(reglin1, lwd=2.0, col="black", lty=1) ## pode-se usar  lty=9
legend("bottomright", bty="n", lwd=c(1.5), pch=c(15), col=c("black"),
legend=c(expression(y=='1,516x'~+~'9,241'~','~R^2=='0,58'^'*')), cex=1)
plot(EI30, LS, pch=16, col="black",
     ylim=c(0,3), xlim=c(0,50), cex.main=1.2, cex.lab=0.9, cex.axis=0.9,
     xlab="",
     ylab=expression(bold(Perda~do~solo~(g~L^{-1}))), font.lab=2,
     main="Liteira", font.main=2)
abline(reglin2, lwd=2.0, col="black", lty=2)
legend("bottomright", bty="n", lwd=c(1.5), pch=c(16), col=c("black"),
legend=c(expression(y=='0,018x'~+~'0,323'~','~R^2=='0,16'^'*')), cex=1)


Aguardo suas consideraçoes.

Att.

Omar Cubas Encinas
Instituto Nacional de Pesquisas da Amazônia / INPA
Aluno de Mestrado em Agricultura do Tropico Úmido / ATU
Manaus - AM - Brasil
(092) 9102-2542





-- 
MSc. Gilson Sánchez Chia

Laboratório de Fisiologia Vegetal

Embrapa Amazônia Ocidental

Fone: (92) 3303-7841
---
Antes de imprimir, pense na sua responsabilidade com o Meio Ambiente.
Please consider the environment before printing this email.
 Renew, Reduce and Recycle. The Planet Knows.
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110708/f6a5ceac/attachment.html>


Mais detalhes sobre a lista de discussão R-br