similar to: Error: cannot take a sample larger than the population

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, (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, 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

Hi, I have a large group of jobs, some of them are running on the same directory. All of them in batch mode. What are the best ways to protect from corrupting the results two or more jobs running on the same directory. One, I would think can be to run each job in a separate directory, collect the results and after remove the directories. But I have thousands of jobs that will run in parallel

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.

> 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, 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, 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 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 all, I tried to follow an online tutorial to run openBUgs but the package BRugs has been removed from R repository. Could someone provide a link for where to download BRugs? Thanks a lot! -- View this message in context: http://www.nabble.com/Where-to-download-BRugs-tp17411410p17411410.html Sent from the R help mailing list archive at Nabble.com.

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, 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, 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

WinXP, Splus7 and R2.3.1. All, I have been using the SSfpl and SSlogis self-starting functions in the nlme library to fit 4PL and restricted 4PL models. I need to adapt these routines to fit the alternative model f(x) = A + (B-A)/(1 + abs(x/EC50)^C) My Question: How do I obtain good starting values for this alternative model? (The pseudo-code found on pages 517 - 520 of "Mixed

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