similar to: How to display y-axis labels in Multcomp plot

Hi there, I wonder if anyone can help me. I'm trying to use bargraph.CI in the Sciplot package when there is a missing combination of the factor levels. Unfortunately the standard errors on the plot do not appear to be correct. Consider an analysis consisting of two factors A and B. When all factor level combinations are present all appears fine: library(sciplot) #all data

fixed @ R-forge. New version should appear on CRAN soon. Thanks for the report! Torsten > > ---------- Forwarded message ---------- > Date: Sat, 3 Sep 2011 23:56:35 +0200 > From: Ulrich Halekoh <Ulrich.Halekoh at agrsci.dk> > To: "r-help at r-project.org" <r-help at r-project.org> > Subject: [R] glht (multcomp): NA's for confidence intervals using

Hej, Calculation of confidence intervals for means based on a model fitted with lmer using the package multcomp - yields results for calpha=adjusted_calpha - NA's for calpha=univariate_calpha Example: library(lme4) library(multcomp) ### Generate data set.seed(8) d<-expand.grid(treat=1:2,block=1:3) e<-rnorm(3) names(e)<-1:3 d$y<-rnorm(nrow(d)) + e[d$block]

Hello: If anyone could guide me with this I would greatly appreciate it. Thanking you in advance for your assistance. Using a 3-level input factor alternative so that a function(below) can compute both a two-sided and one-sided p-values. Making the two-sided test the default. And produce output information about which alternative was tested. Where would I place the ifelse statement?

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

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

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

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

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

It's been a few weeks I'm racking my brains on how to use the function glht the package multcomp to test interactions. Unfortunately, the creator of the package forgot to put a sample in pdf package how to do it. I have looked in several places, but found nothing. If someone for the love of God can help me I'll be extremely grateful. The model is glm. -- View this message in context:

Hi All, I am trying use 'multcomp' for multiple comparisons after my ANOVA analysis. I have used the following code to do ANOVA: dat <- matrix(rnorm(45), nrow=5, ncol=9) f <- gl(3,3,9, label=c("C", "Tl", "T2")) aof <- function(x) { m <- data.frame(f, x); aov(x ~ f, m) } amod <- apply(dat,1,aof) Now, how can I use

Hi all, I'd like to compare all levels of my interaction with each other. I read the pdf 'Additional multcomp Examples' but even though there is an example with an interaction it doesn't work for me when I want to compare within and between groups. Here is an example: #### d.fr<-data.frame(id=rep(1:16,3),treat1=rep(as.factor(LETTERS[1:3]),each=

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

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

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 list members, Can anyone please point to an example of how to use glht(multcomp) with lmer objects? I am trying: summary(glht(lmerObject, linfct = mcp(x = "Tukey"))) as I would for a glm object, but with no luck. Thank you, Andy. [[alternative HTML version deleted]]

i ran this model > model2<-glm(rojos~ageandsex+sector+season+sector:season,quasipoisson) > glht(model2,linfct=mcp(ageandsex="Tukey")) General Linear Hypotheses Multiple Comparisons of Means: Tukey Contrasts Linear Hypotheses: Estimate M - H == 0 0.2898 SUB - H == 0 -0.2261 SUB - M == 0 -0.5159 I tried to do the same changing factor season

I am trying to do Tukey HSD comparisons on a repeated measures expt. I found the following example on r-help and quoted approvingly elsewhere. It is broken. Can anyone please tell me how to get it to work? I am using R 2.4.1. > require(MASS) ## for oats data set > require(nlme) ## for lme() > require(multcomp) ## for multiple comparison stuff > Aov.mod <- aov(Y ~ N + V +