In this post I will show my first try at a commodity index substitute. Regular readers know my frustration with proprietary data as I try to demonstrate various techniques to users who might not have the resources to pay for the data. I hav...

I’m a HUGE fan of Structural Equation Modeling. Heck, I’ve even begun teaching short SEM workshops. (seriously, if you want, email me and we’ll talk). For those of you unfamiliar with the technique, it’s awesome for three main reasons. It’s a method of teasing apart direct and indirect interactions in your data. It allows you

While we were at R/Finance in Chicago over the weekend (more on that later), it was a great pleasure to meet in person Paul Teetor, author of the R Cookbook (published by O'Reilly). If you haven't picked it up yet, I recommend checking it out: if you already have a basic familiarity of the R syntax, but want to...

This was my first year to attend the R/Finance conference that focuses on the use R programming in applied finance. I was unable to get out there until mid-morning on Friday, so I missed Jeff Ryan’s tutorial on Automated Trading … Continue reading →

A friend asked me this morning if there was a way to plot a symbol in R (as a plotting character) representing a half-filled circle. I didn't know, but I figured this out (perhaps it's demonstrated elsewhere -- the ability to use Unicode symbols was added in 2008 or so -- but I didn't stumble across it). First, looking...

Recently, Kevin Ummel gave a short conference about his PhD research at the Sustainable futures – Tyndall Centre researchers’ international conference. His work, named “Optimizing worldwide deployment of wind and solar power technologies via a high resolution, linear programming model”, is built with R. He is using the solaR package for the solar irradiance calculations,

Although the official deadline is long past, I just became aware of this call for a postdoctoral position in Paris starting next Fall. Please feel free to apply: Summary Adaptive Markov Chain Monte Carlo (MCMC) methods are currently a very active field of research. MCMC methods are sampling methods, based on Markov Chains which are

I’m starting on my real R-Scripts now. We got an assesment for some Reports and they needed my graphs and Veerle’s Reports. So the tasks are: Well distribution plot Quality plate control Surface Plot Heatmap CirclePlot and more to come̷...

Are you a non-windows user wanting to share your R package with others who do use windows? There's a handy utility (http://win-builder.r-project.org/)that allows you to upload your .tar.gz file and it checks it, builds it, and sends you a link to downl...

Are you a non-windows user wanting to share your R package with others who do use windows? There's a handy utility (http://win-builder.r-project.org/)that allows you to upload your .tar.gz file and it checks it, builds it, and sends you a link to downl...

(The examples here work with the version of insidefunctor tagged as "v1") I ran into an interesting problem using "inside" functors. Something is wrong in the following code (well, depending on what you thought it should do). > library(insidefunctor) > `%+.%` = fmap(`+`) > x = seq(0, 10, len = 50) > plot(x,...

Hell is new and improved. The R Inferno has been revised. If you don’t know of it, it is a short explanation of a few trouble spots when using the R language. Somehow the short explanation grew to approach book-length. It can be found at the usual place: http://www.burns-stat.com/pages/Tutor/R_inferno.pdf Major improvements An index has been … Continue reading...

The solar irradiance incident outside the earth’s atmosphere is called the extraterrestial or extra-atmospheric irradiance. It is derived from the solar constant only with geometric equations. It can be easily calculated with the calcSol function of the solaR package. With this post I will show an example with some packages from the Spatial task view.

Gompertz Model Visualization # Gomperz growth function gomp <- function(x, a, b, k) a*exp(-b*exp(-k*x)) # Normal model with Gompertz mean function likelihood <- function(weight, age, sigma, a, b, k) { mu <- gomp(age, a, b, k) dnorm(weight, mu, sigma) } # Visualize the model visualize <- function(phi=40, theta=-35) { weight <- seq(0, 250,

I remember a few weeks ago, there was a challenge in the R-help list to make the prime symbol in R graphics. In LaTeX, we simply write $X'$ or $X^\prime$. R has a rough support for math expressions (see demo(plotmath)) and they are certainly unsatisfactory for LaTeX users. In fact we can write native LaTeX