search for: chroms

Greetings, I am new to R, but trying to put in the time to learn. I have read the R manual and several other introductory texts; however, there is nothing like actually putting it into practice. So here is my problem, and its more of a learning exercise for myself than anything else, but I'm stuck and getting extremely frustrated that I can't figure it out. I'm trying to make a

Hi Jim I am trying to prepare a bed file to load as accustom track on the UCSC genome browser. I have a data frame that looks like the one below. > x V1 V2 V3 1 chr1 11255 55 2 chr1 11320 29 3 chr1 11400 45 4 chr2 21680 35 5 chr2 21750 84 6 chr2 21820 29 7 chr2 31890 46 8 chr3 32100 29 9 chr3 52380 29 10 chr3 66450 46 I would like to insert the following 4 lines at the beginning:

...[i-1,6]+1) x<-as.numeric(all.tf7[i,2]) } else if (all.tf7[i,1]!=all.tf7[i-1,1]) { all.tf7[i,6]<-(all.tf7[i-1,6]+1) x<-as.numeric(all.tf7[i,2]) } } #the aim here is to attribute a bin number to each row so that I can then split the dataframe according to those bins. chrom chromStart chromEnd name cumsum bin chr1 10089 10309 ZBTB33 10089 1 chr1 10132 10536 TAF7_(SQ-8) 20221 1 chr1 10133 10362 Pol2-4H8 30354 1 chr1 10148 10418 MafF_(M8194) 40502 1 chr1...

Hello, I am a total beginner with ?R? and found a package ?venn? to create a venn diagram. The problem is, I cannot create the vectors required for the diagram. The manual say: "R> venn(accession, libname, main = "All samples") where accession was a vector containing the codes identifying the RNA sequences, and libname was a vector containing the codes identifying the

.... I'm able to match peaks when they are chromatographed with the same method, but not when there are different methods are used and spectra comes in to play. While searching I found the ALS package which should be usefull for my application, but I couldn't figure it out. I made some dummy chroms with R, which mimic my actual datasets, to play with, but after looking at the manuals of ALS, I'm affraid I can't get the job done. Can someone put me on the right way? Here is my code to generate the dummy chroms, which also plots the 2 chroms and the spectra of the 3 peaks: #2D chromat...

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

Dear list, I wish to plot chromatin precipitation data: I would like to have a rectangles (x:end-start, y:peak) but I do not have an idea how to define x (in terms of qplot syntax) and to choose the correct geom. mydata is a subset of a larger file. > mydata chrom start end peak 1 chr11 5291000 5291926 8 2 chr11 10988025 10988526 7 3 chr11 11767950 11768676 8 4

On Fri, 16 Jan 2009, r-help-request at r-project.org wrote: > Date: Thu, 15 Jan 2009 13:29:03 +0100 > From: Pablo G Goicoechea <pgoikoetxea at neiker.net> > Subject: Re: [R] How to create a chromosome location map by locus ID > To: Sake <tlep.nav.ekas at hccnet.nl> > Cc: r-help at r-project.org > Message-ID: <496F2C0F.3040304 at neiker.net> > Content-Type:

Dear Oscar,   Thanks for your help.It's so nice of you to explain this package to me.   Best Regards,   James LAN 发件人： Oscar Rueda [via R] <ml-node+s789695n4431468h14@n4.nabble.com> 收件人： monkeylan <lanjinchi@yahoo.com.cn> 发送日期： 2012年2月29日, 星期三, 下午 9:21 主题: Re: Bayesian Hidden Markov Models Dear James, The distances are normalized between zero and 1, so in your case all of

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, i have two files (file1.txt and file2.txt) which i would like to merge, based on certain criteria, i.e. it combines data based on matching geneID and exons. i have used the merge option, but it does not give me the desired outcome. merged.txt shows the result i would like. *File1. txt* ** AffyProbe ProbeType Flag GeneSymbol GeneID Exons Chrom Strand Affytart AffyEnd 1

Dear R users, I need to merge two data frames based on both equal and unequal comparisons. The "sqldf" package used to work well , but today, I cannot resolve the following error by reinstallation of the sqldf package. Can anyone suggest a different way to perform this kind of merge function? Thank you, Ding > DMRlog2pbde47DMS <- sqldf("select * from DMR_log2pbde47 as a

Hello all, I'm trying to calculate the standar desviation and I'm using the function sd(x,na.rm=TRUE) and I have this error:? Error in var(x, na.rm = na.rm) : no complete element pairs . Why happen this?, What can I do to solve it?. x is list of three numbers which I have from a table. Thanks so much from Spain Carlos Morales Diego

I have a function that does some plotting. I then add lines to the plot. If executed one line at a time, there is not a problem. If I execute the function, though, I get: Error in ans[[1]] : subscript out of bounds This always occurs after the second lines command, and doesn't happen with all of my data points (some do not have errors). Any ideas? Thanks, Sean

Hi, i have two files (file1.txt and file2.txt) which i would like to merge, based on certain criteria, i.e. it combines data based on matching geneID and exons. i have used the merge option, but it does not give me the desired outcome. merged.txt shows the result i would like. *File1. txt* ** AffyProbe ProbeType Flag GeneSymbol GeneID Exons Chrom Strand Affytart AffyEnd 1

I have a data frame (narrow) with 431 rows and 6 columns containing information on chromosome, position, lod1, lod2, lod3, lod4, looking like this: > narrow chr pos lod1 lod2 lod3 lod4 1 1 3.456 -0.025 -0.003 -0.209 -0.057 2 1 5.697 -0.029 -0.005 -0.200 -0.058 3 1 8.434 -0.049 -0.012 -0.247 -0.092 4 1 9.466 -0.074 -0.025 -0.300 -0.136 5 1 9.706

Hello R Developers I am using DNACopy package http://bioconductor.org/packages/1.9/bioc/html/DNAcopy.html I am not able to figure out how to (if at all possible) to modify output format, namely, I am getting the following: "ID" "chrom" "loc.start" "loc.end" "num.mark" "seg.mean" Is it a way to get also median and standard deviation?

Hello! First time posting here: Here is my code: x <- c(1:22) finaloutput=cidrm=NULL finaldiversityoutput=diversitym=NULL diversityinfo=read.table("Diversity_info.txt", header=T, sep="\t", row.names=NULL) attach(diversityinfo) diversitynr=nrow(diversityinfo) diversitytemp <- matrix(0,nrow=diversitynr,ncol=1) for(j in 1:length(x)) {

Hi, I have a run of 5 graphs that I want to place them under the same page. Everything works fine to place them in a pdf file , or eps file, but when it comes to have a high quality of 300 dpi these graphs are not good. For example I open the eps file with Adobe Illustrator (AI) and it shows that it is a 72dpi graph. If I start with a 72dpi graph AI cannot improve this to 300 dpi. Q: HOW CAN A

New user here. My goal is pull daily averages from a long dataset. I've been working with some code I got from this list from https://stat.ethz.ch/pipermail/r-help/2009-March/191302.html The code how I have been using it is as follows: library(zoo) library(chron) DB<-read.table("/Users/me/Desktop/R/data.csv", sep=",", header=TRUE, as.is =TRUE) z<-zoo(LTER6$temp,