# ------------------------------------------
# Example using the result of the analysis of
# the 6 arrays in the AffySpikeU95Subset data set
# Loading the data
data(AffySpikeU95Subset)
# Defining design and contrast matrix
group<-factor(rep(1:2,each=3))
design<-model.matrix(~group-1)
contrast<-matrix(c(1,-1),1,2)
# Analyzing using plw
model1<-plw(AffySpikeU95Subset,design=design,contrast=contrast,
epsilon=0.01)
## Look at fitted curve for scale parameter
scaleParameterPlot(model1)
Run the code above in your browser using DataLab