similar to: New version of `multcomp' on CRAN

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

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

# R for Windows will not send your bug report automatically. # Please copy the bug report (after finishing it) to # your favorite email program and send it to # # r-bugs at r-project.org # ###################################################### This report is joint from Richard Heiberger <rmh at temple.edu> and Burt Holland <bholland at temple.edu>. Burt Holland is the coauthor

Hi, all: Two examples are shown below. I want to use the multiple comparison of Dunnett. It succeeded in upper case "example 1". However, the lower case "example 2" went wrong. In "example 2", the function pmvt return NaN, so I cannot show this simtest result. Is there any solution? (I changed the variable "maxpts" to a large number in front of the

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

Hi I'm encountering an error/warning when doing multiple comparisons with the package multcomp on a coxme model. My data: I'm looking at the removal of brood from the nest according to three treatments I applied on the brood. The brood and the workers caring about the brood in the nest, belonged to different colonies. Factor: treatment (3 levels: tx,uv,meta) Random effect 1: origin of

Hi, I would like to derive p-values for pair-wise comparison (Tukey's) of effects when the response is a count. I am trying a test case where y ~ Po( lambda(x) ). x has three levels : A, B and C with lambda(x) = 10, 20 and 20 respectively. Hence, p-values for the contrast C - B should distribute uniformally. I have implemented this test case as below but do not get uniform

# Your mailer is set to "none" (default on Windows), # hence we cannot send the bug report directly from R. # Please copy the bug report (after finishing it) to # your favorite email program and send it to # # r-bugs@r-project.org # ###################################################### The multcomp functions work on "lm" objects as anticipated. They do not work on

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

Hello R users, I am having difficulting getting multcomp to run. I have a dataframe attached with a numeric variable q12a and a numeric variable quota (which is really a classification variable). quota has 10 levels and unequal sample sizes. a12a has some missing data. I am interested in doing pairwise testing across the 10 quota groups on q12a. Using the ctest package the following code

Dear all, I was trying to understand how "multcomp" package works by running the examples given in the documentation. However I still don't understand when it comes to multiple comparison set by user (please refer to "Ksub" in the code). Therefore I run 2 other cases along with the original example (case 1), with the expectation I'll get the point from the output. The

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

Dear R-users I analysed the same data set by two different ways; analysis of covariance by using lm and anova functions and multiple comparison by using simtest function in the multcomp library. The output from the analysis of covariance is; > y<-lm(D~Cond+Q1,data=x) > anova(y) Analysis of Variance Table Response: D Df Sum Sq Mean Sq F value Pr(>F) Cond 2

Conditional Inference Procedures in a Permutation Test Framework The `coin' package implements a general framework for conditional inference procedures, commonly known as permutation tests, theoretically derived by Strasser & Weber (1999). The conditional expectation and covariance for a broad class of multivariate linear statistics as well as the corresponding multivariate limiting

Conditional Inference Procedures in a Permutation Test Framework The `coin' package implements a general framework for conditional inference procedures, commonly known as permutation tests, theoretically derived by Strasser & Weber (1999). The conditional expectation and covariance for a broad class of multivariate linear statistics as well as the corresponding multivariate limiting

I used the multcomp package sometime back for doing multiple comparisons. I see that it has been updated and the methods like simint are no longer supported. When I run the program it prompts to me to use glht. How do I get the lower and upper conf int and the pValues using glht? Does anyone have an example? Thanks ../Murli [[alternative HTML version deleted]]

Dear list, I have been studying R and I would like the aid of more experienced to solve the problems of the analysis below: r = gl(3, 8, label = c('r1', 'r2', 'r3')) e = rep(gl(2, 4, label = c('e1', 'e2')), 3) y = c(26.2, 26.0, 25.0, 25.4, 24.8, 24.6, 26.7, 25.2, 25.7, 26.3, 25.1, 26.4, 19.6, 21.1, 19.0, 18.6, 22.8, 19.4, 18.8, 19.2, 19.8, 21.4,

Here it is again, hope this is more clear I am using the following data (only a small subset is given): Habitat Fungus.yield Birch 20.83829053 Birch 22.9718181 Birch 22.28216829 Birch 24.23136797 Birch 22.32147961 Birch 20.30783598 Oak 27.24047258 Oak 29.7730014 Oak 30.12608508 Oak 25.76088669 Oak 30.14750974 Hornbeam 17.05307949 Hornbeam 15.32805111 Hornbeam 18.26920177 Hornbeam 21.30987049

Hello all I am working with the package multcomp but I have problems with the function simtest; the program say that can not find this function, nevertheless I doesn't have any problem with the function glht that it is in the same package. Someone knows what could be the problem? Thank you [[alternative HTML version deleted]]