Hello, I was able to successfully run TSSi program. However, I have several questions. My questions are - how "0.CUFF.16.1_2" is related to "0.CUFF.20.1_4", - how region is determined to each gene, - how this file could be interpreted, - and how nTss is determined? Thank you again for sparing your precious time. Best regards. -- output of sessionInfo(): When I exported data by using the command bleow, x <- segmentizeCounts(counts=counts, start=start, chr=chromosome, region=region, strand=strand) yFit <- normalizeCounts( x, fit=TRUE ) z <- identifyStartSites( yFit ) segmentsRd <- segmentsAsRangedData( z ) export.gff3( segmentsRd, paste(tmpFile, "gff", sep=".") ) I could see some features corresponding to each gene. s0_+_0.CUFF.16.1_2 rtracklayer sequence_feature 85285 87682 . - . chr=s0;region=0.CUFF.20.1_4;nPos=12;nCounts=860;nTss=11 Actually, according to annotation file, 0.CUFF.16.1_2 has starting site at 48426, whereas 0.CUFF.20.1_4 has starting site at 85285. -- Sent via the guest posting facility at bioconductor.org.

Dear Soyeon Cha, When you split the counts into segments (as you did by calling 'segmentizeCounts'), each counts gets a name assigned. By default, these names have the format '<chromosome>_<strand>_<region>, but you can change it with the 'pattern' argument in 'segmentizeCounts'. As an example, 's0_+_0.CUFF.16.1_2' would be chromosome 's0', '+' strand, and region '0.CUFF.16.1_2'. How this relates to the genes in your organism depends on the regions that you defined initially. But since you have the coordinates of each start site, this should be easy to map back, even if you can't infer it from the region/segment name. 'nTSS' counts the number of identified start sites in each segment, as you also can see by inspecting 'z' with 'show'. Best wishes Julian On 08/03/2013 08:18 AM, from a TSSi user [guest] wrote: > > Hello, > I was able to successfully run TSSi program. However, I have several questions. > > My questions are > - how "0.CUFF.16.1_2" is related to "0.CUFF.20.1_4", > - how region is determined to each gene, > - how this file could be interpreted, > - and how nTss is determined? > > Thank you again for sparing your precious time. > Best regards. > > -- output of sessionInfo(): > > When I exported data by using the command bleow, > > x <- segmentizeCounts(counts=counts, start=start, chr=chromosome, region=region, strand=strand) > yFit <- normalizeCounts( x, fit=TRUE ) > z <- identifyStartSites( yFit ) > segmentsRd <- segmentsAsRangedData( z ) > export.gff3( segmentsRd, paste(tmpFile, "gff", sep=".") ) > > I could see some features corresponding to each gene. > s0_+_0.CUFF.16.1_2 rtracklayer sequence_feature 85285 87682 . - . chr=s0;region=0.CUFF.20.1_4;nPos=12;nCounts=860;nTss=11 > > Actually, according to annotation file, 0.CUFF.16.1_2 has starting site at 48426, whereas 0.CUFF.20.1_4 has starting site at 85285. > > > -- > Sent via the guest posting facility at bioconductor.org. > > _______________________________________________ > Bioconductor mailing list > Bioconductor at r-project.org > https://stat.ethz.ch/mailman/listinfo/bioconductor > Search the archives: http://news.gmane.org/gmane.science.biology.informatics.conductor >