User: Walter F. Baumann

Reputation:
10
Status:
New User
Location:
Last seen:
3 weeks, 3 days ago
Joined:
9 months, 3 weeks ago
Email:
w************@gmail.com

Posts by Walter F. Baumann

<prev • 15 results • page 1 of 2 • next >
0
votes
0
answers
133
views
0
answers
Comment: C: Using sluth with Salmon data, error while sleuth_prep
... Hi, first of all thanks for moving this discussion. There is a bootstrap mode (--numBootstraps)in salmon and this added information that is needed by wasabi. I set it on 10. The bootstraps are needed to estimate the technical variance in your sample, which is needed by sleuth (and I guess this is wh ...
written 5 weeks ago by Walter F. Baumann 10
3
votes
0
answers
133
views
0
answers
Using sluth with Salmon data, error while sleuth_prep
... Hi, I want to use sleuth on the data salmon returned. In the first step, because sleuth needs actually kallisto, I used wasabi to generate the abundance.h5 file. However, after using so <- sleuth_prep(s2c, ~ condition, num_cores = 1) (s2c looks like that: ) sample condition 1 sample1 ...
R salmon sleuth salmon rna seq written 5 weeks ago by Walter F. Baumann 10
1
vote
1
answer
111
views
1
answer
Difference between SummarizeOverlaps and HTSeq
... Hi,  I compared the counts per gene of summarizeOverlaps and HTSeq (python). The correlation was ~0.98. Although the correlation is very good, I was surprised that it was not roughly or equal 1, because summarizeOverlaps is according to the documentation designed after the counting modes in HTSeq ( ...
R summarizeoverlaps htseqcounts written 6 weeks ago by Walter F. Baumann 10 • updated 6 weeks ago by thokall90
3
votes
2
answers
337
views
2
answers
SingleCellExperiment function in SingleCellExperiment package returns error (when reproducing example)
... Hi, I want to get into Single Cell RNA-Seq and found already a good tutorial (https://www.bioconductor.org/help/workflows/simpleSingleCell/). There a package called SingleCellExperiment is introduced. However, the package does not seem to work for me. I had problems with the installation (it always ...
R bioconductor single cell singlecellexperiment written 8 weeks ago by Walter F. Baumann 10 • updated 7 weeks ago by lahrmannu0
1
vote
1
answer
114
views
1
answer
Width in IRanges objects on "+" and "-" strand
... Hi,  I have to create some GRanges objects, which are located on the + and - strand (The features represent regions on transcripts.). Here a reproducible example:  GRanges(seqnames = Rle(c("chr1", "chr2", "chr3", "chr4")), ranges = IRanges(10:13, width=3), strand = Rle(c("+", "-" ...
iranges granges grangeslist written 5 months ago by Walter F. Baumann 10 • updated 5 months ago by Michael Lawrence9.8k
3
votes
1
answer
216
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 6 months ago by Walter F. Baumann 10 • updated 6 months ago by Hervé Pagès ♦♦ 13k
0
votes
1
answers
220
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 6 months ago by Walter F. Baumann 10
1
vote
1
answer
220
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 7 months ago by Walter F. Baumann 10 • updated 7 months ago by Martin Morgan ♦♦ 21k
0
votes
3
answers
640
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 9 months ago by Walter F. Baumann 10
0
votes
1
answers
241
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 9 months ago by Walter F. Baumann 10

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: 321 users visited in the last hour