similar to: "Pipelining" programs in R

Hello. I am creating an R package that I'd like to submit to CRAN (OS Windows XP). How do I distinguish among 'public' functions, e.g., those that are intended to be called by users of the package and for which I am providing documentation & examples, and 'private' functions, which are used internally by the 'public' functions, but for which I do not wish to

Hi.... I am assembling an R package (under Windows XP, R v.2.3) and have a very basic question. Running R CMD build does not generate a pdf. R CMD check generates a .dvi, but I cannot figure out how to automatically create the pdf. I thought from reading the "writing r extensions' manual that R CMD build was supposed to generate the pdf. I am having no success using various

Hello... Is it possible to use "density" or another kernel density estimator to identify the mode of a distribution? When I use 'density', the resulting density plot of my data is much cleaner than the original noisy histogram, and I can clearly see the signal that I am interested in. E.g., suppose my data is actually drawn from two or more normal (or other)

This might be of interest: http://math.usu.edu/~minnotte/research/software/modetree.r (I was not able to get to the link, but google has a cached version.) Prof. Marron's SiZer maps may also be of interest, but AFAIK the code is in Matlab only. Andy From: Dan Rabosky > > > Hello... > > Is it possible to use "density" or another kernel density > estimator to

Hi I have a question about the creation of variables within lists in R. I am running simulations and am interested in two parameters, ESM and ESMM (the similarity of these names is important for my question). I do simulations to generate ESMM, then plug these values into a second simulation function to get ESM: x <- list() for (i in 1:nsimulations) { x$ESMM[i] <- do_simulation1()

Hello - I am having a problem with the function 'integrate'. I am running R on OSX (R 2.5.1). I am trying to suppress the error message when 'integrate' attempts to integrate across a parameter set giving a non-finite function value. I'm using it in a MCMC / simulated annealing algorithm, and it is entirely possible that some parameter sets will give non-

Dear all R-helper, if I want to use stdout from other language as my R program input ,which is the best way for design the API,using Pipe function or produce a temporary file,using scan function to read the file ? thanks Michael

Hi, I am trying to figure a way to allocate a string SEXP so that gc() won't ever collect it. Here is a little bit of a background. Suppose I want to write a .Call-callable function that upon each call returns the same value, say mkChar("foo"): SEXP getFoo() { return mkChar("foo"); } The above implementation doesn't take advantage of the fact that

The release schedule has now been set with a release date on April 18. The detailed procedure can be found at http://developer.r-project.org/ (it is not quite there yet, but will appear after a short propagation delay.) The main point for non-developers is that we start making alpha tarballs on March 21 and beta tarballs on April 4. Only very simple and/or critical bugs are fixed in the last

I have used Sweave a lot the latest year, but never really used any long function calls. If I have code which look like this ------------------------------------------------------------- gof <- benthic.flux(ID="Gulf of Finland", meas.conc=conc, bw.conc=bw.conc, time=times,

Hi, for presentation purposes I would like to clear to whole console window (like in a UNIX terminal: 'clear'). Is there such a function? If not, I could image that is not too hard to write such a function. Cheers, Hans

Can anyone point out why this is not working? y<-read.table(pipe(' awk '{ n = $1; sub( ".*" $1 " " ,"") ; while ( n-- ) print }' temp.txt ')) [[alternative HTML version deleted]]

Hello, This code works fine but is so fast I can't see anything but the last plot. for (i in nrow(X)){ plot(as.numeric(d[i,])) } I'd like to view a plot every 500 milliseconds, nrow(X) = 400. How? Thanks -- View this message in context: http://www.nabble.com/Plot-timer-in-a-for-loop-tp17741975p17741975.html Sent from the R help mailing list archive at Nabble.com.

Dear all, I would like to start a process from an R program in such a way that I can both feed input into the process and read the process's output. It seems that in R, I can have a pipe for writing into another process's input or a pipe for reading from another process's output, but not both. Doing both necessitates forking, such that the child can start the external process and

Dear all, I would like to start a process from an R program in such a way that I can both feed input into the process and read the process's output. It seems that in R, I can have a pipe for writing into another process's input or a pipe for reading from another process's output, but not both. Doing both necessitates forking, such that the child can start the external process and

How precise is R numerically? For example I wrote the following function for calculating the volume of the ball inscribed in the unit cube in m dimensions. In order to see what happens in 40 dimensions, I created an output of 24 digits. But how many are precise? Thanks Josef Eschgf?ller Ferrara --------------------------------------- Vol = function (m) {if (m<=1) 1 else Vol(m-2)*pi/(m+m)}

I have been using the fork package on a cluster to spawn jobs. I have noticed that I end up with many "R defunct" (zombie) processes following the use of fork that do not die completely until I close down R. Initially, I thought it may be my code but I ran the example from the latest fork package (see code below) on R.2.3.1 and ended up with the same problem. The code from the

Hi, I am using the following command to print to a file (I omitted the file details): cat( paste( paste(orientation, start, end, names,"\n"), paste(start, end, "exon\n"), sep="")) where "orientation" and "names" are character vectors and "start" and "end" are integer vectors. The problem is that R coerce the integer vectors

Can anyone please help me understand whats happening here? Thanks Tom getAmpRatio<-function(v_amps){ #calculates the amplitude ratios between the 3 largest amps and the rest bigamp<-0 map<-rep(TRUE,length(v_amps)) for(iLoc in 1:3){ bigamp<-bigamp+max(v_amps) map[which.max(v_amps)]<-FALSE v_amps<-v_amps[map]

Hi Dumb question time again, for which I apologise. I have a variable that contains the following numerical text "04010". This is the name to access a list: > as.list(KEGGPATHID2NAME)$"04010" [1] "MAPK signaling pathway" Marvellous! Except I want to do that when "04010" is assigned to a variable called path and I can't figure out how to do it!