## highlight 0.2-5

December 7, 2010
By

I pushed highlight 0.2-5 on CRAN. This release improves the latex renderer and the sweave driver so that multiple lines character strings are properly rendered. This example vignette shows it: \documentclass{report} \begin{document} <<echo=FALSE,results=hide>>= old.op <- options( prompt = " ", continue = " " ) @ <<>>= require( inline ) require( Rcpp ) convolve <- cxxfunction( signature( a = "numeric", b...

Read more »

## R 2.12.1 scheduled for December 16

December 6, 2010
By

The next update to R will be a patch release: R 2.12.1 will be released on December 16, as announced today by the R Core Team. As is typical for a patch release, this version will include some minor bug fixes plus a few new features (from the current build's NEWS file): The DVI/PDF reference manual now includes the...

Read more »

## R 2.12.1 scheduled for December 16

December 6, 2010
By

The next update to R will be a patch release: R 2.12.1 will be released on December 16, as announced today by the R Core Team. As is typical for a patch release, this version will include some minor bug fixes plus a few new features (from the current build's NEWS file): The DVI/PDF reference manual now includes the...

Read more »

## Jeromy Anglim on Reproducible Research and R

December 6, 2010
By

Jeromy Anglim, fellow social scientist and R aficionado from across the globe, gave a great talk to the Melbourne R Users Group last week on the joys of creating reproducible results. A subject near and dear to me, but not one that is given enough attention in research training. Jeromy discusses tools for generating reproducible

Read more »

## What my R code looks and feels like (Vanilla)

December 6, 2010
By

Geoff Robinson discusses how to write simple and reusable R code. He provides several examples with code.

Read more »

## Reproducible Research and R Workflow

December 6, 2010
By

Jeromy Anglim discusses best practices and tools used to support reproducible research in R. He discusses tools for R scripting and version control; and provides examples of SWeave. This presentation was given to the Melbourne R Users Group (MelbURN)...

Read more »

## 3 weak days in a row

December 6, 2010
By

Recently, Trading the odds posted one of many flavors of mean reverting strategies and I decided to get my hands dirty by writing R code and testing it. You can find full description of the strategy by following latter link above. Long story short – if SPY shows lower open, high and close 3 days in

Read more »

## JAGS – Bayesian Analysis

December 6, 2010
By

JAGS is used for Bayesian analysis using MCMC and stands for Just Another Gibbs Sampler.  It is an alternative to WinBUGS and can be accessed through R just like WinBUGS (via R2jags or RJags).  It will work on a Mac unlike WinBUGS.  The ...

Read more »

## Using the "Divide by 4 Rule" to Interpret Logistic Regression Coefficients

December 6, 2010
By

I was recently reading a bit about logistic regression in a book on hierarchical/multilevel modeling when I first learned about the "divide by 4 rule" for quickly interpreting coefficients in a logistic regression model in terms of the predicted probabilities of the outcome. The idea is pretty simple. The logistic curve (predicted probabilities) is steepest at the center where...

Read more »

## Example 8.17: Logistic regression via MCMC

December 6, 2010
By

In examples 8.15 and 8.16 we considered Firth logistic regression and exact logistic regression as ways around the problem of separation, often encountered in logistic regression. (Re-cap: Separation happens when all the observations in a category sha...

Read more »

## Electoral Marimekko Plots

December 6, 2010
By

To be reductive, visual displays of quantitative information might be reasonably categorized on a continuum between “data display” and “statistical graphics.” By statistical graphics, I mean a plot that displays some summary of or relationship amongst several variables, likely having undergone some processing or analysis. This may be as simple as a scatterplot of a … Read more

Read more »

## Bear hunting

December 6, 2010
By

When were there bear and bull markets in US stocks since 1950? Smoothing While we’d really like to estimate the expected return at each point in time, finding bear markets is ambitious enough.  The plan starts by smoothing the daily returns through time, as in Figure 1. Figure 1: Smoothed returns with a 4 year … Continue reading...

Read more »

## Forecasting workshop: Switzerland, June 2011

December 6, 2010
By

I will be running a workshop on Statistical Forecasting: Principles and Practice in Switzerland, 20-22 June 2011. Check out the venue: Waldhotel Doldenhorn, Kandersteg! So if you fancy a trip to the beautiful Swiss Alps next June, read on… Outline Forecasting is required in many situations: deciding whether to build another power generation plant in

Read more »

## R 2.12.0 and Eclipse with StatET installation

December 5, 2010
By

A while back I outlined the setup process for running R under the Eclipse integrated development environment with the help of WalWare.de’s StatET plugin for Eclipse, in Windows 7. This post is an update on that procedure for the newer versions of...

Read more »

## Bayesian adaptive sampling

December 5, 2010
By

In the continuation of my earlier post on computing evidence, I read a very interesting paper by Merlise Clyde, Joyee Ghosh and Michael Littman, to appear in JCGS. It is called  Bayesian adaptive sampling for variable selection and model averaging. The sound idea at the basis of the paper is that, when one is doing

Read more »

## Rethinking ‘loess’ for Binomial-Response Pitch F/X Strike Zone Maps

December 5, 2010
By

So after a long hiatus, I'm back for today. I've been crazy busy with a number of different things--including getting engaged and helping plan out wedding dates and things of that sort--and unfortunately have not kept up here on this blog (or on Fanta...

Read more »

## Pareto plot party!

December 5, 2010
By

A Pareto plot is an enhanced bar chart. It comes in useful for deciding which bars in your bar chart are important. To see this, take a look at some made up DVD sales data. set.seed(1234) dvd_names <- c("Toy Tales 3", "The Dusk Saga: Black Out", "Urban Coitus 2", "Dragon Training for Dummies", "Germination", "Fe

Read more »

## Genetic optimization for Trading Strategies using Rapidminer and R

December 5, 2010
By

That is the second tutorial of Rapidminer and R extension for Trading and the first in Video. In the last example the ROC obtained is not as good as it should be to make money in this business, To improve the strategy we will try to optimize the trading strategy. Different methods of optimization and objective functions...

Read more »

## Genetic optimization for Trading Strategies using Rapidminer and R

December 5, 2010
By

That is the second tutorial of Rapidminer and R extension for Trading and the first in Video. In the last example the ROC obtained is not as good as it should be to make money in this business, To improve the strategy we will try to optimize the trading strategy. Different methods of optimization and objective functions...

Read more »

## GLMM and R issues

December 4, 2010
By

I have been trying to run a Generalized Linear Mixed Model (GLMM) for some count data with repeated measures on sub-sampled sites and fixed effects at the site level with covariates at both the sub-plot and time levels.  Plus there are different n...

Read more »

## Comparison of results

December 4, 2010
By

I am doing a simple comparison of different estimation procedures in dealing with a simple binomial model. Here is where I got started:---------------------------------------------library(INLA)library(npmlreg)library(MCMCglmm)library(DPpackage)data(See...

Read more »

## Comparison of results

December 4, 2010
By

I am doing a simple comparison of different estimation procedures in dealing with a simple binomial model. Here is where I got started:---------------------------------------------library(INLA)library(npmlreg)library(MCMCglmm)library(DPpackage)data(See...

Read more »

## Root finding

December 4, 2010
By

Numerical root finding methods use iteration, producing a sequence of numbers that hopefully converge towards a limits which is a root. In this post, only focus four basic algorithm on root finding, and covers bisection method, fixed point method, Newton-Raphson method, and secant method. Read More: 1886 Words Totally

Read more »

## Root finding

December 4, 2010
By

Numerical root finding methods use iteration, producing a sequence of numbers that hopefully converge towards a limits which is a root. In this post, only focus four basic algorithm on root finding, and covers bisection method, fixed point method, Newton-Raphson method, and secant method. Read More: 1896 Words Totally

Read more »

## Le Monde puzzle [48: resolution]

December 4, 2010
By
$Le Monde puzzle [48: resolution]$

The solution to puzzle 48 given in Le Monde this weekend is rather direct (which makes me wonder why the solution for 6 colours is still unavailable..) Here is a quick version of the solution: Consider one column, 1 say. Since 326=5×65+1, there exists one value c with at least 66 equal to c. Among

Read more »

## A Draft of ProjectTemplate v0.2-1

December 3, 2010
By

I’ve just uploaded a new binary of ProjectTemplate to GitHub. This is a draft version of the next release, v0.2-1, which includes some fairly substantial changes and is backwards incompatible in several ways with previous versions of ProjectTemplate. Foremost of the changes is that most of the logic for load.project() is now built into the

Read more »

## Fun with infochimps: Animated Blog Post Hit Map

December 3, 2010
By

In a few weeks I will be visiting Chicago, and JD Long—the organizer of the local R users group—has graciously invited me to give a presentation. Ostensibly, the presentation will be on my recently released infochimps package, so I thought it was a good time to start actually putting together some examples and documentation for

Read more »

## Some ideas on communicating risks to the general public

December 3, 2010
By

SOME EMPIRICAL BASES FOR CHOOSING CERTAIN RISK REPRESENTATIONS OVER OTHERS This week DSN posts some thoughts (largely inspired by the work of former colleagues Stephanie Kurzenhäuser, Ralph Hertwig, Ulrich Hoffrage, and Gerd Gigerenzer) about communicating risks to the general public, providing references and delicious downloads where possible. Representations to use less often Single-event probabilities as

Read more »

## Because it’s Friday: The 4th Amendment X-ray T-shirt

December 3, 2010
By

Here's a clever way to express your objection to the obtrusiveness of the new backscatter X-ray scanners at airports: A T-shirt with the text of the 4th Amendment printed in metallic ink: Also available in underwear format. Cargo Collective: 4th Amendment Wear

Read more »