Sorry for one more email, but here is a more concise and pertinent "ranef"
example:

# Allow S4 methods to find S3 asreml objects
setOldClass("asreml")

if(!isGeneric("ranef")){
  setGeneric("ranef", function(object, ...) standardGeneric("ranef"))
}

ranef.asreml <- function(obj){
  effects <- obj$coefficients$random
  class(effects) <- c("ranef.asreml",class(effects))
  return(effects)
}

print.ranef.asreml <- function(x, ...){
  print.default(x)
}


Kevin


On Fri, Mar 5, 2010 at 10:47 AM, Terry Therneau <therneau@mayo.edu> wrote:

> 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 the last one wins.  The coxme
> object is an S3 class BTW. (Suggestions to convert it to S4 will not be
> well recieved -- but that's a different and unnecessary thread.)
>
> Odds are that this has a simple solution which I have missed.
>     Terry Therneau
>
>
> tmt935% R --vanilla
>
> R version 2.10.0 (2009-10-26)
> Copyright (C) 2009 The R Foundation for Statistical Computing
> ISBN 3-900051-07-0
>
> R is free software and comes with ABSOLUTELY NO WARRANTY.
> You are welcome to redistribute it under certain conditions.
> Type 'license()' or 'licence()' for distribution details.
>
>  Natural language support but running in an English locale
>
> R is a collaborative project with many contributors.
> Type 'contributors()' for more information and
> 'citation()' on how to cite R or R packages in publications.
>
> Type 'demo()' for some demos, 'help()' for on-line help, or
> 'help.start()' for an HTML browser interface to help.
> Type 'q()' to quit R.
>
> > library(coxme)
> Loading required package: survival
> Loading required package: splines
> Loading required package: bdsmatrix
>
> > fit <- coxme(Surv(time, status) ~ age + (1|ph.ecog), lung)
> > ranef(fit)
> $ph.ecog
> Intercept
> 0.1592346
>
> > library(lme4)
> Loading required package: Matrix
> Loading required package: lattice
>
> Attaching package: 'lme4'
>
>
>        The following object(s) are masked from package:coxme :
>
>         fixef,
>         ranef
>
> > ranef(fit)
> Error in function (classes, fdef, mtable)  :
>  unable to find an inherited method for function "ranef", for signature
> "coxme"
> > q()
>
> ______________________________________________
> R-devel@r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-devel
>



-- 
Kevin Wright

	[[alternative HTML version deleted]]

