User: Walter F. Baumann

Reputation:
10
Status:
New User
Location:
Last seen:
1 year, 5 months ago
Joined:
2 years, 1 month ago
Email:
w************@gmail.com

Posts by Walter F. Baumann

<prev • 15 results • page 1 of 2 • next >
0
votes
0
answers
740
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 17 months ago by Walter F. Baumann 10
3
votes
0
answers
740
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 17 months ago by Walter F. Baumann 10
1
vote
1
answer
447
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 17 months ago by Walter F. Baumann 10 • updated 17 months ago by thokall140
3
votes
2
answers
1.7k
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 18 months ago by Walter F. Baumann 10 • updated 18 months ago by lahrmannu0
1
vote
1
answer
300
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 21 months ago by Walter F. Baumann 10 • updated 21 months ago by Michael Lawrence10k
3
votes
1
answer
555
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 23 months ago by Walter F. Baumann 10 • updated 23 months ago by Hervé Pagès ♦♦ 13k
0
votes
1
answer
443
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 23 months ago by Walter F. Baumann 10
1
vote
1
answer
443
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 23 months ago by Walter F. Baumann 10 • updated 23 months ago by Martin Morgan ♦♦ 23k
0
votes
3
answers
2.5k
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 2.1 years ago by Walter F. Baumann 10
0
votes
1
answer
1.1k
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 2.1 years ago by Walter F. Baumann 10

Latest awards to Walter F. Baumann


Help
Access

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