similar to: anova.mlm again

Hi, I'm encountering a very odd problem with calls to anova.mlm() from within a function. Consider the following code (data.n is a matrix of numeric values): mlmfit <- lm(data.n ~ 1) mlmfit0 <- lm(data.n ~ 0) print(mlmfit) anova(mlmfit,mlmfit0,test="Spherical") If I run it just like this from the console, it works just fine. If, however, I call it from within a function,

Hi, yet another anova.mlm problem - it doesn't seem to end. This time, I have a setup with a few within-subject factors and a between-subject factor (SGROUP). Consider the most simple case with only one within-factor (apo): > mlmfit0 <- lm(data.n ~ 0 + SGROUP) > mlmfit1 <- lm(data.n ~ 1 + SGROUP) > anova(mlmfit1,mlmfit0,test="Spherical",M=~hemi,X=~1) Analysis of

I've rolled up R-2.1.1.tar.gz a short while ago. This is a maintenance release containing mainly bugfixes. See the full list of changes below. You can get it from http://cran.r-project.org/src/base/R-2/R-2.1.1.tar.gz (give it some time to arrive there) or wait for it to be mirrored at a CRAN site nearer to you. If you're *really* impatient, www.biostat.ku.dk/~pd/R-release should work

I've rolled up R-2.1.1.tar.gz a short while ago. This is a maintenance release containing mainly bugfixes. See the full list of changes below. You can get it from http://cran.r-project.org/src/base/R-2/R-2.1.1.tar.gz (give it some time to arrive there) or wait for it to be mirrored at a CRAN site nearer to you. If you're *really* impatient, www.biostat.ku.dk/~pd/R-release should work

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

I would like to perform multiple comparisons or post-hoc testing on the independent variable in an Anova.mlm object generated by the Anova function of the car package. I have defined a multivariate linear model and subsequently performed a repeated measures ANOVA as per the instructions in section #3 of the following comprehensive tutorial on the subject from the Gribble lab at UWO:

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

Dear R-devel'ers Below is an outline for a set of routines to improve support for multivariate linear models and "classical" repeated measurements analysis. Nothing has been coded yet, so everything is subject to change as loose ideas get confronted by the harsh realities of programming. Comments are welcome. They might even influence the implementation... -pd General

Hi, I'm using anova.mlm sphericity tests/corrections, and I'm getting different values than SAS. In order to be able to use these values for publications, I'd need to know more about the SAS bug that is mentioned in the Reference Manual. - What exactly causes the different values? - Is it just a slight difference, or can I expect significant differences in H-F/G-G epsilons 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, 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

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

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

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

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'm still trying to figure out that GLM procedure in SAS. Let's start with the simple example: PROC GLM; MODEL col1 col3 col5 col7 col9 col11 col13 col15 col17 col19 col21 col23 =/nouni; repeated roi 6, ord 2/nom mean; TITLE 'ABDERUS lat ACC 300-500'; That's the same setup that I had in my last email. I have three factors: facSubj,facCond and facRoi. I had this pretty

I forgot to attach the notes.txt with the objdump. On Thu, Jul 23, 2015 at 3:55 PM, Martell Malone <martellmalone at gmail.com> wrote: > Hi again rui, :) > > I've got all the patches into llvm and clang for supporting mingw-w64 via > compiler-rt and now we are able to build a full mingw-w64 toolchain without > gcc :) > With great help from yaron and rnk. > >

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

Dear list, I?m maintainsing a package containing only datasets (152): http://dutangc.free.fr/pub/RRepos/web/CASdatasets-index.html <http://dutangc.free.fr/pub/RRepos/web/CASdatasets-index.html> When R CMD checking the package, I get the following NOTE * checking data for non-ASCII characters ... NOTE Note: found 4 marked UTF-8 strings I wonder how to find which dataset(s) (all recorded