similar to: lme4 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

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

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.

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

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"

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,

I am intersted in time-space clustering of local industry. I made the point file, polygons and time table then run 'stkhat' function in Splancs, but it generates only time k function. I noticed that it has only one ks value. ... Can anybody help me?. thanks in advance. countiespt <- read.table('d:/dissertation/dailynew/jun10/countiespt_ok.txt', sep=",")

Dear R users, Just to clarify. I am not offering to pay someone to do my Dissertation. These 4-5 commands on Kalman Filter would be only a tiny part of my 10,000 words dissertation. A part that even after trying for a few days, I am still stuck on. I am offering ?50, just to say thanks. Regards -- View this message in context:

Dear All, With kind help from several friends on the list, I am getting close. Now here are something interesting I just realized: for random effects, lmer reports standard deviation instead of standard error! Is there a hidden option that tells lmer to report standard error of random effects, like most other multilevel or mixed modeling software, so that we can say something like "randome

I''m using a plugin which adds a method to ActiveRecord::Base http://similetimelinerailshelper.googlecode.com/svn/trunk/simile_timeline/ This adds an "acts_as" type declaration to the Models, via the file in vendor/plugins/simile_timeline/lib/simile_timeline.rb Rails is working with these declarations, and they are being used successfully in ''regular

Hi!, I'm really hoping someone out there will be able to help me. I recently started my MSc dissertation on Population Projection Matrices, which has been going well until now. I am trying to set-up a general script that does a pairwise comparison of all elements in my matrices. So for example, given that I have the following matrix S: > S [,1] [,2] [,3] [1,]

Hello, I am trying to add R code to my dissertation and the length of a line of text that Sweave generates is longer than the right margin in my dissertation. I have found no way to tell Sweave how to set the line length. I am using an iMac with tiger ver. 10.4.10, and R: ver. 2.5.1. Here is the R code that is longer than my right margin: input to Sweave:

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 *

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

Hey Everyone, I wrote up a library for my dissertation work, and I'd love to include it as an appendix. I have found ways to export the Rd files as html, txt and pdf/ps, but not a raw tex export. I've had to resort to dumping it into text files, and just using verbatim in my LaTeX code. Does anybody have any recommendations or experiences with dealing with this? Please CC me on

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

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've just finished my Master's dissertation of Political Science about FOSS. In the third chapter i used of examples of Samba community. That is the reason why i want to share this with you and you might find the text at the following addresses : in portuguese : http://bit.ly/2Egkbv in french : http://bit.ly/2VJLPp It was quite difficult for me to talk about technology and

Hi, can anyone tell me how to nest two fixed factors using glmer in lme4? I have a split-plot design with two fixed factors - A (whole plot factor) and B (subplot factor), both with two levels. I want to do GLMM as I also want to include different plots as a random factor. But I am interested on the effect of A a B and their interaction on the response variable. I tried

Actually, I re-read the post and think it needs clarification. We may both be right. If the question is "I am building a model and want to know if I should retain this random effect?" (or something like that) then the LRT should be used to compare the fitted model against another model. This would be accomplished via anova(). In other multilevel programs, the variance components are