From: David Reitter <david.reitter_at_gmail.com>

Date: Tue 31 Jan 2006 - 09:01:33 EST

https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html Received on Tue Jan 31 09:15:23 2006

Date: Tue 31 Jan 2006 - 09:01:33 EST

I have used the same call to "predict" with "glm" models without problems. I assume I'm doing something wrong, but I have no idea what it is. If someone could help me (even by telling me how to trace this properly - is there an interactive tracer/debugger I can use?), that'd be fantastic.

Here's what I'm doing:

> summary(model)

...

Random effects:

Formula: ~log(distance) | target.utt

...

Fixed effects: primed ~ log(distance) * role * source - log
(distance):source

...

> x=10:500*0.1

> new <- data.frame(distance=x, role="r", source="m" )

> yp = predict(model, newdata=new, type="response", level=0)

Error in X %*% fixef(object) : non-conformable arguments

> traceback()

4: predict.lme(object, newdata, level = level, na.action = na.action) 3: predict(object, newdata, level = level, na.action = na.action) 2: predict.glmmPQL(model, newdata = new, type = "response", level = 0) 1: predict(model, newdata = new, type = "response", level = 0) ______________________________________________R-help@stat.math.ethz.ch mailing list

https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html Received on Tue Jan 31 09:15:23 2006

*
This archive was generated by hypermail 2.1.8
: Fri 03 Mar 2006 - 03:42:15 EST
*