This function initialise raw gene expression values in a matrix. Return either a matrix of (1) continuous values or (2) binary values.
Note that kmeans clustering as binarisation only works well if the data has a bimodal distribution.
character. Specify whether max expression value is the lowest (as in qPCR), or the highest (as in RNAseq and microarray). Option: 'low', 'high'. Default to 'high'.
uni_thre
numerical. Speficy threshold for unimodality test. Default to 0.2.
scale
logical. Whether to scale the data to a range of 0-1. Default to T.
discretised
logical. Whether to return discretised data. Default to F.