..... qn <- QuantileNormalization(csBC, typesToUpdate="pm") csN <- process(qn, verbose=verbose)
Then how to get the probeset data with labels(with probeset_id as row names, like a data frame)? And for FIRMA, firma <- FirmaModel(plmEx) fit(firma, verbose=verbose) fs <- getFirmaScores(firma) exFirma <- extractDataFrame(fs,addNames=TRUE) exFirma[,6:ncol(exFirma)]<-log2(exFirma[,6:ncol(exFirma)]) .... Suppose I have 100 samples which can be classified as 2 groups(50 samples each group). After this step, how should I know which exon are differentially spliced in the two groups? Thanks! --~--~---------~--~----~------------~-------~--~----~ When reporting problems on aroma.affymetrix, make sure 1) to run the latest version of the package, 2) to report the output of sessionInfo() and traceback(), and 3) to post a complete code example. You received this message because you are subscribed to the Google Groups "aroma.affymetrix" group. To post to this group, send email to [email protected] To unsubscribe from this group, send email to [EMAIL PROTECTED] For more options, visit this group at http://groups.google.com/group/aroma-affymetrix?hl=en -~----------~----~----~----~------~----~------~--~---
