[R] masking a regular lat/lon grid to extract map boundaries

2008-10-07 Thread claudia tebaldi
Dear R-helpers, I have lat/lon coordinates of regularly spaced grid points, about 4Km apart, covering the entire US continental region. I would like to mask this rectangular grid in order to extract all and only the grid points within a specific region. Today I want to extract Montana, say, from

Re: [R] How to validate model?

2008-10-07 Thread Ajay ohri
This is an approach Run the model variables on hold out sample. Check and compare ROC curves between build and validation datasets. Check for changes in parameter estimates (co efficients of variables) p value and signs. Check for binning (response versus deciles of individual variables).

Re: [R] How to validate model?

2008-10-07 Thread Frank E Harrell Jr
Ajay ohri wrote: This is an approach Run the model variables on hold out sample. Check and compare ROC curves between build and validation datasets. Check for changes in parameter estimates (co efficients of variables) p value and signs. Check for binning (response versus deciles of

Re: [R] masking a regular lat/lon grid to extract map boundaries

2008-10-07 Thread Duncan Murdoch
On 07/10/2008 6:08 PM, claudia tebaldi wrote: Dear R-helpers, I have lat/lon coordinates of regularly spaced grid points, about 4Km apart, covering the entire US continental region. I would like to mask this rectangular grid in order to extract all and only the grid points within a specific

Re: [R] R and computer heat

2008-10-07 Thread Alexandre Aguiar
Em Terça 07 Out 2008, Dr Eberhard Lisse escreveu: Why don't you hack a little script that looks at system load, temperature and date/time and writes it somewhere. Then Done! Columns are time, % of cpu usage and cpu temp in Celsius degrees. The R script used was complete in about 2 min 10 sec

[R] Fortune?

2008-10-07 Thread Rolf Turner
Ben Bolker wrote (to the r-sig-mixed-models list): I still don't know what to do about the compromise between how statistics should be done and how journal editors seem to insist it should be done ... This seems to me to be worthy of being included as a fortune.

[R] follow up on [Rd] NAMESPACE methods guidance, please ( http://tolstoy.newcastle.edu.au/R/e4/devel/08/06/1901.html )

2008-10-07 Thread Tao Shi
This is a follow-up on the discussion originally posted on the R-devel list ( http://tolstoy.newcastle.edu.au/R/e4/devel/08/06/1901.html ), as I have encountered the exact same issue mentioned in Martin's email. Here is a simplified version of my problem:

Re: [R] How to validate model?

2008-10-07 Thread Pedro.Rodriguez
Hi, Yes, from my humble opinion, it doesnt make any sense to use the (2-class) ROC curve for a rating system. For example, if the classifier predicts 100% for all the defaulted exposures and 0% for the good clients, then even though we have a perfect classifier we have a bad rating system.

[R] help to debug C codes called from R

2008-10-07 Thread Droit Arnaud
Hello everybody, I have a package with a C codes called from R. I want to debug the C functions to check variables values and to include some breakpoints in the C codes. I am wondering if anyone knows of any tools to easily help debug in the R environment. Thank you very much for your help.

[R] Histogram colours in lattice.

2008-10-07 Thread Rolf Turner
I am trying to do a histogram lattice plot and I would like the histogram to be filled with a different colour in each panel. Note: I want every bar in each histogram to be the same colour, but that there should be different colours *between* histograms. Can't seem to get this to work. I

Re: [R] Histogram colours in lattice.

2008-10-07 Thread hadley wickham
On Tue, Oct 7, 2008 at 8:54 PM, Rolf Turner [EMAIL PROTECTED] wrote: I am trying to do a histogram lattice plot and I would like the histogram to be filled with a different colour in each panel. Note: I want every bar in each histogram to be the same colour, but that there should be

Re: [R] ordering problem

2008-10-07 Thread Moshe Olshansky
Try AA - apply(A,1,function(x) paste(x,collapse=)) and work with AA. --- On Tue, 30/9/08, Jose Luis Aznarte M. [EMAIL PROTECTED] wrote: From: Jose Luis Aznarte M. [EMAIL PROTECTED] Subject: [R] ordering problem To: [EMAIL PROTECTED] Received: Tuesday, 30 September, 2008, 8:43 PM Hi

Re: [R] masking a regular lat/lon grid to extract map boundaries

2008-10-07 Thread Greg Snow
Another approach that will probably be more work up front, but may be better in the long run if you are going to be doing a lot of this is to use the maptools and sp packages. Create a spatialpointsdataframe from your data, read in a us map (there are shapefiles that you can get from the US

Re: [R] How to validate model?

2008-10-07 Thread Frank E Harrell Jr
Ajay ohri wrote: the purpose of validating indirect measures such as ROC curves. Biggest Purpose- It is useful while in more marketing /sales meeting context ;) That is far from clear. It seems that ROC curves are being used to impress non-statisticians more than for shedding light on

[R] genoud nonlinear least squares optimisation

2008-10-07 Thread Rhys Whitley
Hello, I am trying to optimise a nonlinear model to derive 'best-fit' parameter esimates using the genoud function. I have been using the genetic algorithm - gafit - in order to do this, but I am getting parameter estimates that do not always reach the global minimum. I am very keen to apply

Re: [R] R and Multi threading

2008-10-07 Thread Prof Brian Ripley
On Tue, 7 Oct 2008, pejpm wrote: I will preface this message by saying that I am not an R developer and no very little about R...but here is my situation: One of my users has developed a model for analysing commodity prices. At the moment when he runs this model on his daily data set it takes

[R] question from Braun/Murdoch book

2008-10-07 Thread Erin Hodgess
Hi R People: I am looking at the Braun/Murdoch book, A First Course in Statistical Programming in R, and I have a question about a function there. It's on page 52, Example 4.5; the sieve of Erastosthenes. There is a line: primes - c() Is there a difference between using that and primes - NULL

<    1   2