similar to: changes in the structure of mer objects?

Maybe I do not explain well what I would like to do. I do not want to change the labels of the axis, but the scale. What I want is a general procedure for changing the scale. Its like using a logarithmic scale on a plot. Labels are the same, but the increases of x along the x-axis are defined by a known monotone and continuous function. Florent Bresson ----- Message d'origine ---- De :

Dear members, I would like to switch from nlme to lme4 and try to translate some of my models that worked fine with lme. I have problems with the pdMat classes. Below a toy dataset with a fixed effect F and a random effect R. I gave also 2 similar lme models. The one containing pdLogChol (lme1) is easy to translate (as it is an explicit notation of the default model) The more parsimonious

Hi, I' trying to learn how to use lme for Linear Mixed Models and I have a problem when I have to include more than one random effect in my model. I know that this could be a stupid question to ask, but I'm not able to solve it by myself... One example: if my model is response = operator + block + day with operator and block as fixed effects and day as random effect, I use res.lme

Hi, I have a lattice xyplot that contains panels according to FactorA, and curves for the 2 levels of Factor B within a panel. I try to add text in the panels of a lattice graph. I suppose I have to write a custom function (panel.txt). What I really would like is to adapt the text in the panel according to the levels of FactorA. In the manuals, I find examples for the strips using which.given

Dear all, I'm trying to store/extract the mean& standard error of the fixed effects parameter and the variance of the random effects parameter from "lmer" procedure from mlmre4 package developed by bates n pinheiro. while storing fixed effects parameter is straight forward, the same is not true for storing the variance parameter of the random effects. kindly help me ~prabhu

Could anybody explain this difference in the function order() between R under windows and R under Linux ? It seems that under Linux the order in of character dsitinguishes between caps and normal letters and sorts them starting with the capitals (see below). How can I avoid this (I mean: get the same results as I get under Windows)? Thanks Joris Under Windows R 1.6.2 > mstring <-

Dear all, OS: x86_64-suse-linux 9.2 CPU: Intel(R) Xeon(TM) CPU 3.20GHz R-version: R-2.1.0 I've started using a new Linux server, upgraded at the same time to R-2.1.0 (see above) and have problems with some elementary analysis that ran without a problem on my previous configuration. anova.glm gives the following error: Error in La.chol2inv(x, size) : lapack routines cannot be loaded This

Hi, I have been having a tedious issue with lmer models with lots of factors and lots of levels. In order to get the basic information at the beginning of the print out I also have to generate these enormous tables as well. Is there a method command to leave off all of the effects and correlations? Or, do I have to go to string commands?

After loading the package "Matrix" (version 0.995-12), using the summary function with an lme (package nlme version 3.1-75) object results in an error message saying Fehler in dim(x) : kein Slot des Namens "Dim" f?r dieses Objekt der Klasse "correlation" (translated: 'Error in dim(x) : no slot of the name "Dim" for this object of class

I'm using FC4 and R 2.3.1 to fit a mixed effects logistic regression. The response is 0/1 and both the response and the age are the same for each pair of observations for each subject (some observations are not paired). For example: id response age 1 0 30 1 0 30 2 1 55 2 1 55 3 0 37 4 1 52 5 0 39 5 0 39 etc. I get the

LS, I'm estimating multilevel regression models, using the lme-function from the nlme-package. Let's say that I estimated a model and stored it inside the object named 'model'. The summary of that model is shown below: Using summary(model)$tTable , I receive the following output: > summary(model)$tTable Value Std.Error DF t-value

Full_Name: Din Chen Version: 2.4.0 OS: Windows XP Submission from: (NULL) (66.17.122.18) I just updated the R.2.4.0. and got the error message for random effect modelling, which was working on R.2.3.1. library(lme4) mmod <- lmer(bright ~ 1+(1|operator), pulp) summary(mmod) Then when I tried to extract the residuals and random effect using: resid(mmodr) ranef(mmodr) I got error message:

Can anybody explain me what is going wrong with the rbind command or how I could avoid the problem? I have a script with loops which time has a dataframe as output (called T2sbdata). After each loop I would like to append that dataframe to the results I obtained from the previous loops (called T2sbdataALL) and use for that purpose rbind. Below a part of the loop and the result I get. The

Hello, Could somebody please tell me what am I doing wrong in following? I try extract coefficients (using arm-package) from the lmer frunction, but I get the following warning: a<-data.frame(coef(res)) Error in as.data.frame.default(x[[i]], optional = TRUE, stringsAsFactors = stringsAsFactors) : cannot coerce class "coef.mer" into a data.fram I think I have done it before

Hi All. I need a bit of help extracting the residual error variance from the VarCorr structure from lmer. #Here's a 2-way random effects model lmer.1 <- lmer(rating ~ (1|person)+(1|rater), data = dat) #Get the structure vc.fit <- VarCorr(lmer.1) #results in..... $person 1 x 1 Matrix of class "dpoMatrix" (Intercept) (Intercept) 0.7755392 $rater 1 x 1 Matrix

All, How does one extract the level-1 variance from a model fit via lmer()? In the code below the level-2 variance component may be obtained via subscripting, but what about the level-1 variance, viz., the 3.215072 term? (actually this term squared) Didn't see anything in the archives on this. Cheers, David > fm <- lmer( dv ~ time.num*drug + (1 | Patient.new), data=dat.new )

I am having a problem extracting from "mer" objects.    I have constructed my problem using existing datasets.   Using the following commands:   require(lme4) fm1 <- lmer(Yield ~ 1 + (1 | Batch), Dyestuff) fixef(fm1) I get the following error message: "Error in UseMethod("fixef") : no applicable method for "fixef""   I know that "fixef" is in

Trying to compile the package Matrix_0.9975-11.tar.gz with newest R-2.5.0 alpha (2007-03-31 r40986) on FreeBSD 7.0-CURRENT (i386) I get the following error: ----- R CMD INSTALL Matrix_0.9975-11.tar.gz * Installing to library '/usr/local/lib/R/library' * Installing *source* package 'Matrix' ... ** libs ** arch - "Makefile", line 10: Missing dependency operator

Hi helpeRs, I've been working through some examples in Pinhiero & Bates( 2000) trying to understand how to translate to the new Lme4 syntax but without much luck. Below is what I think I should do, but either the answers don't come out the same or I get errors. In the Oxide problems I'm particularly interested in obtaining the levels coeficients but this options no longer seems

Dear All, I have a data as follows: ID <- 1:100 Y <- rnorm(100) X <- rnorm(100) type <- as.factor(rep(1:3,100,time=1)) df <- as.data.frame(cbind(ID, Y,X,type)) I want to plot Y versus X by specifying the pch to be as follows: Subjects having type = 1 must be plotted with a character pch = "A" Subjects having type = 2 must be plotted