Fixed bug in code - sampling theta from correct full conditional