User: Yang Liao

Yang Liao130
Reputation:
130
Status:
Trusted
Location:
Australia
Last seen:
6 days, 21 hours ago
Joined:
6 years, 2 months ago
Email:
l***@wehi.edu.au

Posts by Yang Liao

<prev • 21 results • page 2 of 3 • next >
0
296
views
0
... From the error message in your post, it is a "general format error" encountered by featureCounts. In other words, featureCounts couldn't correctly parse the current alignment record in the BAM file. I generated a BAM file myself using rsem-tbam2bgam but featureCounts ran correctly on it. Can you ...
written 8 months ago by Yang Liao130 • updated 8 months ago by Gordon Smyth38k
1
415
views
1
... Thanks for the details. I noticed that you ran propmapped() on the BAM files, then gave the return value from propmapped() directly to featureCounts(). As you may have noticed, the return value from propmapped() is a data.frame object, and featureCounts cannot take a data.frame as the argument for i ...
written 8 months ago by Yang Liao130
2
379
views
2
... Hi Gwenael, The "GZIP ERROR:-2" message sounds like an old bug in featureCounts. Which version of featureCounts did you use? If it was not the latest version (v1.6.3), you can try our latest version, which is available on   https://sourceforge.net/projects/subread/files/subread-1.6.3/ since this ...
written 10 months ago by Yang Liao130
1
714
views
1
... I'm facing the same problem and I found that the makeTxDbFromGFF() function already dropped the exons if they span multiple chromosomes and/or multiple strands. I built a GTF file that contains only four exons from one gene, but from different strands of chr1. The makeTxDbFromGFF() function gave wa ...
written 15 months ago by Yang Liao130
3
713
views
3
... Thanks for the detailed explanation! Yes, I can now reproduce the error by adding "-O 0" to the compiler parameters. It is a little surprising to see that an optimisation option can change the legitimate behaviour of the compiler. I have removed this keyword from my codes.  ...
written 16 months ago by Yang Liao130
3
713
views
3
... I understand that the "inline" keyword instructs the compiler to embed the function into all its callers, and modern compilers are smarter to automatically make the decision, leaving this keyword obsolete. However, it should be harmless when you build a shared object from the codes including inline ...
written 16 months ago by Yang Liao130
1
670
views
1
... Thanks! I've downloaded the data, and found that featureCounts had a bug in the routine for fixing the BAM format. It did not keep the read group names when fixing the BAM files from Picard or STAR.  This bug has been fixed. The next released version should work on your BAM files. Cheers, Yang ...
written 2.0 years ago by Yang Liao130
1
670
views
1
... It still works when I converted the SAM file to BAM by using Picard (SortSam). Is it possible if you can send us the BAM file you're using?  Cheers, Yang # /picard-tools-2.9.4/bin/SortSam I=user-RG.bam o=user-RG-Picard.bam SO=coordinate > x<-featureCounts(c('user-RG.bam','user-RG-Picar ...
written 2.0 years ago by Yang Liao130
1
670
views
1
... Hi Ge Tan, I used Rsubread_1.26.1 to assign the lines you provided, and the read groups were correctly treated by featureCounts: x<-featureCounts(c('user-RG.bam','user-RG.sam') , annot.inbuilt = "hg19", byReadGroup = T) NCBI RefSeq annotation for hg19 (build 37.2) is used. //================= ...
written 2.0 years ago by Yang Liao130
1
2.4k
views
1
... It looks that the current directory has no enough space, or some reads in the SAM file are extremely long (e.g., longer than 1200 bps).  ...
written 4.7 years ago by Yang Liao130

Latest awards to Yang Liao

Scholar 15 months ago, created an answer that has been accepted. For A: Error with featureCounts