From: Paul Johnson <pauljohn32_at_gmail.com>

Date: Wed, 27 Apr 2011 10:55:09 -0500

dat <- data.frame(v1,v2,v3)

rm(v1,v2,v3)

model<- lm(log(v1)~log(v2)+v3, data=dat)

Date: Wed, 27 Apr 2011 10:55:09 -0500

Hi. Comments below

On Wed, Apr 27, 2011 at 2:32 AM, agent dunham <crosspide_at_hotmail.com> wrote:

> Hi, thanks, I think I've changed the previous as you told me but I'm having

*> this error, what does it mean?
**>
**>
**> model<- lm(log(v1)~log(v2)+v3, data=dat)
**>
**> newax<- expand.grid(
**> v2 = seq(min(log(dat$v2)), max(log(dat$v2)), length=100),
**> v3= seq(min(dat$v3), max(dat$v3), length=100))
**>
**> fit <- predict(model,newax)
**>
**> graph <- regr2.plot(x=seq(min(log(dat$v2)), max(log(dat$v2)), length=100),
**> y=seq(min(dat$v3), max(dat$v3), length=100),
**> z= matrix(fit, 100, 100), main="Least-squares",
**> resid.plot= TRUE, plot.base.points= TRUE, expand=0.5,
**> ticktype="detailed", theta=-45)
**>
**>
**> Error en cbind(x, y, z, 1) %*% pmat : argumentos no compatibles
**>
**> Thanks again, user_at_host.com
**>
*

I've struggled with these 3d things before. You should supply us the full testable program to fix. Here I've fiddled around and can get a plot from persp or regr2.plot, but I understand it is not informative. But the plot does run, and maybe you will see how to fix.

Generally, I would suggest you do the separate work in separate steps, not all jammed together inside regr2.plot. If you do each bit separately, then you can inspect the result and be sure it is good.

See:

##PJ 2011-04-27

## testing to address question in r-help on persp and regr2.plot

v1 <- rnorm(100,m=100,s=10) v2 <- rnorm(100, m=1000, s=100) v3 <- rnorm(100)

dat <- data.frame(v1,v2,v3)

rm(v1,v2,v3)

model<- lm(log(v1)~log(v2)+v3, data=dat)

rv2 <- range(dat$v2)

plv2 <- seq(rv2[1], rv2[2], length.out=100)

rv3 <- range(dat$v3)

plv3 <- seq(rv3[1], rv3[2], length.out=100)

newax<- expand.grid(

v2 = plv2,

v3= plv3)

fit <- predict(model,newax)

zmat <- matrix(fit,100,100)

persp(x=plv2, y=plv3, z=zmat)

require(HH)

graph <- regr2.plot(x= plv2,

y=plv3, z=zmat, main="Least-squares", resid.plot= TRUE, plot.base.points= TRUE, expand=0.5,ticktype="detailed", theta=-45)

> ______________________________________________

*> 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.
**>
*

-- Paul E. Johnson Professor, Political Science 1541 Lilac Lane, Room 504 University of Kansas ______________________________________________ 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.Received on Wed 27 Apr 2011 - 16:02:05 GMT

Archive maintained by Robert King, hosted by
the discipline of
statistics at the
University of Newcastle,
Australia.

Archive generated by hypermail 2.2.0, at Wed 27 Apr 2011 - 17:20:34 GMT.

*
Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help.
Please read the posting
guide before posting to the list.
*