similar to: Lmer, P-values and mixed logistic regression

I'm fitting a GLMM to some questionnaire data. The structure is J individuals, nested within I areas, all of whom answer the same K (ordinal) questions. The model I'm using is based on so-called continuation ratios, so that it can be fitted using the lme4 package. The lmer function fits the model just fine, but using mcmcsamp to judge the variability of the parameter estimates produces

(To the list moderator: I just subscribed to the list. Apologies for not having done so longer before trying to post.) Hi all, I am currently using lmer to analyze data from an experiment with a single fixed factor (treatment, 6 levels) and a single random factor (block). I've been trying to follow the online guidance for estimating p-values for parameter estimates on these and other

Dear all, I am running a mixed model using lmer. In order to obtain CI of individual coefficients I use mcmcsamp. However, I need advice which values that are most appropriate to present in result section of a paper. I have not used mixed models and lmer so much before so my question is probably very naive. However, to avoid to much problems with journal editors and referees addicted to

Dear R users, I have now tried out several options of obtaining p-values for (quasi)poisson lmer models, including Markov-chain Monte Carlo sampling and single-term deletions with subsequent chi-square tests (although I am aware that the latter may be problematic). However, I encountered several problems that can be classified as (1) the quasipoisson lmer model does not give p-values when

Hi, I have fitted a generalized linear mixed effects model using lmer (library lme4), and the family = quasibinomial. I have tried to obtain a MCMC sample, but on calling mcmcsamp(model1, 1000) I get the following error which I don't understand at all: Error in .local(object, n, verbose, ...) : Update not yet written traceback() delivers: 4: .Call(mer_MCMCsamp, ans, object) 3:

Dear R users I am trying to obtain p-values for (quasi)poisson lmer models, including Markov-chain Monte Carlo sampling and the command summary. > > My problems is that p values derived from both these methods are totally different. My question is (1) there a bug in my code and > (2) How can I proceed, left with these uncertainties in the estimations of > the p-values? > > Below

Hi, I've got a linear mixed model created using lmer: A6mlm <- lmer(Score ~ division + (1|school), data=Age6m) (To those of you to whom this model looks familiar, thanks for your patience with this & my other questions.) Anyway, I was trying this to look at the significance of my fixed effects: A6post <- mcmcsamp(A6mlm, 50000) library(coda) HPDinterval(A6post) ..but I got this

I've been trying to keep track with lmer, and now I have a couple of questions with the latest version of Matrix (0.995-4). I fit 2 very similar models, and the results are severely rounded in one case and rounded not at all in the other. > y <- 1:10 > group <- rep (c(1,2), c(5,5)) > M1 <- lmer (y ~ 1 + (1 | group)) > coef(M1) $group (Intercept) 1 3.1 2

Dear all, I am currently analyzing count data from a hierarchical design, and I?ve tried to follow the suggestions for a correct estimation of p-values as discusssed at R-Wiki (http://wiki.r-project.org/rwiki/doku.php?id=guides:lmer-tests&s=lme%20and%20aov). However, I have the problem that my model only consists of parameters with just 1 d.f. (intercepts, slopes), so that the

Hi, I would appreciate help with the following model <<1>>= gunload <- read.table(hh('datasets/gunload.dat'), header = T) gunload$method <- factor(gunload$method, labels = c('new', 'old')) gunload$physique <- factor(gunload$group, labels = c('slight', 'average', 'heavy')) gunload$team9 <- factor(rep(1:9, each = 2)) @ This

Hello, I am a chemical student and I make use of 'lme/lmer function' to handle experiments in split-plot structures. I know about the mcmcsamp and I think that it's very promissory. I would like knowing "the concept behind" of the mcmcsamp function. I do not want the C code of the MCMCSAMP function. I would like to get the "pseudo-algorithm" to understanding that

Thanks to some help by Doug Bates (and the updated version of the Matrix package), I've refined my question about fitting nested and non-nested factors in lmer(). I can get it to work in linear regression but it crashes in logistic regression. Here's my example: # set up the predictors n.age <- 4 n.edu <- 4 n.rep <- 100 n.state <- 50 n <- n.age*n.edu*n.rep age.id

Hi, I am having problems trying to assess the significance of random terms in a generalized linear mixed model using lme4 package. The model describes bird species richness R along roads (offset by log length of road log_length) as a function of fixed effects Shrub (%shrub cover) and Width (width of road), and random effect Site (nested within Site Cluster). >From reading answers to previous

I am helping my wife do some statistical analysis. She is a biologist, and she has performed some measurements on various genotypes of mice. My background is in applied mathematics and engineering, and I have a fairly good statistics background, but I am by no means a PhD level expert in statistical methods. We have used the lmer package to fit various models for the various experiments that she

How can I get prediction intervals from a mixed-effects model? Consider the following example: library(nlme) fm3 <- lme(distance ~ age*Sex, data = Orthodont, random = ~ 1) df3.1 <- with(Orthodont, data.frame(age=seq(5, 20, 5), Subject=rep(Subject[1], 4), Sex=rep(Sex[1], 4))) predict(fm3, df3.1, interval='prediction') # M01 M01

Hello, I’m trying to perform a REML analysis using the lmer() function (lme4 package). Well, it seems to work well, except that I’m not getting any p-value (see example below). Can someone tell me what I did wrong? Thanks for your help, Amélie > library(gdata) > dive <- read.xls("C:/Documents and Settings/Amelie/My Documents/Postdoc/CE 2005-2006/divebydive.xls",

Hello, I've just located the illuminating explanation by Douglas Bates on degrees of freedom in mixed models. The take-home message appears to be: don't trust the p-values from lme. Questions: Should I give up hypothesis testing for fixed effects terms in mixed models? Has my time spent reading Pinheiro & Bates been in vain? Is there a publication on this issue? Thanks, Dan Bebber

Hi folks, This is really minor but to someone not familiar with the various tentacles of the lmer package it could be really annoying. I was trying to plot the posterior density of the fixed effect parameters of a lmer model, > hr.mcmc = mcmcsamp(hr.lmer, n=50000) > densityplot(hr.mcmc, plot.points=F) There is this error, "Error in densityplot(hr.mcmc, plot.points = F) : no

Hi, I would like to obtain the dispersion parameter for a quasipoisson model for later use in calculating QAIC values for model comparison.Can anyone suggest a method of how to go about doing this? The idea I have now is that I could use the residual deviance divided by the residual degrees of freedom to obtain the dispersion parameter. The residual deviance is available in the summary

I've got a problem with "mcmcsamp" used with glmer objects produced with "lmer" from the lme4 package. When calling mcmcsamp, I get the error Error in if (var(y) == 0) { : missing value where TRUE/FALSE needed This does not occur with all models, but I can't find anything wrong with the dataset. If the error is in my data, can someone tell me what I am looking