cdf <- AffymetrixCdfFile$byChipType("HuEx-1_0-st-v2") print(cdf) cdf <- AffymetrixCdfFile$byChipType("HG-U133_Plus_2") cs <- AffymetrixCelSet$byName("celfiles", cdf=cdf) # rawdata in cwd. HuEx-1_0-st-v2 in celfiles # this should match cdf file names too print(cs) bc <- RmaBackgroundCorrection(cs) #backgroun csBC <- process(bc,verbose=verbose) qn <- QuantileNormalization(csBC, typesToUpdate="pm") #quantile norm csN <- process(qn, verbose=verbose) getCdf(csN) ##get CDF infor plmTr <- ExonRmaPlm(csN, mergeGroups=TRUE) #entrire transcript #for exon-by-exon mergeGroup make FALSE print(plmTr)