KF: simulate always started at mean.  If we give it distribution then we should start randomly from distribution.
KF: Distributions has a logpdf function.  Changed code to use logpdf() instead of log(pdf())