similar to: AOV/LME

I am doing an analysis and would like to use lme() and the multcomp package to do multiple comparisons. My design is a within subjects design with three crossed fixed factors (every participant sees every combination of three fixed factors A,B,C). Of course, I can use aov() to analyze this with an error term (leaving out the obvious bits): y ~ A*B*C+Error(Subject/(A*B*C)) I'd also like

Dear R-Community! The example "oats" in MASS (2nd edition, 10.3, p.309) is calculated for aov and lme without interaction term and the results are the same. But I have problems to reproduce the example aov with interaction in MASS (10.2, p.301) with lme. Here the script: library(MASS) library(nlme) options(contrasts = c("contr.treatment", "contr.poly")) # aov: Y ~

(A repost of this request with a bit more detail) Hi, All. I'd like to calculate effect sizes for aov or lme and seem to have a bit of a problem. partial-eta squared would be my first choice, but I'm open to suggestions. I have a completely within design with 2 conditions (condition and palette). Here is the aov version: > fit.aov <- (aov(correct ~ cond * palette +

Hi I searched in the list and only found questions without answers e.g. http://finzi.psych.upenn.edu/R/Rhelp02a/archive/19955.html : Is there a way to get the same results with lme as with aov with Error()? Can anybody reproduce the following results with lme: id<-c(1,1,1,2,2,2,3,3,3,4,4,4,5,5,5,1,1,1,2,2,2,3,3,3,4,4,4,5,5,5,1,1,1,2,2,2,3,3,3,4,4,4,5,5,5)

I'm working with a nested model (mixed). I have four factors: Patients, Tissue, sex, and tissue_stage. Totally I have 10 patients, for each patient, there are 2 tissues (Cancer vs. Normal). I think Tissue and sex are fixed. Patient is nested in sex,Tissue is nested in patient, and tissue_stage is nested in Tissue. I tried aov and lme as the following, > aov(gene ~ tissue + gender +

First of all, I apologize for asking a question that has appeared recurrently in this mailing list. However, I have googled for it, have looked at the mailing list archives, and also looked at Pinheiro & Bates book (although not very thoroughly, I must confess), to no avail. Here is the question: I am trying to obtain with lme or lmer the same exact numerical results (p-values) that I obtain

I need to specify error terms properly in a mixed-effects anova model. I know you can add error terms in aov using Error and can specify random factors in lme but I am not sure how these get treated. When making the calculations for fixed and random factors, are the correct error terms used and how can you get aov or lme to use different error terms for fixed and random effects? I'm

Sorry for reposting, but even after extensive search I still did not find any answers. using: summary(aov(pointErrorAbs~noOfSegments*turnAngle+Error(subj/(noOfSegments+turnAngle)), data=anovaAllData )) with subj being a random factor and noOfSegments and turnAngle being fixed factors, I get the following results: ---------------------------------------------- Error: subj Df Sum

Dear R user, I am trying to reproduce the results in Montgomery D.C (2001, chap 13, example 13-1). Briefly, there are three suppliers, four batches nested within suppliers and three determinations of purity (response variable) on each batch. It is a two stage nested design, where suppliers are fixed and batches are random. y_ijk=mu+tau_i+beta_j(nested in tau_i)+epsilon_ijk Here are the

I have 2 questions on ANOVA with 1 between subjects factor and 2 within factors. 1. I am confused on how to do the analysis with aov because I have seen two examples on the web with different solutions. a) Jon Baron (http://www.psych.upenn.edu/~baron/rpsych/rpsych.html) does 6.8.5 Example 5: Stevens pp. 468 - 474 (one between, two within) between: gp within: drug, dose aov(effect ~ gp * drug *

I am analyzing data from an experiment with two factors: Carbon (+/-) and O3 (+/-), with 4 replicates of each treatment, and 4 harvests over a year. The treatments are assigned in a block design to individual Rings. I have approaches this as a repeated measures design. Fixed factors are Carbon, O3 and Harvest, with Ring assigned as a random variable. I have performed repeated measures analysis

Hi, we have a split-plot experiment in which we measured the yield of crop fields. The factors we studied were: B : 3 blocks I : 2 main plots for presence of Irrigation V : 2 plots for Varieties N : 3 levels of Nitrogen Each block contains two plots (irrigated or not) . Each plot is divided into two secondary parcels for the two varieties. Each of these parcels is divided into three subplots

Hi all, I have collected response time data from 178 participants ('sub') for each combination of 4 within-Ss factors ('con','int','tone','cue'). Additionally, I have recorded the gender of each participant, so this forms a between-Ss factor ('gender'). Normally this would be analyzed using aov:

I have used lme() on data from a between-within subjects experiment. The correct ANOVA table is known because this is a textbook example (Experimental Design by Roger Kirk Chapter 12: Split-Plot Factorial Design). The lme() F-values differ from the known results. Please help me understand why. d<-read.table("kirkspf2.dat",header=TRUE) for(j in 1:4) d[,j] <- factor(d[,j]) ### Make

I am trying to understand better an analysis mean RT in various conditions in a within subjects design with the overall mean RT / subject as one of the factors. LME seems to be the right way to do this. using something like m<- lme(rt~ a *b *subjectRT, random= ~1|subject) and then anova(m,type = "marginal"). My understanding is that lme is an easy interface for dummy coding

Dear R Users, I want to write a function that applies to the dataframe and variables that were used in a previous call to lm or aov. In order to do this, I need to write a function that applies to the output of lm or aov, and yields the names of the dataframe and variables that were used in the lm or aov analysis. For example, suppose that I give the command: aov.out <- aov( Rt ~

Dear listmembers, I would like to define a class with a slot that takes either an object of class aov or NULL. I have been reading "S4 Classes in 15 pages more or less" and "Lecture: S4 classes and methods" #First I tried with list and NULL setClass(listOrNULL") setIs("list", "listOrNULL") setIs("NULL", "listOrNULL") #doesn't

Hi, I have a question of aov. e.g. aov.ex = aov(x~y) summary(aov.ex) The aov summary will print to the screen. How can I extract the aov result, in particular the values of Pr(>F) and F value into a vector so that I can use them for other use? Thanks. -- Waverley @ Palo Alto

Hello, I have a simple question about using the aov function syntax (ie. * + or :) for the interaction of 2 factors. I have read the help files, and researched other sites, and have included my source files. My goal is to measure the signifigance of the interaction between population and condition (aka. population:condition). I can't seem to figure it out. 1. In the first example the

Dear Madame, Dear Sir, I am able to obtain the coefficients from a 'summary' of 'lm', but NOT from a 'summary' of 'aov'. The following example shows my steps. ## Initialize rm(list = ls()) # remove (almost) everything in the working environment utils::data(npk, package="MASS") # get data model <- yield ~ block + N*P*K ## Using lm npk.lm <-