similar to: Conflicts using Rcmdr, nlme and lme4

Hello all! R2.0.1, W2k I posted this question to the list last Sunday without getting any replies on the list. I got two off the list though, suggesting me to plot "manually" as a second step, from estimable() or intervals() objects respectively. As this was not really what I wished for, I take the risk to upset somebody with a trivial question, and re-post it (just a little

Hi, I cannot seem to get lme4 to work. I have installed the lme4 and Matrix package with apt-get. and both can be found in /usr/lib/R/site-library. When I tried an example for lmer, R could not find the function lmer(), > library(lme4) Attaching package: 'lme4' The following object(s) are masked from package:nlme : getCovariateFormula getResponseFormula

Hello all! System: R2.0.1, W2k. All packages apdated with update.packages(). I use the default graphic device for my system (no active choice). I?m trying to fine-tune scatterplot() graphs (package car), for publication. I want to control plotting symbols and colours to match plots from other, less flexible, systems. I also need to make all text and symbols bigger, to enable printing of the

Hello all! I?ve been using Rcmdr for some time, as a quick way of producing graphics and basic statistics. I run R1.8.1, OS W2000. Two days ago the dataset loader stopped working. Normally, the button <No active dataset> is clickable to give you the opportunity to choose dataset to load in the Rcmdr context. Clicking on the button now produces this: Rcmdr Version 0.9-3 Error in parse(file,

I'm trying to use the lmeSplines package together with lme4. Below is (1) an example of lmeSplines together with nlme (2) an attempt to use lmeSplines with lme4 (3) then a comparison of the random effects from the two different methods. (1) require(lmeSplines) data(smSplineEx1) dat <- smSplineEx1 dat.lo <- loess(y~time, data=dat) plot(dat.lo) dat$all <- rep(1,nrow(dat)) times20

Hi, I would like to use the effect() function (actually a slightly modified version of it) on the output of the lmer() function in the lme4 package. But the effects package requires the nlme pacvkage, which is incompatible with lme4. Workaround? ______________________________________________ Professor Michael Kubovy University of Virginia Department of Psychology for mail add: for FedEx or

Hello all, I posted a similar question recently but I suspect that it was not well enough formulated to trigger any answers. So I try again: Is there a way to produce boxplots (or something similar) that uses the estimated fixed effects of an lme{nlme} object? When I want to know the mean, I use estimable{gmodels} together with an appropriate matrix. I did look for a similar tool to estimable in

Hi, I ran the following model using nlme: model2<-lme(log(malrat1)~I(year-1982),random=~1|Continent/Country,data=wbmal10) I'm trying to run a Poisson GlMM to avoid the above transformation but I don't know how to specify the model using lmer in the lme4 library: model3<-lmer((malrat1)~I(year-1982) + ??,data=wbmal10,family=poisson) How can I introduce a random factor of the

Hello: I have an two factorial random block design. It's a ecology experiment. My two factors are, guild removal and enfa removal. Both are two levels, 0 (no removal), 1 (removal). I have 5 blocks. But within each block, it's unbalanced at plot level because I have 5 plots instead of 4 in each block. Within each block, I have 1 plot with only guild removal, 1 plot with only enfa removal,

Dear Sir/madam, I use lmer to extract model in your package of lme4. It seems works well. But the problem is when I use anova/summary the extracted model, no p-value is shown at all. In previous version(nlme), I mainly use p-value to judge which term is significant or not, and then make a decision to keep this term or not. Does it means that sth wrong with my installation of package/R? or you use

I am modeling the effects of an environmental variable (X) on fish recruitment (Y) for several, ecologically related species (i) using a mixed-effects model. The linear relationship between X and Y includes a fixed effect that is common across all species and random effects that vary by species. In the lmer() notation from the package lme4: Model<- lmer(Y ~ X +(X|i)) Because the residuals of

Hi all, I have a problem with package installing in Windows, on my PC machine. The end goal is to be able to use the lme() function. Here's what I did so far: > install.packages("lme4") Warning in install.packages("lme4") : argument 'lib' is missing: using 'C:\Users\Angela\Documents/R/win-library/2.9' --- Please select a CRAN mirror for use in this

Dear list, I am trying to fit some mixed models using packages lme4 and nlme. I did the model selection using lmer but I suspect that I may have some autocorrelation going on in my data so I would like to have a look using the handy correlation structures available in nlme. The problem is that I cannot translate my lmer model to lme: mod1<- lmer(y~x + (1|a:b) + (1|b:c), data=mydata)

Hi List,   I’ m analysing the selectivity of resting site use by forest carnivores through mixed modelling techniques and I wonder which will be the best r package to deal with several aspects simultaneously: -          binomial variable response; -          possible spatial and/or temporal correlation; I have tried nlme (lme function) and lme4 (lmer function) packages, however I realize that the

Dear R listers, I know I'm breaking the rules by asking a "homework" related question-- I hope you'll forgive me. I am a social psychology graduate student, and the only one in my department who uses R. I successfully completed my multiple regression and structural equation modeling courses using R (John Fox's car and sem packages were a big help, as was his book).

Hello all, R2.2.1, W2k I posted a similar question last week after having problem with a group update (using update.packages()), without getting any answer. Now I have updated all packages separately, and found that the problem comes from nlme: Everything comes home, the md5 checksums are ok, but then the old version can??t be removed from my computer, leaving a "nlme" library with

Thanks a lot to Demitris for a prompt answer some minutes ago on another tread (see below). To avoid excess mails on the list, I move onto next question: I have another small plotting problem that confuses me. I want to plot results from a field trial series, using the numbers of the trials as symbols in the plot. pch = as.character(trial_no) works fine, but truncates the trial number to the

This question comes up periodically, probably enough to give it a proper thread and maybe point to this thread for reference (similar to the 'conservative anova' thread not too long ago). Moving from lme syntax, which is the function found in the nlme package, to lmer syntax (found in lme4) is not too difficult. It is probably useful to first explain what the differences are between the

I'm still relatively new to R, so my apologies if this is covered somewhere. I've been running some mixed-effect models in R using lme{nlme}, but read in Faraway's recent book, Extending the Linear Model with R, that lmer in package lme4 is a much improved version. I tried using this approach, but the output for the fixed effects doesn't report a p-value or the degrees of freedom

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