Hi all, I know pnorm() is able to calculate a sequence of probabilities by providing a vector of means, and a single number for location and standard deviation. For example, pnorm(0.5,c(0.5,2),1) will produce [1] 0.5000000 0.0668072. However, I wonder if its multivariate counterpart mvtnorm can do the same thing? Namely, if I provide a sequence of (vector-)means and a single vector for boundaries, and a single correlation matrix, can I get a sequence of probabilities? Many thanks! Jianing -- Life is a Markov Chain, your future will be independent with yesterday. Life is also a super-martingale, once you lose it, you are not expected to get it back. [[alternative HTML version deleted]]