Exemplo n.º 1
0
arma::rowvec rmvlx(arma::mat R, arma::rowvec Z){
	Rcpp::RNGScope scope;
	int m = R.n_rows;
	double e  = Rf_rexp(1);
	arma::rowvec ans(m);
	ans = sqrt(e) * rmvnormx(R, Z);
	return ans;
}
Exemplo n.º 2
0
Type rexp(Type rate)
{
  return Rf_rexp(asDouble(rate));
}