similar to: 3D plot

HI, this is my problem I want to subset this file df, using only unique df$exon printing the line once even if df$exon appear several times: unique(df$exon) will show me the unique exons If I try to print only the unique exon lines with df[unique(df$exon),] -this doesn't print only the unique ones :( could you help? thanks Nat exon size chr start

Hello I am am new to R and any help with the following would be appreciated: I have a list (example attached) and I would like to create a new list which is a filtered version of this list. I.e I would like a list that only contains elements with this value: Chr 10 : 21853562 - 21855482 Any pointers/tips would be great. Thanks! Nikol -- Bioinformatician/Computer Associate Cambridge Institute

Dearl all, I have to fit a function y = f(x1, x2) to data experiemntal data describing the measured behavior of y. x1 and x2 are the independent variables. Could you suggest me wich R package can I use for this purpose? Thanks, Paola. -- *Paola Lecca, PhD* *The Microsoft Research - University of Trento* *Centre for Computational and Systems Biology* *Piazza Manci 17 38123 Povo/Trento,

Hi R experts. how can I average two sets of columns? dd <- data.frame(b = c("chr2", "chr1", "chrY", "chr13", "chrX"),         w=11:15, x = 1:5, y = c(8, 3, 9, 9,7),        z = c(1, 1, 1, 2, 8)) expected result for mean(w,x) and mean(y,z) is 1  chr2 6  4.5 2  chr1 7  2 3  chrY 8  5 4 chr13 9  5.5 5  chrX 10 7.5 Yu [[alternative

i am having problem running my own data. yesterday it was working just fine. today it is not. this is the code i was using as an example to follow. this code ALSO worked just fine yesterday, and is no longer working at all. i suspect it is a problem with either my computer or the software, at this point. if THIS won't even run.... something is wrong. i can assure you this isn't

Dearl all, I have setup samba2.2.7 with cups1.1.14. In win2k client, the browse the printer share and install the windows 2000 printer successfuly. When I try to print some to it, nothing will be printed out. I check the cups's error_log found the following message: "E [11/Dec/2002:20:35:43 +0800] print_job: Unsupported format 'application/octet-stream'!" Also, I check

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

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, 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: How do I select different variables from a list of data frames. I have a list of 13 that looks like below. Each data frame has more variables than I need. How do I go through the list and select the variables that I need. In the example below, I need to get the variables "a", and "q10" and "q14" to be returned to two separate data frames. Thank you. Yours, Simon

Hola a todos: Estoy intentando visualizar unos resultados de un programa en R de forma dinámica. Dichos resultados se pueden representar, básicamente, con un plano que gira sobre uno de sus ejes: para dar sensación de movimiento he utilizado la función rgl.clear(), ya que no encuentro los equivalentes en rgl de los parámetros add o new. El problema es que al limpiar la pantalla con rgl.clear la

Hello: Forgive me, this is surely a simple question but I can't figure it out, having consulted the help archives and "Data Manipulation With R" (Spector). I have a list of 11 data frames with one common variable in each (prov). I'd like to use lapply to go through and recode one particular level of that common variable. I can get the recode to work, but it only returns the

Ok so I am trying to add a subtitle to my xYplot using Hmisc and I have some dummy code of x<-seq(1,10,1) y<-rev(seq(1,10,1)) ci<-y*.10 ciupper<-y+ci cilower<-y-ci this code works fine: xYplot(Cbind(y,ciupper,cilower)~x, type=c("b"),plot.points = TRUE, method='filled bands', main='main title') but when I add sub=. at the end and use this

I have a multiple text files, separated by a single space, that I need to combine into a single data.frame. Below is the track I'm on using list.files to capture the names of the files and then lapply with read.table. But I run into problems making a usable dataframe out of the data. #Creating example data in similar format to data I have sub <- rep(1,10) trial <- seq(1,10,1) size

I have been using lattice xyplot and am quite pleased, and I can use the type=c("b","g") to have it print gridlines into the page, yet if I want to have a line plot with points on it, how do I get the xYplot to print gridlines (I use Hmisc xYplot because of its bands method which allows plotting of confidence intervals). Any suggestions? I have looked at the panel functions but

I'd like to use rgl.surface (or some other function if more appropriate) to create a horizontal and vertical transparent grey slice (plane) running through both the x and y origins and extending across the z axis, i.e. the 3-d equivalent of the normal 2-d coordinate axes we are all familiar with. The examples for rgl.surface are a bit more complex than what I need and I am having trouble

In the "filled bands" part of xYplot of the Hmisc package, is there a way to have multiple bands with multiple lines? or does it just allow one for now? So I had an example bit ago had a made up line and CI, now if I wanted to make a second line with a CI filled in can I put them on the same plot? x<-seq(1,10,1) y<-seq(1,10,1) ci<-y*.10 ciupper<-y+ci

Can anyone tell me how to remove several numbers for a sequence. For example: xx<- c(1,5,7,10) yy<-seq(1,10,1) how do I get take xx away from yy to get the new sequence 2,3,4,6,8,9 Many thanks in advance -- View this message in context: http://r.789695.n4.nabble.com/Remove-several-numbers-from-a-sequence-tp4640630.html Sent from the R help mailing list archive at Nabble.com.

mboost contains a blackboost method to build tree-based boosting models. I tried to write my own "cost-sensitive" ada family. But obviously my understanding to implement ngradient, loss, and offset functions is not right. I would greatly appreciate if anyone can help me out, or show me how to write a cost-sensitive family, thanks! Follows are some families I wrote ngradient <-

Hi! I have problems with labeling x axis while plotting time series data. I have 40 monthly measurement. One period lasts 4 months. I'd like to have 40 ticks on x axis (10 larger, the rest smaller) and labels just at the beginning of each period, just like in the image <http://r.789695.n4.nabble.com/file/n4660465/2221.jpg> My code leaves x axis empty: > data <-