Dear all,
I have a question regarding control probes in the mogene10sttranscriptcluster array. I wish to remove the control probes (not duplicated) from a dataset downloaded from GEO (getGEO) which is based on the Affymetrix Mouse Gene 1.0 ST Array (MoGene-1_0-st).
Any suggestions with a few lines of code would be appreciated.
Regards,
Francois
1 answer
Hey,
You wish to remove them from the normalised dataset within R, I presume? If the data has been processed by the affy package, then this should work:
#Get the control probe indices
ControlProbes <- grep("^AFFX", featureNames(NormalisedData))
ControlProbes
#Remove the control probes
NormalisedData <- NormalisedData[-ControlProbes,]
#Fit the linear model
fit <- lmFit(NormalisedData, design=DesignMatrix)
You should download the library information for your chip type from Affymetrix / Thermo-Fisher, though, which would give a more detailed description of control probes and their IDs on your chip type. Usually, they begin with AFFX, though.
Let me know if this works.
Log in to answer this question.