similar to: simple sphericity test using bartlett and cov?

I haven't been able to find an answer on this that's direct, only implied. In several places I have read that when people asked for sphericity tests they were guided toward lme or mlm models. But, there is no direct indication that the lme method is not subject to the sphericity assumption. In fact, it seems like it should be. Its just a linear model that handles random and

Hello and thanks in advance I am running a glm in R the code is as follows with residual diagnostic code below model4<-glm(Biomass~(Treatment+Time+Site)^2, data=bobB, family=quasi(link="log", variance="mu")) par(mfrow=c(2,2)) plot(model2) to test the effect of grazing exclusion of feral horses for a Phd with following factors: Treatment - 3 levels which are grazed

Hello Dear R user, I want to conduct a Principal components analysis and I need to run two tests to check whether I can do it or not. I found how to run the KMO test, however i cannot find an R fonction for the Bartlett's test of sphericity. Does somebody know if it exists? Thanks for your help! Thibault [[alternative HTML version deleted]]

Apologies for the obvious, but just to clarify: there is no reason to "justify" a PCA -- it's just an eigen decomposition of a matrix and is therefore "justified" by linear algebra. If one wants to determine whether some subset of the eigenvectors = principal components suffice to "represent" the data in some sense, then that is where distributional

I believe that the following function is a good function for eyeballing sphericity. But, I can't find functions for tests anywhere (g-g or h-f). Anyone who wants to can feel free to use it. It is bad that R does not provide any tests for sphericity AND cannot do a MANOVA on repeated measures (although lme is an alternative add on package). One of these solutions must be provided in the

I have 3 questions (below). Background: I am teaching an introductory statistics course in which we are covering (among other things) repeated measures anova. This time around teaching it, we are using R for all of our computations. We are starting by covering the univariate approach to repeated measures anova. Doing a basic repeated measures anova (univariate approach) using aov() seems

I'm an experimental psychologist and when I run ANOVA analysis in SPSS, I normally ask for a test of non-sphericity (Box's M-test). I also ask for output of the corrections for non-sphericity, such as Greenhouse-Geisser and Huhn-Feldt. These tests and correction factors are commonly used in the journals for experimental and other psychology reports. I have been switching from SPSS to R

I have been trying to use lavaan (version 0.4-7) for a simple path model, but the program seems to be computing far less degrees of freedom for my model then it should have. I have 7 variables, which should give (7)(8)/2 = 28 covariances, and hence 28 DF. The model seems to only think I have 13 DF. The code to reproduce the problem is below. Have I done something wrong, or is this something I

Hello everyone, I am trying to do within subjects repeated measures anova followed by the test of sphericity (sample dataset below). I am able to get either mixed model or linear model anova and TukeyHSD, but have no luck with Repeated-Measures Assuming Sphericity or Separate Sphericity Tests. I am trying to follow example from "car" package, but it seems that I am not getting something

[corrected dataset below] Hello everyone, I am trying to do within subjects repeated measures anova followed by the test of sphericity (sample dataset below). I am able to get either mixed model or linear model anova and TukeyHSD, but have no luck with Repeated-Measures Assuming Sphericity or Separate Sphericity Tests. I am trying to follow example from "car" package, but it seems

Hi, I have a dataset which has around 138 variables and 30,000 cases. I am trying to calculate a mahalanobis distance matrix for them and my procedure is like this: Suppose my data is stored in mymatrix > S<-cov(mymatrix) # this is fine > D<-sapply(1:nrow(mymatrix), function(i) mahalanobis(mymatrix, mymatrix[i,], S)) Error in solve.default(cov, ...) : system is computationally

Hi, Is there a way to do a sphericity test in R for repeated measures ANOVA (using aov or lme)? I can't find anything about it in the help. Thanks Orou [[alternative HTML version deleted]]

Hi R-users Does anybody know if there are functions or packages in R for sphericity tests like Bartlett's or Mauchly's? Thanks! __________________________________________________________ Eng. Agr., M.Sc. Eduardo Dutra de Armas __________________________________________________________ Centro de Energia Nuclear na Agricultura (CENA/USP) Laboratório de Ecotoxicologia Av.Centenário 303,

Dear Prof. Baron, thankyou for your pdf paper on using R, at http://www.psych.upenn.edu/~baron/rpsych.pdf I am currently interested in the sphericity notes on pp. 45-47. The equation on p. 46 suggests the following R calculation: D <- k^2 * ( mean(diag(S)) - mean(S) )^2 while the notes in the pdf indicate the following: D <- k^2 * ( mean(S) - mean(diag(S)) )^2 Can you please consider

Is there a function that calculates the power of a repeated-measure ANOVA design, e.g., 2 groups, 4 within-subject factors, an average 0.40 correlation between the 4 within factors, etc. I don't think I can use power.anova.test() because it does not consider corr=0.40. I am hoping that someone has already implemented the method by Muller & Barton (1989, JASA 84: 549-555) on how to

Hi, I'm getting to grips with using R to perform repeat measures analysis of variance and was wondering is there a potted method to obtain Huynh and Feldt or other adjustments for sphericity violations? Thanks Dominic Dr Dominic J Barraclough Dept Neurobiology and Anatomy, University of Rochester, 601 Elmwood Av, Box 603, Rochester NY 14642 tel (716) 2753627 dominic_barraclough at

Does anyone have any suggestions for how to test for violation of the sphericity assumption and apply the Huynh-Feldt correction to a repeated-measure model? My model follows the format below: Rep.anova=aov(Y~X1*X2+Error(Subject/(X1*X2),data=DataFile)) I've already looked at the explanations aov() and Error() in Baron's “Notes on the use of R for psychology experiments and

With respect to calculating the epsilon index of sphericity for ANOVA, discussed on pp. 45-47 of: http://www.psych.upenn.edu/~baron/rpsych.pdf It notes that epsilon is not required for a repeated measures design with only k=2 levels, as the minimum value of epsilon (e) is given by: e = 1/(k-1) so for k=2, we have e = 1 (ie, no correction of the F test df; see p. 46). These notes apply to a

Dear helpers, I have a 2x2 mixed design with two "groups" (between-subjects) and two presentation-types (within-subjects). The difference between groups is in the order of manipulations: group.CD having first a block of present.type.C and then a block of present.type.D, each block containing 31 trials. group.DC having first a block of present.type.D and then a block of present.type.C,

why does teh following code not work? a<-rnorm(50) dim(a)<-c(10,5) mycov<-function(mat,i,j){ sum((mat[,i]-mean(mat[,i]))*(mat[,j]-mean(mat[,j]))) } covmat<-function(X){ outer(1:dim(X)[2],1:dim(X)[2],function(v1,v2){mycov(X,v1,v2)})} from what i know covmat(a) should give the covariance matrix of a, but it gives a matrix with constant