This book, “Bayesian Cognitive Modeling: A Practical Course,” by Michael Lee and E. J. Wagenmakers, has a bunch of examples of Stan models with discrete parameters—mixture models of various sorts—with Stan code written by Martin...

This book, “Bayesian Cognitive Modeling: A Practical Course,” by Michael Lee and E. J. Wagenmakers, has a bunch of examples of Stan models with discrete parameters—mixture models of various sorts—with Stan code written by Martin...

There’s something for everyone! I had a lot of fun just copying the titles to make this list, as I’d already forgotten about a lot of this stuff. Here are the scheduled posts, in order through 31 Dec: Fitting models with discrete parameters in Stan How to use lasso etc. in political science? An unconvincing […] The post On deck through the rest of 2015 appeared first on Statistical Modeling,…

You've had a long day. You've implemented a custom algorithm in the SAS/IML language. But before you go home, you want to generate some matrices and test your program. If you are like me, you prefer a short statement—one line would be best. However, you also want the flexibility to […] The post Ten "one-liners" that create test matrices for statistical programmers appeared first on The DO Loop.

I’ve been banging my head against this problem, on and off, for a couple months now. It’s an EP-like algorithm that a collaborator and I came up with for integrating external aggregate data into a Bayesian analysis. My colleague tried a simpler version on an example and it worked fine, then I’ve been playing around […] The post The Final Bug, or, Please please please please please work this time!…

Check this out this "Retraction Watch" post, forwarded by a reader:http://retractionwatch.com/2014/11/11/overly-honest-references-should-we-cite-the-crappy-gabor-paper-here/Really funny. Except that it's a little close to home. I suspect th...

Attention conservation notice: A ponderous, scholastic joke, which could only hope to be amusing to those who combine a geeky enthusiasm for over-written horror stories from the early 20th century with nerdy enthusiasm for truly ancient books. I wish...

Attention conservation notice: A ponderous, scholastic joke, which could only hope to be amusing to those who combine a geeky enthusiasm for over-written horror stories from the early 20th century with nerdy enthusiasm for truly ancient books. I wish...

Attention conservation notice: Publicity for an upcoming academic talk, of interest only if (1) you care about quantifying uncertainty in statistics, and (2) will be in Pittsburgh on Monday. I am late in publicizing this, but hope it will help drum u...

Attention conservation notice: Publicity for an upcoming academic talk, of interest only if (1) you care about quantifying uncertainty in statistics, and (2) will be in Pittsburgh on Monday. I am late in publicizing this, but hope it will help drum u...

In some of my recent public talks (for example: here and here) I have mentioned a desire for “a deeper theory of fitting and testing.” I thought I would expand on what I meant by this. In this note I am going to cover a lot of different topics to try and suggest some perspective. … Continue reading A deeper theory of testing

OK, explain to me this email: God day, How are you? My name is **. I came across your contact email at the University of Cyprus, Department of Economics. I seek for a private Economics teacher for my Daughter. I would like to know if you would be available for job. If you would be […] The post Annals of Spam appeared first on Statistical Modeling, Causal Inference, and Social…

Tuesday, at the end of my 5-hour crash course on machine learning for actuaries, Pierre asked me an interesting question about computational time of different techniques. I’ve been presenting the philosophy of various algorithm, but I forgot to mention computational time. I wanted to try several classification algorithms on the dataset used to illustrate the techniques > rm(list=ls()) > myocarde=read.table( "http://freakonometrics.free.fr/myocarde.csv", head=TRUE,sep=";") > levels(myocarde$PRONO)=c("Death","Survival") But the dataset is rather small,…

Inspired by X’s blog’s Le Monde puzzle entries, I have a little Stan coding puzzle for everyone (though you can solve the probabilty part of the coding problem without actually knowing Stan). This almost (heavy emphasis on “almost” there) makes me wish I was writing exams. Puzzle #1: Inferring Ability from Streaks Suppose a player […] The post Stan Puzzle #1: Inferring Ability from Streaks appeared first on Statistical Modeling,…

Statisticians often need to integrate some function with respect to the multivariate normal (Gaussian) distribution, for example, to compute the standard error of a statistic, or the likelihood function in of a mixed effects model. In many (most?) useful cases, these integrals are intractable, and must be approximated using computational methods. Monte-Carlo integration is one … Continue reading Notes on Multivariate Gaussian Quadrature (with R Code) →

For the problem I reported yesterday, in which my R package was working fine locally but failing on Travis, the key solution is to run update.packages(ask=FALSE) locally, and maybe even update.packages(ask=FALSE, type="source") to be sure to grab the source of packages for which binaries are not yet available. I now know to do that. In […]

“The samples were collected in privacy, using passive drool procedures, and frozen immediately.” Anna Dreber sends along a paper, “Assessing the Robustness of Power Posing: No Effect on Hormones and Risk Tolerance in a Large Sample of Men and Women,” which she published in Psychological Science with coauthors Eva Ranehill, Magnus Johannesson, Susanne Leiberg, Sunhae […] The post Low-power pose appeared first on Statistical Modeling, Causal Inference, and Social Science.

Jane Frazier spoke at our research team meeting today on “Reproducibility in computational research”. We had a very stimulating and lively discussion about the issues involved. One interesting idea was that reproducibility is on a scale, and we can all aim to move further along the scale towards making our own research more reproducible. For […]

Hmmmm, coverage for travel delay, that might not be so bad. This is Amtrak, after all. Let’s click through to the fine print: They’ll cover me for a departure delay of six or more hours, huh? Nice try in your attempt to scam me out of $8.50. It didn’t work this time, but, hey, why […] The post Amtrak is evil appeared first on Statistical Modeling, Causal Inference, and Social…

I've released BayesFactor 0.9.12-2 to CRAN; it should be available on all platforms now. The changes include:Added feature allowing fine-tuning of priors on a per-effect basis: see new argument rscaleEffects of lmBF, anovaBF, and generalTestBFFixed bug...

Somehow when my code stops working, my first (and second, and third) reaction is to blame everything except my own code. (“It’s not me, it’s you.”) And almost always, it’s my own code that’s the problem (hence the title of this post). I spent the day trying to resolve a bug in my early-in-development R […]