Hi, I am wondering if there is a way to compute a FIRMA score for each transcript. Currently I only have FIRMA score for each probeset or group. I did as follows:
1. plmTr <- ExonRmaPlm(csN, mergeGroups=TRUE) 2. fit(plmTr) 3. firma<-FirmaModel(plmTr) 4.fit(firma) 5.fs<-getFirmaScores(firma) Or with the FIRMA score of each probeset, find out which transcripts are differentially expressed. Now I am trying to use limma to solve my problems and want some directions if anyone knows. Thanks, Libing --~--~---------~--~----~------------~-------~--~----~ 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 aroma-affymetrix@googlegroups.com To unsubscribe from this group, send email to aroma-affymetrix-unsubscr...@googlegroups.com For more options, visit this group at http://groups.google.com/group/aroma-affymetrix?hl=en -~----------~----~----~----~------~----~------~--~---