similar to: truncated warning messages

Hi, I was wondering I'm going about this in the correct way. I need to test if there are coding sequences or exons in hg19 which match a string of 100bp "D" i.e. [A,G or T]. However I'm getting a strange result. I get a hit on chr7, using the 100bp search however when I search with 60bp sequence of "D" I don't get any hits. library("BSgenome")

Hi, I am currently building an R package and I am facing a peculiar problem where some of the functions does not work within the package. However, if I source the script the function works. For example, in a method for parallelization of analysis on each chromosome simultaneously I am receiving error at the following position of the code: # this profile the information chromosome wise and

Hello, I am using R 2.11.0. I have a curious problem where I get a warning in R CMD check which is seemingly not relevant to my Rd file. The warning says : * checking Rd \usage sections ... WARNING Bad \usage lines found in documentation object 'enrichmentCalc': <unescaped bksl>S4method{enrichmentCalc}{GenomeDataList, BSgenome}(rs, organism, seqLen=NULL, ...) <unescaped

Hello, I am using R 2.11.0. I have a curious problem where I get a warning in R CMD check which is seemingly not relevant to my Rd file. The warning says : * checking Rd \usage sections ... WARNING Bad \usage lines found in documentation object 'enrichmentCalc': <unescaped bksl>S4method{enrichmentCalc}{GenomeDataList, BSgenome}(rs, organism, seqLen=NULL, ...) <unescaped

I'm in the process of converting some S3 methods to S4 methods. I have this function : setGeneric("enrichmentCalc", function(rs, organism, seqLen, ...){standardGeneric("enrichmentCalc")}) setMethod("enrichmentCalc", c("GenomeDataList", "BSgenome"), function(rs, organism, seqLen, ...) { ... ... ... })

Hello, While doing some enrichment tests using chisq.test() with simulated p-values, I noticed some strange behaviour. The computed p-value was extremely small, so I decided to dig a little deeper and debug chisq.test(). I noticed then that the simulated statistics returned by the following call tmp <- .Call(C_chisq_sim, sr, sc, B, E) were all the same, very small numbers. This, at first,

Hello, Everyone, how to use "loop" to make the process automatic and fast? When compute each sample, the script type in R almost the same, just the input and output file's name is changed(chr1 change to chr2, chr3,chr4...). The first sample's script like this: >chr1=MEDIPS.readAlignedSeqences(BSgenome="hg19", file="chr1",numrows= )

Hi all, I am using the package ChIPpeakAnno, and I have a problem with the function getAllPeakSequence. This is related to object oriented programming I think, I have the following message: > peaksWithSequences <- getAllPeakSequence(peakList, upstream = 100, > downstream = 100, genome = Hsapiens) Error in validObject(.Object) : invalid class "GRanges" object: superclass

Dear all, I have data.frame object in R. I want to export it in tab-delimited file with several lines of header initiated with comment sign (#). I do not know how to do that in R. Could you please give helps on this problem? Thanks in advance. Best, Jian-Feng, ################################################################## The lines I want to write in the header lines look like, with words

Hi, Found in the man page for update.packages: repos: character vector, the base URL(s) of the repositories to use, i.e., the URL of the CRAN master such as ?"http://cran.r-project.org"? or its Statlib mirror, ?"http://lib.stat.cmu.edu/R/CRAN"?.Can be ?NULL? to install from local files (?.tar.gz? for source packages). - Using

Over the last several days, I've been working hard to get all of the Fedora R packages rebuilt against R 4.0 in rawhide (in the F33-R-4 side tag). With the exception of R-biomaRt, R-BSgenome, R-GenomicAlignments, and R-rtracklayer, I believe everything is built and updated to the latest versions. And of those packages, they're all ready to go when Fedora infrastructure is working reliably

BioInformatics Scientist Job Code: 10-TR25 Location: Cambridge, MA Description We are seeking a highly motivated, independent bioinformatics scientist to join a group of scientists, analysts and programmers to develop tools and methodologies for large-scale gene expression data analysis. The group supports a variety of research projects in target and drug discovery as well as biomarker

On Tue, 9 Jun 2020 at 04:42, Tom Callaway <tcallawa at redhat.com> wrote: > > Over the last several days, I've been working hard to get all of the Fedora > R packages rebuilt against R 4.0 in rawhide (in the F33-R-4 side tag). With > the exception of R-biomaRt, R-BSgenome, R-GenomicAlignments, and > R-rtracklayer, I believe everything is built and updated to the latest

On Tuesday, 9 June 2020 03.40.52 WEST Tom Callaway wrote: > Over the last several days, I've been working hard to get all of the Fedora > R packages rebuilt against R 4.0 in rawhide (in the F33-R-4 side tag). With > the exception of R-biomaRt, R-BSgenome, R-GenomicAlignments, and > R-rtracklayer, I believe everything is built and updated to the latest > versions. And of those

FYI, I plan on implementing this for F31 if no issues arise. ---------- Forwarded message --------- From: Ben Cotton <bcotton at redhat.com> Date: Tue, 2 Jul 2019 at 10:55 Subject: Fedora 31 System-Wide change proposal: Automatic R runtime dependencies To: <devel-announce at lists.fedoraproject.org>, Development discussions related to Fedora <devel at lists.fedoraproject.org>

R-BiocFileCache is now branched for f32 (finally). You should be able to build it if/when the PDC comes back up. Lotta random outages right now. Tom On Mon, Jul 6, 2020 at 10:40 AM Jos? Ab?lio Matos <jamatos at fc.up.pt> wrote: > On Tuesday, 9 June 2020 03.40.52 WEST Tom Callaway wrote: > > Over the last several days, I've been working hard to get all of the > Fedora >

Hello, I would like to plot specific SNPs with their exact locations on a chromosome. Based on my genotyping results I would like to separate these SNPs in three different categories: 1, 2 and 3 and use different colours to represent these categories. The script below generates the sample data. I can plot these with the image function using the following: val <- 1:3 samp <- sample(val,

Need way to resize an existing graphics window. This should be applicable across platforms (as part of a package). Context: function1() draws main plot (I'm using grid), function2() adds smaller plot above main plot, but this one can sometimes overflow the original graphics window area. Thanks, Sigal

Hi, I want to analyse DNA sequence data (mtDNA) in R as in calculate Fst, Heterozygosity and such summary statistics. Package Adagenet converts the DNA sequence into retaining only retaining the polymorphic sites and then calcuates Fst.. but is there any other way to do this? I mean analyse the DNA sequence as it is.. and calculate the statistics? Thanks! [[alternative HTML version deleted]]

Hi All, I have an CBS segmentation algorithm output for 10 tumor samples each from 2 different tumors. Now, I am in an urgent need to assign gene (followed by all genes present) that belong to a particular segment after I removed all the CNVs from segment data. The format of the data is: Sample Chromosome Start End Num_Probes Segment_Mean Sample1A-TA 1 51598 76187 15