Normal CDF implementation in stan

Currently the likelihood calculations use the stan `Phi()` function for CDF terms but `normal_lpdf()` for PDF terms. 
Investigate whether there are computational (speed, stability) or consistency gains to be made by using `Phi_approx()` or `normal_(l)cdf()` instead.