use optim() to compute proposal distribution