similar to: FW: fairly simple file I/O

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

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

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.

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

I have a multidimensional array "a", for example, > a , , 1 [,1] [,2] [1,] 1 3 [2,] 2 4 , , 2 [,1] [,2] [1,] 5 7 [2,] 6 8 So a is 2x2x2. I have another array "b", for example, > b [,1] [,2] [1,] 9 11 [2,] 10 12 So b is 2x2. I want to "tack" b onto a so that b becomes a new "plane", so to speak,

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

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 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" }

I'm about to buy a new computer and give my old computer to my Mom. Because I will giving my old computer to my mom, I will be loosing my Windows XP license. I can either get a nice computer with a Windows Lincense or I can buy and really honking mega awesome computer without a Windows License. Naturally, I would rather buy a "honking mega awesome computer" instaed of a "nice

My grateful thanks to whoever can guide me in implementing this... I have a need to place calls via Asterisk Manager Protocol to a legacy PBX and twiddle its MWI lights. By some means I get notified that an MWI light has to be changed. Via AMP: 1. Send an ORIGINATE command honking up an incoming port on the PBX via a TA804. a. I presume that the TA804 should be a SIP channel. yes/no? b.

Thanks in advance for your patience and assistance. Since a system administrator's work is most successful when nobody knows you've done anything....(ugly that...you get neither raises nor equipment funding as a result), I prefer to learn from other's experiences rather than fail in the same pits. We are considering a switch from UW to Dovecot, since our users persist in using

I am writing a FLAC encoder/decoder and when I use the flac_static.lib that was compiled in debug mode in MS Visual Studio 2003, the seek_absolute function works fine. However, if I compile flac_static in release mode in the same environment, the seek_absolute function no longer works. No parts of my code change, the only difference is whether I link in the release lib instead of the debug lib.

I selected the defaults in partitioning my drive and df -h shows: Filesystem Size Used Avail Use% Mounted on /dev/mapper/VolGroup00-LogVol00 18G 2.5G 15G 15% / /dev/hda1 99M 12M 82M 13% /boot none 126M 0 126M 0% /dev/shm Now I am getting ready to install Scalix CE and it is telling me a whole story about what my