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)