User: Walter F. Baumann
Walter F. Baumann • 10
- Reputation:
- 10
- Status:
- New User
- Location:
- Last seen:
- 11 months, 3 weeks ago
- Joined:
- 1 year, 8 months ago
- Email:
- w************@gmail.com
Posts by Walter F. Baumann
<prev
• 15 results •
page 1 of 2 •
next >
0
votes
0
answers
527
views
0
answers
... 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 12 months ago by
Walter F. Baumann • 10
3
votes
0
answers
527
views
0
answers
... 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 ...
written 12 months ago by
Walter F. Baumann • 10
1
vote
1
answer
320
views
1
answer
... 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 ( ...
written 12 months ago by
Walter F. Baumann • 10
• updated 12 months ago by
thokall • 120
3
votes
2
answers
1.3k
views
2
answers
... 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 ...
written 13 months ago by
Walter F. Baumann • 10
• updated 12 months ago by
lahrmannu • 0
1
vote
1
answer
242
views
1
answer
... 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("+", "-" ...
written 16 months ago by
Walter F. Baumann • 10
• updated 16 months ago by
Michael Lawrence ♦ 10k
3
votes
1
answer
424
views
1
answer
... 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 ...
written 17 months ago by
Walter F. Baumann • 10
• updated 17 months ago by
Hervé Pagès ♦♦ 13k
0
votes
1
answers
370
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 17 months ago by
Walter F. Baumann • 10
1
vote
1
answer
370
views
1
answer
... 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 ...
written 18 months ago by
Walter F. Baumann • 10
• updated 18 months ago by
Martin Morgan ♦♦ 22k
0
votes
3
answers
1.8k
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 20 months ago by
Walter F. Baumann • 10
0
votes
1
answers
854
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 20 months ago by
Walter F. Baumann • 10
Latest awards to Walter F. Baumann
No awards yet. Soon to come :-)
Use of this site constitutes acceptance of our User
Agreement
and Privacy
Policy.
Powered by Biostar
version 2.2.0
Traffic: 233 users visited in the last hour