Move Elogπ and ElogdetΛ out of ρnk into expectations()
Ween off threestats(), use a stats() implementation, almost there...