similar to: vectorization with subset?

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 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

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'm working with a genomic data-set with ~31k end-points and have performed an F-test across 5 groups for each end-point. The QA measurments on the individual micro-arrays all look good. One of the first things I do in my work-flow is take a look at the p-valued distribution. it is my understanding that, if the findings are due to chance alone, the p-value distribution should be uniform. In

Hi, i have a fairly large amount of genomic data. I have created a dataframe which has "Reference" as one column and "Variation" as another. I want to plot a ROC curve based on these 2 columns. I have serached the R manual but I could not understand. Can anybody help me with the R code for plotting ROC curve. Thnx ashu6886 -- View this message in context:

*Postdoctoral Researcher - Bioinformatics/Statistics University of Oxford, UK *An exciting opportunity has arisen for a Postdoctoral Researcher in Bioinformatics/Statistics to join the Department of Oncology at the University of Oxford. The postholder will work under the supervision of Dr Francesca Buffa and Prof Adrian Harris, and will work closely with staff in the Molecular Oncology

The following citEntry includes a url with %3A and other encodings citEntry(entry="article", title = "Software for Computing and Annotating Genomic Ranges", author = personList( as.person("Michael Lawrence" )), year = 2013, journal = "{PLoS} Computational Biology", volume = "9", issue =

Hi, I am very new to R and wanted to know if there is a package that, given very long nucleotide sequences, searches and identifies short (7-10nt) motifs.. I would like to look for enrichment of certain motifs in genomic sequences. I tried using MEME (not an R package, I know), but the online version only allows sequences up to MAX 60000 nucleotides, and that's too short for my needs..

Hi all, is there a function (or variable) I can use in a custom init script that identifies the init script name? i.e. I'm porting some init scripts from gentoo, where the $SVCNAME variable identifies the init script name within the script itself... d /* Davide Cittaro Cogentech - Consortium for Genomic Technologies via adamello, 16 20139 Milano Italy tel.: +39(02)574303007 e-mail:

Hi all, I have two data frames, two columns each, 1000s of rows. Each row represents a segment of the genome where a deletion has occurred. First column is start position of the deletion in genomic distance, second is end position. So, e.g., first 3 rows of data frame A is: 1003 1023 5932 6120 12348 12689 first 3 rows of data frame B is: 852 5305 1010 1015 8500 9500 10000 13000 I want to merge

Dear all, I am still fairly new to R and try to analyze large tables of data generated from genomic experiment. Currently, I am trying to plot pair of experiments coming from different file, trying to look at the behavior of individual feature in pair of experiment. My problem is that I have independent list from different source and I would like to plot the pair of value using a common key. As

Hi all, I have a matrix with 10000 rows and 10 columns. The last columns contains another identifiers but the values are not uniques so that I want to generate another matrix with rows with unique values in the last column. If I did tmp<-unique(my_mat$col10) this will give me 8560 unique entries so the ideal matrix will be 8560X10 columns now then. I tried sub_mat<-my_mat[tmp,] but

Hi, I don't know much about the RGL package, and I have read the documentation and tried some parameters, with no luck... I would like to generate a movie from a 3D object (code below), where the vortex A is closer to the observer, and then the object rotates and the B vortex gets closer. I would like to capture this movie to a file. By the way, I am not being able to insert unicode text

I have run into a problem loading a just saved R object using R-devel. I have been saving and loading this particular type of R object for a long while and never ran into this problem. I save, then immediately reload (to test save) and get "ReadItem: unnknown type 65". This error is reproducible after logout from server and restart of emacs and R. Below is my output and

I have run into a problem loading a just saved R object using R-devel. I have been saving and loading this particular type of R object for a long while and never ran into this problem. I save, then immediately reload (to test save) and get "ReadItem: unnknown type 65". This error is reproducible after logout from server and restart of emacs and R. Below is my output and

Hi list. I would like to know how do I insert a North arrow, stored as a png file in my computer, in a map? I found lots of post asking similar things, one of them mentioned the pixmap package. The map was done using map() and shapefiles (the code is below). I’m using the pixmap () and addlogo() functions. Well I can import the png with pixmap() function (I guess, once there’s no error message),

Hi all, I'm new to Centos, just moved here from Gentoo Linux. I have to install a server for bioinformatics purposes and I see that default yum repositories do not include any bioinformatics software (i.e. ncbi-toolkit, blat, and others). I'm googling a bit but I can't find a valuable solution: which is (or which are) the best repository I should add to have a satisfying list

Hello, in a three dimensional coordinate system, I'd like to find all my experimental data points that fall within an ellipsoid around a fixed coordinate. The fixed point is defined by (x.coord.point, y.coord.point, z.coord.point). The coordinates of the ellipsoid are given by the three vectors x,y,z. In a previous version of my code, I simply used a box instead of an ellipsoid to sort

Hello, I'm wondering how I can merge two featuresets into one. My dataset is two sets of microarray data and it looks like followings: > rawData $v1 TilingFeatureSet (storageMode: lockedEnvironment) assayData: 2197815 features, 59 samples element names: channel1, channel2 protocolData rowNames: LT290677RU_D1_2011-02-16 LT286300LU_D1_2010-07-24 ... LT003990RU_D1_2010-11-04 (59

Hi all, I have 6 datasets(dataframes Assem_ContigsLen7 through all_ContigsLen12) containing 3 columns (contig_id, contig_length, read_count). Each dataset is composed of 3 types of contigs (assemblies of genomic fragments), 1- all Bacterial fragments, 2 - all Viral fragments, 3 - mixed fragments. I identified the type of contig through a merge with another table with just contig_id and