similar to: Confidence interval on parameters from optim function

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 updated my R version from 2.7 to 2.8 (patch) and since then I can not open HDF5 files that I saved previously with the command hdf5save from the hdf5 library. As anyone had that problem yet? Any help would be appreciated. Thanks Emmanuel ------------------------------------------- Dr. Emmanuel Devred Bedford Institute of Oceanography, 1 Challenger Drive, Dartmouth, Nova

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

Dear all, I have a problem with the function Julian, may be a bug in the function ? Here is a vector of character, which represents dates (May 18 to May 20 2000): > amj <- c("2000-05-18","2000-05-18","2000-05-18","2000-05-19","2000-05-19" > ,"2000-05-19", "2000-05-19", "2000-05-20", "2000-05-20",

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

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

Dear all, I'm trying to merge 2 dataframes, but I'm not being entirely successful and I can't understand why. Dataframe x1 State_prov Shape_name bob2009 bob 2010 bob2011 Nova Scotia Annapolis 0 0 1 Nova Scotia Antigonish 0 0 0 Nova Scotia Gly NA NA

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

Hello, I'm playing around with cluster analysis, and am looking for methods to select the number of clusters. I am aware of methods based on a 'pseudo F' or a 'pseudo T^2'. Are there packages in R that will generate these statistics, and/or other statistics to aid in cluster number selection? Thanks, John. --

Hello, I'm working with a very large dataset (250,000+ lines in its' current form) that includes presence only data on various species (which is nested within different sites and sampling dates). I need to convert this into a dataset with presence/absence for each species. For example, I would like to expand "My current data" to "Desired data": My current data

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

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

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

i want to estimate parameters with maximum likelihood method with contraints (contant numbers). for example sum(Ai)=0 and sum(Bi)=0 i have done it without the constraints but i realised that i have to use the contraints. Without constraints(just a part-not complete): skellamreg_LL=function(parameters,z,design) { n=length(z); mu=parameters[1]; H=parameters[2]; Apar=parameters[3:10];

Hello All, I'm looking to test a variable in a logit model (glm(..., binomial(link="logit"))) for exogeneity (endogeneity). At this point I am planning to try implementing Jeffery Grogger's "A Simple Test for Exogeneity in Probit, Logit, and Poisson Regression Models", Economic Letters, 1990. To do this, I need to be able to do an instrumental variables NLS

Dear R-users, I would like to know if there is any way to constraints optimized parameters using the function lm, glm or others that are written in the form: Lm( formula, data ...) As I understand, formula are of the type y ~ X1 +X2+ ... Xi (where Y, X1, X2 ..Xi are vectors). In my case I would like the estimates of this linear combination computed with lm to be positives. I haven't found

I got a problem. Samba is in place and working on my linux machine but when I try and use Windows to map a drive it asks for a password and doesn't know the one I enter. How can I enter a new one, or find the one I currently have in place? --- Ted Gervais 1-902-679-2253 Coldbrook, Nova Scotia Canada.

Hi. I'm not having any luck installing the binary for the new Macintosh OS called "Tiger". I get error messages as below. Can anyone offer me some advice? (PS: I checked the bug list first but a search on "Tiger" turned up nothing so I hope it's OK to post here.) Error in dyn.load(x, as.logical(local), as.logical(now)) : unable to load shared

I've been using (and loving) R for quite a while now, but I have to admit that something simple is still stumping me. The question is how I can control the box within which a plot is drawn, in cases where I'm controlling the aspect ratio with the "asp" argument. The problem comes up in pdf() and png() output formats, but also with interactive windows, and so I'll