plot_convergence function
First column is deviance
plot_convergence finished
stat_IRWLS partly works