A multidimensional "which" function

September 16, 2011
By
A multidimensional "which" function

The well-known which function accepts a logical vector and returns the indices where its value equals TRUE. Actually, which also accepts matrices or multidimensional arrays. Internally, R uses a single index to run through such two- or higher-dimension...

Read more »

Soil-Landscape Block Diagrams in SoilWeb

September 16, 2011
By
Soil-Landscape Block Diagrams in SoilWeb

Users of our Google Earth interface to USDA-NCSS soils information will now see links to soil-landscape block diagrams listed within map unit descriptions. Automated Linking to NCSS Block Diagrams read more

Read more »

Always put comments in your code!

September 16, 2011
By
Always put comments in your code!

I have a paper which I wrote some years ago, which has not been finished, and which should be accompanied by an R package. So far nothing special, but at that time, I was only at the beginning of my affair with R, and so I made several mistakes (OK – I did also some

Read more »

Soil Series Query for SoilWeb

September 16, 2011
By
Soil Series Query for SoilWeb

A map depicting the spatial distribution of a given soil series can be very useful when working on a new soil survey, updating an old one, or searching for specific soil characteristics. We have recently added a soil series query facility to SoilWeb, w...

Read more »

Simulation studies in R – Using all cores and other tips

September 16, 2011
By

After working more seriously with simulations I noticed some updates were necessary to my previous setup. Most notably are the following three: It is very handy to explicitly call the different scenarios instead of using nested loops Storing intermediate results in single files obliviates the need to rerun an almost finished but crashed analysis and

Read more »

Beeswarm Plot with ggplot2

September 16, 2011
By
Beeswarm Plot with ggplot2

A colleague showed me results of his study project with beeswarm plots made by GraphPad. I was wondering if it could be implemented in R and more specifically with ggplot2. There is a R package allowing to draw such graphs, the beeswarm package (beeswa...

Read more »

Performance with ggplot2

September 16, 2011
By
Performance with ggplot2

Now after Reporting Good Enough to Share, let’s use ggplot2 and PerformanceAnalytics to turn this into this From TimelyPortfolio I have been notified that the colors aren’t great.  How does everyone like this? R code (click to download)...

Read more »

Statistics and Data Analysis in Python with pandas and statsmodels

September 16, 2011
By

Wes McKinney is a prominent figure in the scientific Python community, and has made tremendous contributions to several core statistical computing libraries in that language. This month, Wes will be speaking specifically about two packages he has crea...

Read more »

Datasets to Practice Your Data Mining

September 16, 2011
By
Datasets to Practice Your Data Mining

There are many datasets available online for free for research use. Some of them are listed below. - The R Datasets Package: There are around 90 datasets available in the package. Most of them are small and easy to feed … Continue reading →

Read more »

Project Euler: problem 2

September 16, 2011
By
Project Euler: problem 2

Each new term in the Fibonacci sequence is generated by adding the previous two terms. By starting with 1 and 2, the first 10 terms will be:1, 2, 3, 5, 8, 13, 21, 34, 55, 89, ...By considering the terms in the Fibonacci sequence whose values do not exc...

Read more »

How Lloyd’s of London uses R for Insurance

September 15, 2011
By
How Lloyd’s of London uses R for Insurance

Lloyd's is the world's leading specialist insurance market, and is often the first to insure new, unusual or complex risks. So it's no surprise that Lloyd's is one of the many companies that use R and its advanced capabilities for data analysis to help manage its insurance risks. At the useR! conference last month, Lloyd's analysts Markus Gesmann, Viren...

Read more »

Using Google Spreadsheets with R: an update

September 15, 2011
By

Prompted by a rush of visitors from Andrew Gelman's blog, I went back and updated the details of my post from 2009 on reading data from Google Spreadsheets into R. Since then, Google had switched to using a secure (https) connection for Google Docs, which required some tweaks to the code. If you haven't seen it before, it's a...

Read more »

Correlations among US Stocks: Is it really time to fire your adviser?

September 15, 2011
By
Correlations among US Stocks: Is it really time to fire your adviser?

Note: This post is NOT financial advice!  This is just a fun way to explore some of the capabilities R has for importing and manipulating data. The Financial Times says it's time to "Fire your Adviser" because correlations among US stocks ar...

Read more »

Reporting Good Enough to Share

September 15, 2011
By
Reporting Good Enough to Share

Sorry to all my faithful readers for my absence recently. I started a new job at a new firm, so my blogging has moved down the priority list but only temporarily. I am still committed to documenting my thoughts, especially finance and R thoughts as dis...

Read more »

littler 0.1.4

September 15, 2011
By

Matthias Klose, the tireless force behind the Debian / Ubuntu gcc, python, and what have you packages, sent me a minimal patch to let littler build when the ld linker uses the --as-needed option (as Ubuntu builds now do): all it took was a little reor...

Read more »

Project Euler: problem 1

September 15, 2011
By
Project Euler: problem 1

To be fairly honest (assuming there are degrees of honesty), I do know a little about math and programming but I don't know much math or any programming. I've loved math for a long time, but started to learn and understand fairly recently. So during th...

Read more »

Recent Updates in the aqp (Algorithms for Quantitative Pedology) Package for R

September 14, 2011
By

New version of our 'aqp' package for quantitative soils investigations, available on CRAN (version 0.99-5) and R-Forge (0.99-8). Some of the major changes are listed below: -------------------------- aqp 0.99-8 (2011-09-14) -------------------------- ...

Read more »

R Fork Bomb

September 14, 2011
By
R Fork Bomb

So maybe I’m a strange guy, but I think fork bombs are really funny.  What’s a fork bomb?  The basic premise is that you spawn a process that spawns a process that spawns a process…, ad infinitum. The most beautiful example of a fork bomb, and really one of the most beautiful lines of code

Read more »

Shortest paths to/from nodes of a certain type

September 14, 2011
By
Shortest paths to/from nodes of a certain type

Elijah asked the following via SOCNET mailing list: I was wondering if anyone knew of a script or tool which would give me the network distance of nodes to a particular class of nodes.  I think of this as an Erdos number, except instead of getting the distance to one node, I want the distance

Read more »

Revolution Analytics Fall Webinar Series

September 14, 2011
By

We've lined up what we think is an amazing series of R-related webinars over the next couple of months. These free 30-60 minute webinars will cover a wide range of topics: big-data analysis in R with the RevoScaleR package, Hadoop and Netezza; introductions to R for SAS users and for R users new to Revolution R; and applications of...

Read more »

Example 9.5: New stuff in SAS 9.3– proc FMM

September 13, 2011
By
Example 9.5: New stuff in SAS 9.3– proc FMM

Finite mixture models (FMMs) can be used in settings where some unmeasured classification separates the observed data into groups with different exposure/outcome relationships. One familiar example of this is a zero-inflated model, where some observat...

Read more »

How to program MapReduce jobs in Hadoop with R

September 13, 2011
By

MapReduce is a powerful programming framework for efficiently processing very large amounts of data stored in the Hadoop distributed filesystem. But while several programming frameworks for Hadoop exist, few are tuned to the needs of data analysts who typically work in the R environment as opposed to general-purpose languages like Java. That's why the dev team at Revolution Analytics...

Read more »

More sas7bdat progress

September 13, 2011
By

The development version of the read.sas7bdat function (in the sas7bdat package) now reads field labels and formats. In addition, errors of the type "found <x> <type> subheaders where 1 expected" are now a thing of the past. These improvements are largely due to work by Clint Cummins. The function also works on some files generated

Read more »

Backtesting a Simple Stock Trading Strategy

September 13, 2011
By
Backtesting a Simple Stock Trading Strategy

Note: This post is NOT financial advice!  This is just a fun way to explore some of the capabilities R has for importing and manipulating data.   I recently read a post on ETF Prophet that explored an interesting stock trading strategy in Ex...

Read more »

Speed up recursion in R 600-fold with Rcpp

September 12, 2011
By

Rcpp package co-author Dirk Eddelbuettel provides another case study in speeding up R code by rewriting repeatedly-called R code as inline C++ functions, using the classic Fibonacci recursion algorithm as an example. The speed gains here are impressive -- over 600x compared to native recursive R code -- but you could also improve performance by using a more efficient,...

Read more »

Why you should care about reproducible research

September 12, 2011
By

This week's Economist has an in-depth article on the consequences of failures reproducible research, adding more detail to the report in the New York Times in July. Errors in data analysis by researchers at Duke University led to patients in clinical trials being assigned the wrong drug: Dr Potti and his colleagues had mislabelled the cell lines they used...

Read more »

Testing and significance

September 12, 2011
By
Testing and significance

Julien Cornebise pointed me to this Guardian article that itself summarises the findings of a Nature Neuroscience article I cannot access. The core of the paper is that a large portion of comparative studies conclude to a significant difference between protocols when one protocol result is significantly different from zero and the other one(s) is(are)

Read more »

Forbush events

September 12, 2011
By
Forbush events

As noted here there is a new paper linking Forbush events with changes in DTR. Simply, during a Forbush event  cosmic rays are modulated ( the flux reaching the earth decreases. The theory goes something like this. If GCRs play a role in cloud formation, then when they decrease you should be able to detect an

Read more »

RQuantLib 0.3.8

September 12, 2011
By

A bug-fix release RQuantLib 0.3.8 is now on CRAN and in Debian. RQuantLib combines (some of) the quantitative analytics of QuantLib with the R statistical computing environment and language. Thanks to Helmut Heiming who noticed a side-effec t f...

Read more »