On Sat, May 18, 2013 at 7:05 AM, Stephen Milborrow <milbo@sonic.net> wrote:

> Paul Johnson <pauljohn32@gmail.com> wrote:
>>
>> m1 <- lm(log(y) ~ log(x), data = dat)
>>
>> termplot shows log(y) on the vertical.  What if I want y on the vertical?
>>
>
> plotmo in the plotmo package has an inverse.func argument,
> so something like the following might work for you?
>
> Thanks, I will read plotmo. I did not hear of that one before.

It looks like we have been working on same problem. Take at look at my
package "rockchalk" and run the examples for "plotSlopes" and
"plotCurves."  Exact same ideas you are thinking about.

I don't think it will help in this case because it still relies on the user
to know about  "exp" as the inverse of log. When users do predict on glm,
it "just knows" how to put predictions on the response or link scales, and
I am aiming for something automagical like that.




> library(MASS)
> library(plotmo)
> log.brain <- log(Animals$brain)
> log.body  <- log(Animals$body)
> m2 <- lm(log.brain ~ log.body)
>
> myplot <- function(...)
>    plotmo(m2, do.par=F, nrug=-1, col.resp=2, pch=20,        se=1, main="",
>  xlab="log(body)", ...)
>
> par(mfrow = c(2, 2))
> myplot(ylab="log(brain)")
> myplot(ylab="brain", inverse.func=exp)
> termplot(m2, se=T, rug=T) #  for comparison
>
>
> --
Paul E. Johnson
Professor, Political Science      Assoc. Director
1541 Lilac Lane, Room 504      Center for Research Methods
University of Kansas                 University of Kansas
http://pj.freefaculty.org               http://quant.ku.edu

	[[alternative HTML version deleted]]

