similar to: setting par(srt) according to plot aspect ratio

I'm able to do this as follows, but am wondering if anyone knows a simpler way which still avoids explicit loops? > (mystring <- letters[1:5]) [1] "a" "b" "c" "d" "e" > unlist(sapply(mystring[-length(mystring)], + function(x) paste(x,mystring[(grep(x,mystring)+1):length(mystring)],sep=""))) a1 a2 a3

Hi STATISTICA has a function called "Two-way joining" (see http://www.statsoft.com/TEXTBOOK/stcluan.html#twotwo) and the reference material states that this is based on the method as published by Hartigan (found this paper: http://www.jstor.org/pss/2284710 through wikipedia). What is the relationship (if any) between the "heatmap" function in R and this technique? Is there an

Dear R users, I have various vector geometry operations to perform on 3-D coordinate data located on multiple (500+) csv files. The code I have written for the calculations works just fine. I have written a 'for' loop to automate the task of extracting the coordinates from the files and perform the analyses. The loop works reasonable well, but if the number of csv files is greater than

For example, subtracting 1:2 from the rows of a two-column matrix: > t(apply(matrix(1:6,ncol=2),MARGIN=1,function(y) y - 1:2)) [,1] [,2] [1,] 0 2 [2,] 1 3 [3,] 2 4 > sweep(matrix(1:6,ncol=2),MARGIN=2,1:2,FUN="-") [,1] [,2] [1,] 0 2 [2,] 1 3 [3,] 2 4 Is there a logic to this difference, or is it just a quirk of the history of these

The order() function allows you to specify multiple vectors, which are used successively to break ties. If I want to use many vectors to break ties (say, 25 or more), that are columns of a matrix or elements of a list, does anyone know a shortcut to do this without passing 25 arguments to order()? -- Levi Waldron post-doctoral fellow Jurisica Lab, Ontario Cancer Institute Division of Signaling

Can someone recommend a more sophisticated way to annotate heatmaps than the ColSideColors argument of heatmap and heatmap.2? In particular, I would like to be able to annotate columns with more than one piece of information, like in Figure 1 of the article at http://www3.interscience.wiley.com/cgi-bin/fulltext/117905619/HTMLSTART / doi:10.1002/hep.22256. Some example data and a heatmap:

Hi, I'm trying to figure out how to safely make sure that I close the same graphics device that I opened earlier in a thread (and not one opened by a parallel thread). In a *single-thread* environment, one can do the following to open and close a device: makePlot <- function(i) { filename <- sprintf("foo%d.png", i); png(filename); idx <- dev.cur();

Dear list, I'm using the brew package to generate a report containing various plots. I wrote a function that creates a plot in png and pdf formats, and outputs a suitable text string to insert the file in the final document using the asciidoc syntax, <% tmp <- 1 makePlot = function(p, name=paste("tmp",tmp,sep=""), width=300) {

Hi, I'm dealing with a lattice plot inserted into a tk widget and would like to know when a user has clicked on the plot area of a plot (i.e. inside the axes). For example, library(tkrplot) library(lattice) tt <- tktoplevel() makePlot <- function() print(xyplot(1 ~ 1)) printCoords <- function(x, y) print(c(x, y)) img <- tkrplot(tt, makePlot) tkbind(img, "<1>",

Hi, I have a question about filterMicroRna in AgiMicroRna package function for filtering probes in Agilent microRNA dataset. >ddPROC = filterMicroRna(ddNORM.micro, dd, control = TRUE, IsGeneDetected = TRUE, wellaboveNEG = FALSE, limIsGeneDetected = 75, limNEG = 25, makePLOT = TRUE, target.micro, verbose = TRUE) If in a dataset there are two or more sample replicates and in the step of

Hello, I am running into this error: /Maximum number of connections from user+IP exceeded (mail_max_userip_connections=10)/ The suggested solution in hundreds of support requests on this mailing list and throughout the internet is to increase the number of maximum userip connections. But this is not curing the problem, it is just postponing it to the moment when the new limit is reached. When i

I'm trying to help with http://bugs.codeweavers.com/show_bug.cgi?id=434 and came across my thread from Feb, that I must have missed. >> ----------------------------------- >> Non-Working (wine-2002-02-20 cvs) >> trace:win:SetCapture Calling Event_Capture >> trace:win:EVENT_Capture (0x20021) Window #!! -Maybe? >> trace:win:EVENT_Capture About to Return from Event

In the following, can anyone tell me why length(eee) returns 9? I was expecting 15398, and when I try to add this vector to a data frame with that many rows, it fails complaining that the vector is of length 9. In what I thought was an identical situation with a related dataset, the same code worked as expected. > length(fff) [1] 15398 > str(fff) int [1:15398] 20010102 20010102

Dear Friends, Yesterday I posed a question to the list concerning the possibility of doing animation examples in R. Here is an example S-Plus (4.5 for Windows ) that I wrote to illustrate my problem. If I try this in R (comment out the graphsheet and guilocator calls) I don't see my results until after the function has iterated through the 100 trials. I would like to plot each iteration.

Hi everyone, I am desperately looking for a method to interpolate strings within an SQL statement as follows: I get a lot of rows out of a database (in my example POSITION_to_ZIPCODE Database with holds records for German ZIP Code <--> Gauss-Krueger Coordinate System ) and want this to be selected and computed individually row by row as follows: library(RODBC) channel <-

can i have some instruction on how to draw this type of plot in R : http://sitmo.com/doc/Image:Ll_ou_conpron.png#filelinks ? Regards, Send instant messages to your online friends http://uk.messenger.yahoo.com [[alternative HTML version deleted]]

Dear All, I am not very into object-oriented programming, but I would like to learn the ropes for some R applications. Quoting from the online R language definition (paragraph 5.1) > Consider the following simple example. A point in two-dimensional > Euclidean space can be specified by its Cartesian (x-y) or polar > (r-theta) coordinates. Hence, to store information about the location

Full_Name: Naoki Takebayashi Version: 0.65.0+R-release.diff (Oct 6, 1999) OS: Linux/Alpha Submission from: (NULL) (129.79.224.171) This will fix the "problem 2 (crash in fft)" in Bug ID #277 On Linux/Alpha, make check failed because R could not handle the following example in base-Ex.R ##___ Examples ___: # The Old Faithful geyser data data(faithful) : : ## Missing values: x <-

We are looking at adding to our storage. We would like ~20TB-30 TB. we have ~ 200 nodes (1100 cores) to feed data to using nfs, and we are looking for high reliability, good performance (up to at least 350 MBytes /second over 10 GigE connection) and large capacity. For the X45xx (aka thumper) capacity and performanance seem to be there (we have 3 now) However, for system upgrades , maintenance

I've wrote some code to simulate a random walk in 2 dimensions on a lattice. Basically I want to add something in to make it plot it point by point so you can see what is going on. Heres my code for the random walk in 2d RW2D<-function(N) { i<-0 xdir<-0 ydir<-0 xpos<-vector() xpos[1]<-xdir ypos<-vector() ypos[1]<-ydir for (i in 1:N-1)