 Dear R-helpers,

My question is related to how to impose constraints when when sampling from a distribution.

For example, suppose I'm sampling a vector from a multivariate normal distribution

vbeta <- 100*diag(2)
mbeta <- c(1,1)
ans <- beta <- c(rmvnorm(1,mbeta,vbeta))

ans will thus be a vector with two elements. 

My question is how do I place a restriction on one of the elements of the ans vector. For example, suppose my goal is to reject any draw for ans[2] below -1 and above 1, how should I implement that in the above setting?

Thanks in advance for your help.

Kind Regards

Anup



      
	[[alternative HTML version deleted]]

