Add nsim argument to vi_permute()

[bgreenwell]

The nsim argument (which will default to 1) enables multiple replications of shuffling each feature to reduce Monte Carlo error and help stabilize results from run to run (and maybe even produce standard deviations?).

imp <- vi_permute(fit, nsim = 10, ...)