Search Archives

Dear R-Help, When I submit an update to one of my packages, I decided to try to avoid having to fix errors that sometimes occur in CRAN's reverse-dependency checks by performing the same checks ahead of time...

I don't understand why this is a bug in usage. Is it because the 2nd argument is not named? I get the same behavior if I do name it: ===== [R version 2.6.2 (2008-02-08), Windows XP...

Package developers list: I recently added an 'as.mcmc.list' method for objects in one of my packages. Now, in doing an --as-cran check, I get this message: * checking S3 generic/method consistency ... NOTE Found the following apparent S3...

Regardless of what kind of model is used, a natural cubic spline has two end conditions -- that the second derivative is zero at each end. If there are no interior knots, then that forces it to be a straight line...

> Hello, > > I have a question related to mixed effect modeling and how to do multiple comparisons. > > We have a longitudinal study with different groups and many dependent variables such as of brain cortical volume in different areas, etc. > > I am...

An easy way to get the results you want based on the adjusted variance-covariance matrix is: library(lsmeans) lsmeans(fm1, "Machines") The adjustments and d.f. are made using pbkrtest's 'vcovAdj' and 'Lb_ddf' functions. Russ Lenth -----Original...

I received a request that I remove the 'plyr' package from the Imports for my package, because plyr is retired. Indeed, the README file for plyr states: > plyr is retired: this means only changes necessary to keep it on CRAN...

Just use lsm.both <- lsmeans(mTCFd, c("temp", "generation")) cld(lsm.both) -----Original Message----- Message: 1 Date: Wed, 27 Jul 2022 12:40:32 +0000 (UTC) From: Alexandre Santos <alexandresantosbr at yahoo.com.br> To: "r-sig-mixed-models at...

Thanks. This helped me learn a few things... 1. buildVignettes() operates on an installed package. So in RStudio, I do "build and reload" and it copies the vignette sources to the right place. Unfortunately, it does not also copy the...

CIs for the slopes in the fixed part may be obtained easily via the lsmeans package: fm <- lme(Y ~ X * fixfactor, random=~X|randomfactor) require("lsmeans") lstrends(fm, "fixfactor", var="X") You may also obtain pairwise comparisons by calling 'pairs...

Bob, I'm not talking about using svd as a generic method. I am talking about a method FOR svd results, e.g. an S3 method like foo.svd(), for which there already exist other methods, say foo.default and...

I think we tend to highly overrate the idea of trying to engineer factor codings in order to get estimates of interesting contrasts among the regression coefficients. It gets confusing when there are multiple factors, and often, only some of...

Hmmm. I guess try adding ', nesting = NULL' to the call for EMM. Sent from my iPad > On Mar 9, 2022, at 2:11 PM, Timothy MacKenzie <fswfswt at gmail.com> wrote: > > ?Thanks so much Ben and Russ. > > Russ, I get...

Dear package developers: My 'lsmeans' package provides, among other things, 'as.mcmc' and 'as.mcmc.list' methods for its 'ref.grid' objects. Those methods allow obtaining posterior predictions, or averages thereof, from certain models fitted using MCMC methods. My question...

Dear developers, I maintain a package, emmeans, that provides post-hoc analyses for a variety of models. These models are supported by writing S3 methods for recover_data and emm_grid. Quite a few models are supported internally in emmeans...

Thanks, Luke. That's more efficient for sure. I'd just have to rename about a zillion existing internal methods -- but only once! Additional advantage I see is ability to control WHICH methods I'd allow to be overridden. Russ...

Full_Name: Russell Lenth Version: 2.6.2 OS: Windows XP Pro Submission from: (NULL) (128.255.132.36) The poly() function allows a higher-degree polynomial than it should, when raw=FALSE. For example, consider 5 distinct 'x' values...

This would be straightforward using the emmeans package. Just do: library(emmeans) EMM <- emmeans(m2, ~ item_type * time) EMM # to list the estimates Then, note the order of those four estimates, and set up the desired contrasts manually: CON <- list...

Thanks to all who responded. I am pleased to say that with your help, I have managed to work around this problem by dynamically registering the methods. My file zzz.R contains code to register various methods having generics in...

All, I'm going with the previous approach, as modified by Duncan, primarily because I want to keep recover_data.lm and a few others. So if I fit a model 'm3' using stanreg::rstanarm: class(m3) [1] "stanreg" "glm...