## Create heat maps with PROC SGPLOT

August 17, 2015
By

When SAS 9.4m3 was released last month (including SAS/STAT and SAS/IML 14.1), I was happy to see that a HEATMAP statement had been added to the SGPLOT procedure. Although heat maps in the SAS/IML language have been available for several releases, you previously had to use the Graph Template Language […] The post Create heat maps with PROC SGPLOT appeared first on The DO Loop.

## R, Python, and SAS: Getting Started with Linear Regression

August 17, 2015
By

Consider the linear regression model, $$y_i=f_i(\boldsymbol{x}|\boldsymbol{\beta})+\varepsilon_i,$$ where $y_i$ is the response or the dependent variable at the $i$th case, $i=1,\cdots, N$ and the predictor or the independent variable is the $\boldsymbol{x}$ term defined in the mean function $f_i(\boldsymbol{x}|\boldsymbol{\beta})$. For simplicity, consider the following simple linear regression (SLR) model, $$y_i=\beta_0+\beta_1x_i+\varepsilon_i.$$ To obtain the (best) estimate of $\beta_0$ and $\beta_1$, we solve for the least residual sum of…

## A case in which metric data are better analyzed by an ordinal model

August 17, 2015
By

Here we consider some data that might have been smoothly distributed over a metric scale, but ended up being concentrated on only a few values. The usual treatment of the data as normally or t-distributed is not appropriate, and instead the data are bi...

## Rockin the tabloids

August 16, 2015
By

Rick Gerkin points me to this opinion piece from a couple years ago by biologist Randy Schekman, titled “How journals like Nature, Cell and Science are damaging science” and subtitled “The incentives offered by top journals distort science, just as big bonuses distort banking.” Here’s Schekman: The prevailing structures of personal reputation and career advancement […] The post Rockin the tabloids appeared first on Statistical Modeling, Causal Inference, and Social…

## Seattle histogram

August 16, 2015
By

Filed under: pictures, R, Statistics, Travel Tagged: histogram, sculpture, Seattle, Washington Convention Center

## Manski on Uncertainty in Official Statistics

August 16, 2015
By

Chuck Manski has a fascinating forthcoming Journal of Economic Literature piece, "Communicating Uncertainty in Official Economic Statistics: An Appraisal Fifty Years After Morgenstern."Manski's subtitle refers to Oskar Morgenstern's book, On ...

## The Rise of the Robots (Advisors…)

August 15, 2015
By

The Asset Management industry is on the verge of a major change. Over the last couple of years Robots Advisors (RA) have emerged as new players. The term itself is hard to define as it encompasses a large variety of services. Some are designed to help traditional advisers to better allocate their clients money and some are […]

## Calling R from Scala sbt projects using rscala

August 15, 2015
By

Overview In the previous post I showed how the rscala package (which has replaced the jvmr package) can be used to call Scala code from within R. In this post I will show how to call R from Scala code. I have previously described how to do this using jvmr. This post is really just … Continue reading Calling R from Scala sbt projects using rscala

## Calling R from Scala sbt projects using rscala

August 15, 2015
By

Overview In the previous post I showed how the rscala package (which has replaced the jvmr package) can be used to call Scala code from within R. In this post I will show how to call R from Scala code. I have previously described how to do this using jvmr. This post is really just … Continue reading Calling R from Scala sbt projects using rscala

## Calling Scala code from R using rscala

August 15, 2015
By

Introduction In a previous post I looked at how to call Scala code from R using a CRAN package called jvmr. This package now seems to have been replaced by a new package called rscala. Like the old package, it requires a pre-existing Java installation. Unlike the old package, however, it no longer depends on … Continue reading Calling Scala code from R using rscala

## Calling Scala code from R using rscala

August 15, 2015
By

Introduction In a previous post I looked at how to call Scala code from R using a CRAN package called jvmr. This package now seems to have been replaced by a new package called rscala. Like the old package, it requires a pre-existing Java installation. Unlike the old package, however, it no longer depends on … Continue reading Calling Scala code from R using rscala

## Why couldn’t Breaking Bad find Mexican Mexicans?

August 15, 2015
By

Watching “Breaking Bad” . . . I’m told on good authority that may of the actors playing Mexicans are not actually Mexican; some of them can barely speak Spanish at all. Whassup with that? How hard is it to find a Mexican actor in LA...

