Admin: Hervé Pagès

gravatar for Hervé Pagès
Hervé Pagès ♦♦ 13k
Reputation:
13,010
Status:
Trusted
Location:
United States
Last seen:
1 week ago
Joined:
12 years, 9 months ago
Email:
h*****@fhcrc.org

Posts by Hervé Pagès

<prev • 1,585 results • page 2 of 159 • next >
0
votes
0
answers
71
views
0
answers
Comment: C: Error: the "phase" metadata column must contain 0, 1, or 2, for all the CDS feat
... Would you mind providing a reproducible example and your `sessionInfo()`? Thanks! ...
written 5 weeks ago by Hervé Pagès ♦♦ 13k
0
votes
1
answers
192
views
1
answers
Answer: A: New Version of R cannot coerce "GRanges" to "RangesList"
... Hi, In BioC 3.7, the RangesList class was replaced with the IntegerRangesList class. So: library(GenomicRanges) gr <- GRanges(Rle(c("A", "B"), 3), IRanges(9:4, 10)) gr # GRanges object with 6 ranges and 0 metadata columns: # seqnames ranges strand # <Rle> <IRanges> ...
written 9 weeks ago by Hervé Pagès ♦♦ 13k
0
votes
1
answers
84
views
1
answers
Answer: A: How does GenomicFeature calculate 3'UTR?
... Hi, UCSC and Ensembl always consider the stop codon to be part of the CDS. This means that it's not considered to be part of the 3' UTR. For consistency with UCSC and Ensembl, makeTxDbFromGFF() imports the stop codon as part of the CDS. The makeTxDbFromGRanges() function, which makeTxDbFromGFF() is ...
written 11 weeks ago by Hervé Pagès ♦♦ 13k
2
votes
0
answers
109
views
0
answers
News: New Bioconductor package for dbSNP151
... Hi all, A new SNPlocs package is now available with the content of dbSNP build 151:   http://bioconductor.org/packages/SNPlocs.Hsapiens.dbSNP151.GRCh38 SNPlocs packages contain SNP locations and alleles for various dbSNP builds. SNPs can be looked up by chromosome, range, or rs id, using snpsBySe ...
bsgenome news snplocs.hsapiens.dbsnp151.grch38 written 11 weeks ago by Hervé Pagès ♦♦ 13k
0
votes
2
answers
161
views
2
answers
Comment: C: nearest without overlapping
... Note that the relist/unlist and split/unsplit combos are elegant, powerful, and reliable idioms, but the split/unlist combo is a risky one. It only works if the splitting factor is sorted which is indeed the case here because the class of Hits object hits is actually SortedByQueryHits. And because ...
written 11 weeks ago by Hervé Pagès ♦♦ 13k
0
votes
2
answers
161
views
2
answers
Comment: C: nearest without overlapping
... I see. Yes finding the nearest-but-not-overlapping ranges requires a little bit more work. One way to go is to take advantage of the fact that precede() and nearest() both exclude overlapping ranges: ### A version of distanceToNearest() that finds all the ### nearest-but-not-overlapping ranges in ...
written 11 weeks ago by Hervé Pagès ♦♦ 13k
0
votes
1
answers
90
views
1
answers
Answer: A: extract codon position via BSgenome
... Hi, In other words you want to skip matches that overlap with a previous match. Note that this is what regular expressions do:  > gregexpr("TAT", "ATATATATGCAT") [[1]] [1] 2 6 attr(,"match.length") [1] 3 3 attr(,"index.type") [1] "chars" attr(,"useBytes") [1] TRUE But in your case it seems th ...
written 11 weeks ago by Hervé Pagès ♦♦ 13k
0
votes
2
answers
161
views
2
answers
Answer: A: nearest without overlapping
... Hi, Don't remove anything from subject **before** you call distanceToNearest() or you will loose hits that you want to keep: library(GenomicRanges) x <- GRanges(c("chr1:2-4", "chr1:5-8")) subject <- GRanges("chr1:7-10") distanceToNearest(x, subject) # Hits object with 2 hits and 1 metadata ...
written 11 weeks ago by Hervé Pagès ♦♦ 13k
0
votes
1
answers
165
views
1
answers
Comment: C: Rhtslib fails to find liblzma.dylib on Mac OSX
... That seems to have addressed the issue. I consequently bumped the version of all packages that have Rhtslib in their LinkingTo field (e.g. csaw, bamsignals, deepSNV, etc...), but only in release (RELEASE_3_7 branch) for now. See for example commit 6ea86c79cfbc2b9d765ed31b35364b5eb39c7671 in the RELE ...
written 3 months ago by Hervé Pagès ♦♦ 13k
0
votes
1
answers
129
views
1
answers
Answer: A: Non-Redundant Edges in Connected Components
... Hi Dario, Here is a fast solution based on an unexported utility defined in the S4Vectors package: connectedCompEdges <- function(concomp) { hits <- S4Vectors:::makeAllGroupInnerHits(lengths(concomp), hit.type=1L) nodes <- unlist(concomp, use.names=FALSE) from <- nodes[que ...
written 3 months ago by Hervé Pagès ♦♦ 13k

Latest awards to Hervé Pagès

Teacher 6 weeks ago, created an answer with at least 3 up-votes. For A: Extracting Coordinates of startcodon from Grangeslist
Scholar 3 months ago, created an answer that has been accepted. For A: coverage() output SimpleRleList cannot be converted to GRanges
Scholar 3 months ago, created an answer that has been accepted. For A: Extracting Coordinates of startcodon from Grangeslist
Scholar 3 months ago, created an answer that has been accepted. For A: Feature request: syntax for dropping unused seqlevels in GRanges objects
Scholar 4 months ago, created an answer that has been accepted. For A: coverage() output SimpleRleList cannot be converted to GRanges
Teacher 4 months ago, created an answer with at least 3 up-votes. For A: Extracting Coordinates of startcodon from Grangeslist
Teacher 14 months ago, created an answer with at least 3 up-votes. For A: Extracting Coordinates of startcodon from Grangeslist
Teacher 14 months ago, created an answer with at least 3 up-votes. For C: How to create tx2gene data.frame when there's no TxDb object for the organism yo
Commentator 14 months ago, created a comment with at least 3 up-votes. For C: Error in GRanges: could not find symbol "recursive" in environment of the generi
Scholar 14 months ago, created an answer that has been accepted. For A: Extracting Coordinates of startcodon from Grangeslist
Teacher 14 months ago, created an answer with at least 3 up-votes. For A: overlapsAny for elements in GRangesList
Teacher 14 months ago, created an answer with at least 3 up-votes. For A: How do I merge a list of GRanges?
Teacher 14 months ago, created an answer with at least 3 up-votes. For A: about setdiff() and meta-information in GRanges
Teacher 14 months ago, created an answer with at least 3 up-votes. For A: extract subsequence(s) from DNAStringSet with getSeq
Scholar 14 months ago, created an answer that has been accepted. For A: coverage() output SimpleRleList cannot be converted to GRanges
Popular Question 14 months ago, created a question with more than 1,000 views. For unique, duplicated, order, etc... on a GRanges object
Scholar 14 months ago, created an answer that has been accepted. For A: makeTxDbFromUCSC("hg38", "refGene") gives "not supported" error
Scholar 14 months ago, created an answer that has been accepted. For C: How to create tx2gene data.frame when there's no TxDb object for the organism yo
Scholar 14 months ago, created an answer that has been accepted. For A: about setdiff() and meta-information in GRanges
Scholar 14 months ago, created an answer that has been accepted. For A: extract subsequence(s) from DNAStringSet with getSeq
Scholar 14 months ago, created an answer that has been accepted. For A: Feature request: syntax for dropping unused seqlevels in GRanges objects
Scholar 14 months ago, created an answer that has been accepted. For A: What's the best way to retrieve the sizes (lengths) of chromosomes for a referen
Teacher 14 months ago, created an answer with at least 3 up-votes. For A: Error when running "locateVariants" function in "VariantAnnotation" package
Scholar 14 months ago, created an answer that has been accepted. For A: A bug of "seqlevels<-" on TxDb
Scholar 14 months ago, created an answer that has been accepted. For A: subsetByOverlaps with Multimapping

Help
Access

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