Petr Savicky-2, thanks you. It seems the code works very well. I use rate=1, so "common" is between 0 to 1. Thanks again. -- View this message in context: http://r.789695.n4.nabble.com/How-to-generate-bivariate-exponential-distribution-tp3596944p3602539.html Sent from the R help mailing list archive at Nabble.com.