similar to: How to protect two jobs running on the same directory at the same time not to corrupt each other results:

Hi, Is there a smart way in the R graphs to create a line that is broken in intervals based on the indicator given below. following is a small test graph Location,indicator,otherinfo 1.2,1,2.2 2.5,1,2.5 3.7,1,2.3 20.1,2,4.3 22.5,2,5.2 25.0,2,3.4 27.3,2,2.2 35.1,3,3.4 37.0,3,7.2 38.0,3,6.1 40.1,3,5.4 52.9,3,3.3 Right now in the plot the line is continuous, but I would like to have it broken

Hi, If I use x<-1:10 y<-rnorm(10,0,1) ### pdf(file="c:\\aldi\\test.pdf") plot(x,y) segments(x,y,x+2,y+2) segments(x,y,x+0.5,y+0.5,col=3) ### dev.off() ### q() Is there a way that I can imbed in the graph plot for each point defined by x to x+0.5 and y to y+0.5 (and colored in green) a different hyperlink? For example point one (together with the green tail) will have the

Hi, I have a data.frame with many variables for which I am performing the mean by subgroup, for a pair of variables at a time, where one of them for each pair defines the subgroup. The subgroups in the x$cm1 are 0, 1 and 2. x ph1 cm1 0.2345 2 1.2222 1 2.0033 0 0.0000 2 1.0033 1 0.2345 0 1.2222 2 2.0033 0 0.0000 1 1.0033 2 > meanbygroup <- as.vector(with(x, by(x$ph1, x$cm1, mean)))

Hi, I am trying to install the genetics package on a server with Linux, Fedora. I installed it in a PC and worked fine. In a server since I am not used with R, I am not sure what do I need to change so genetics pkg can see some package dependencies: Any suggestion is appreciated, Aldi Note: genetics expects gregmisc and mvtnorm to be installed already. gregmisc creates gdata etc dependencies.

Hi, In Splus7 this statement xlrmN1 <- sample(c(0,1,2),400 ,prob=c(0.02 ,0.93 ,0.05 )) worked fine, but in R the interpreter reports that the length of the vector to chose c(0,1,2) is shorter than the size of many times I want to be selected from the vector c(0,1,2). Any good reason? See below the error. > xlrmN1 <- sample(c(0,1,2),400 ,prob=c(0.02 ,0.93 ,0.05 )) Error in

Hi, Using R 2.15.1 on Windows 7. Have installed both versions 32 and 64bit. In both of them among others I have installed a package rgenoud When I open R gui of 32bit and write library(rgenoud) it responds by showing a functional rgenoud version 5.7-8. The same it does on Rgui 64bit. Now I am working in cygwin (v. 1.12.4.0) with xwin. Normally before when I had installed a package, I only had

> Date: Mon, 20 Feb 2006 16:43:55 -0600 > From: Aldi Kraja <aldi at wustl.edu> > > Hi, > Using package gclus in R, I have created some graphs that show the > trends within subgroups of data and correlations among 9 variables (v1-v9). > Being interested for more details on these data I have produced also the > var-covar matrices. > Question: From a pair of two

Hi, Using package gclus in R, I have created some graphs that show the trends within subgroups of data and correlations among 9 variables (v1-v9). Being interested for more details on these data I have produced also the var-covar matrices. Question: From a pair of two subsets of data (with 9 variables each, I have two var-covar matrices for each subgroup, that differ for a treatment on one

Hi, Based on par function, I can split the screen into two parts left and right. I wish x occupies the half left screen, and all plants occupy half right screen, which happens right now. But I wish the right screen, to be split in 3 or more vertical parts where each pair of the same type of plant, are together in its own block of boxplot, because each plant has its own unit of measure.

Hi, I am working with some extremely small p-values and I want to capture the corresponding quantiles. I see the help file it says: 'qnorm' is based on Wichura's algorithm AS 241 which provides precise results up to about 16 digits. What happen after the 16th digits? If I am running R in a server 64-bit, can that improve the chances that beyond 16th digits to still have

Hi, I have a data.frame of this kind: x obs movie earned rating 1 P1 3.2 xx 2 P1 4.2 xx 3 P1 5.2 xx 4 P1 6.2 xx 5 P2 3.5 xx 6 P2 6.5 xx 7 P2 7.5 xx 8 P2 4.5 xx 9 P2 4.5 xx 10 P3 4.8 x1 11 P4 7.3 x2 12 P4 3.2 x2 13 P4 3.3 x2 I want to overlay the

Hi, I am running a large number of jobs (thousands) in parallel (linux OS 64bit), R version 2.14.1 (2011-12-22), Platform: x86_64-redhat-linux-gnu (64-bit). Up to yesterday everything ran fine with jobs in several blocks (block1, block2 etc) of submission. They are sent to an LSF platform to handle the parallel submission. Today I see that only one of the blocks (the 19) has not finished

Hi Test made in: R in windows Vista OS, R version 2.8.1 From FAQ: http://cran.r-project.org/doc/FAQ/R-FAQ.html#How-do-I-convert-factors-to-numeric_003f "It may happen that when reading numeric data into R (usually, when reading in a file), they come in as factors. If |f| is such a factor object, you can use as.numeric(as.character(f)) to get the numbers back." 1: Why it may happen?

Hi, (R version 2.15.0) I am running a pgm with 1 response (earlier standardized Y) and 44 independent vars (Xi) from the same data =a2: When I run the 'lm' function on single Xi at a time, the beta coefficient for let's say X1 is = -0.08 (se=0.03256) But when I run the same Y with 44 Xi-s with the 'step' function (because I left direction parameter empty, I assume a backward

Hi, I was wondering if there is a way to figure out why in SAS random beta coefficients are 0 vs. in R the beta-s are non zero. The variables of the data are nidl, time, and sub (for subject). Time and nidl are continuous variables. I am applying random coefficients model. Any input is greatly appreciated, Thanks, Aldi 1. mixed model in SAS: ====================== ods output SolutionR =

Hi, I have a vector Y1 <-c(8, 11, 7, 5, 6, 3, 6, 3, 3) and an index iy <-c(c(1, 2),c(1 2), c(1, 2, 3, 4), c(2, 3, 5), c(4), c(5, 6, 7), c(7, 8, 9)) how can I produce the mean, or the sum of the elements specified in the index iy from the vector Y1? expecting something like this for the sum: Y2 19 19 31 24 5 15 12 I thought lapply function may perform this, but does not work:

Hi, If I have an R object UUU, where the second element is U2, based on "g" column of my.table my.table of UUU is: mmm ggg gindex map Info aaa123 U1 1 1 1 aaa124 U1 1 2 1 bbb1378 U2 2 1 1 bbb8888 U2 2 2 0 bbb1389 U2 2 3

Dear list member, I think that I have detected a strange behavior of the save() command: > year <- "2000" > assign(paste0("Var_", year), list(A=10, B=20)) > get(paste0("Var_", year)) $A [1] 10 $B [1] 20 # At this point all is ok, I have created a list of name Var_2000 > save(paste0("Var_", year), file=paste0("Var_", year,

Hola, ¿Como estas?, la siguiente tal vez sea una pregunta simple. Quiero crear nombres de variables pero que el nombre lleve el valor de un objeto. Ej: tengo estos dos valores x1<-10 x2<-20 Quiero crear las variables: var_10 <- var_20 <- Pero usando los objetos x1 y x2 como referencia. Solo como referencia, en stata esto lo hacia de la siguiente forma local x1 = 10 local x2 =

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