search for: aov

Displaying 20 results from an estimated 334 matches for "aov".

2007 Jan 15
model simplification in repeated anova
Hi, I''ve done a repeated ANOVA looking as follows: model<-aov(y~a*b*c+Error(d)) Now I want to get rid of some or all of the interactions. But neither AIC nor anova work. What is the appropriate function to test my model?
2008 May 28
Tukey HSD (or other post hoc tests) following repeated measures ANOVA
...that even I could understand... ;-) Now, I want to calculate a post-hoc test following up a within-subjects ANOVA. The dv is reaction time (RT), there is a 3-level Condition factor (Cond; within-subjects), a number of subjects (Subj), and the dataframe is called WMU3C. The model is > RT.aov <- aov(RT~Cond + Error(Subj/Cond), WMU3C) I understand that TukeyHSD only works with an aov object, but that RT.aov is an aovlist object. > class(RT.aov) [1] "aovlist" "listof" I''ve tried to work around it using the "maiz" example in the MMC do...
2003 Oct 15
aov and non-categorical variables
It is unclear to me how aov() handles non-categorical variables. I mean it works and produces results that I would expect, but I was under impression that ANOVA is only defined for categorical variables. In addition, help(aov) says that it "call to ''lm'' for each stratum", which I presume means...
2002 Jul 05
balance in AoV (was aov() and NaN)
<ripley at> wrote: > Hint 2: in the absence of balance, ...., and lme can do that Would it be possible to make aov like wrappers to the various special lm variants allowing for a uniform syntax for anova? aov(resp~f1*f2+Error(S/(f1*f2))) ## uses lm aov.lme(resp~f1*f2+Error(S/(f1*f2))) ## uses lme aov.rlm(resp~f1*f2+Error(S/(f1*f2))) ## uses rlm ... I''d do it, but I do not nearly understa...
2002 Sep 06
Using Anova Sums of Squares
Dear all, I''m trying to access the Sums of Squares resulting from a summary(aov(....)) so I can use them in a function. Is there an easy way to access these? The structure of a summary.aov object is something like this: > str(summary(fhc.aov)) List of 5 $ Error: PPNR :List of 1 ..$ :Classes anova and `data.frame'': 1 obs. of 5 variables: .. ..$ Df...
2013 Mar 06
aov() and anova() making faulty F-tests
Dear useRs, I''ve just encountered a serious problem involving the F-test being carried out in aov() and anova(). In the provided example, aov() is not making the correct F-test for an hypothesis involving the expected mean square (EMS) of a factor divided by the EMS of another factor (i.e., instead of the error EMS). Here is the example: Expected Mean Square df...
2002 Jul 05
aov() and NaN
hi, another anova question !! how does aov() treat missing values ? since it might depend on my experimental design: i have a within subject design with two factors ... i use aov() like this: aov(performance~type*block+Error(subj/(type+block)),data=anovaData) now it might happen that i have a couple of NaN data points in my result...
2006 Aug 03
between-within anova: aov and lme
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 ( does 6.8.5 Example 5: Stevens pp. 468 - 474 (one between, two within) between: gp within: drug, dose aov(effect ~ gp * drug * dose + Error(...
2010 Sep 29
Problems extracting p-value from summary(aov(...))
Hi, I did a aov and used summary to obtain the p-value. I tried many ways to extract the p-value from the summary result but failed. Among others I tried the following: > test.summary <- > summary(aov(data[,1]~time.points+Error(subject/time.points))) > test.summary Error: subject Df...
2011 Dec 11
multiple comparison of interaction of ANCOVA
...5.55 2010 C"), header = TRUE) > closeAllConnections() The experiments is designed without replication, thus I have to use ANCOVA or linear mixed effect model to analyze the data. In the model, variable year is coded as a continuous variable, and Trt is factor variable. > mydata.aov <- aov(y~Trt*year, mydata) > anova(mydata.aov) Analysis of Variance Table Response: y Df Sum Sq Mean Sq F value Pr(>F) Trt 2 140.106 70.053 197.9581 3.639e-08 *** year 1 0.610 0.610 1.7246 0.221600 Trt:year 2 8.804 4.402 12.4387 0.002567 ** Re...
2007 Oct 07
Question about aov
Hello R gurus, I am a beginner with R. I am doing an ANCOVA analysis using ''aov,'' and need some help understanding how ''aov'' works. I have a dataset (taken from looking at hypnotic induction. The variable ''X'' is a measure of how susceptible the subject is to being hypno...
2012 Oct 08
aov() usage
Hi R-listers, I am wondering if the function aov() in plyr is appropriate for two different types of tests: 1) > summary(aov(EDI ~ VegIndex, AND 2) > summary(aov(HSuccess ~ VegIndex + Aeventexhumed + VegIndex:Aeventexhumed, the later inclusive of an interaction of the two explanatory variabl...
2002 Dec 09
Problem with differences between S+ and R in parsing output tables with $
...t doesn''t seem to retrieve the appropriate number from a table. Note that some of the cosmetic differences in the two tables have to be dealt with by the parser. Any ideas what''s happening? REX # Begin R ------------------------------------------------------------- # > summary(aov(model)) Df Sum Sq Mean Sq F value Pr(>F) fa 1 11.1 11.1 0.9327 0.3345 Residuals 698 8279.3 11.9 # R: ... trying to parse the table gives a NULL for the probability of F... > (summary(aov(model))$"Pr(>F)")[1] NULL # End R -----...
2008 Dec 16
model.tables error from aov
...Mean :13.40 18 : 3 3rd Qu.:16.00 2 : 3 Max. :59.00 (Other):40 Note that the TrainErrs column is the only numeric column, and I forced everything else to be a factor. (Is that correct?) I then do the following: <- aov(TrainErrs ~ (idio*Type) + Error(Subject/Type) + (idio), So, idio is between-subjects and Type is within-subjects. This is based on examples I''ve found elsewhere. summary( This seems to work fine: Error: Subject Df Sum Sq Mean Sq...
2008 Mar 10
question for aov and kruskal
Hi R users! I have the following problem: how appropriate is my aov model under the violation of anova assumptions? Example: a<-c(1,1,1,1,1,1,1,1,1,1,2,2,2,3,3,3,3,3,3,3) b<-c(101,1010,200,300,400, 202, 121, 234, 55,555,66,76,88,34,239, 30, 40, 50,50,60) z<-data.frame(a, b) fligner.test(z$b, factor(z$a)) aov(z$b~factor(z$a))->ll TukeyHSD(ll) Now from...
2006 Nov 03
ANOVA in Randomized-complete blocks design
...quite confident how to represent this model in R, see below) The ANOVA table from S&R looks like this: MSB series 3 0.021 0.007 10.23 ** MSA genotypes 2 0.010 0.005 6.97 * MSE error 6 0.004 0.001 In R, I am using the following model (is this correct?) aov.ob = aov( genotype ~ genotype + series ) However, my results are Df Sum Sq Mean Sq F value Pr(>F) series 3 0.0135867 0.0045289 6.6360 0.02469 * genotype 2 0.0056732 0.0028366 4.1563 0.07367 . Residuals 6 0.0040948 0.0006825 What am I doing wrong? Regards, Jan...
2008 Mar 15
Hi all, I apologize for what might be a silly question. I am interested in doing a one way anova. This is not too hard in and of itself, either with anova, aov or oneway.test . However, I need to 1) get pvalues, 2) do a posthoc analysis with Tukey HSD, 3) and have (sometimes) an unbalanced design. I just can''t seem to put all the pieces together. Any suggestions? Thanks in advance, Dan. -- Daniel C. Jupiter, Ph.D. Postdoctoral Research Ass...
2004 Jan 16
error message in plot(aov-object) -- repost
Hi all, I posted this question several days ago, but did not get any answer until now. Since I still have no clue about the source of this error message, I repost a description of the problem including some code: A student at our institute fitted an aov model, and got the following error message: > plot(p.aov) Hit <Return> to see next plot: Hit <Return> to see next plot: Error in plot.window(xlim, ylim, log, asp, ...) : need finite ylim values The second plot was not produced. I then tried to produced the plot manuall...
2000 Jul 05
Tukey.aov with split-plot designs
I am using R 1.1 with Redhat 6.2 and RW 1.001 with Win98 (the upkey doesn't work on my IBM either as has been previously reported by others). The function aov doesn't return either the residuals or the residual degrees of freedom for split-plot designs. If you use the following code from Baron and Li's "Notes on the use of R for psycology experiments and questionnaires" data <- c(30,14,24,38,26,28,18,20,34,28,16,10,18,20,14,34,22,...
2008 Sep 13
moving from aov() to lmer()
Hello, I''ve used this command to analyse changes in brain volume: mod1<-aov(Volume~Sex*Lobe*Tissue+Error(Subject/(Lobe*Tissue)),data.vslt) I''m comparing males/females. For every subject I have 8 volume measurements (4 different brain lobes and 2 different tissues (grey/white matter)). As aov() provides only type I anovas, I would like to use lmer() w...