Today I did some experiments to learn how to combine R and LaTeX to create reproducible research reports. Here are my first results:First Demo ...

Today I did some experiments to learn how to combine R and LaTeX to create reproducible research reports. Here are my first results:First Demo ...

The jobs Within the New Zealand Ministry of Business, Innovation and Employment, the Sector Trends team has recently secured resourcing for additional analysts on a range of statistical programmes. That’s the team that I usually manage, although for the next few months I’m doing a stint on a similar team, different topics. The formal details and position descriptions...

Last week, in our mathematical statistics course, we’ve seen the law of large numbers (that was proven in the probability course), claiming that given a collection of i.i.d. random variables, with To visualize that convergence, we can use > m=100 > mean_samples=function(n=10){ + X=matrix(rnorm(n*m),nrow=m,ncol=n) + return(apply(X,1,mean)) + } > B=matrix(NA,100,20) > for(i in 1:20){ + B=mean_samples(i*10) + } > colnames(B)=as.character(seq(10,200,by=10)) > boxplot(B) It is...

The Problem I have been attempting to use RMarkdown rather than LaTeX to produce R examples. One issue that has slowed my conversion has been my struggles with how to reference figures and tables. The examples below illustrate how I have been using captioner to solve this problem. The Solution Foundation For a simple example, I load FSAdata for the RuffeSLRH92...

My previous post covered the basics of logistic regression. We must now examine the model to understand how well it fits the data and generalizes to other observations. The evaluation process involves the assessment of three distinct areas – goodness of fit, tests of individual predictors, and validation of predicted values – in order to

“…likelihood inference is in a fundamental way more complicated than the classical method of moments.” Carlos Amendola, Mathias Drton, and Bernd Sturmfels arXived a paper this Friday on “maximum likelihood estimates for Gaussian mixtures are transcendental”. By which they mean that trying to solve the five likelihood equations for a two-component Gaussian mixture does not