User: tangming2005

gravatar for tangming2005
tangming200590
Reputation:
90
Status:
Trusted
Location:
United States
Last seen:
2 months ago
Joined:
4 years, 1 month ago
Email:
t***********@gmail.com

Posts by tangming2005

<prev • 114 results • page 2 of 12 • next >
2
votes
1
answer
466
views
1
answer
POU5F1(OCT4) /ENTREZID: 5460 is not in the TxDb.Hsapiens.UCSC.hg19.knownGene
... Hi, I was trying to get the coordinates of some genes, and was surprised that POU5F1 gene is not in the txdb. Thank you for looking into this. Tommy library(TxDb.Hsapiens.UCSC.hg19.knownGene) txdb<- TxDb.Hsapiens.UCSC.hg19.knownGene hg19.genes<- genes(txdb) library(AnnotationDbi) libra ...
txdb.hsapiens.ucsc.hg19.knowngene written 13 months ago by tangming200590 • updated 12 months ago by jma199130
0
votes
1
answers
763
views
1
answers
Comment: C: reads2pairs from TEQC error:Error in reads2pairs(reads) : read pair IDs do not s
... samtools view -F 256 will remove some reads, but not all not uniquely mapped reads. Have you tried TEQC on any BWA-MEM aligned bam files? from Heng Li: https://www.biostars.org/p/59281/#59303 Eland is probably the first short read aligner. It reports a tag, which can be '[UR][0-2]|NM', for each r ...
written 14 months ago by tangming200590
0
votes
1
answers
763
views
1
answers
Comment: C: reads2pairs from TEQC error:Error in reads2pairs(reads) : read pair IDs do not s
... I was using 3.1.0. A side question is how TEQC determine the not uniquely mapped reads for BWA-MEM aligned bam files. I read here: https://www.biostars.org/p/76893/   https://www.biostars.org/p/167892/   BWA-mem added XA tag to a read if it is mapped to multiple places. some people also use  sam ...
written 14 months ago by tangming200590
0
votes
1
answers
763
views
1
answers
Comment: C: reads2pairs from TEQC error:Error in reads2pairs(reads) : read pair IDs do not s
... I actually used bwa mem to map to get the bam files, so bwa mem is keeping the multi-mapped reads I guess.  You mentioned that `get.reads` only read in the primary alignments? but why it is still there?  Maybe I need to remove those reads in bam before feeding into TEQC? Thank you! Tommy ...
written 14 months ago by tangming200590
0
votes
1
answer
763
views
1
answer
reads2pairs from TEQC error:Error in reads2pairs(reads) : read pair IDs do not seem to be unique
... Hi, I was running QC for my whole exome data (paired-end sequencing) and got this error: > readpairs <- reads2pairs(reads) Error in reads2pairs(reads) : read pair IDs do not seem to be unique sessionInfo() R version 3.2.3 (2015-12-10) Platform: x86_64-pc-linux-gnu (64-bit) Running under: ...
teqc written 14 months ago by tangming200590
0
votes
2
answers
870
views
2
answers
Comment: C: TCGAbiolinks Error in gene.location$ensembl_gene_id : $ operator is invalid f
... using 2.3.22 gave me this   > LUSC_rna_data <- GDCprepare(query_rna_LUSC.hg38) |========================================================================================================| 100%    1 MB |============================================================================================ ...
written 18 months ago by tangming200590
0
votes
2
answers
870
views
2
answers
Comment: C: TCGAbiolinks Error in gene.location$ensembl_gene_id : $ operator is invalid f
... yes, that worked if I install psych first. ...
written 18 months ago by tangming200590
0
votes
2
answers
870
views
2
answers
Comment: C: TCGAbiolinks Error in gene.location$ensembl_gene_id : $ operator is invalid f
... installing from github gave me this error. thanks for looking into it > devtools::install_github(repo = "BioinformaticsFMRP/TCGAbiolinks") Downloading GitHub repo BioinformaticsFMRP/TCGAbiolinks@master from URL https://api.github.com/repos/BioinformaticsFMRP/TCGAbiolinks/zipball/master Installin ...
written 18 months ago by tangming200590
0
votes
2
answers
870
views
2
answers
Comment: C: TCGAbiolinks Error in gene.location$ensembl_gene_id : $ operator is invalid f
... Now it is randomly use genome annotations... it is human data set, but using elephant genes...   > LUAD_rna_data <- GDCprepare(query_rna_LUAD.hg38) |================================================================================| 100%    1 MB                 |  33% |======================= ...
written 18 months ago by tangming200590
0
votes
2
answers
870
views
2
answers
Comment: C: TCGAbiolinks Error in gene.location$ensembl_gene_id : $ operator is invalid f
... thanks, there is another error showed up.   > LUAD_rna_data <- GDCprepare(query_rna_LUAD.hg38) |================================================================================| 100%    1 MB                 |  18% |============================================================================ ...
written 18 months ago by tangming200590

Latest awards to tangming2005

Popular Question 5 months ago, created a question with more than 1,000 views. For IdeogramTrack from Gviz with chromosome name "4" rather "chr4"
Popular Question 17 months ago, created a question with more than 1,000 views. For How to bin a GRanges to small bins
Popular Question 17 months ago, created a question with more than 1,000 views. For DESeq2 for ChIP-seq differential peaks
Great Question 17 months ago, created a question with more than 5,000 views. For DESeq2 for ChIP-seq differential peaks
Centurion 17 months ago, created 100 posts.
Popular Question 17 months ago, created a question with more than 1,000 views. For Gviz GeneRegionTrack draw only one gene in the range
Popular Question 17 months ago, created a question with more than 1,000 views. For annotate GRanges with upstream or downstream nearest gene
Popular Question 2.5 years ago, created a question with more than 1,000 views. For DESeq2 for ChIP-seq differential peaks
Popular Question 2.5 years ago, created a question with more than 1,000 views. For annotate GRanges with upstream or downstream nearest gene
Supporter 2.5 years ago, voted at least 25 times.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.2.0
Traffic: 141 users visited in the last hour