Re: [R] How to plot 3-D surface graph from lmer mixed models?

2009-08-15 Thread willow1980
Hi, First of all, I thank Professor Murdoch and Dr Sarkar for providing important information. Now, this problem has been basically solved. The package akima helps to fultil the objective originally resorting to expand.grid. Now, I can use contour or perspective 3-D plot. Hopefully, this

[R] How to plot 3-D surface graph from lmer mixed models?

2009-08-13 Thread willow1980
Dear R users, I have a problem in plotting 3 dimensional graph using mixed models. My model is sur_prop ~ afr_c+I(afr_c^2)+I(afr_c^3)+byear_c+I(byear_c^2)+I(byear_c^3)+I(byear_c^4)+(1|Studyparish)+afr_c:byear_c

Re: [R] How to plot 3-D surface graph from lmer mixed models?

2009-08-13 Thread Duncan Murdoch
willow1980 wrote: Dear R users, I have a problem in plotting 3 dimensional graph using mixed models. My model is sur_prop ~ afr_c+I(afr_c^2)+I(afr_c^3)+byear_c+I(byear_c^2)+I(byear_c^3)+I(byear_c^4)+(1|Studyparish)+afr_c:byear_c

Re: [R] How to plot 3-D surface graph from lmer mixed models?

2009-08-13 Thread willow1980
Dear Professor Murdoch, That is exactly the difficulty for me. I don't know how to make a prediction with lmer using expand.grid; at the moment, I can use “mo...@x%*%fixef(model)” to get predicted values for existing observational data, but not data by expand.grid. Actually, if I know this, I can

Re: [R] How to plot 3-D surface graph from lmer mixed models?

2009-08-13 Thread Deepayan Sarkar
On Thu, Aug 13, 2009 at 6:34 AM, willow1980jianghua@shef.ac.uk wrote: Dear Professor Murdoch, That is exactly the difficulty for me. I don't know how to make a prediction with lmer using expand.grid; at the moment, I can use “mo...@x%*%fixef(model)” to get predicted values for existing

Re: [R] How to plot 3-D surface graph from lmer mixed models?

2009-08-13 Thread Duncan Murdoch
willow1980 wrote: Dear Professor Murdoch, That is exactly the difficulty for me. I don't know how to make a prediction with lmer using expand.grid; at the moment, I can use “mo...@x%*%fixef(model)” to get predicted values for existing observational data, but not data by expand.grid. Actually, if