This function is used to perform multiple imputation for missing data given the joint distribution.
Imputation(data, theta, psi, method = "MaxProb")- incomplete dataset
- vector of probability for each component
- specific probability for each variable in each component
- methods for imputation, including "Sampling" and "MaxProb". In "Sampling" method, sample missing values from the conditional distribution randomly. In "MaxProb" method, impute missing values with maximal probability from the conditional distribution.
ImputedData - dataset has been imputated.