Search Archives

An embedded and charset-unspecified text was scrubbed... Name: not available URL: <https://stat.ethz.ch/pipermail/r-sig-mixed-models/attachments/20140626/bb3da97c/attachment.pl>

An embedded and charset-unspecified text was scrubbed... Name: not available URL: <https://stat.ethz.ch/pipermail/r-sig-mixed-models/attachments/20140712/ba6ff3d2/attachment.pl>

An embedded and charset-unspecified text was scrubbed... Name: not available URL: <https://stat.ethz.ch/pipermail/r-sig-mixed-models/attachments/20130912/edbcee06/attachment.pl>

An embedded and charset-unspecified text was scrubbed... Name: not available URL: <https://stat.ethz.ch/pipermail/r-sig-mixed-models/attachments/20130918/0122309d/attachment.pl>

An embedded and charset-unspecified text was scrubbed... Name: not available URL: <https://stat.ethz.ch/pipermail/r-sig-mixed-models/attachments/20140724/18f62544/attachment.pl>

An embedded and charset-unspecified text was scrubbed... Name: not available URL: <https://stat.ethz.ch/pipermail/r-sig-mixed-models/attachments/20140724/19e2663f/attachment.pl>

For predictions, you might take a look at the lsmeans package and see if it provides what you need. It's pretty simple to use when you want to predict on a rectangular grid. Russ Lenth

Oops -- I meant R version 2.5.1, not 1.5.1. My apologies. -- Russell V. Lenth, Professor Department of Statistics & Actuarial Science (319)335-0814 FAX (319)335-3017 The University of Iowa russell-lenth at uiowa.edu Iowa...

From: Andrew Simmons akwsimmo at gmail.com<mailto:akwsimmo at gmail.com> > Doesn't apply already work on any dimensional arrays using argument MARGIN? Egads, yes it does! I can't explain why I didn't know that... Russ [[alternative...

OK. Well, as Noah has pointed out, I can just use apply to do what I need, and get plyr out of the picture. But for the broader question, Jeff is saying that there really are 700 packages that are...

Dear R-Devel, I noticed that if we run base::svd(x), we obtain an object of class "list". Shouldn't there be an "svd" class, in case someone (e.g., me) wants to write methods for singular value decompositions...

Dear R-devel, Consider the following clip (in R version 3.6.0, Windows): > requireNamespace("ggplot2", quietly = TRUE) Registered S3 methods overwritten by 'ggplot2': method from [.quosures rlang c.quosures rlang print.quosures rlang It seems to me that if...

Doug, With all respect, you are the one who knows exactly how you want to change the Wikipedia article, and I think you should be the one to repair it. You can't expect someone else to do it the...

Hmmmm, I thought of an approach -- a kind of manual dispatch technique. My generic is recover_data <- function(object, ...) { rd <- try(getS3method("recover_data", class(object)[1], envir = .GlobalEnv, silent = TRUE)) if (!inherits(rd, "try-error")) rd(object, ...) else UseMethod...

I made an unwise choice in version 1.2 of emmeans, by putting brms in Imports. That has the effect of expanding the dependencies in emmeans to over 100 other packages. I have rectified that in the latest github version...

Thanks, Neonira. Always something new to learn! Russ From: neonira Arinoem <neonira at gmail.com> Sent: Saturday, August 10, 2019 2:22 PM To: Lenth, Russell V <russell-lenth at uiowa.edu> Cc: I?aki Ucar <iucar at fedoraproject.org...

The emmeans package handles this quite easily. You can do something like this: model <- lm(log(y) ~ treatment + ..., ...) # where treatment is a factor library(emmeans) emm <- emmeans(model, "treatment", type = "response") emm # shows the estimated geometric means contrast(emm, "pairwise...

If I understand this correctly, I think you can use the 'trend' argument in lsmeans to get what you want. See the man page for lsmeans and the examples near the end. Russ Lenth Sent from my iPad > On Oct...

Emmeans 1.2.1 is now on its way to CRAN and it?s mirrors. It should be a lot easier to install. Russ > On May 21, 2018, at 8:47 AM, Lenth, Russell V <russell-lenth at uiowa.edu...

The lsmeans package can do the comparisons with control for each gene using something like this: lsmeans(fitted.model, trt.vs.ctrl ~ trt | gene, ref = 1) ... assuming that the control is the 1st lvel of treatment. lsmeans does not do...