diff --git a/stats.py b/stats.py
index <HASH>..<HASH> 100644
--- a/stats.py
+++ b/stats.py
@@ -79,13 +79,14 @@ def sample_mn(Sigma,M,K):
     right = np.linalg.cholesky(K)
     return M + left.dot(np.random.normal(size=M.shape)).dot(right.T)
 
-def sample_mniw(dof,kappa,lmbda,M,K):
+def sample_mniw(dof,lmbda,M,K):
     Sigma = sample_invwishart(lmbda,dof)
-    return sample_mn(Sigma,M,K), Sigma # TODO TODO divide by kappa here?
+    return sample_mn(Sigma,M,K), Sigma
 
 def sample_invwishart(lmbda,dof):
     # TODO make a version that returns the cholesky
     # TODO allow passing in chol/cholinv of matrix parameter lmbda
+    # TODO lowmem! memoize!
     n = lmbda.shape[0]
     chol = np.linalg.cholesky(lmbda)
added comments about a new sample_invwishart function needed