similar to: Crosscorr.plot

Hello, R community, I have been using the lmer and mcmcsamp functions in R with some difficulty. I do not believe this is my code or data, however, because my attempts to use the sample code and 'sleepstudy' data provided with the lme4 packaged (and used on several R-Wiki pages) do not return the same results as those indicated in the help pages. For instance: > sessionInfo() R

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

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

Dear All, I am working with linear mixed-effects models using the lme4 package in R. I created a model with the lmer function including some main effects, a two-way interaction and a random effect. Now I am searching how I could incorporate an interaction between the random effect and one of the fixed effects. I tried to express the interaction in:

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

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

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, using 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 is

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

The object is a mcmc sample from lmer. I am using R v2.7.1. Please let me know what additional information I can provide, hopefully I am just making a simple mistake. Thanks in advance! > data(ratdrink, package = 'faraway') > rd.er <- lmer(wt ~ weeks*treat + (1 | subject), data = ratdrink) > rd.mc <- mcmcsamp(rd.er, 10000) > library(coda) Loading required package:

Using ranef() (package nlme, version 3.1-75) with an 'lme' object I can obtain random effects for intercept and slope of a certain level (say: 1) - this corresponds to (say level 1) "residuals" in MLWin. Maybe I'm mistaken here, but the results are identical. However, if I try to get the standardized random effects adding the paramter "standard=T" to the

hello I read "Practical Time Series" (Gareth Janacek; 2001) and they presented e.g the smoothing functions msmooth(x,k) or the bivariate function crosscorr(x,y,k), but both didn't work on my machine. I only load the ts library, is another library necessary or did this function change since 2001? Is there a more recent and detailed manual for ts? thanks, cheers Martin -- Martin

I have written a function called mcsamp() that is a wrapper that runs mcmcsamp() and automatically monitors convergence and structures the inferences into vectors and arrays as appropriate. But I have run into a very little problem, which is that mcmcsamp() shortens the variable names. For example: > set.seed (1) > group <- rep (1:5,10) > a <- rnorm (5,-3,3) > y <-

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

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

Hi, Can anyone tell me why I am not getting the correct intervals for fixed effect terms for the following generalized linear mixed model from HPDinterval: > sessionInfo() R version 2.4.1 (2006-12-18) i386-pc-mingw32 locale: LC_COLLATE=English_United States.1252;LC_CTYPE=English_United States.1252;LC_MONETARY=English_United States.1252;LC_NUMERIC=C;LC_TIME=English_United States.1252

Dear all,   I did ANOVA for many response variables (Var1, Var2, ....Var75000), and i got the result of p-value like below. Now, I want to select those predictors, which have pvalue less than or equal to 0.05 for each response variable. For example, X1, X2, X3, X4, X5 and X6 in case of Var1, and similarly, X1, X2.......X5 in case of Var2, only X1 in case of Var3 and none of the predictors in case

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

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

Dette er en melding med flere deler i MIME-format. --=_alternative 004613C000257091_= Content-Type: text/plain; charset="US-ASCII" And some more informastion I forgot. R does not crash if I write out the formula: set.seed(123) x1 <- runif(1000) x2 <- runif(1000) x3 <- runif(1000) x4 <- runif(1000) x5 <- runif(1000) x6 <- runif(1000) x7 <- runif(1000) x8 <-