search for: multicomp

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

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", not the subje...

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

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

...y variance analysis and have now the following problem: I''d like to do some multiple Comparisons between different factor levels (two factors, each with some 5 or 6 levels), because there are some significant interactions between these two factors. In S-plus 5, there is the function "multicomp", but in R, there ist no matchable funktion, as far as I could find out. So I ask you: What would you do in similar situations? Thanks very much for your answer. Regards -- Ruth Meili <meili at stat.math.ethz.ch> Seminar fuer Statistik, LEO D6 ETH (Federal Inst. Technology) 8092 Zuri...

...hy if I try "lmer(sp~cla+(1|cla:plotti),..." or "lmer(sp~cla+(1|plotti:cla),...." I obtain the same results? 3) why the anova summary don't say if differences in classes are significance (or not significance)? 4) I'd like to perform a post-hoc test with the package "multicomp" but the lmer function give me a lmer object (and this kind of object is not read by the "multicomp" package). How could I perform my analysis in a different way? Thank you a lot for your help! Giovanni >I'd like to perform a glm analysis with a hierarchically nested desig...

[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 -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[...

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

Hello, Is there an option I can easily set to put labels next to the points in a simple plot? The closest I got to it was to use the option pch, but this replaces the point symbol with the character specified, and it is not what I want. Thanks -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read

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

...ise comparisons with : mc<- simint(ren~ID,type="Tukey",data=donnees) I cannot get a result in a reasonable time (after 2 hours the computer was still working and I interrupted the process). Can anybody tell me if there is some capacity limitation for simint (as well as for simtest)? Multicomp of Splus handles the problem adequately in a short time, but Splus is practically not accessible for student training... Kind regards, Patrick Giraudoux [[alternative HTML version deleted]]

...e 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 how do I work around it? Is there a more "fail-safe" general method of writing complex outputs to a table? The story here is I want to be able to call certain values from the TukeyHSD results to produce custom summaries. Using R version...

...riability is homogenous or not. Is it possible to do the same in R with glm? If not, what is the equivalent for testing these assumptions? Also, is it valid to use Tukey's post hocs to test main effects if necessary? If not, is there another method other than contrasts, something to do with multicomp? In a previous model where Site was random I used the following lme model model2<-lme(Anncov~Treat*Time, random=~1|Site, data=bob,control=list(opt="optim")) And checked residuals using: print(shapiro.test(model1$residuals)) print(qqnorm(model1)) #resids vs fitted plot(...

In the model: lm.1 <- lm(variable ~ BLOC + TIL * YEAR , data=selvanera) I found TIL*YEAR interaction significant. Then I am trying to compare means of the different levels of TIL inside every YEAR using: mc.2 <- glht(lm.1, linfct = mcp(TIL*YEAR="Tukey")) summary(mc.2, test = univariate()) but it does not work. There is any way of doing this, like the SLICE option in

...riability is homogenous or not. Is it possible to do the same in R with glm? If not, what is the equivalent for testing these assumptions? Also, is it valid to use Tukey's post hocs to test main effects if necessary? If not, is there another method other than contrasts, something to do with multicomp? In a previous model where Site was random I used the following lme model model2<-lme(Anncov~Treat*Time, random=~1|Site, data=bob,control=list(opt="optim")) And checked residuals using: print(shapiro.test(model1$residuals)) print(qqnorm(model1)) #resids vs fitted plot(...

..., are not the overall effect but are a comparison within a year (2002). I can still get my overall effect (using anova.lme) but how do I calculate the estimates (with P-Values if possible) for each seperate overall treatment comparison (not within year). I tried 'glht' (package 'multicomp') but this only works if there are no interactions present, otherwise again it gives a comparison for one particular year. Very grateful for any assistance, Mark

I swear I've scoured the help files and several texts before posting what feels like a dumb newbie question. How can I draw two kernel density plots in the same frame ? I have similar variables in two separate data frames, and I would like to show their two histograms/densities in a single picture. Same units, scale, range for both, so I'm simply trying to draw one and then add the

...out[[nm]] <- array(c(center, center - width, center + width), - c(length(width), 3), dnames) + out[[nm]] <- array(c(center, center - width, center + width,pvals), + c(length(width), 4), dnames) } class(out) <- c("multicomp", "TukeyHSD") attr(out, "orig.call") <- x$call diff -ru r-devel/R/src/library/stats/man/TukeyHSD.Rd r-local/R/src/library/stats/man/TukeyHSD.Rd --- r-devel/R/src/library/stats/man/TukeyHSD.Rd 2005-07-14 23:14:39.380653872 -0300 +++ r-local/R/src/library/stats/man/Tu...

Hi to all R lovers, I am thinking of investigating and perhaps improve the anova "aov" function in R. However before I do so I just like to see if anyone knows whether R is capable of multiple comparison tests (eg "Standard Difference of the Means, SED", "Fisher's Least Significant Difference" ...etc)? If it is, how can I do it? Thanks. Ko-Kang Wang

Dear List, are there post-hoc tests like Scheffe, LSD, etc. available after ANOVA test is performed with significant F-statistic? I have tried help.search("Scheffe"), but "No documentation found" (and I have most of packages installed). Probably there are such tests in R, and I am just searching badly... My second question is: Which test/method I should use for ANOVA-like