Wrote script for running the PMF algorithm for different values of K, storing the RMSE and variational lower bound. Will run it on the real data now