search for: amod

...39; 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 'glht' for the above amod. I know that I cannot use simply glht(amod, linfct = mcp(f = "Dunnett")). Also, if I want to use Dunnett for comparing C vs (T1 and T2), how can I specify this in the glht function. Thanks in advance. Regar...

...model nor term. Is it possible to build together a model and term component myself, so that glht will work for repeated measurements ANOVAs? If so, how would I do that? Best regards, Michael Zehetleitner For example (senseless, but it makes the point): --------- snip ------------- > amod <- aov(breaks ~ tension, data = warpbreaks) > amod.err <- aov(breaks ~ tension + Error(wool/tension), data = warpbreaks) > > class(amod) [1] "aov" "lm" > class(amod.err[["wool:tension"]]) [1] "aov" "lm" > > names(amod...

I am not sure whether there is a bug. When I tested the example given for "glht" in the help, I entered the following error: Running commands: amod <- aov(minutes ~ blanket, data = recovery) rht <- glht(amod, linfct = mcp(blanket = "Dunnett"), alternative = "less") Errors are: Error in try(coef.(model)) : could not find function "coef." Error in modelparm.default(model, ...) : no 'coe...

...he point from the output. The three cases only differ in the index I gave for "Ksub". ### Script start#### library(mvtnorm) library(multcomp) data("cholesterol") summary(cholesterol) ### adjusted p-values for all-pairwise comparisons in a one-way layout ### set up ANOVA model amod <- aov(response ~ trt, data = cholesterol) ### set up multiple comparisons object for all-pair comparisons cht <- glht(amod, linfct = mcp(trt = "Tukey")) summary(cht, test = adjusted("Westfall")) ### use only a subset of all pairwise hypotheses K <- contrMat(table(ch...

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

Hi, I want to export the result of glht in R into a LaTeX table, such as that result: Linear Hypotheses: Estimate Std. Error z value Pr(>|z|) Group1 - Group2 == 0 -0.14007 0.01589 -8.813 <0.001 "***" Group1 - Group3 == 0 -0.09396 0.01575 -5.965 <0.001 *** --- Signif. codes: 0 ?***? 0.001 ?**? 0.01 ?*? 0.05

I am not sure whether there is a bug. When I tested the example given for "glht" in the help, I entered the following error: Running commands: amod <- aov(minutes ~ blanket, data = recovery) rht <- glht(amod, linfct = mcp(blanket = "Dunnett"), alternative = "less") Errors are: Error in try(coef.(model)) : could not find function "coef." Error in modelparm.default(model, ...) : no 'coef...

Dear *, I would like to obtain for each factor of my anova model the "response variable vs factor" plot with means and 95% Tukey HSD intervals. I would appreciate any information on how to do that. Cheers -------------------------------------------------------------------- Max MANFRIN Tel.: +32 (0)2 650 3168 IRIDIA - CoDE, CP 194/6

..."Coeff") ) ssB <- selfStart( model = function(x, Coeff, A) { paste(x, Coeff, A) }, initial = function(mCall, data, LHS) { x <- eval(mCall[["x"]], data, parent.frame()) A <- eval(mCall[["A"]], data, parent.frame()) Amod <- paste(A, "mod in B") getInitial(y ~ ssA(x, Coeff, Amod), data) }, parameters = c("Coeff") ) getInitial(y ~ ssA("this", "that", "other"), data= data.frame(x=c("test"))) getInitial(y ~ ssB("this", &...

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

hey, how can i extract the values from the CI's when i use following code for a tukey test? the output shows three CI's and i know it should work with 'names but i don't really know how... thanks library(multcomp) data1$soort<-as.factor(data1$soort) amod<-aov(waarde~soort,data=data1) g<-glht(amod, linfct=mcp(soort = "Tukey")) confint(g) -- View this message in context: http://www.nabble.com/Tukey-in-R%2C-extracting-values-tp16609570p16609570.html Sent from the R help mailing list archive at Nabble.com.

I was running 1.9.18p10 on a test DEC Unix 4.0d box here and I decided to give 2.0beta3 a try. For some reason, I cannot connect to any of the server's shares or even get a list of it's services (via Start->Run->\\server) from Win95 clients. It asks for a password and no matter what I put in there I get "The password is not correct. Try again." I configured Beta3 a la:

...hese times are normally kept in GMT */ time_t mtime; time_t atime; @@ -125,11 +125,11 @@ int blocksize=20; int tarhandle; -static void writetarheader(int f, char *aname, int size, time_t mtime, +static void writetarheader(int f, char *aname, SMB_BIG_UINT size, time_t mtime, char *amode, unsigned char ftype); static void do_atar(char *rname,char *lname,file_info *finfo1); static void do_tar(file_info *finfo); -static void oct_it(long value, int ndgs, char *p); +static void oct_it(SMB_BIG_UINT value, int ndgs, char *p); static void fixtarname(char *tptr, char *fp, int l); stat...

Hello, I am trying to determine the most appropriate way to run post-hoc comparisons on my lme model. I had originally planned to use Tukey HSD method as I am interested in all possible comparisons between my treatment levels. TukeyHSD, however, does not work with lme. The only other code that I was able to find, and which also seems to be widely used, is glht specified with Tukey:

...p a Driver Domain for a NIC. I use ubuntu Jaunty Dom0 and the same kernel for the Driver Domain as well. I have successfully pass-through-ed the PCI device to the DomU ( Driver Domain). But i am not sure as to what i will have to do next to get the IDD functional ? -- -- With Regards , Amod Jaltade Good, Fast, Cheap: Pick any two (you can''t have all three). -- RFC 1925 _______________________________________________ Xen-users mailing list Xen-users@lists.xensource.com http://lists.xensource.com/xen-users

Hi All, I have performed a 3-way ANOVA analysis for my experimental data using aov function. My simple R funtion for this is: 3aof <- function(x){ m <- data.frame(R,S,T, x); anova(aov(x ~ R+S+T+R*S+R*T+S*T+R*S*T, m) ) } Now, I am getting P values for all the main and interactions effects. If I want to perform postdoc test on one of my main effects, say T, what method I should use (i have

Hi list, It appears there is a problem with the OpenMPI I/O library on CentOS 6.2 & 6.3 (package openmpi-1.5.4-1.el6.x86_64). When I compile the attached program it ends up in the error path since MPI_File_open returns 16. The corresponding (unhelpful) message is: MPI_ERR_OTHER: known error not in list I couldn't find any pointers on the net and the same program works with OpenMPI

Hi everyone, I am fairly new to R, and I am aware that others have had this problem before, but I have failed to solve the problem from previous replies I found in the archives. As this is such a standard procedure in psychological science, there must be an elegant solution to this...I think. I would much appreciate a solution that even I could understand... ;-) Now, I want to calculate a

...lp at stat.math.ethz.ch> Message-ID: <929F1767385B5847B3D584E35576A1B90355EEEB at MBX03.stf.nus.edu.sg> Content-Type: text/plain I am not sure whether there is a bug. When I tested the example given for "glht" in the help, I entered the following error: Running commands: amod <- aov(minutes ~ blanket, data = recovery) rht <- glht(amod, linfct = mcp(blanket = "Dunnett"), alternative = "less") Errors are: Error in try(coef.(model)) : could not find function "coef." Error in modelparm.default(model, ...) : no 'coef...

On 22.08.2018 17:52, Ryan Taylor wrote: > This is probably going to effect on other backends and break llvm-lit > for them? Very likely, yes. Can you take a look at how big the fallout is? This might give us a hint about what other frontends might expect, and who needs to be involved in the discussion (if one is needed). Cheers, Nicolai > > On Wed, Aug 22, 2018 at 11:41 AM