similar to: multcomp and lme (followup)

Dear list, As a reply to my recent mail: > simint and TukeyHSD work for aov objects. > Can someone point me to similar functions for lme objects? Douglas Bates wrote There aren't multiple comparison methods for lme objects because it is not clear how to do multiple comparisons for these. I don't think the theory of multiple comparisons extends easily to lme models. One could

I suppose that multcomp in R and multicomp in S-Plus are related and it appears that it is possible to use multicomp with lme in S-Plus given the following correspondence on s-news sally.rodriguez at philips.com 12:57 p.m. 24/04/03 -0400 7 [S] LME summary and multicomp.default() Is it possible to use multicomp with lme in R and if so what is the syntax from a simple readily available

Hi, the following example aiming at a class containing a logical slot segfaults under R-1.9.0 when `gctorture(on = TRUE)' is used: Code code (dummy.c): #include <Rdefines.h> SEXP foo() { SEXP ans; PROTECT(ans = NEW_OBJECT(MAKE_CLASS("test"))); SET_SLOT(ans, install("lgl"), allocVector(LGLSXP, 1)); LOGICAL(GET_SLOT(ans,

I have run the following logistic regression model: options(contrasts=c("contr.treatment", "contr.poly")) m <- glm(wolf.cross ~ null.cross + feature, family = "binomial") where: wolf.cross = likelihood of wolves crossing a linear feature null.cross = proportion of random paths that crossed a linear feature feature = CATEGORY of linear feature with 5 levels:

New package `multcomp' for general multiple comparisons written by Frank Bretz, Torsten Hothorn and Peter Westfall We've uploaded the package `multcomp' to CRAN. The R package allows for multiple comparisons of k groups in general linear models. We use the unifying representations of multiple contrast tests, which include all common multiple comparison procedures, such as the

New package `multcomp' for general multiple comparisons written by Frank Bretz, Torsten Hothorn and Peter Westfall We've uploaded the package `multcomp' to CRAN. The R package allows for multiple comparisons of k groups in general linear models. We use the unifying representations of multiple contrast tests, which include all common multiple comparison procedures, such as the

Hello all, Quick question: I want to do posthoc contrasts for a linear mixed effects model. However, when trying to use the csimtest function in multcomp package I receive an error message saying it cannot find the function, even after installing and loading package multcomp. Any pointers would be greatly appreciated Daniel

Dear R-help list, I would like to perform multiple comparisons for lme. Can you report to me if my way to is correct or not? Please, note that I am not nor a statistician nor a mathematician, so, some understandings are sometimes quite hard for me. According to the previous helps on the topic in R-help list May 2003 (please, see Torsten Hothorn advices) and books such as Venables &

Hi, yesterday I installed the current R-devel and now I get the following warning message for program codes which made no problems under version R-1.2.3: Example: R> read.table("test.dat") Warning message: incomplete final line in: readLines(con, n, ok) where `test.dat' looks like peters at artemis:~ > cat test.dat 1 2 3peters at artemis:~ > i.e. without a newline

The package ipred is uploaded to CRAN. The main focus of the package is the calculation of improved predictors in classification tasks. Misclassification error can be improved by bootstrap aggregated classification trees and/or the framework of indirect classification. Furthermore, a unified interface for the estimation of misclassification error completes the features of ipred. We try to make

The package ipred is uploaded to CRAN. The main focus of the package is the calculation of improved predictors in classification tasks. Misclassification error can be improved by bootstrap aggregated classification trees and/or the framework of indirect classification. Furthermore, a unified interface for the estimation of misclassification error completes the features of ipred. We try to make

Hi dere, I switched over from Windows to Linux and I'm loving it. I do have a few problems with Gmod. 1. It won't work. I downloaded Wine after I downloaded Gmod, would that mess it up? Do I have to redownload it? 2. I also downloaded it from Steam, does that affect it? kthnx :D

Hi group, I'm trying to do a Tukey test to compare the means of a factor ("treatment") with three levels in an lme model that also contains the factors "site" and "time": model = response ~ treatment * (site + time) When I enter this model in csimtest, it takes all but the main factor "treatment" as covariables, not as factors (see below). Is it

Is there a simple say of extracting the shape parameter from gnlr? The return given is the actual function rather than the value. I have looked at all the values returned by names(gnlr). ie gmod<-gnlr(data,....) gmod$shape returns the function definition rather than the value found. Simple case I am using it is for censored normal data (rather simpler than survival4 it seems). The value

Hi, I use glht to make multcomp, using Tukey, from a glm model. It is possible to get a more simplified output of result? Somethink like ordering by letters. Thanks Ronaldo -- Human kind cannot bear very much reality. -- T. S. Eliot, "Four Quartets: Burnt Norton" -- > Prof. Ronaldo Reis J?nior | .''`. UNIMONTES/Depto. Biologia Geral/Lab. de Ecologia | : :' :

Dear r-helpers, I have a question about multiple testing. Here an example that puzzles me: All matrixes and contrast vectors are presented in treatment contrasts. 1. example: library(multcomp) n<-60; sigma<-20 # n = sample size per group # sigma standard deviation of the residuals cov1 <- matrix(c(3/4,-1/2,-1/2,-1/2,1,0,-1/2,0,1), nrow = 3, ncol=3, byrow=TRUE, dimnames =

This is only partly a question about R, as I am not quite sure about the underlying statistical theory either. I have fitted a non-linear mixed-effects model with nlme. In the fixed part of the model I have a factor with three levels as explanatory variable. I would like to use Tukey HSD or a similar test to test for differences between these three levels. I have two grouping factors:

Ubuntu 10.4 32bit I've tried Ubuntu proprietary drivers, and latest Nvidia drivers 256.53 or something. (gpu = evga 8600 gts 256mb) Wine versions 1.2 and up all have my issue When I run Steam games windowed or fullscreened, they run amazing. Until I get a skype/pidgin/rythmbox/email notification popup, at which point the steam app drops dramatically in performance. (1 different frame every 10

Dear R-listers I'm a new R-user who needs some help with a test that I want to do. I have done a field experiment: four treatments (cont, x, y and xy) at three sites (A, B and C), the response is count data (0 - 15). I've done a Poisson regression: >glm(response~as.factor(treatment)*as.factor(site), family=quasipoisson, offset(max.response), data=dat) The "offset" is the

Dear R-helpers: My question is how do I efficient and valid correct for multiple tests in a repeated measurement design: Suppose we measure at two distinct visits with repeated subjects a treatment difference on the same variable. The treatment differences are assessed with a mixed model and adjusted by two methods for multiple tests: # 1. Method: Adjustment with library(multcomp)