Hi. Comments below On Wed, Apr 27, 2011 at 2:32 AM, agent dunham <crossp...@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, u...@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@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@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.