On 10/10/2013 03:52 PM, Rebecca Stirnemann wrote:
Thanks Jim for helping,
Your sample data actually looks like my dataset. The one I put up looks
strange for some reason so please ignore that.
I have three landusenumb variables 1 2 and 3. is rep (1,2,3) correct?
When I run the following code I am getting:
> mod1 <- glmer(frat ~ flandusenumb + ground.cover_lo + (1|fsite)
,family = binomial, data= mao1)
>
> #Calculate predicted values
> newdata1 <- data.frame(ground.cover_lo = c(25,50,100), flandusenumb =
rep(1,2,3))
> pred34 <- predict(mod1,newdata=newdata1,type="response")
Error in UseMethod("predict") :
no applicable method for 'predict' applied to an object of class "mer"
Can you see what I am doing wrong?
What I am aiming for is a graph which looks like this.
Thanks
Rebecca
Okay, so you aren't getting anything to plot. You will have to send your
plot image directly to me, as R-help scrubs most images from emails.
I assume that "landusenumb" is the occasion of measurement (i.e. you
have three repeated measurements). Perhaps you could calculate the three
values from the initial model, remembering that the logarithm of the
odds of predation is your response.
Jim
______________________________________________
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.