search for: amode

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

Dear all, glht (from the multcomp package) needs a term and a model component in it's fitted model. In fitted models from e.g. repeated measurements ANOVAs I do not find neither 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

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

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

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

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

Hi All, Another nls related problem (for background, I'm migrating a complicated modelling package from S-plus to R). Below I've reduced this to the minimum necessary to demonstrate my problem (I think); the real situation is more complicated. Two similar selfStart functions, ssA and ssB. The 'initial' function for ssB modifies its arguments a little and then calls getInital

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

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); stati...

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:

Hi, I am trying to setup 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

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

Dear friends, are there enough R users in Cyprus to form a club? jason Dr. Iasonas Lamprianou Department of Education The University of Manchester Oxford Road, Manchester M13 9PL, UK Tel. 0044 161 275 3485 iasonas.lamprianou at manchester.ac.uk ----- Original Message ---- From: "r-help-request at r-project.org" <r-help-request at r-project.org> To: r-help at r-project.org

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