Control the Generation of Random Numbers for Sampling from the Proposal Distribution¶
Description¶
The functions mvt.sampler
and mvnorm.sampler
can be used to fine-tune the
construction and sampling from the proposal distribution in the importance sampling
algorithm used in the MCEM algorithm to fit latpos
models.
Usage¶
mvt.sampler(df,fix.seed=TRUE) # Multivariate Student's t-distribution
mvnorm.sampler(fix.seed=TRUE) # Multivariate normal distribution
Arguments¶
df
-
a number, the degrees of freedom of the multivariate Student’s t-distribution.
fix.seed
-
a logical value; if TRUE the seed of the random number generated is reset in each MCEM iteration to make sure that approximate marginal log-likelihood surface is continuous while the simulation sample size is constant.
Value¶
Both functions return a list that contains a function to generate random numbers and to compute the log-density function of the distribution from which random numbers are generated.