Compulsory exercise: Results related to Gamma model

Some of you might have got some results that looks strange for the models in problem 3 (all effects seems to be zero).

However, this seems mainly to be due to that the effects become small (but still significant). By using the option digits=5 either in the print or summary commands on the object you get, you will be able to see the results.  Further, if your object is fit.Bayes, you may get a better plot of the random effects by using the command

plot(fit5.Bayes,pars="b")

Geir