User: gundalav

gravatar for gundalav
gundalav290
Reputation:
290
Status:
Trusted
Location:
La La Land
Last seen:
2 months, 2 weeks ago
Joined:
6 years, 10 months ago
Email:
g*******@gmail.com

about me

Posts by gundalav

<prev • 50 results • page 1 of 5 • next >
1
vote
1
answer
135
views
1
answer
How to obtain distinct/uniqe rows from GenomicRanges object
... I have the following [GenomicRanges][1] object created with this: library(GenomicRanges) gr <- GRanges(seqnames = "chr1", strand = c("+", "-","-", "+"),ranges = IRanges(start = c(1,3,3,5), width = 3)) gr That looks like this: GRanges object with 4 ranges and 0 metadata columns ...
R bioconductor genomicrange written 12 weeks ago by gundalav290 • updated 12 weeks ago by zx87548.2k
0
votes
1
answer
316
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 15 months ago by gundalav290 • updated 4 months ago by kathka0
0
votes
2
answers
1.6k
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 16 months ago by gundalav290
6
votes
2
answers
1.6k
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 16 months ago by gundalav290 • updated 15 months ago by finswimmer12k
2
votes
1
answer
1.2k
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 21 months ago by gundalav290 • updated 21 months ago by Sean Davis25k
0
votes
1
answer
847
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 2.0 years ago by gundalav290
4
votes
1
answer
847
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 2.0 years ago by gundalav290
0
votes
2
answers
2.0k
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.4 years ago by gundalav290
11
votes
2
answers
2.0k
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.4 years ago by gundalav290 • updated 2.4 years ago by Benn7.7k
0
votes
3
answers
4.0k
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.6 years ago by gundalav290

Latest awards to gundalav

Popular Question 7 months ago, created a question with more than 1,000 views. For Non-Immune cell type gene expression database
Popular Question 7 months ago, created a question with more than 1,000 views. For Mirna Disease Annotated Corpus
Popular Question 7 months ago, created a question with more than 1,000 views. For Visualizing Up/Down Regulation Of Genes Based On Fold Change
Popular Question 15 months ago, created a question with more than 1,000 views. For Non-Immune cell type gene expression database
Popular Question 15 months ago, created a question with more than 1,000 views. For Mirna Disease Annotated Corpus
Popular Question 15 months ago, created a question with more than 1,000 views. For Visualizing Up/Down Regulation Of Genes Based On Fold Change
Popular Question 15 months ago, created a question with more than 1,000 views. For Your Favorite Mirna-Gene Concordance Database Resource
Popular Question 15 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 16 months ago, created a question with more than 1,000 views. For Non-Immune cell type gene expression database
Popular Question 16 months ago, created a question with more than 1,000 views. For Non-Immune cell type gene expression database
Popular Question 21 months ago, created a question with more than 1,000 views. For Non-Immune cell type gene expression database
Popular Question 2.0 years ago, created a question with more than 1,000 views. For Non-Immune cell type gene expression database
Popular Question 2.0 years ago, created a question with more than 1,000 views. For Visualizing Up/Down Regulation Of Genes Based On Fold Change
Popular Question 2.4 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.6 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.6 years ago, created a question with more than 1,000 views. For Best Way To Measure Mirna-Seq Abundance
Great Question 3.0 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.0 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 4.3 years ago, created a question with more than 5,000 views. For Creating Barplot With Horizontal Marker And Asterisks For Biological Paper
Popular Question 4.3 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 4.3 years ago, created a question with more than 1,000 views. For Best Way To Measure Mirna-Seq Abundance
Popular Question 4.3 years ago, created a question with more than 1,000 views. For How To Embed Ucsc Genome Browser To Web Page
Popular Question 4.3 years ago, created a question with more than 1,000 views. For How To Remove Noise In Microarray Data?
Popular Question 4.3 years ago, created a question with more than 1,000 views. For Best Practice For Extracting Gene Names From Pubmed Abstracts
Popular Question 4.3 years ago, created a question with more than 1,000 views. For Immune System Specific Gene Ontology Source

Help
Access

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