User: Leonard Goldstein

Reputation:
70
Status:
Trusted
Location:
United States
Last seen:
3 weeks, 4 days ago
Joined:
3 years, 1 month ago
Email:
g*******@gene.com

Posts by Leonard Goldstein

<prev • 39 results • page 1 of 4 • next >
0
votes
4
answers
128
views
4
answers
Answer: A: Cannot extract TxFeatures for a defined GRanges object
... exonsBy seems to be dropping a lot of genes, including MYC: > tx <- makeTxDbFromBiomart() > txs_by_gene <- transcriptsBy(tx, "gene") > exs_by_gene <- exonsBy(tx, "gene") > length(txs_by_gene) [1] 63967 > length(exs_by_gene) [1] 36751 > subsetByOverlaps(txs_by_gene, GRang ...
written 25 days ago by Leonard Goldstein70
0
votes
4
answers
128
views
4
answers
Answer: A: Cannot extract TxFeatures for a defined GRanges object
... It looks like MYC is missing from the TxDb object you created, so this is a question about the makeTxDbFromBiomart function in the GenomicFeatures package. Maybe repost your question and add GenomicFeatures as one of the tags, and hopefully someone can help. Alternatively you can try using one of th ...
written 25 days ago by Leonard Goldstein70
0
votes
1
answers
45
views
1
answers
Answer: A: SGFeatureCounts function in SGSeq package failed to generate splice graphs
... Hi Magda, Can you please provide more details -- what commands were used to generate the SGFeatureCounts object and what package versions were used i.e. output of sessionInfo(). For larger analyses I strongly recommend not to use wrapper functions like analyzeFeatures() and instead run individual a ...
written 26 days ago by Leonard Goldstein70
0
votes
1
answers
63
views
1
answers
Answer: A: SGSeq error on getBamInfo
... The error message should list the problematic files but it returns an empty string. Looks like something went wrong that was not caught by the internal checks (probably unrelated to the XS tag). Are your BAM file paths correct and are index files present? You can try running the following and see if ...
written 6 weeks ago by Leonard Goldstein70
0
votes
3
answers
111
views
3
answers
Answer: A: R session crash with SGSeq::convertToTxFeatures
... Hi, since I can't reproduce the problem it is difficult to troubleshoot. As an alternative to TxDb objects, the convertToTxFeatures function also accepts a GRangesList of exons grouped by transcript. So you could try something like this > tx <- exonsBy(hsaEnsembl, "tx", use.names = TRUE) > ...
written 7 weeks ago by Leonard Goldstein70
0
votes
3
answers
111
views
3
answers
Answer: A: R session crash with SGSeq::convertToTxFeatures
... I could run this without any issues. If the problem persists maybe you can try reinstalling R and Bioconductor?  > library(GenomicFeatures) > library(SGSeq) > hsaEnsembl <- makeTxDbFromBiomart(dataset="hsapiens_gene_ensembl") > txf <- convertToTxFeatures(hsaEnsembl) > > se ...
written 8 weeks ago by Leonard Goldstein70
0
votes
1
answers
337
views
1
answers
Comment: C: SGSeq: moving toward diffex from SGSeq analysis
... Hi Sylvain,  Yes variantID() and eventID() are accessor functions for metadata columns in the SGVariant/Counts objects. You can see all available columns by typing mcols(sgv.counts). To extract gene names you can simply use geneNames(sgv.counts). Note this will return an empty CharacterList unless ...
written 11 weeks ago by Leonard Goldstein70
0
votes
1
answers
337
views
1
answers
Answer: A: SGSeq: moving toward diffex from SGSeq analysis
... Hi Sylvain,  Not sure I completely follow your question. You seem to be confusing SGFeatureCounts and SGVariantCounts objects. The idea behind section 11 of the vignette is to test for differential splice variant usage within a splice event, rather than differential exon usage within a gene. This a ...
written 12 weeks ago by Leonard Goldstein70
0
votes
1
answers
554
views
1
answers
Comment: C: Isoforms detection using bioconductor
... Based on the error message it looks like you are working with less than 4 samples, but you try to make plots for columns 1-4 of the SGFeatureCounts object. Regarding out-of-bound plots, the vignette includes several examples on how to specify what genomic regions to include in the plots. Please have ...
written 3 months ago by Leonard Goldstein70
0
votes
1
answers
554
views
1
answers
Comment: C: Isoforms detection using bioconductor
... The input for getBamInfo is a data.frame with columns sample_name and file_bam (see example below). Please read the vignette or help pages. > si <- data.frame(sample_name = c("a", "b"), file_bam = c("file_a.bam", "file_b.bam"), stringsAsFactors = FALSE) > si <- getBamInfo(si) ...
written 3 months ago by Leonard Goldstein70

Latest awards to Leonard Goldstein

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