User: gundalav

gravatar for gundalav
gundalav260
Reputation:
260
Status:
Trusted
Location:
La La Land
Last seen:
3 months, 4 weeks ago
Joined:
6 years, 6 months ago
Email:
g*******@gmail.com

about me

Posts by gundalav

<prev • 49 results • page 1 of 5 • next >
0
votes
1
answer
244
views
1
answer
How to remove the squished track in Give AlignmentTrack()
... I have the following code: library(Gviz) library(Mus.musculus) library(data.table) library(OrganismDbi) library(TxDb.Mmusculus.UCSC.mm10.knownGene) txdb <- TxDb.Mmusculus.UCSC.mm10.knownGene # Define bam file ----- adipose_bam <- "myfile.bam" ...
rna-seq written 11 months ago by gundalav260 • updated 8 days ago by kathka0
0
votes
2
answers
1.2k
views
2
answers
Comment: C: How to add RG tags in a single cell BAM file where each group indicate barcode
... Yes I have. BTW what I mean by RG tag is the one indicated in every read not as `@RG` header. ...
written 12 months ago by gundalav260
6
votes
2
answers
1.2k
views
2
answers
How to add RG tag into the optional field in a BAM file
... I have a BAM file that looks like this: ![enter image description here][1] As you notice that the barcodes are included as part of read names. Then I'm trying to use a [tool called chromVAR][2] that require `RG` tag to be included in the BAM file. RG tags are used to distinguish reads from diffe ...
single-cell bam samtools sequencing written 12 months ago by gundalav260 • updated 11 months ago by finswimmer11k
1
vote
1
answer
972
views
1
answer
How to show all columns after performing GenomicRanges::subsetByOverlaps
... I have the following two [GenomicRanges][1] objects. First one `gr1` is this: library(GenomicRanges) set.seed(1) gr1 <- GRanges( seqnames=Rle(c("chr1", "chr2", "chr1", "chr3"), c(1, 3, 2, 4)), ranges=IRanges(1:10, width=10:1, names=head(letters,10)), ...
R bioconductor genomicrange gwas written 17 months ago by gundalav260 • updated 17 months ago by Sean Davis25k
0
votes
1
answer
731
views
1
answers
Comment: C: Can or should we perform feature counting with `bedtools intersect -c`?
... @ATPoint Thanks for your reply. Speed is less of a concern at the moment, but rather the **correctness** of Method1. I tried using `cut -f1,2,6,7,8,9`. But the Method1 is totally uncorrelated (different) from Method2 and Method3. I need to be convinced why Method1 is incorrect. Your enlightenment wi ...
written 20 months ago by gundalav260
3
votes
1
answer
731
views
5 follow
1
answer
Can or should we perform feature counting with `bedtools intersect -c`?
... I have the following BAM file and a annotation file. The BAM is a paired end alignment. What I want to do is to calculate reads (as fragment) count in I can do it with the following three methods. - **Method1:** [bedtools intersect][1] - **Method2:** [bedtools multicov][2] - **Method3:** [fea ...
bedtools featurecounts rna-seq written 20 months ago by gundalav260
0
votes
2
answers
1.6k
views
2
answers
Comment: C: How to recover treated/control count from DESeq2 output
... Thanks. `airway@assays$data$counts` is that normalized or not? If not how can I get the normalized one from S4? ...
written 2.1 years ago by gundalav260
11
votes
2
answers
1.6k
views
10 follow
2
answers
How to recover treated/control count from DESeq2 output
... I have the following DESeq2 code and the corresponding results suppressMessages(library(DESeq2)) suppressMessages(library(airway)) data(airway) airway_se <- airway airway_dds <- DESeqDataSet(airway_se, design = ~cell + dex) deseq <- DESeq(airway_dds) #> ...
rna-seq dseq2 written 2.1 years ago by gundalav260 • updated 2.1 years ago by Benn6.8k
0
votes
3
answers
3.6k
views
3
answers
Comment: C: How to calculate RPKM from featureCounts output
... @b.nota i.e. 100 + 30 = 130 (sum of 7th column of featureCounts output)?? ...
written 2.2 years ago by gundalav260
0
votes
3
answers
3.6k
views
3
answers
Comment: C: How to calculate RPKM from featureCounts output
... @b.nota in other words the `totalReads` in the calculation *is not the total assigned reads to the genes*? So I have to directly count from the BAM file? Would prefer to recode cause later need to calculate TPM as well. ...
written 2.2 years ago by gundalav260

Latest awards to gundalav

Popular Question 11 months ago, created a question with more than 1,000 views. For Non-Immune cell type gene expression database
Popular Question 11 months ago, created a question with more than 1,000 views. For Mirna Disease Annotated Corpus
Popular Question 11 months ago, created a question with more than 1,000 views. For Visualizing Up/Down Regulation Of Genes Based On Fold Change
Popular Question 11 months ago, created a question with more than 1,000 views. For Your Favorite Mirna-Gene Concordance Database Resource
Popular Question 11 months ago, created a question with more than 1,000 views. For Why mapping rate is low after adapter trimming paired-end data?
Popular Question 12 months ago, created a question with more than 1,000 views. For Non-Immune cell type gene expression database
Popular Question 12 months ago, created a question with more than 1,000 views. For Non-Immune cell type gene expression database
Popular Question 17 months ago, created a question with more than 1,000 views. For Non-Immune cell type gene expression database
Popular Question 20 months ago, created a question with more than 1,000 views. For Non-Immune cell type gene expression database
Popular Question 20 months ago, created a question with more than 1,000 views. For Visualizing Up/Down Regulation Of Genes Based On Fold Change
Popular Question 2.1 years ago, created a question with more than 1,000 views. For Visualizing Up/Down Regulation Of Genes Based On Fold Change
Epic Question 2.2 years ago, created a question with more than 10,000 views. For Creating Barplot With Horizontal Marker And Asterisks For Biological Paper
Popular Question 2.2 years ago, created a question with more than 1,000 views. For Best Way To Measure Mirna-Seq Abundance
Great Question 2.7 years ago, created a question with more than 5,000 views. For Creating Barplot With Horizontal Marker And Asterisks For Biological Paper
Popular Question 2.7 years ago, created a question with more than 1,000 views. For How To Make Cytoscape Export Network In Cluster Format (Non-Sif)
Great Question 3.9 years ago, created a question with more than 5,000 views. For Creating Barplot With Horizontal Marker And Asterisks For Biological Paper
Popular Question 3.9 years ago, created a question with more than 1,000 views. For How To Identify Tissue Specific Genes From Mixed Samples In Microarray Data
Popular Question 3.9 years ago, created a question with more than 1,000 views. For Best Way To Measure Mirna-Seq Abundance
Popular Question 3.9 years ago, created a question with more than 1,000 views. For How To Embed Ucsc Genome Browser To Web Page
Popular Question 3.9 years ago, created a question with more than 1,000 views. For How To Remove Noise In Microarray Data?
Popular Question 3.9 years ago, created a question with more than 1,000 views. For Best Practice For Extracting Gene Names From Pubmed Abstracts
Popular Question 3.9 years ago, created a question with more than 1,000 views. For Immune System Specific Gene Ontology Source
Popular Question 3.9 years ago, created a question with more than 1,000 views. For Visualizing Up/Down Regulation Of Genes Based On Fold Change
Popular Question 3.9 years ago, created a question with more than 1,000 views. For Ways To Do Gene Ontology Analysis For Microarray Data
Student 3.9 years ago, asked a question with at least 3 up-votes. For Creating Barplot With Horizontal Marker And Asterisks For Biological Paper

Help
Access

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