similar to: How to use summary.mer inside a package?

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

Hi, I am currently preparing a new version of my package papeR. When I run R CMD check using the development version of R I get the following note: Package in Depends field not imported from: ?nlme?, ?lme4?, ?survival? These packages needs to imported from for the case when this namespace is loaded but not attached. I now have problems to fix this issue. It is easy to get rid of two of the

Hi all I have problems in accessing a mer object called model.01 from a workspace that was created with R 2.8.1 and saved with save into an .RData file (on Windows XP or Ubuntu 8.10, don't remember anymore). Now I want to open it in R 2.9.0 on Ubuntu 8.10. I use # load workspace load("name.RData") which seems to work: ls() # all objects in there [1] "all"

Hello, I am trying to do contrasts after applying a binomial mixed effect model with the function lmer. I have to extract the fix effect values, but as I write fixef(model), I get this error message: Error in UseMethod("fixef") : no method for 'fixef' with objects of class "mer" Has anybody some ideas why? And how can I then post-hoc testing my data if I cannot do

Dear all, I try to run the example of lmer and get the following error message. > library(lme4) > example(lmer) lmer> (fm1 <- lmer(Reaction ~ Days + (Days | Subject), sleepstudy)) [[1]] Error in get(x, envir, mode, inherits) : variable "as.dpoMatrix" was not found This error message is similar to what I get with other models. It looks like the mer class has a slightly

Dear R-users, I am using the most updated package of lme4 (lme4_0.999375-2). I have a data set consisting of ~900 observations at ~440 independent survey sites. There are an unbalanced number of surveys at each site. I am attempting to develop several models evaluating the presence/absence of a species (PRES) at these random sites (SITE) using a number of predictor variables. The

I am making the question clear. Please help. > Dear R experts > > I was using kinship package to fit mixed model with kinship matrix. > The package looks like lme4, but I could find a way to extract p-value > out of it. I need to extract is as I need to analyse large number of > variables (> 10000). > > Please help me: > > require(kinship) > > #Generating

The coxme package has a ranef() method, as does lme4. I'm having trouble getting them to play together, as shown below. (The particular model in the example isn't defensible, but uses a standard data set.) The problem is that most of the time only one of lme4 or coxme will be loaded, so each needs to define the basic ranef function as well as a method for it. But when loaded together

What is the easiest way to extract a value from a 'mer' object from glmer? The first I need is the trtpair Std.Dev. which in this case is 0.17542? I've managed to get the fixed effects numbers from summary(fednmaout)@coeffs but no luck with the Random effects stuff. TIA Jim > fednmaout Linear mixed model fit by REML Formula: lor ~ as.factor(t2) + as.factor(t3) + (1 |

Hello, I hope that you can bare with me. I am new to models, but I think I have a pretty godd understanding of how to run them now, including how to use AICc and Anova. The issue is that I have many factors that I wish to compare so doing each one at a time would take forever. I came across the MuMIn package and I was so excited, however I am getting an error message and i don't know why.

Dear all, For a project I need to calculate the conditional AIC of a mixed effects model. Luckily, I found a reference in the R help forum for a function to be used: CAIC <- function(model) { sigma <- attr(VarCorr(model), 'sc') observed <- attr(model, 'y') predicted <- fitted(model) cond.loglik <- sum(dnorm(observed,

The object is a mcmc sample from lmer. I am using R v2.7.1. Please let me know what additional information I can provide, hopefully I am just making a simple mistake. Thanks in advance! > data(ratdrink, package = 'faraway') > rd.er <- lmer(wt ~ weeks*treat + (1 | subject), data = ratdrink) > rd.mc <- mcmcsamp(rd.er, 10000) > library(coda) Loading required package:

Hi All, I am trying to shift from running mixed models in SAS using PROC MIXED to using lme4 package in R. In trying to match the coefficients of R output to that of SAS output, I came across this problem. The dataset I am using is this one: http://support.sas.com/documentation/cdl/en/statug/63033/HTML/default/viewer.htm#statug_mixed_sect034.htm If I run the following code: proc mixed data=rc

On 16/02/2015 8:20 AM, Therneau, Terry M., Ph.D. wrote: > >> > I'm testing out a new version of coxme and R CMD check fails with "could not find function >> > ranef" (or random.effects or fixef, or anything from nlme). The NAMESPACE file has the >> > line below >> >> > importFrom(nlme, ranef, random.effects, fixef,

Hello All, I would like to retrieve some of the results from the lmer(...) function in library lme4. If I run a model, say fm.1 <- lmer(y ~ 1 + (1 | x), data = dog) and try names(fm.1), I get NULL. Is there anyway to retrieve the information? Thanks

I'm testing out a new version of coxme and R CMD check fails with "could not find function ranef" (or random.effects or fixef, or anything from nlme). The NAMESPACE file has the line below importFrom(nlme, ranef, random.effects, fixef, fixed.effects, VarCorr) and nlme is declared in the DESCRIPTION file as an import. I feel that I must be staring at some obvious (but

Dear R help list, I have done a lot of searching but have not been able to find an answer to my problem. I apologize in advance if this has been asked before. I am applying a mixed model to my data using lmer. I will use sample data to illustrate my question: >library(lme4) >library(arm) >data("HR", package = "SASmixed") > str(HR) 'data.frame': 120 obs.

I have a function (see below) that does some bootstrapping (I am happy to expand offline why I could use existing functions.) I put my results into and empty matrix and add a row of results with each iteration. My problem is i am a new user to R and I don't understand data frames, matrices, elements, and vectors well. What I would like is to have a data frame I can manipulate outside of the

Hi, I am using lmer (from the package lme4) to predict a binary response variable (REL) from a bunch of fixed effects and two random effects (Speaker_ID and NPhead_lemma): fit <- lmer(REL ~ SPEAKER_GENDER + log(SPEECHRATE) + SQSPEECHRATE + ..... + (1|Speaker_ID) + (1|NPhead_lemma), family="binomial", data=data.lmer, method="Laplace", model=T, x=T) I

Dear r-helpers, Spencer Graves and Manual Morales proposed the following methods to simulate p-values in lme4: ************preliminary************ require(lme4) require(MASS) summary(glm(y ~ lbase*trt + lage + V4, family = poisson, data = epil), cor = FALSE) epil2 <- epil[epil$period == 1, ] epil2["period"] <- rep(0, 59); epil2["y"] <- epil2["base"]