Monthly Archives: September 2012

Unit root, or not ? is it a big deal ?

September 10, 2012
By
Unit root, or not ? is it a big deal ?

Consider a time series, generated using set.seed(1) E=rnorm(240) X=rep(NA,240) rho=0.8 X=0 for(t in 2:240){X=rho*X+E} The idea is to assume that an autoregressive model can be considered, but we don't know the value of the parameter. ...

Read more »

Extending Gold time series

September 10, 2012
By
Extending Gold time series

While back-testing trading strategies I want all assets to have long history. Unfortunately, sometimes there is no tradeable stock or ETF with sufficient history. For example, I might use GLD as a proxy for Gold allocation, but GLD is only began trading in November of 2004. We can extend the GLD’s historical returns with its

Read more »

PDF and CDF for normal distributions with R

September 10, 2012
By
PDF and CDF for normal distributions with R

Below, we give the R code to plot the PDF and the CDF for normal distributions. We wish to get charts quite similar to the ones read on Wikipedia (Normal Distribution). The resulting charts are shown at the bottom. Notice that … Continue reading →

Read more »

igraph 0.6 issues: Changed numbering of Vertices

September 10, 2012
By

I have tried one of my previous scripts with an updated igraph version and I got an interesting (pretty much unexpected) error:At type_indexededgelist.c:269 : invalid (odd) length of edges vector, Invalid edge vectorThe problem is that it was a well-te...

Read more »

R-bloggers submission

September 10, 2012
By

I might be the 393rd blogger on the R-bloggers.com :-)  http://www.r-bloggers.com/Wohoo!

Read more »

Item Response Theory: Developing Your Intuition

September 10, 2012
By
Item Response Theory: Developing Your Intuition

Suppose that you accepted my argument from the last two posts on halo effects and bifactor models.  As you might recall, I argued that when respondents complete rating scales, they predominating rely on their generalized impression with a more minor role played by the specific features that the ratings were written to measure.  Consequently, we...

Read more »

Example 10.1: Read a file byte by byte

September 10, 2012
By
Example 10.1: Read a file byte by byte

More and more makers of electronic devices use standard storage media to record data. Sometimes this is central to the device's function, as in a camera, so that the data must be easy to recover. Other times, it's effectively incidental, and the device maker may not provide easy access to the stored data....

Read more »

Population health management with RevoScaleR

September 10, 2012
By

This guest post is by Douglas McNair MD PhD, Engineering Fellow & President, Cerner Math Inc. -- ed. RevoScaleR scaling big-data modeling performance for real-time health data analysis at Cerner The size of data sets is increasing much more rapidly than the speed of cores, of RAM, and of disk drives. This is particularly true of electronic health records...

Read more »

INLA functions (continued)

September 10, 2012
By
INLA functions (continued)

I have polished up one of the two functions I've thought of implementing for INLA and it's now available in the development version of the R package. So if you've got INLA installed in your R version (see how you can do it here, if you don't), you...

Read more »

Progress bars in R using winProgressBar

September 10, 2012
By
Progress bars in R using winProgressBar

Using progress bars in R scripts can provide valuable timing feedback during development and additional polish to final products. winProgressBar and setWinProgressBar are the primary functions for creating progress bars in R. Progress bars, and progress indicators in general, are relatively uncommon in R programming. This makes sense, as they can add bloat and, being The post Progress...

Read more »