OR1Expt.R

Our Perl scripts. {{#security: Sven,Goldfinger}}

1. ------------------------- OR1 dataset analysis ------------------------------ #

sourceDir <- "/Volumes/HD2/Clinton/Data" filename <- dir(sourceDir, pattern="OR1.+\.txt")

dset <- read.table(file.path(sourceDir, filename), sep="\t", header=TRUE)

dim(dset) summary(dset)

1. Convert names

colnames(dset)[4:5] <- c("AM","RE")

library(lattice)

1. Formulas
2. y ~ x
3. y ~ x | z
4. y ~ x | (z * w)

1. help(xyplot), help(panel.superpose)

xyplot( RE ~ AM, data=dset) xyplot( RE ~ Timepoint| AM * Treatment, data=dset, panel.superpose=dset$Amount, groups=Amount)

1. Anova example
2. help(aov), help(summary)

aovobj <- aov(RE ~ AM*Treatment, data=dset) summary(aovobj)