similar to: Whole genome searching of 100bp "D" sequence

Hi, I was installing hundreds of packages on a machine with a single call to install.packages() and after a long time the call to install.packages() finally returned with the following warnings and errors: Warning messages: 1: packages ?hgu133aprobe?, ?hgu95av2.db?, ?BSgenome.Celegans.UCSC.ce2?, ?BSgenome.Mmusculus.UCSC.mm10?, ?BSgenome.Dmelanogaster.UCSC.dm3.masked?,

I do have a bunch of genes ( nearly ~50000) from the whole genome, which read in genomic ranges A range(gene) can be seem as an observation has three columns chromosome, start and end, like that seqnames start end width strand gene1 chr1 1 5 5 + gene2 chr1 10 15 6 + gene3 chr1 12 17 6 + gene4 chr1 20 25 6 + gene5

I do have a bunch of genes ( nearly ~50000) from the whole genome, which read in genomic ranges A range(gene) can be seem as an observation has three columns chromosome, start and end, like that seqnames start end width strand gene1 chr1 1 5 5 + gene2 chr1 10 15 6 + gene3 chr1 12 17 6 + gene4 chr1 20 25 6 + gene5

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,

Hi, I am trying to use 'GridR' package for the first time, and I'm running into a strange error from grid.check: > grid.check(gridFun) Error in exists(add) : invalid first argument After playing around in recover mode, I see that this because the variable 'add' created by grid.check is blank: Browse[1]> split[[1]][k] [1] " : " Browse[1]> add [1]

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'm trying to use coxph() function to fit a very simple Cox proportional hazards regression model (only one covariate) but the parameter space is restricted to an open set (0, 1). Can I still obtain a valid estimate by using coxph function in this scenario? If yes, how? Any suggestion would be greatly appreciated. Thanks!!! [[alternative HTML version deleted]]

Hello, I have an input file: http://r.789695.n4.nabble.com/file/n3700031/testOut.txt testOut.txt where col 1 is chromosome, column2 is start of region, column 3 is end of region, column 4 and 5 is base position, column 6 is total reads, column 7 is methylation data, and column 8 is the strand. I would like a summary output file such as:

Hi, is there any package/function in Bioconductor that can do this: if given the chromosome positions of a fragment, find out all genes within, and with the information about which strand is the sense strand. And vice versa. Thanks a lot. ----- Appreciate your time & attention! -- View this message in context: http://www.nabble.com/list-genes-w-n-a-genomic-fragment-tp18331452p18331452.html

I've got two tables.... first one(table1): ID chrom start end Ex1 2 152 180 Ex2 10 2000 2220 Ex3 15 3000 4000 second one ( table2): chrom location name 2 160 Alv 2 190 GNN 2 100

Great R people, I have noticed a strange behaviour in read.delim() and friends in the R 2.7.0 version. I will describe you the problem and also the solution I already found, just to be sure it is an expected behaviour and also to tell people, who may experience the same difficulty, a way to overcome it. And also to see if it is a proper behaviour or maybe a correction is needed. Here is the

I have a dataframe in the general format: chr1 0.5 chr1 0 chr1 0.75 chr2 0 chr2 0 chr3 1 chr3 1 chr3 0.5 chr7 0.75 chr9 1 chr9 1 chr22 0.5 chr22 0.5 where the first column is the chromosome location and the second column is some value. What I'd like to do is have a histogram created for each chr location (i.e. a separate histogram for chr1, chr2, chr3, chr7, chr9, and chr22). I am just

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

Genome Institute of Singapore (GIS) Biostatistics Group The Genome Institute of Singapore ( http://www.gis.a-star.edu.sg <http://www.gis.a-star.edu.sg> ) is looking for statisticians who are interested in a wide range of biomedical probems. We are especially interested in hearing from senior researchers interested in leading their own group. The institute is well-funded and engaged

We are seeking a talented and driven postdoctoral fellow to join the laboratory of Jake Michaelson, PhD, assistant professor in the department of psychiatry at the University of Iowa. The Michaelson lab investigates the effect of genetic variation on the development and function of the brain, particularly in the context of neurodevelopmental conditions such as autism and language impairment. Our

Hello, I am using R to look at whole-genome gene expression data. This means about 27,000 genes, each with a vector of numbers reflecting expression at different tissues and times. I need to do an all against all co-expression calculation (basically, just calculate Pearson's r for every gene-gene pair). I try to store the result of such a thing in a 27000x27000 matrix, but r seems not to like

Dear all, I am currently analyzing a set of arrays hybe on the lattest affy Drosophila 2.0 GeneChip. I am trying to run simple annaffy analysis but can¹t find what is the name of the annotation file I need to use. Here is the output of the AffyBatch object I am using: > expData AffyBatch object size of arrays=732x732 features (16749 kb) cdf=Drosophila_2 (18952 affyids) number of samples=4

Hi, I have a set of genes and its chromosomal physical position in a text file. I want to view those genes in the chromosome using R package GenePlotter. Could any one please tell how to view this. Thanks in advance. Yours sincerely, S.Mahalakshmi [[alternative HTML version deleted]]

I've been looking for how to change a certain percentage of values in a data frame, but I've been struggling to find information in R. For example: #################example data############## > data V1 V2 V3 V4 V5 V6 V7 1 chr1 500 500 CHH 0 0.5 + 2 chr1 550 550 CHH 0 0.0 + 3 chr2 700 700 CHH 0 0.0 + 4 chr2 1000 1000 CHH 0 0.0 + 5 chr3

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