OR1Expt.R

Template:Scheme {{#security: Sven,Goldfinger}}

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)

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