make dummy variables and covariates

probabilities for S, Psi and p