this function generates missing data in a complete data matrix
insertMVs(original, mean.THR, sd.THR, MNAR.rate)
complete data matrix containing all measurements
- parameters of the threshold distribution which controls the MVs rate (mean.THR should be initially set such that the result of the initial thresholding, in terms of no. of NAs, equals the desired total missing data rate) - example: if one wants to generate 30 mean.THR can be set as follows: mean.THR = quantile(pepExprsData, probs = 0.3) - sd.THR is usually set to a small value (e.g. 0.1)
percentage of MVs which are missing not at random
A list that contains the original complete data matrix, the data matrix with missing data and the percentage of missing data