R/isTrMatrix.R

In tbuder/CellTrans: Quantifying stochastic cell state transitions

#' Check for transition matrix.
#'
#' Verification whether a matrix is a stochastic matrix of a Markov process or not, i.e. with non-negative entries and row sums equal to one.
#' @param A The matrix that is checked to be a transition matrix.
#' @keywords transition matrix, Matrix process
#' @export

isTrMatrix=function(A) {
  sumisOne=TRUE
  for (i in 1:nrow(A) ) { if (abs(sum(A[i,])-1)>1E-8) {sumisOne=FALSE}
  }
  
  return((all(A>=0)) & (sumisOne))
}