search for: baayen

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 been...

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, withMCMC = TRUE) aovlmer.fnc(latinsquare.lmer,...

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 &url=http%3A%2F%2Fwww.ualberta.ca%2F~baayen%2Fpublications%2FbaayenDavidsonB ates.pdf&ei=FhqTUoXuJKKV7Abds4GYBA&usg=AFQjCNFst7GT7mBX7w9lXItJTtELJSKWJg&si g2=KGA5MHxOvEGwDxf-Gcqi6g...

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 * Z_nsize * frebase_ln + (1|word) data = x) (abreviado) = print (resultado, c=F) Fixed effects: Estimate Std. Error t value (Intercept) 6.640496 0.034490 192....

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. So question number one is : Has pvals been replaced? by pvals.fnc in the package la...

...e-test (proposed by Kruskal and Wish, 1981), to determine optimal number of dimensions. However, I am trying to find references to this two-dimensional gold-standard. Can anyone point me to authors which explicitly states that two-dimensions are typical and easiest to represent graphically? In Baayen's book (2008) one can find this statement. Are there more? Thanks, PM

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 analysis for my dissertation, I woul...

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 analysis for my dissertation, I woul...

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

...escribió: 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 * Z_nsize * frebase_ln + (1|word) data = x)  (abreviado) = print (resultado, c=F) Fixed effects:                    Estimate Std. Error t value (Intercept)        6.640496  ...

...*t*-score, but I'm not sure where I could get the standard error of the differences between coefficients (my understanding is that the standard error of the coefficients I would get from ranef(test)$Subject is not going to be exactly the same as the standard error reported for the fixed effect [Baayen, 2008:247]). Thank you very much for your feedback, Steve Politzer-Ahles -- Stephen Politzer-Ahles University of Kansas Linguistics Department http://www.linguistics.ku.edu/ [[alternative HTML version deleted]]

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, what is the appropriate syntax to get the random error correct when performing repeated measures anova with 'lme'. let's say i have 3 independent variables, with 'aov', i would write something like: aov(dep_var~(indep_var1*indep_var2*indep_var3) + Error(subject/(indep_var1*indep_var2*indep_var3)). With 'lme' however, i can't find the right formula. i tried

...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 in the book I am using (Baayen, 2008) either have continuous variables or factors with only two levels. What I find particularly confusing are the interaction terms in the output. The output doesn't present the full interaction (3 X 3) as I would expect with an ANOVA. Instead, it only presents an interaction term for one...

...ores 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 * Z_nsize * frebase_ln + > > (1|word) data = x) (abreviado) > > > > = print (resultado, c=F) > > > > >...

...specification for random effects: contradictory reults 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 &url=http%3A%2F%2Fwww.ualberta.ca%2F~baayen%2Fpublications%2FbaayenDavidsonB ates.pdf&ei=FhqTUoXuJKKV7Abds4GYBA&usg=AFQjCNFst7GT7mBX7w9lXItJTtELJSKWJg&si g2=KGA5MHxOvEGwDxf-Gcqi6g...

...xed-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 Baayen (2008). Analyzing Linguistic Data: A Practical Introduction to Statistics Using Ror on the net did not work out. Here is what I got:   > prec0_va2.lmer Generalized linear mixed model fit by the Laplace approximation Formula: prec_0 ~ (verb + agent)^2 + (1 | subject)    Data: risuvane1_binom_t...

Hola, ¿qué tal? Victor Peinado (http://nlp.uned.es/~victor/) y yo estamos organizando un grupo de lectura. Junto con quienes se nos sumen, vamos a ir leyendo el libro "Analyzing Linguistic Data: A practical introduction to Statistics" (http://www.ualberta.ca/~baayen/publications/baayenCUPstats.pdf), que trata de: - R (instalación, gráficos, etc.) - Métodos estadísticos con R (modelos lineales, clústering, clasificación, modelos mixtos) - Lingüística (que es el contexto en el que se aplica lo anterior). Está particularmente indicado para quienes tengan interé...

...ly reporting the estimates, t-values and p-values. However, since in this model there are two coefficients associated with the facA x facC interaction, I believe I should report the F-statistic in this case as this tells us whether the interaction overall is significant (e.g. as in Section 7.2.2 in Baayen's textbook). However, since the anova table is calculated stepwise, how do I decide whether facB should appear before facC in the model specification (the F values depends on the order)? Also, the contrast coefficients for facA, facB and facC in the model above are dependent on what the referen...

...to find out if there is a breakpoint after which the effect of Age of Acquisition (AOA) of the second language changes. I.e. if the slope of AOA before the breakpoint (at a certain AOA) is different from the slope past the breakpoint. For this purpose, I can use the breakpoint analysis proposed by Baayen (2008) using lm or lmer. Using that approach I clearly find a breakpoint which improves the fit of the model significantly over the simpler model (on the basis of model comparison). Given that I can model a non-linearity in a gam (from mgcv) using s(...), I'd prefer to show the presence of a b...