[R] plotting log regression

Andrews, Chris chrisaa at med.umich.edu
Mon Dec 17 14:14:28 CET 2012


Your example seems strange because a line fits on the x-y scale; not on the log(x)-log(y) scale.  Anyway, here is my example.  You can build on it for more general data.

x <- exp(1:10)
y <- exp(10:1 + rnorm(10))


logmod <- lm(log(y)~log(x))
logypred <- predict(logmod)


plot(y~x)
lines(exp(logypred)~x, col=2)


plot(y~x, log="xy")
lines(exp(logypred)~x, col=2)

Chris


-----Original Message-----
From: Rui Barradas [mailto:ruipbarradas at sapo.pt] 
Sent: Friday, December 14, 2012 12:06 PM
To: mrkooper
Cc: r-help at r-project.org
Subject: Re: [R] plotting log regression

Hello,

How can you expect to see the fit line if you are ploting x and y values, not their logarithms?
And your definitions of x and y are wrong, they should use c().

x <- c(1,2,3,4,5)
y <- c(6,7,8,9,10)

plot(log(x), log(y))

fit <- lm(log(y) ~ log(x))  # Same as glm
abline(fit)


Please read R-intro, in the doc folder of your installation of R.

Hope this helps,

Rui Barradas
Em 14-12-2012 13:58, mrkooper escreveu:
> I want to plot the regression line of a log regression into a plot 
> with my normal, nonlog, data.
>
> for example
>
> x <- (1,2,3,4,5)
> y <- (6,7,8,9,10)
>
> plot (x,y)
>
> I tried a log-regression by
>
> a <- glm (log(y) ~ log(x))
>
> and then i tried to insert the answer to my graph, where the standard 
> values are shown:
> abline (a, col="red")
>
> unfortunately it does not work....
>
>
>
>
> --
> View this message in context: 
> http://r.789695.n4.nabble.com/plotting-log-regression-tp4653087.html
> Sent from the R help mailing list archive at Nabble.com.
>
> ______________________________________________
> R-help at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide 
> http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.


**********************************************************
Electronic Mail is not secure, may not be read every day, and should not be used for urgent or sensitive issues 




More information about the R-help mailing list