## ShinyStan v2.0.0

August 14, 2015
By

For those of you not familiar with ShinyStan, it is a graphical user interface for exploring Stan models (and more generally MCMC output from any software). For context, here’s the post on this blog first introducing ShinyStan (formerly shinyStan) from earlier this year. ShinyStan v2.0.0 released ShinyStan v2.0.0 is now available on CRAN. This is […] The post ShinyStan v2.0.0 appeared first on Statistical Modeling, Causal Inference, and Social Science.

## Harry S. Truman, Jesus H. Christ, Roy G. Biv

August 14, 2015
By

Are there any others? The post Harry S. Truman, Jesus H. Christ, Roy G. Biv appeared first on Statistical Modeling, Causal Inference, and Social Science.

## STAN trailer [PG+53]

August 13, 2015
By

[Heading off to mountainous areas with no Internet or phone connection, I posted a series of entries for the following week, starting with this brilliant trailer of Michael:] Filed under: Kids, R, Statistics, University life Tagged: Andrew Gelman, Hami...

## STAN trailer [PG+53]

August 13, 2015
By

[Heading off to mountainous areas with no Internet or phone connection, I posted a series of entries for the following week, starting with this brilliant trailer of Michael:] Filed under: Kids, R, Statistics, University life Tagged: Andrew Gelman, Hami...

## Hey—Don’t trust anything coming from the Tri-Valley Center for Human Potential!

August 13, 2015
By

Shravan sends along this article by Douglas Peters and Stephen Ceci, who report: We selected 12 already published research articles by investigators from prestigious and highly productive American psychology departments, one article from each of 12 highly regarded and widely read American psychology journals with high rejection rates (80%) and nonblind refereeing practices. With fictitious […] The post Hey—Don’t trust anything coming from the Tri-Valley Center for Human Potential! appeared…

## Odds and Probability: Commonly Misused Terms in Statistics – An Illustrative Example in Baseball

Yesterday, all 15 home teams in Major League Baseball won on the same day – the first such occurrence in history.  CTV News published an article written by Mike Fitzpatrick from The Associated Press that reported on this event.  The article states, “Viewing every game as a 50-50 proposition independent of all others, STATS figured the […]

## Odds and Probability: Commonly Misused Terms in Statistics – An Illustrative Example in Baseball

Yesterday, all 15 home teams in Major League Baseball won on the same day – the first such occurrence in history.  CTV News published an article written by Mike Fitzpatrick from The Associated Press that reported on this event.  The article states, “Viewing every game as a 50-50 proposition independent of all others, STATS figured the […]

## Classic Data Visualizations

August 12, 2015
By

My thanks to Veronica Johnson at Investech.com for drawing my attention a recent piece of theirs relating to Classic Data Visualizations.As they say:"A single data visualization graphic can be priceless. It can save you hours of research. They’re eas...

## UCLA Statistics 2015 Commencement Address

August 12, 2015
By

I was asked to speak at the UCLA Department of Statistics Commencement Ceremony this past June. As one of the first graduates of that department back in 2003, I was tremendously honored to be invited to speak to the graduates. When I arrived I was just shocked at how much the department had grown. When

## Correlation is not a measure of reproducibility

August 12, 2015
By

Biologists make wide use of correlation as a measure of reproducibility. Specifically, they quantify reproducibility with the correlation between measurements obtained from replicated experiments. For example, the ENCODE data standards document states A typical R2 (Pearson) correlation of gene expression (RPKM) between two biological replicates, for RNAs that are detected in both samples using RPKM or read counts, should

## Reprint of “Observational Studies” by William Cochran followed by comments by current researchers in observational studies

August 12, 2015
By

Dylan Small organized this discussion in the new journal, Observational Studies. Cochran’s 1972 article is followed by comments from: Norman Breslow Thomas Cook David Cox & Nanny Wermuth Stephen Fienberg Joseph Gastwirth & Barry Graubard Andrew Gelman Ben Hansen & Adam Sales Miguel Hernan Jennifer Hill Judea Pearl Paul Rosenbaum Donald Rubin Herbert Smith Mark […] The post Reprint of “Observational Studies” by William Cochran followed by comments by current…