Ween off returning rnk for lowerbound, by computing E[log p(Z|π)] 10.72 and E[log q(Z)] 10.75 in rnk()