similar to: axis name editing

Dear R guru: If I got a variable aaa<- "up.6.11(16)" how can I extract 16 out of the bracket? I could use substr, e.g. substr(aaa, start=1, stop=2) [1] "up" But it needs start and stop, what if my start or stop is not fixed, I just want the number inside the bracket, how can I achieve this? Many thanks yan

Dear R user How to make the program stop at the spot where the error occurs? I mean inside the iterations,. Many thanks yan ********************************************************************** This email and any files transmitted with it are confide...{{dropped:10}}

Dear all R users, Suppose, Dear all R users, Suppose, pauto.cor = pacf(lh, plot=F) max.lag = max(pauto.cor$lag) min.lag = min(pauto.cor$lag) centre = (max.lag - min.lag)/2 pauto.cor = pauto.cor$acf pauto.cor = pauto.cor[-1] par(mar=c(3,0,1,1)) barplot(pauto.cor, axes=F,xlim=c(max(pauto.cor), min(pauto.cor)), space=0, col="green4",border="green",horiz=T) #This plots

Dear brainy R users, I need to output a matrix, with two colors , meaning some elements using different color I normally use write.table(table.m, file="table file name.csv", sep=","), how could I pass the index for the color ? Many thanks yan ********************************************************************** This email and any files transmitted with it are

Dear R gurus, If I got a vector with string characters like "abcd_efgh_12ab3_dfsfd", how could I extract "12ab3", which is the characters after second underscore and before the third underscore? Tons of thanks yan ********************************************************************** This email and any files transmitted with it are confide...{{dropped:10}}

Dear R users I need to solve the finding all cliques in a graph problem, is there a R package implementing Bron-Kerbosch algorithm? Many thanks YAn ********************************************************************** This email and any files transmitted with it are confide...{{dropped:10}}

Hi I am trying to save some plots in a postscript file. When I generate the plots in the main window, they appear correctly - their orientation is landscape (i.e., horizontal). However, when I open the .ps file with GSview, the whole page appears vertically, and the plot appears horizontally, which means that the plot is only partially visible (example here

Hi, I am using the function heatmap(stats) to draw a microarray heatmap, columns are samples and rows are gene features. I did a 2D clustering during the heatmap drawing. The features and samples indeed cluster into several blocks both vertically and horizontally. I can get the index of re-ordered rows and columns after the heatmap drawing by typing the the return variable of the heatmap

Dear R users, I need to add 0 in front of a series of numbers, e.g. 1->001, 19->019, Is there a fast way of doing that? Many thanks yan [[alternative HTML version deleted]]

I am using hclust and plot to produce dendrograms. Using my input data I am able to complete an analysis and obtain a vertical plot. I want to be able to plot the dendrogram horizontally.I am using version 2.6 of R and have updated my packages recently. Using the sample script for dendrograms I can produce a horizontal plot using the instruction horiz = TRUE in plot(). When I use the same

Dear All I'm trying to sort a matrix using function order, Some thing really odd: e.g. abc<-cbind(c(1,6,2),c(2,5,3),c(3,2,1))## matrix I want to sort if I do abc[ order(abc[,3]), increasing = TRUE] the result is correct [,1] [,2] [,3] [1,] 2 3 1 [2,] 6 5 2 [3,] 1 2 3 But if I want to sort in decresing order: abc[ order(abc[,3]), decreasing = TRUE] the

Hi, I'm new to R, and trying to map stuff I knew how to do in Matlab to R. For instances, in matlab, sum(x) gave the column sums of matrix(/vector) x. In R it just gives the sum of all the elements, regardless of their form (vector or matrix). Is there a simple way to return a vector of column sums for a matrix? Also, in Matlab binary operators exist for concatenating a matrix onto the end or

Dear all, I would like to draw a dendrogram and mark some parts/branches (by using "segments") including their labels. If I draw it without specifying the length of x axix, I am able to do that (as in My dendrogram 1 of the following codes). However, if I want to specify the x axix, I am not able to draw marking line (by using "segments") including labels (as in My dendrogram

Hi, I've just upgraded to 1.9.0 and one of my Sweave files that produces a number of barplots in a standard manner now produces them in a different way. I have made a couple of small changes to my code to get the back the output I was getting before upgrading and now (mostly out of curiosity) would like to understand what has changed. I *think* I've tracked it down to tapply() and/or

Dear all, I'm just wondering how to submit a package to cran? I followed the instruction, using anonymous as username and my email address as password, but it didn't connect. Any hints? Yan [[alternative HTML version deleted]]

Dear All, I'm wondering if there is a R function could give me all the combinations of the grouping/cluster result, given the number of the groups. e.g. 3 objects: x1 x2 x3, number of groups is 2 so the result will be group1:x1,x2; group2: x3 group1: x1;group2: x2,x3 group1: x1,x3;group2: x2 many thanks yan

Hi, I have a question about the heatmap dendrogram in R. I loaded my data matrix in command heatmap(), and it gave me a heatmap accordingly, and the rows (genes) also were clustered accordingly. But now I don't want the genes clustered in that way, I have a new order of these genes, and want the rows of the heatmap are drawed as the order I give. Could anybody help me on that? Thank you!

Hi, Using library(lattice), is there any way to tell xyplot to plot panels top to bottom, then left to right (i.e. panels are appended vertically, then horizontally). as.table changes the plot direction from left-to-right then top-to-bottom, to right-to-left then bottom- to-top, but that's not quite what I want to do. Thanks Yan

Dear R users, I'm trying to add a regression line on my boxplots (something like:boxplot(c(1:3),c(4:6),c(5:8))) But I can't see it. Please help !!! It's not a April fool's joke!!! Yan [[alternative HTML version deleted]]

Hi, I'm building a graph (barplot) in which the X axis label disappears. I tried to use the option mgp of par() and I could not get the desired result. Note that want the axis labels horizontally. caes = c(37,20,19,16,75,103) names(caes) = c("Pinscher", "Pastor \n Alem?o", "Poodle", "Rottweiller", "SRD", "Outros") caess =