A correction for estimating error when using the Local Pooled Error
statistical test
Carl Murie and Robert Nadon
You can access the ajdustedLPE.r script below.
adjustedLPE.r
A data set taken from the Affymetrix Latin Square HGU133A experiment can be used to test the script.
HGU133A_RMA.txt
Documentation is located within the script but a simple execution of the code can be done with the following lines of code:
source("adjustedLPE.r")
data <- read.table("HGU133A_RMA.txt")
lpeResults <- runLPE(data)