Coloring Sudokus
Someday you will find me caught beneath the landslide (Champagne Supernova, Oasis) I recently read a book called Snowflake Seashell Star: Colouring Adventures in Numberland by Alex Bellos and Edmund Harris...continue reading.
Category: Rstat
Statistical methods accounting for human coding habits in databases
This is some work I did in 2013 that never led to funding to take it further. I did a … Morecontinue reading.
Does conspecific negative density dependence in forests really correlate with species rarity and latitude?
By Lisa Hülsmann and Florian Hartig (U Regensburg), with short comments by Ryan Chisholm (NUS) and Matteo Detto (Princeton). By Betteridge’s law, the answer to this question is of course no. Or better:...continue reading.
Safe Disposal of Unexploded WWII Bombs
Abstract Unexploded WWII bombs are ticking threats despite being dropped more than 70 years ago. In this post we explain how statistical methods are used to plan the search and...continue reading.
A Comparative Review of the RKWard GUI for R
Introduction RKWard is a free and open source Graphical User Interface for the R software, one that supports beginners looking to point-and-click their way through analyses, as well as advanced...continue reading.
PubMed retractions report has moved
A brief message for anyone who uses my PubMed retractions report. It’s no longer available at RPubs; instead, you will find it here at Github. Github pages hosting is great,...continue reading.
SAS Implementation of ZAGA Models
In the previous post https://statcompute.wordpress.com/2017/09/17/model-non-negative-numeric-outcomes-with-zeros/, I gave a brief introduction about the ZAGA (Zero-Adjusted Gamma) model that provides us a very flexible approach to model non-negative numeric responses. Today, I...continue reading.
Elegant regression results tables and plots in R: the finalfit package
The finafit package brings together the day-to-day functions we use to generate final results tables and plots when modelling. I spent many years repeatedly manually copying results from R analyses...continue reading.
Scaling of MCMC with dimension (experiments)
Hi all, In this post, I’ll go through numerical experiments illustrating the scaling of some MCMC algorithms with respect to the dimension. I will focus on a simple setting, to...continue reading.
The Pleasing Ratio Project
Music is a world within itself, with a language we all understand (Sir Duke, Stevie Wonder) This serious man on the left is Gustav Theodor Fechner, a German philosopher, physicist...continue reading.
Model averaging in ecology: a review of Bayesian, information-theoretic and tactical approaches for predictive inference
This (co-)guest post by Carsten F. Dormann & Florian Hartig summarizes a comprehensive review on model averaging for predictive inference, just published in Ecological Monographs.continue reading.
Mimicking SQLDF with MonetDBLite
Like many useRs, I am also a big fan of the sqldf package developed by Grothendieck, which uses SQL statement for data frame manipulations with SQLite embedded database as the...continue reading.
An R Tutorial: Visual Representation of Complex Multivariate Relationships Using the R qgraph Package, Part Two Repost
This is a repost of the original article that was posted as an embedded PDF file. Douglas M. Wiig April 8, 2018 Abstract This article is part of my...continue reading.
A macro to execute PROC TTEST for multiple binary grouping variables in SAS (and sorting t-test statistics by their absolute values)
In SAS, you can perform PROC TTEST for multiple numeric variables in the same procedure. Here is an example using the built-in data set SASHELP.BASEBALL; I will compare the number...continue reading.
Do you really need a multilevel model? A preview of powerlmm 0.4.0
In this post I will show some of the new simulation features that will be available in powerlmm 0.4.0. You can already install the dev version from GitHub. # GitHub...continue reading.
MLE with General Optimization Functions in R
In my previous post (https://statcompute.wordpress.com/2018/02/25/mle-in-r/), it is shown how to estimate the MLE based on the log likelihood function with the general-purpose optimization algorithm, e.g. optim(), and that the optimizer...continue reading.
Career Panel at the 2018 Canadian Statistics Student Conference – McGill University, Montreal, Quebec
I will speak on the career-advice panel at the 2018 Canadian Statistics Student Conference. It will be held on Saturday, June 2, at McGill University. If you will attend this...continue reading.
Yes, statistical errors are slowing down scientific progress!
Over at dynamic ecology, Jeremy Fox argues that Technical statistical mistakes are overrated; ecologists (especially students) worry too much about them. Individually and collectively, technical statistical mistakes hardly ever appreciably...continue reading.
Bayesian Inference with Backfitting MCMC
Previous posts in this series on MCMC samplers for Bayesian inference (in order of publication): Bayesian Simple Linear Regression with Gibbs Sampling in R Blocked Gibbs Sampling in R for...continue reading.
Read Random Rows from A Huge CSV File
Given R data frames stored in the memory, sometimes it is beneficial to sample and examine the data in a large-size csv file before importing into the data frame. To...continue reading.
