Semester page for STK4021 - Autumn 2022
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
Nov. 2, 2022 10:40 AM
The last part about comparing waic (or loo) is a bit tricky and is ok if you skip this point. For those interested, note that the pointwise object gives elppd. In order to get lppd you need to add p_waic (see eq (7.13) in the book). Then you can adjust the individual lppd values according to the transformation before you subtract p_waic again in order to get the elppd value on the right scale.
Oct. 27, 2022 11:26 AM
Syllabus
Student representatives
Exercises
Resources