similar to: (hopefully) simple array op

I would like to do something like the following: Fancyhist<-function(x,...) { # first, process x into xprocess somehow, then ... if (is.null(breaks)) { # yes, I know this is wrong # define the histogram breaks somehow, then call hist: hist(xprocess,breaks=breaks,...) } else { # use breaks give in calling argument hist(xprocess,...) } } But, those of you who know R better

I have a four-dimensional array. I have a three-dimensional array I want to add to it. My particular problem is conceptually very simple because the new array is just going to go into the end of the 4D array; i.e. will be added into the final index. I think I have seen an R command that can do this in one shot, but for the life of me, I can't find it again. Am I making this up, or is there

Am doing some vector plots with the arrows() function. Works well. But what I need to do is supply an arrow for scaling for the reader. I need to plot an arrow of some known magnitude somewhere on the page (preferably outside the bounds of the plot, so that it can be seen clearly) with some text underneath it that says, for instance, "10 kg-m/sec". Any ideas? Thanks.

OK, I needed to plot a set of vectors on top of a contour plot. I figured out a way to do this. I create a panel function that calls "larrows()" with arguments constructed from my vector data. Then, when I go to do the contour plot, I call contourplot() with the "panel" argument set to point to my newly created panel function. So far, so good. Now, I need to do a

I would like to plot some lat-lon data in a filled contour, and then overlay a map of the globe on top. I am trying something like this: filled.contour(lons, lats, glb.data, plot.axes={axis(1);axis(2);map(projection='rectangular',parameters=0,add =T)} ) The filled contour plots fine, and the map shows up, but only as a very tiny little black rectangle in a small area of the

OK, my apologies. I am sure this is a question that has been answered before. But I have looked all over the web and can't find an answer for it. I promise, wasting your time and bandwidth is my last resort. So here goes: I have an ASCII file formatted like so: Label 1.1 Time 1 Label 1.2 Array of data from time 1 Label 2.1 Time 2 Label 2.2 Array of data from time 2 Label

I want to do a series of contour plots, 4 in all. The data is coming from a data frame named "nd.frame", which has elements "xdf", "ydf", "zdf", and "pndt". I am treating "pndt" as a factor, and it has four levels. I make a call to the lattice graphics routine "contourplot" like so:

Sorry if I am cluttering up this list with too many rank beginner questions. But I *am* a rank beginner, and I am at my wit's end with this one: I have a data frame df, with components x,y,z. This command to make a contour plot with lattice graphics works just fine: > contourplot(z~x*y,data=df,cuts=20) Then I try > contourplot(z~x*y,data=df,cuts=20,region=T) And it dies on me with

Am trying to install package JGR and, by necessity, JavaGD on a few Linux platforms. R rev. 2.14.0. Had some success on one platform, but am running into a problem on a new platform. I installed rJava, JavaGD, iplots, and JGR. Things went - I thought - without a hitch. Now, when I try to use JGR, I start R and issue the command "library(JGR)" and get an error message: Error in

Am trying to install 2.14.0 on a CentOS system. The configure script apparently can't find something to do with the readline library: Configure: error: --with-readline=yes (default) and headers/libs are not available Yes, I know I need readline 4.2 or later. "yum" assures me that rev 5.1 is already installed. And I know I need to provide configuration flags if it is in

I just want to integrate a 3D data set along one dimension to obtain a 2D data set. Something like: (given array "d" with dim nx,ny,nz ...) data_int<-array(dim=c(nx,ny)) for (n in 1:ny) { for (m in 1:nx) { data_int[m,n]<-sum(d[m,n,]) } } The thing is, given R's facility with integers, it seems that I should be able to obtain data_int without the

I am trying to install R on a pretty up-to-date CentOS system. I have tried installing 2.14.0 and 2.13.2. In both cases, the configure step fails, and does not produce a Makefile. So, of course, I can't issue "make". The errors that I can find in config.log are, first, several missing files: ac_nonexistent.h, minix/config.h., readline/history.h, and readline/readline.h. Then,

Here it is: read.inputs<-function(infile) { for (counter in 1:length(input.list)) { seek(infile,where=0,origin='start') newline<-readLines(infile,n=1) while(length(newline)>0) { if(!is.na(grep(as.character(input.list[[counter]][1]),newline)[1])) { break } newline<-readLines(infile,n=1) } if (length(newline)>0)

I would like to open an existing netCDF file and add a variable to it. I am using the ncdf package. This test code gives the idea of what I am trying to do: library(ncdf) print('here we go') print('first, construct netCDF file') t<-dim.def.ncdf('t','',1:1,unlim=T,create_dimvar=F) print('defining first vars')

I am using readLines to read a fairly large ASCII file. readLines reads a fixed number of lines, then other R code processes the data, then readLines reads the same number of lines again, then other R code processes the data, then .... Sort of like: conn<-file('filename','r') for (chunk in 1:100000) { Lines<-readLines(conn,n=25) # process "Lines" }

Dear R Help Community, I have a question and an answer (based on reading this forum and online research), but I though I should share both since probably there's a much better way to go about my solution. My question is specifically about how to best visualise multiple response contingency tables. What I mean by 'multiple response' is that the total number of responses per row of a

For a data frame, we can add an additional row or column easily. For example, we can add an additional row of zero and an additional row of column as follows. Is there an easy and similar way to add zeros in each dimension? For example, for array(1:12, dim=c(2,2,3))? Thanks for your help!! Hanna > x <- as.data.frame(matrix(1:20,4,5))> x[5,] <- 0> x[,6] <- 0> x V1 V2

Dear R-users, we are using library(snow) for computation on a linux cluster with RMPI. We have a problem with clusterEvalQ: after launching clusterEvalQ it seems loading the required library on each node but if we type a function belonging to the loaded package R doesn't find it. > library(snow) # making cluster with 3 nodes > cl <- makeCluster(3, type = "MPI") Loading

Hi everyone I've been playing with do.call() but I'm having problems understanding it. I have a list of "n" elements, each one of which is "d" dimensional [actually an n-by-n-by ... by-n array]. Neither n nor d is known in advance. I want to bind the elements together in a higher-dimensional array. Toy example follows with d=n=3. f <-

I'm trying to massage some data from Matlab into R. The matlab file has a "struct array" which when imported into R using the R.matlab package, becomes an R list with 3+ dimensions, the first of which corresponds to the structure fields, with corresponding row names, and the second and third+ dimensions correspond to the dimensions of the original struct array (as matlab