User: Walter F. Baumann

Reputation:
0
Status:
New User
Location:
Last seen:
1 week, 1 day ago
Joined:
4 months ago
Email:
w************@gmail.com

Posts by Walter F. Baumann

<prev • 10 results • page 1 of 1 • next >
2
votes
1
answer
79
views
1
answer
Extracting Coordinates of startcodon from Grangeslist
... Hi,  I work with ribosome profiling data and I want to check the coverage of reads on the  start codon.  My start is with an GRAngesList object I obtained from cdsBy(txdb, by=„tx).  In the end I want to obtain a GrangesList object with transcript id as names, and start of startcodon and stop of s ...
biomart grangeslist bioconductor ribosome profiling cds written 29 days ago by Walter F. Baumann 0 • updated 28 days ago by Hervé Pagès ♦♦ 12k
0
votes
1
answers
106
views
1
answers
Comment: C: Problems with summarizeOverlaps
... Thanks. I tried to apply the given examples of the Rsamtools and the GenomicAlignments manual on my case, but so far I was not successful. (This is a single-end RNA-Seq experiment).  myalignment <- BamFile(file=„/pathto/Aligned.bam", index="/pathto/Aligned.bam.bai", ...
written 4 weeks ago by Walter F. Baumann 0
1
vote
1
answer
106
views
1
answer
Problems with summarizeOverlaps
... Hi, I want to create a table with counts per region (column), per transcript or gene (row). First thing I did was to define the regions: txdb <- makeTxDbFromGFF("/pathto/gencode.v19.annotation.gtf", dataSource="Gencode", organism="Homo sapiens", format="gtf") txdb_5utr_transcript <- five ...
rnaseq bioconductor summarizeoverlaps written 6 weeks ago by Walter F. Baumann 0 • updated 6 weeks ago by Martin Morgan ♦♦ 20k
0
votes
3
answers
246
views
3
answers
Answer: A: DESeq2 Following RSEM
... The DESeq2 manuals says: "As input, the DESeq2 package expects count data as obtained", which means just handles whole numbers. I doubt that "round" is a good approach here, because you change the counts. Where does the e.g. 0.64 count from Sample1 A1BG come from? ...
written 3 months ago by Walter F. Baumann 0
0
votes
1
answers
127
views
1
answers
Comment: C: Error using biocLite
... It should be biocLite("org.Hs.eg.db") (with quotation marks)... However, I just tried it and it works as expected, but I use R 3.3.1.  ...
written 3 months ago by Walter F. Baumann 0
0
votes
0
answers
174
views
0
answers
Finding candidate binding sites for RNA binding proteins in R
...   Hi, I need to scan mRNAs for binding motifs of RBPs (RNA binding proteins) in human.  After some research I came across this tutorial on how to find candidate binding sites of known transcription factors. So far I did not saw an implemented tool to directly work with RNA or RBPs. I downloaded ...
0
votes
1
answers
253
views
1
answers
Comment: C: Troubles with new version of GenomicFeatures(?)
... Thank you for the useful link, but it unfortunately did not solve my problem. I reinstalled the dependencies of GenomicFeatures: biocLite(c("BiocGenerics", "S4Vectors", "IRanges", "GenomeInfoDb", "GenomicRanges", "AnnotationDbi", "GenomicFeatures")) I use the same versions of packages on my other ...
written 3 months ago by Walter F. Baumann 0
3
votes
1
answer
253
views
1
answer
Troubles with new version of GenomicFeatures(?)
... Hi, I have the following command which always worked well, until now.      txdb <- makeTxDbFromGFF("/path/to/gencode.v19.annotation.gtf", dataSource="Gencode", organism="Homo sapiens", format="gtf") Now I get the following output:      Import genomic features from the file as a GRanges object ...
1
vote
1
answer
129
views
1
answer
summarizeOverlaps counting exons and 5UTR reads
... Hey, I try to compare counts from 2 RNA seq data. I am interested in the 5'UTR and want to compare this number with the total number of reads per transcript. Here my code:      txdb is the gtf file from Gencode (v19) "Comprehensive gene annotation"      txdb_5utr_transcript <- fiveUTRsByTransc ...
rnaseq summarizeoverlaps written 3 months ago by Walter F. Baumann 0 • updated 3 months ago by Hervé Pagès ♦♦ 12k
3
votes
3
answers
171
views
3
answers
Comparing two RNAseq experiments using summarizeOverlaps
... Hi,  I have two RNA-seq experiments. I first checked them to see whether they are single or pairedEnd using: bamtools view -c -F 1 myfile.bam, bamtools view -c -f 1 myfile2.bam There I saw that all reads are singleEnd. In the next step I uploaded my data in R, and counted the reads per feature by ...
rnaseq R bioconductor summarizeoverlaps rna-seq written 4 months ago by Walter F. Baumann 0 • updated 3 months ago by Hervé Pagès ♦♦ 12k

Latest awards to Walter F. Baumann

No awards yet. Soon to come :-)

Help
Access

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