Added sigma parameter (noise sd)
Runs on simdata -- not sure if working