similar to: lmer t value for 3 levels of fixed factor

Hello, I just joined this list today, so am worried about proper protocol, but would like to post a question about lme4. In Baayen, Davidson, and Bates (2008), Mixed-effects modeling with crossed random effects for subjects and items, the authors describe steps for a Latin Square Design (p. 402) in which they compare 3 levels of the experimental conditions. I am considering replicating this

Dear folks (or Dear Professor Bates), I'm quite confused as to the current status of some of the available functions applicable to lmer objects. Following the examples in Baayen, Davidson, Bates (2006), my plan is to run mcmcsamp on a random effect model created by lmer in package lme4, then use the (perhaps outdated) pvals to estimate p-value. But then I couldn't find pvals anywhere.

Hola a todos, quería preguntaros un medio para obtener los valores p usando lmer. He tratado con pvals.fnc, que es lo que me habían recomendado, pero por algún motivo no está ya disponible etc. Ésta es la función que tengo, pero da las "t", sin los valores p. Aunque Baayen indica que valores por encima de 2 son significativos necesito saber las p. resultado = lmer(rt_ln ~ (fre_ln *

Hi Everyone, I am using R 2.10.1. lmer function works properly, however pvals.fnc () does not despite the fact that I uploaded: - library(lme4) - library(coda) - library(languageR) This is the error message I get pvals.fnc(lexdec3.lmerE2, nsim=10000)$fixed Error in pvals.fnc(lexdec3.lmerE2, nsim = 10000) : MCMC sampling is not yet implemented in lme4_0.999375 for models with random

Hi there, I am looking to fit a two-way random effects model to an *unblalanced* layout, y_ijk = mu + a_i + b_j + eps_ijk, i=1,...,R, j=1,...,C, k=1,...,K_ij. I am interested first of all in estimates for the variance components, sigsq_a, sigsq_b and sigsq_error. In the balanced case, there are simple (MM, MLE) estimates for these; In the unbalanced setup,

Hi All, I was wondering if someone could help me to solve this issue with lmer. In order to understand the best mixed effects model to fit my data, I compared the following options according to the procedures specified in many papers (i.e. Baayen <http://www.google.it/url?sa=t&rct=j&q=&esrc=s&source=web&cd=1&ved=0CDsQFjAA

Hi all, I'm planning to work through the book "Analyzing Linguistic Data" by R.H. Baayen, which is an introduction to R used for, well, what the title says. ;-) On the first page of the book, Baayen says that in order to work with the book, R needs to download and install a number of packages from CRAN. The problem is that one of these packages, "Design", has apparently

Hola Manuel lo he tratado de hacer pero me sale Error: unexpected string constante in: "anova(a,as,test=Chisq") no tengo ni idea de por qué... Me resulta alucinante no poder contar ya con pvals.fnc. ¿Será imposible hacerse con ello? Saludos, Miguel -------------------------------------------- El vie, 13/6/14, Manuel Azcárate <mazcarategarcia en gmail.com> escribió:

Dear R list, I have a recurring problem with the languageR package, specifically the aovlmer.fnc function. When I try to run the following code (from R. H. Baayen's textbook): # Example 1: library(languageR) latinsquare.lmer <- lmer(RT ~ SOA + (1 | Word) + (1 | Subject), data = latinsquare) x <- pvals.fnc(latinsquare.lmer,

Hi All. I'm trying to run a simple model from Baayan, Davidson, & Bates and getting a confusing error message. Any ideas what I'm doing wrong here? # Here's the data..... Subj <- factor(rep(1:3,each=6)) Item <- factor(rep(1:3,6)) SOA <- factor(rep(0:1,3,each=3)) RT <- c(466,520,502,475,494,490,516,566,577,491,544,526,484,529,539,470,511,528) priming

Hi,   I am trying to use the lmer function from the lme4 package in R 2.8.0. to fit a generalized mixed-effects model for a dependent variable with a binomial distribution (for more info on my experiment, look below). However, I encounter a major problem: How is it possible to find the general test statistic and see the relative importance of the predictors? The methods which I found described in

Existe discusión sobre el uso de los p-valores en modelos mixtos. Como se ha dicho antes, para mi lo más adecuado es comparar modelos mediante la función anova. Por Internet se puede encontrar un buen libro de Douglas Bates y en español, busca modelos mixtos con R de Luis Cayuela, enfocado hacia ecología, pero está muy bien El 13/06/2014 14:00, "Jorge I Velez"

Hello. I am new to R and new to linear mixed effects modeling. I am trying to model some data which has two factors. Each factor has three levels rather than continuous data. Specifically, we measured speech at Test 1, Test 2 and Test 3. We also had three groups of subjects: RepTP, RepNTP and NoRepNTP. I am having a really hard time interpreting this data since all the examples I have seen

Dear all, I am performing a binomial glmm analysis using the lmer function in the lme4 package (last release, just downloaded). I am using the "Laplace method". However, I am not sure about what I should do to test for the significance of fixed effects in the binomial case: Is it correct to test a full model against a model from which I remove the fixed effect I want to test

Hola a todos, Ajusté un modelo lineal mixto en el cual tengo una interacción significativa entre dos facotres, cada uno con 3 niveles (osea un diseño 3x3). Esta es la tabla de medias             P1                P 2             P3 M1      220.66       311.85        260.80 M2      348.57        89.33         191.71 M3      103.57      381.37         511.62 Con el fin de interpretar

Hello everyone, I am testing a model in which I have a two-level factor (let's call it First [1, 2]) nested under a four-level factor (let's call it Second [A, B, C, D]). I have used the following model to get coefficients representing whether, for each level of Second, there is a significant difference (in the outcome variable, Latency) between the levels of First: test <- lmer(

Dear Thierry, thank you for the quick reply. I have only one question about the approach you proposed. As you suggested, imagine that the model we end up after the model selection procedure is: mod2.1 <- lmer(dT_purs ~ T + Z + (1 +T+Z| subject), data =x, REML= FALSE) According to the common procedures specified in many manuals and recent papers, if I want to compute the p_values relative to

Hi all, I am currently doing a consistency analysis using an lmer model and trying to use parametric bootstrapping for the confidence intervals. My model is like this: model<-lmer(y~A+B+(1|C/D)+(1|E),binomial) where E is the individual level for consistency analysis, A-D are other fixed and random effects that I have to control for. Following Nakagawa and Scheilzeth I can work out the

I have installed lme4 library, but when I try something with lmer() function, I receive error message. On the other hand, I can use lme() function from the same library. Are those two the very same function or not? I am a bit confused. I am using: $platform: "i386-pc-linux-gnu" $arch: "i386" $os: "linux-gnu" $system: "i386, linux-gnu" $major: "2"

Hello lmer() - users, A call to the lmer() function causes my installation of R (2.11.1 on Mac OS X 10.5.8) to crash and I am trying to figure out the problem. I have a data set with longitudinal data of four subsequent performance measures of 1133 individuals nested in 88 groups. The data is in long format. I hypothesize a performance increase for each individual over time and intend to