Beta to M

BetaToM.R

BetaToM <- function(beta){
  matrix <- log2(beta) - log2(1-beta)
  impute.matrix(matrix)
}

MToBeta <- function(mmatrix){
  matrix <- 2^mmatrix/(2^mmatrix+1)
  impute.matrix(matrix)
}

impute.matrix <- function (matrix) {
	missing <- which(is.na(matrix) | !is.finite(matrix), arr.ind=TRUE)
	if (nrow(missing)!=0){
		for (j in 1:nrow(missing)){
			mean <- mean(matrix[missing[j,1],][is.finite(matrix[missing[j,1],])], na.rm=TRUE)
			matrix[missing[j,1],missing[j,2]] <- mean
		}
	}
    	matrix
}