similar to: Getting Sphericity Tests for Within Subject Repeated Measure Anova (using "car" package) (Adjusted Dataset)

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

Hello everyone, I have tried to look for this everywhere and so far have no luck. I have a plotMeans(DV, IV1, IV2, error.bars="se") graph that plots my data (DV-continuous, IVs are factors, IV1 - two levels, IV2-four levels). I am trying to increase a scale of my y-axis (to be consistent with my other graphs), but unfortunately nothing works with "plotMeans" function, which

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,

Hello - I am attempting to add interaction terms to a model using the gamm function in the mgcv package. I stripped it down to only two terms to simplify it. I have no trouble with individual terms, continuous or factors, nor with interacting factor terms; I cannot, however, figure out how to add interactions between continuous*factor terms. I applied a smoother to my continuous term because it

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

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

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

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

In the following code I have created the posterior density for a Bayesian survival model with four parameters. However, when I try to use the adapt function to perform integration in four dimensions (on my old version of R I get an error message saying that I have applied a non-function, although the function does work when I type kernel2(param0, theta0), or on the newer version of R the computer

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

Data set up as one observation/subject looks like (with a total of 10 subjects) Two treatments: shoe type with 3 categories and region with 8 categories ==> 24 "treatment" columns Subject PHallux PMidToes PLatToe PMTH1 PMidMTH PLatMTH PMidfoot PRearfoot LHallux LMidToes LLatToe LMTH1 LMidMTH LLatMTH LMidfoot LRearfoot DHallux DMidToes DLatToe DMTH1 DMidMTH DLatMTH

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 all, is there anyone who can help me extracting p-values from an Anova object from the car library? I can't seem to locate the p-values using str(result) or str(summary(result)) in the example below > A <- factor( rep(1:2,each=3) ) > B <- factor( rep(1:3,times=2) ) > idata <- data.frame(A,B) > fit <- lm( cbind(a1_b1,a1_b2,a1_b3,a2_b1,a2_b2,a2_b3) ? sex,

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

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

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

Olá pessoal, estou realizando uma ANOVA com medidas repetidas e estou utilizando a função "Anova" do pacote "car". Medi o biovolume de algas a cada dois dias durante 10 dias (no banco de dados abaixo só coloquei até o 4° dia). Tenho 2 tratamentos ("c","t") e o experimento foi realizado em tréplicas ("A","B","C"). > Pa2

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