Which of the following two charts (both created with R) to you prefer? This dotchart: Or this bar chart? Andrew Gelman (who, incidentally, is speaking at the October NYC UseR meeting) prefers the dotchart prefers a line plot (update: see Gelman's comment, below), but personally I think the bar chart ... [Read more...]

[Read more...]

If you use blogger or even wordpress you've probably found that it's complicated to post code snippets with spacing preserved and syntax highlighting (especially for R code). I've discovered a few workarounds that involve hacking the blogger HTML templ... [Read more...]

Background: In some cases it is necessary to conduct the same analysis multiple times on either the same or different data. In such circumstances it is worth writing a function to simplify the code. In this example the location quotient provides a simple calculation easily written in to a function. ... [Read more...]

Since posting last night, there's been several more aftershocks, including several big ones. The total now stands (at 7:00pm) at 304 aftershocks, 70 greater than 4 on the Richter scale, 29 greater than 4.5, and 10 greater than 5. A couple of the fives ... [Read more...]

In the recent GDAT class, confidence intervals (CI) for performance data were discussed. Their generalization to confidence bands (CB) for scalability projections using the USL model also came up informally. I showed a prototype plot but it was an ugly hack. Later requests from GDAT attendees to apply CBs to ... [Read more...]

Wasnt that hard. green is land, blue is sea. new data awaits you. old data taken down For Zeke: Green is landcoast cells, blue is seacoast.  The cells are AREA weighted per cce request but not Fractiionally weighted. That is if a cell is 90% land and 10% water, this weighting is ... [Read more...]

Having purchased the September edition of La Recherche because of its (disappointing!) coverage on black matter, I came by a short coverage on an Intel circuit producing “truly random” numbers… I already discussed this issue in an earlier post, namely that there is no reason physical generators are “more” random ... [Read more...]

This is an updated version of my Making Maps with R tutorial. I think the code is lot simpler and it also includes some data for you to play around with. Background: Spatial data are becoming increasingly common, as are the tools available in R to process it. Of course ... [Read more...]

An interesting thing that I found while doing the analysis for my previous post on the Canterbury earthquakes, was the difficulty of constraining figures less than a negative number. For example:abcd__ a[1] 1 2 3 4 5 6 7 8 9 10__ b[1] 4 5 6__ c[1] TRUE TRUE TRUE TRUE TRUE TRUE FALSE FALSE FALSE FALSE__ d[1] -1 -2 -3 -4 ... [Read more...]

Some of you may know this R reshape package already, I have started to play with it after the post Handling Large CSV Files in R. It is really an excellent one worthing a new post to introduce formally.What is reshape package? reshape: Flexibly reshape data, Reshape lets you ... [Read more...]

4:35 am on Saturday morning, I was woken up by a 7.1 magnitude earthquake, the epicentre being approximately 30 km west of where I live in Christchurch, New Zealand. The earthquake was one of the largest that New Zealand has experienced, and has been the most destructive one since the earthquake that struck ... [Read more...]

UPDATE: cce asked for a chart of coastal locations broken out by their SST component and their Land componement. That’s a little tricky, but once you figure out the mask logic its just addition and such.There are TWO masks: One mask contains ALL zeros, except for coast cells ... [Read more...]

Having missed the earliest entry by Radford last month, due to disconnection in Yosemite, I was stunned to read his three entries of the past month about R performances being significantly modify when changing brackets with curly brackets! I (obviously!) checked on my own machine and found indeed the changes ... [Read more...]

Replaces 4.0, minor changes made here and there. totally replaces 4.0 as the testing move forward. see 4.1 Run: DownloadAll. Run: DataSetUp. SSTTest reruns the SST analsyis and compares at HADSST2 Land1900Present: reruns the old version of program with a bench against CRU land figures. Next version will integrate more of raster ... [Read more...]

[Read more...]

UPDATE: Ron Broberg has a more definitive explanation of the difference which indicates that 5sig issue is not the main cause of the difference. See his exposition here CRU, it appears, trims out station data when it lies outside 5 sigma. Well, for certain years where there was actually record cold ... [Read more...]

UPDATE: Ron Broberg has a more definitive explanation of the difference which indicates that 5sig issue is not the main cause of the difference. See his exposition here A short update. I’m in the process of integration the Land Analysis and the SST analysis into one application. The principle ... [Read more...]

[Read more...]

Introduction I’ve just uploaded a new package to CRAN based on a stimulus generation algorithm that I use for my experiments on vision. The FourierDescriptors package provides methods for creating, manipulating and visualizing Fourier descriptors, which are a representational scheme used to describe closed planar contours. The canonical reference ... [Read more...]