similar to: speed issue with 'multcomp' ?

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

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

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

Dear useRs, `multcomp' version 0.991-1 will be shortly available from CRAN near you. Nearly all functionality contained in the package has been re-implemented from scratch. The focus of the package has been extended to general linear hypotheses in arbitrary parametric models and the most important function to check out is `glht()'. Multiple comparison of means procedures (for example

Dear useRs, `multcomp' version 0.991-1 will be shortly available from CRAN near you. Nearly all functionality contained in the package has been re-implemented from scratch. The focus of the package has been extended to general linear hypotheses in arbitrary parametric models and the most important function to check out is `glht()'. Multiple comparison of means procedures (for example

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

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

# 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

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

# 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

I am trying to use multcomp to do a Tukey posthoc on growth increments among genetic crosstypes. #Fixed effect model m1 <- lm(inc ~ 0 + Age+ Crosstype + Sex, data = Data.age) summary(m1) RESULTS of the model: summary(m1) Call: lm(formula = inc ~ 0 + Age + Crosstype + Sex, data = Data.age) Residuals: Min 1Q Median 3Q Max -0.87180 -0.34002 -0.02702 0.27710 2.17820

Hi! I use simtest fonction of multcomp package to compile a Dunnett's test. I have 10 treatments and one control group, so i create a matrix with: m<-matrix(0,10,11) m[1,1]<--1 m[1,2]<-1 m[2,1]<--1 m[2,3]<-1 m[3,1]<--1 m[3,4]<-1 m[4,1]<--1 m[4,5]<-1 m[5,1]<--1 m[5,6]<-1 m[6,1]<--1 m[6,7]<-1 m[7,1]<--1 m[7,8]<-1 m[8,1]<--1 m[8,9]<-1

Dear R community I currently try to get post hoc multiple comparisons with the package multcomp from a cox mixed-effects model, where the survival is explained by two variables (cover with levels: nocover and cover; treatment with levels: tx, uv, meta), whose interaction is significant. I read Hothorn, T. 2011: Additional multcomp Examples and there is an example involving a two-way ANOVA with

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

Hi, I have data with binomial response variable (survival) and 2 categorical independent variables (site and treatment) (see below).? I have run a binomial GLM and found that both IVs and the interaction are significant.? Now I want to do a post-hoc test for all pairwise comparisons to see which treatment groups differ.? I tried the glht function in the multcomp package, but I get surprising

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

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