similar to: Multiple comparisons on Anova.mlm object

Hi, I wish to compute multivariate test statistics for a within-subjects repeated measures design with anova.mlm. This works great if I only have two factors, but I don't know how to compute interactions with more than two factors. I suspect, I have to create a new "grouping" factor and then test with this factor to get these interactions (as it is hinted in R News 2007/2), but

Hi, Am newish to stats and R, so I certainly appreciate any help. Basically I have 50 inidividuals whom I have 6 photos each of their optic nerve head. I want to check that the orientation of the nerve head is consistent, ie the 6 replicates show minimal or preferably no rotation differences. I'll draw an arbitrary line between some blood vessels (same reference in each set of replicates) and

I have design with two repeated-measures factor, and no grouping factor. I can analyze the dataset successfully in other software, including my legacy DOS version BMDP, and R's 'aov' function. I would like to use 'Anova' in 'car' in order to obtain the sphericity tests and the H-F corrected p-values. I do not believe the data are truly deficient in rank. I

Hi, sorry if this is explained somewhere but I didn't find anything. How can I use "Anova" from the car package to test a modell without between subject's factors? Suppose I have the following data mat.1 mat.2 mat.3 di ex 1 85 85 88 1 1 2 90 92 93 1 1 3 97 97 94 1 1 4 80 82 83 1 1 5 91 92 91 1 1 6 83 83

For some time I have been looking for a convenient way of performing post-hoc analysis to Repeated Measures ANOVA, that would be acceptable if sphericity is violated (i.e. leaving aside post-hoc to lme models). The best solution I found was John Fox's proposal to similar requests in R-help: http://tolstoy.newcastle.edu.au/R/e2/help/07/09/26518.html

[hope this doesn't get posted twice, my first mail bounced] Hi, again, this is about the call bug in anova.mlm. I originally reported it as PR#7898 and I suggested a fix at PR#7904. (None of these message were forwarded to R-Devel, instead I received a bounce telling me that my provider's SPF settings are incorrect. They are, though, so there seems to be a problem with the

Hello all, I am getting the following error: Error in linear.hypothesis.mlm(mod, hyp.matrix.1, SSPE = SSPE, V = V, : The error SSP matrix is apparently of deficient rank = 7 < 11 After running: mod.ok <- lm(as.matrix(dat[,-1]) ~ DC, data=dat) (av.ok <- Anova(mod.ok, idata=idata, idesign=~week)) Although if I jitter the data in "dat", the function seems to work. What

Env: R 2.8.1, Win Xp, Eclipse/StatET In a .Rd file, I have an example containing the lines: # calculate Y M, using polynomial contrasts trends <- as.matrix(VocabGrowth) %*% poly(8:11, degree=3) colnames(trends)<- c("Linear", "Quad", "Cubic") [At the risk of a long message, I'll append the complete .Rd file at the end of this message, in case this was

Dear R folk, First off I want to thank those of you who responded with comments for my R quick and dirty stats tutorial. They've been quite helpful, and I'm in the process of revising them. When it comes to repeated measures MANOVA, I'm in a bit of a bind, however. I'm beginning to see that all of the documentation is written for psychologists, who have a slightly

Hello everyone, I've conducted a Type III repeated-measures ANOVA using Anova() from the car package, based on the suggestions at http://blog.gribblelab.org/2009/03/09/repeated-measures-anova-using-r/(option 3) and http://languagescience.umd.edu/wiki/EEG#ERP_ANOVA_in_R. My ANOVA has two factors: Condition (3 levels) and Region (6 levels) and their interaction. Below is code to run the Anova

Dear R users, I use the excelent Anova function of the library car because the easy way to get sphericity correction. Unless I use the scan function. I have not been able to access the values ​​of sum squares and degrees of freedom for each effect in the univariate summary table. Example of the car library for Anova function: library(car) phase <- factor(rep(c("pretest",

Hi, I have encountered this error when attempting a One-way Repeated-measure ANOVA with my data. I have read the "Anova in car: SSPE apparently deficient rank" thread by I'm not sure the within-subject interaction has more degrees of freedom than subjects in my case. I have prepared the following testing script: rm(list = ls())

To whom it may concern I made some analysis with R using the command Anova. However, I found some problmes with the output obtained by selecting type II o type III sum of squares. Briefly, I have to do a 2x3 mixed model anova, wherein the first factor is a between factor and the second factor is a within factor. I use the command Anova in the list below, because I want to obtain also the sum

Dear list, I would like to run an ANCOVA using car::Anova with repeated measures factors, but I can't figure out how to do it. My (between-subjects) covariate always interacts with my within-subject factors. As far as I understand ANCOVA, covariates usually do not interact with the effects of interest but are simply additive (or am I wrong here?). More specifically, I can add a covariate as

When I fit a multivariate linear model, and the formula is defined outside the call to lm(), the method summary.mlm() fails. This works well: > y <- matrix(rnorm(20),nrow=10) > x <- matrix(rnorm(10)) > mod1 <- lm(y~x) > summary(mod1) ... But this does not: > f <- y~x > mod2 <- lm(f) > summary(mod2) Error en object$call$formula[[2L]] <- object$terms[[2L]]

Dear list member, My question is related to input file format to an Anova from car package. Here is an example of what I did: My file format is like this (and I dislike the idea that I will need to recode it): Hormone day Block Treatment Plant Diameter High N.Leaves SH 23 1 1 1 3.19 25.3 2 SH 23 1 1 2 3.42 5.5 1 SH 23 1 2 1 2.19 5.2 2 SH 23 1 2 2 2.17 7.6 2 CH 23 1 1 1 3.64 6.5 2 CH 23 1 1 2

Hello dear R members. I have been learning the Anova syntax in order to perform an SS type III Anova with repeated measures designs (thank you Prof. John Fox!) And another question came up: where/what are the (between/within) residuals for my model? ############ Play code: phase <- factor(rep(c("pretest", "posttest", "followup"), c(5, 5, 5)),

Hello, I've just fitted a model with multi-responses, and I get an object of class "lm" "mlm". My problem is that as soon as I invoke the predict method for a dataframe "newdata", the methods runs and give me back prediction at the fitting points but not for newdata. Does someone has an explanation for this behavior, and some ideas to make predict.mlm work

It seems that confint.default returns an empty data.frame for objects of class mlm. For example: ``` nobs <- 20 set.seed(1234) # some fake data datf <- data.frame(x1=rnorm(nobs),x2=runif(nobs),y1=rnorm(nobs),y2=rnorm(nobs)) fitm <- lm(cbind(y1,y2) ~ x1 + x2,data=datf) confint(fitm) # returns: 2.5 % 97.5 % ``` I have seen proposed workarounds on stackoverflow and elsewhere, but

Full_Name: Yves Rosseel Version: 2.2.1 OS: i686-pc-linux-gnu Submission from: (NULL) (157.193.116.152) Dear developers, For the single-model case, the anova.mlm() function does not seem to handle multi-parameter predictors (eg factors) correctly. A toy example illustrates the problem: Y <- cbind(rnorm(100),rnorm(100),rnorm(100)) A <- factor(rep(c(1,2,3,4), each=25)) fit <- lm(Y ~ A)