Monday, March 4, 2013

Microarray data analysis

Transcripts whose probe sets with detection above background P-value ,0.05 in at least 2 out of 3
replicates at the 0 min time point were considered expressed
and used for subsequent analyses

### 1, extract probe level intensities
apt-cel-extract -o ./GSE21236/0000_pre/raw_probe_signal.txt ./GSE21236/RawData/*.CEL -c /Home/liang/data/Array/Affymetrix/Mouse/MoGene-1_0-st-v1.r4.clf -p /Home/liang/data/Array/Affymetrix/Mouse/MoGene-1_0-st-v1.r4.pgf  --b /Home/liang/data/Array/Affymetrix/Mouse/MoGene-1_0-st-v1.r4.bgp

### 2, calculate dabg p-value
apt-probeset-summarize -a dabg -c /Home/liang/data/Array/Affymetrix/Mouse/MoGene-1_0-st-v1.r4.clf -p /Home/liang/data/Array/Affymetrix/Mouse/MoGene-1_0-st-v1.r4.pgf -b /Home/liang/data/Array/Affymetrix/Mouse/MoGene-1_0-st-v1.r4.bgp -o ./GSE21236/0000_pre/ ./GSE21236/RawData/*.CEL

No comments:

Post a Comment