similar to: multcomp and lme

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 just realized that in the call to `csimint' the argument `asympt=TRUE' is missing since we need to compute the confidence intervals for a glm based on the normal approximation. Torsten --------------------------------------------------------------------- library(multcomp) set.seed(290875) # a factor at three levels group <- factor(c(rep(1,10), rep(2, 10), rep(3,10))) # Williams

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:

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:

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, 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,

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-listers, I am trying to compute simultaneous confidence intervals with simint from the package multcomp. 230 measures (abundance) have been taken in 23 sites (factor) of a data.frame (donnees: a file can be sent on request, saved with save(donnees,file="donnees")). I would like to get all pairwise comparisons with : mc<- simint(ren~ID,type="Tukey",data=donnees) I

I am using R for a class which uses S. Often I can not find functions that in S but not in R, example, multicomp(). What is your advise. Thanks. -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "body",

Dear R people, I have heard rumors of a function called (I believe) multicomp in S+ for doing multiple comparisons for linear models. It doesn't seem to be in the version of R I am using, 0.90.0. Is there something equivalent to this in R, or does it simply not exist? I see there is are more recent versions of R. For my platform, SuSE Linux, the most recent rpm available appears to be

Hi R users: I have de following data frame (called "Sx") Descripcion Nitratos Cont85g 72.40 Cont85g 100.50 Cont85g 138.30 Cont80g 178.33 Cont80g 79.01 Cont80g 74.16 Cont75g 23.70 Cont75g 15.80 Cont75g 16.20 Patron80g

Is there a fonction for multiple comparison tests (similar to "multicomp" in Splus) in a package of R? Thanks in advance for any hint... Cheers, Patrick Giraudoux University of Franche-Comté Department of Environmental Biology EA3184 af. INRA F-25030 Besançon Cedex tel.: +33 381 665 745 fax.: +33 381 665 797 http://lbe.univ-fcomte.fr [[alternative HTML version deleted]]

[Moderator: This was erronously sent to R-announce, and filtered fortunately ] Hi, This is my first time to use R. I'm wondering which function I can use to do multiple comparison. I have an lm object and want to do multiple comparison based on this object. Splus has a function multicomp() but I could not find a similar one in R. Thanks in advance. Really appreciate it. Julie

Dear list members, First of all thank you for your helpful advices. After your answeres to my firt mail I studied a lot (R-News n?5) and I tried to perform my analysis: First, to fit a GLM with a nested design I decided to use the function "lmer" in package "lme4" as suggested by Spencer Graves and Filippo Piro. I remember to you that my data were: land use classes, 3 levels

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

Hello, I'm a totally newbie to R and I'm taking a class using S+. In the class we use the multcomp command which takes a aov object and calculates confidence intervals for all pairwise differences by the Fisher least significant differences method. How can I do this in R. Thank you for taking the time with such a basic question. I've been looking on the net for a few days and I

I'm having trouble finding an R equivalent to the S-Plus "multicomp" function, which does post-hoc comparisons of treatments means in ANOVAs. Am I missing something obvious? Thanks, Peter ------------------------ Peter Adler, PhD Dept. Ecology, Evolution and Marine Biology University of California Santa Barbara, CA 93106 tel: (805) 893-7416

So I've come across a few cases where complex outputs from functions will not write to tables. The most recent case involves the TukeyHSD function in the stats package. If I save the TukeyHSD call and print it, that obviously goes fine, but when I try writing to a table, I get an error message that says "cannot coerce class \multicomp\" into dataframe. What does this mean, and

Hi, I've got a problem building a contrast matrix for the Dunnet contrast in package multcopm. The following works fine: > summary(simtest(adiff ~ trial)) Simultaneous tests: Dunnett contrasts Data: adiff by trial Contrast matrix: trial1 trial2 trial3 trial4 trial5 trial2-trial1 -1 1 0 0 0 trial3-trial1 -1 0 1 0 0

Hi I have some experimental data where I have counts of the number of insects collected to different trap types rotated through 5 different location (variable -location), 4 different chemical attractants [A, B, C, D] were applied to the traps (variable - semio) and all were trialled at two different CO2 release rates [1, 2] (variable CO2) I also have a selection of continuous variables