[This article was first published on Rmetrics blogs, and kindly contributed to R-bloggers]. (You can report issue about the content on this page here) Want to share your [Read more...]

Last week (October 1-3) MPK Analytics attended the annual Cyber Summit in Banff. The theme for this year was; “Leading the Way in the Age of Big Data“. As might […] [Read more...]

[This article was first published on Yet Another Blog in Statistical Computing » S+/R, and kindly contributed to R-bloggers]. (You can report issue about the content on this page here) Want to share your content on R-bloggers? click here if you have a blog, or here if you don't. ################################################# ## FIT A MULTIVARIATE ADAPTIVE REGRESSION ## ## SPLINES MODEL (MARS) USING MDA PACKAGE ## ## DEVELOPED BY HASTIE AND TIBSHIRANI ## ################################################# # LOAD LIBRARIES AND DATA library(MASS); library(mda); data(Boston); # FIT AN ADDITIVE MARS MODEL mars.fit <- mars(Boston[, -14], Boston[14], degree = 1, prune = TRUE, forward.step = TRUE) # SHOW CUT POINTS OF MARS cuts <- mars.fit$cuts[mars.fit$selected.terms, ]; dimnames(cuts) <- list(NULL, names(Boston)[-14]); print(cuts); factor <- mars.fit$factor[mars.fit$selected.terms, ]; dimnames(factor) <- list(NULL, names(Boston)[-14]); print(factor); # EXAMINE THE FITTED FUNCTION BETWEEN EACH IV AND DV par(mfrow = c(3, 5), mar=c(2, 2, 2, 2), pty="s") for (i in 1:13) { xp <- matrix(sapply(Boston[1:13], mean), nrow(Boston), ncol(Boston) - 1, byrow = TRUE); xr <- sapply(Boston, range); xp[, i] <- seq(xr[1, i], xr[2, i], len=nrow(Boston)); xf <- predict(mars.fit, xp); plot(xp[, i], xf, xlab = names(Boston)[i], ylab = "", type = "l"); } To leave a comment for the author, please follow the link [...] [Read more...]

The Barron’s article Still Too Pricey by Andrew Bary looks at the share price of the Facebook and based on the P/E ration valuation metrics concludes that even at the current prices, stock is overvalued. I want to show how to do this type of fundamental analysis using ... [Read more...]

So...I'm back with Your Shape: Fitness Evolved 2012 for XBox Kinect. Why? Because I want to loose some weight and get back in shape of course -;)The reason I stop playing the game is simple...I'm lazy...but this time, I have come back with a goal...bur... [Read more...]

[Read more...]

[This article was first published on Rmetrics blogs, and kindly contributed to R-bloggers]. (You can report issue about the content on this page here) Want to share your [Read more...]

[This article was first published on Rmetrics blogs, and kindly contributed to R-bloggers]. (You can report issue about the content on this page here) Want to share your [Read more...]

One Paragraph Summary Always explore your data visually. Whatever specific hypothesis you have when you go out to collect data is likely to be worse than any of the hypotheses you’ll form after looking at just a few simple visualizations of that data. The most effective hypothesis testing framework ... [Read more...]

I gave a short talk today to the [Davis R Users’ Group] about ggplot. This what I presented. Additional resources at the bottom of this post ggplot is an R package for data exploration and producing plots. It produces fantastic-looking graphics and allows one to slice and dice one’s ... [Read more...]

As we noted last month, the new Themes feature in ggplot2 helps you customize the design of R charts to your liking. Now, R user Jeffrey Arnold has built on this feature to create standardized themes to make R graphics looks like those from major publications and other software systems. ... [Read more...]

BSMAP is an aligner for bisulfite sequencing reads. It outputs aligned reads as well as methylation ratios per base (via methratio.py script). The methylation ratios can be read into R via methylKit package and regular methylKit analysis can ... [Read more...]

[Read more...]

Henry John-Alder told me once that in a marathon, twice as runners cross the line at 2h 59m than at 3h 00m. He pointed out that this anomaly in the distribution of finishers per minute (roughly normal shaped) is due … Continue reading → [Read more...]

This Gist is mostly for my future self, as a reminder of how to find distances between each row in two different matrices. To create a distance matrix from a single matrix, the function dist(), from the stats package is sufficient. There are times, ho... [Read more...]

Open your sources.list file in geditsudo gedit /etc/apt/sources.listand add the following line:deb http://cran.cnr.berkeley.edu/bin/linux/ubuntu/ precise/Note that you don't have to use that mirror. You may use any mirror from the list here : http://cran.r-project.org/mirrors.htmlAdd ... [Read more...]

I have previously described and back-tested the Permanent Portfolio strategy based on the series of posts at the GestaltU blog. Today I want to show how we can improve the Permanent Portfolio strategy perfromance using following simple tools: Volatility targeting Risk allocation Tactical market filter First, let’s load the ... [Read more...]

Another bug-fix release of Armadillo, now at version 3.4.3 whike the 3.4.* stabilizes, and with it a version 0.3.4.3 of RcppArmadillo, our wrapper for R and Armadillo. The new version is already on CRAN as of earlier today. Once again no R level or i... [Read more...]

Release 0.2.6 of RProtoBuf arrived on CRAN earlier this morning. RProtoBuf provides GNU R bindings for the Google Protobuf data encoding library used and released by Google. This release was once more driven largely by Murray whom we have now add... [Read more...]

This is another of those confluence style posts, where a handful of things I’ve read in quick succession seem to phase lock in my mind: (brought to mind in part via @downes a week or so ago: How to Synch 32 Metronomes) The first was a post by Alan Levine ... [Read more...]