similar to: Summary: A speed improvement challenge

Dear All, I am going through a worked example provided by Harrell, Lee and Mark (1996, Stats in Medicine, 15, 361-387). I know that the code provided is for S-PLUS and R but the languages don't differ enough for this to be a problem. I am using the Hmisc and Design libraries and have used the following code (as shown in the example provided in the referenced paper): '%in%' <-

Buenos días: Después de más de 20 años en SAS, ahora igual me toca trabajar en R. Y la transición es muy dura: la "lógica" de los dos lenguajes es, para mi, totalmente diferente. Disculpen si lo que pregunto es una obviedad, pero llevo todo el día con esto y no logro encontrar la respuesta adecuada con google. Estoy intentado crear una función que realice una descriptiva simple de

Hi all, In checking my R codes, I encountered the following problem. Is there a way to fix this? I tried to specify options(digits=). I did not fix the problem. Thanks so much for your help! Hanna > cdf(pmass)[2,2]==pcum[2,2][1] FALSE> cdf(pmass)[2,2][1] 0.9999758> pcum[2,2][1] 0.9999758 [[alternative HTML version deleted]]

Hi, When writing formulae in Rd files I use \leq to get a less-than-or-equal-to sign, as this is what the LaTeX references I consulted suggested. This works correctly for the PDF output (Rcmd Rd2dvi.sh --pdf), however when the Rd file is converted to text or html this command is converted to "<=q". A workaround is to use \le instead of \leq which works with my LaTeX distribution

Hello. I have another rect question. If I want to make a lot of rectangular, without borders, and without any spaces between them, how could I do it? I.e. > plot(-1:1., -1:1., type = "n") > rect(0,-1,1,1,col="red",lwd=0,border=F,xpd=T) > rect(-1,-1,0,1,col="blue",lwd=0,border=F,xpd=T) will produce two rectangular, but there is space between them. How

Hi, Check the FAQ 7.31 https://cran.rstudio.com/doc/FAQ/R-FAQ.html#Why-doesn_0027t-R-think-these-numbers-are-equal_003f And read the posting guide too... https://www.r-project.org/posting-guide.html HTH, Ivan -- Dr. Ivan Calandra TraCEr, Laboratory for Traceology and Controlled Experiments MONREPOS Archaeological Research Centre and Museum for Human Behavioural Evolution Schloss Monrepos 56567

Hello. I am porting Splus library to R, which use S "graphics.h" and "device.h", the library makes a number of high resolution plots using these. Does anyone have an experience with similar porting? I don't really want to rewrite all C code, so any recommendation would be helpful. Thanks. --------------------------- Vadim Kutsyy vadim at kutsyy.com

Another question "graphapp.h" in R1.1 defines number of plotting function. is there an example of how to use them? Thank. --------------------------- Vadim Kutsyy vadim at kutsyy.com http://www.kutsyy.com -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info",

Another question "graphapp.h" in R1.1 defines number of plotting function. is there an example of how to use them? Thank. --------------------------- Vadim Kutsyy vadim at kutsyy.com http://www.kutsyy.com -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info",

Is there a quick way to make boxplots groups by two variables? By that I mean, that if x axes have values ("A","B","C"), than at each value there would be a few boxplots each for a value of second variable (say ("1","2","3")). Thank You. -- Vadim Kutsyy http://www.kutsyy.com vadim at kutsyy.com The University of

Is there a quick way to make boxplots groups by two variables? By that I mean, that if x axes have values ("A","B","C"), than at each value there would be a few boxplots each for a value of second variable (say ("1","2","3")). Thank You. -- Vadim Kutsyy http://www.kutsyy.com vadim at kutsyy.com The University of

I am writing a function to work with class I am defining. I have a question about using nargs() inside of parentheses function. nargs() shows the same for supplying 1 argument, or no arguments at all. Here is a small example: > "[.myclass"<-function(x,...) print(nargs()-1) > x<-c(1,2,3) > class(x)<-"myclass" > x[] [1] 1 > x[1] [1] 1 > x[1,2] [1] 2

Dear Rich, I think that it's generally a bad idea to give statistical (as opposed to simply technical) advice by email without knowing the context of the research. I think that you'd do well to seek help from a statistician, and not just do what I suggest below. Interpolating the data only makes sense if there's no random component to the response (mag in your data). Otherwise, it

Hello, I am trying to create subsets of grouped data (by area size), and use the area size as part of the output name. The code below works for area (xout) 1 and 50, the other files are given NA for an area. A simple example: xout <- c(1,5,10,25,50,100) for(i in xout) { print(paste("Areal_Ppt_",xout[i],"sqmi.txt", sep="")) } [1] "Areal_Ppt_1sqmi.txt"

Hello, I'm fairly new to R and I can't work out how to produce a double inequality like (LaTeX) $0 \leq x \leq 1$ in the legend of a graph. If I try > legend(50, 0.1, legend = c(expression(0 <= x <= 1), c(2 <= x <= 3)), pch = c(1,1), col = c(2, 3)) then I get an error message "unexpected '<=' in ...". I've checked the help files for plotmath and

Here is a complete 104-line native code compiler for a tiny subset of OCaml that is expressive enough to compile an external Fibonacci program: type expr = | Int of int | Var of string | BinOp of [ `Add | `Sub | `Leq ] * expr * expr | If of expr * expr * expr | Apply of expr * expr type defn = | LetRec of string * string * expr open Camlp4.PreCast;; let expr = Gram.Entry.mk

Hello all, I'm struggling to grasp how I might use lapply() instead of looping to run approx() on a list consisting of multiple arrays - each of equal dimension. But simpler than that, I haven't been able to successfully apply approx() to an array, unless I loop through the third dimension and extract the matrix. See example code below... Any suggestions will be gratefully received. Thanks

Greetings, Sorry, the last message was sent by mistake! Here it is again: I encounter a strange problem computing some numerical integrals on [0,oo). Define $$ M_j(x)=exp(-jax) $$ where $a=0.08$. We want to compute the $L^2([0,\infty))$-inner products $$ A_{ij}:=(M_i,M_j)=\int_0^\infty M_i(x)M_j(x)dx $$ Analytically we have $$ A_{ij}=1/(a(i+j)). $$ In the code below we compute the matrix

Hello everybody, Data is myd <- data.frame(id1=rep(c("a","b","c"),each=3),id2=rep(1:3,3),val=rnorm(9)) I want to get a cumulative sum over each of id1. trying aggregate does not work myd$pcum <- aggregate(myd[,c("val")],list(orig=myd$id1),cumsum) Please suggest a solution. In real the dataframe is huge so looping with for and subsetting is not a

On Sat, Feb 14, 2009 at 00:17, <ashrafi@ucdavis.edu> wrote: Hi, I was trying to read ~400 chips in an affybatch and I got the same message. Could you find a remedy for that. My server has 128 GB of RAM. However, R halted ever before it uses the memory. I have been able to load upto 250 CEL files but this time I wanted to test what would happen if I want to normalize 400 chips. Thanks