[BioC] What does maanova contrast?

Loren Engrav engrav at u.washington.edu
Fri Jul 10 18:37:51 CEST 2009


Am using the maanova package 1.14.0

Model includes two breeds  (fixed breed factor) and 5 time points (random
time factor)

The three commands are

> BDeset_gcrma_NoCon25LRmadata.raw.WithAffyID081109 <-
read.madata("ABeset_gcrmaNoCon_25PsLRs081107.txt",
designfile="DesignV4LR081109.txt", header=TRUE, spotflag=FALSE, intensity=2,
probeid=1)

> BGanova.Br_T.mix081109 <-
fitmaanova(BDeset_gcrma_NoCon25LRmadata.raw.WithAffyID081109,
formula=~Breed+Time, random= ~Time, covariate = ~1, verbose=TRUE)

> ftest.Breed.mix_1000perm081109 <- matest
(BDeset_gcrma_NoCon25LRmadata.raw.WithAffyID081109, BGanova.Br_T.mix081109, term
= "Breed", n.perm = 1000, critical = .9, test.type = c("ftest"), shuffle.method
= c("sample"), MME.method = c("REML"), test.method = c(1,1), pval.pool = TRUE,
verbose = TRUE)

With two breeds and five time points there are 9 possible contrasts.  Does
the Ftest in maanova study all 9 breed contrasts or just the 5 breed
contrasts at the specific time points?

Thank you

-- 
Loren Engrav, MD
Univ Washington
Seattle



More information about the Bioconductor mailing list