EBSeq/R/GetNg.R

   1 GetNg<- function(IsoformName, GeneName){
   2         GeneNg = tapply(IsoformName, GeneName, length)
   3         IsoformNg = GeneNg[GeneName]
   4         names(IsoformNg) = IsoformName
   5         GeneNgTrun=GeneNg
   6         GeneNgTrun[GeneNgTrun>3]=3
   7         IsoformNgTrun=IsoformNg
   8         IsoformNgTrun[IsoformNgTrun>3]=3
   9         out=list( GeneNg=GeneNg, GeneNgTrun=GeneNgTrun, IsoformNg=IsoformNg, IsoformNgTrun=IsoformNgTrun)
  10         }