similar to: HDF5 and R 2.8.0patched

Dear R-users, I haven't found a way in the searchable archive to overplot a contour (lines) over a surface. I have a (n,m) matrix that represents sea surface temperature that I have plotted using image.plot(), filled.contour() or image(). I would like to overplot this image with some contour lines of mixed layer depth values(same size matrix). How can I do this? Any help is appreciated,

Hi everyone, I have two questions: I would like to get confidence intervals on the coefficients derived from the optim() function. I apply optim() to a given function f > res <- optim(c(0.08,0.04,1.),f,NULL,method="L-BFGS-B",lower=c(0.,0.,0.)) And I would like to get the p-value and confidence intervals associated with > res$par My second question deals with error message. I

Had centos 4.1 set up and running with X. Had to replace the monitor to an Acer LCD monitor (model AL1714). When I try to reconfigure X, "system-config-display" seems ok going through the first stage of detecting the graphics card but then returns a floating "Input not supported" message to the screen and nothing happens afterward. Any hints on solving this issue will be

Hi, I have been working with R for the last year and using the UKFSST package to look at satellite tag track data and SST information. Fpr those not familiar, the package uses the positions estimated by the satellite tags themselves and the associated SST data from servers (in this case, from the University of Hawaii and the NOAA Coastwatch website) for the time preiod to estimate a reasonable

Hi, I am having problems using .First() with the Mac OS X gui version of R 2.7.2 (under OS X 10.5.5); .First is ignored on startup. .First() does work in the command-line version of R under Mac OS X but my specific task requires the gui version. I can find no reference to this issue in the Mac OS X FAQ or in other manuals, does anyone have a solution for this? Thanks, Ian -- Ian Jonsen

I have Samba 2.2.3a running on Tru64 Unix v4.0f (?) I have Win NT4.0 (sp6a) as our local PC OS. I have created a Samba password file for users. I've tried several permutations, but currently I have security = user encrypt passwords = yes unix password sync = yes username map = /../../../usermap.txt I also have standard password encryption on the NT side, i.e. not plaintext. If I sync

Hi there, I am currently working on something that uses hdf5 library. I think hdf5 is a great data format, I've used it somewhat extensively in python via PyTables. I was looking for something similar to that in R. The closest I can get is this library: hdf5. While it does not work the same way as PyTables did, but it's good enough to let them exchange data via hdf5 file. There is just 1

I've installed hdf5 [library & includes] here locally, last week. [see ?hdf5save in R] There is a remaining compilation / linking problem for me: The hdf5 library uses "compress2" a routine defined in Gnu's zlib, but not available in the zlib that comes with X11R6.3 (which we have installed). Now unfortunately, the final link statement for R (0.65) [for me, with hdf5

HDF5 is a high quality data storage format from NCSA (http://hdf.ncsa.uiuc.edu/index.html) Marcus Daniels of the Swarm Development Group has made available an R library called hdf5 which allows loading and saving of data in hdf5 format. The source for his library is here: ftp://ftp.swarm.org/pub/swarm/src/testing/hdf5_1.4.7.tar.gz I prefer to use RPMs to maintain software across many linux

Hello everybody, I'm using hdf5 files to store results from intermediate calculations. These are usually part of a list, called "res". As I want the hdf-files to contain all the members of res in its top "directory", I used to do attach(res) do.call("hdf5save", args=c(fileout=file.path(dir, ofile), as.list(names(res)))) detach(res) which did what I

Hello, I'm having difficulty passing an object name to a lapply function. Can somebody tell me the trick to make this work? #Works T13702 <- TRACKDATA[["13702.xls"]][["data"]] min(unlist(lapply(list(T13702), function(x) mdy.date(x[1, 2], x[1, 1], x[1, 3])))) 16553 #Works d<-2 assign(paste("T",substr(names(TRACKDATA)[d],1,(nchar(names(TRACKDATA)[d]

You can test if the problem is accumulation in memory registers, which is certainly what this sounds like. Just do a loop over a reasonably small number of iterations and store or print the time between each iteration. If memory accumulation it will run optimally for the first few iterations, after which the time will increase noticeably (essentially exponentially, hence ultimately freezes up). If

Hello, Might anyone know how to conduct Wald-type F-tests of the fixed effects estimated by glmmPQL? I see this implemented in SAS (GLIMMIX), and have seen it recommended in user group discussions, but haven't come across any code to accomplish it. I understand the anova function treats a glmmPQL fit as an lme fit, with the test assumptions based on maximum likelihood, which is inappropriate

Hello, I'm working with a package that uses download.file in functions to extract information from remote databases. My current environment is Windows XP Pro SP3, R 2.7. A full extraction can be a great deal of data, so the download is accomplished in generally manageable packets, such that a single download will result in many files, which are written to a directory. It is not uncommon for a

I tried to use calcMin with a function that uses a number of ... arguments (all args from resid on) besides the vector of parameters being fit. Same idea as optim, nlm, nlminb for which this form of ... syntax works. But with calcMin I get an error regarding unused arguments. No partial matches to previous arguments that I can see. Anybody know the reason or fix for this?

Hi. When I try to use rsync 2.5.7 from a Windows 2000 machine with CygWin and ZoneAlarm Pro 4.5.538, I get this error: Read from remote host myweb.com: Connection reset by peer rsync: read error: Connection reset by peer rsync error: error in rsync protocol data stream (code 12) at /home/lapo/packaging/tmp/rsync-2.5.7/io.c(177) And this is the command that I am executing: rsync -vaz

I am trying to get contents of a REST response: getURL("http://localhost/myweb-app/rest-ws") This is a web application (myweb-app) which is providing a REST web service (rest-ws)... Unfortunately, the HTTP status sent back is 404. If I request the url using Chrome/IE, I get a HTTP status 200 OK. In Opera the request does not succeed either. I am using 2.15.1 (Win7, 64Bit) and just

You could calculate the confidence interval of the correlation at your desired df: http://davidmlane.com/hyperstat/B8544.html The below code takes as arguments the observed correlation, N, and alpha, calculates the confidence interval and checks whether this includes 0. cor.test2=function(r,n,a=.05){ phi=function(x){ log((1+x)/(1-x))/2 } inv.phi=function(x){

--=-=-= Included below are: 1) Diffs to to remove HDF5 from R 2) Diffs to pass thru --with-* options to a module's configure script during "R INSTALL". (To reiterate a plea for this capability...) This applies only to Unix, the Windows (mingw32 in my case) build approach of passing in make variables like *-CFLAGS and *-LDFLAGS works fine. The only hitch

Hi all, Looking to calculate the expected mean and variance of a ratio distribution where the source distributions are gaussian with known parameters and sample values are correlated. I see (from wikipedia: http://en.wikipedia.org/wiki/ Ratio_distribution#Gaussian_ratio_distribution) that this calculation is quite involved, so I'm hoping that someone has already coded a function to