R/make_Z_fortran.r

In multimix: Fit Mixture Models Using the Expectation Maximisation (EM) Algorithm

#'  Read Z from FORTRAN output. Make into R matrix
#'
#' The FORTRAN version of Multimix produces two output files:
#' GENERAL.OUT and GROUPS.OUT. The latter mainly contains the Z matrix.
#'
#' This function facilitates the obtaining of Multimix R output given
#' Multimix FORTRAN output.
#'
#' @param gr.out string containing  a file name.
#'
#' @return a matrix containing a \eqn{Z}{Z} matrix.
#' @importFrom utils read.table
#' @export
#' @author Murray Jorgensen
#' @examples
#' Z <- make_Z_fortran(system.file('extdata', 'GROUPS-BP-Multimixf90.OUT', 
#'                     package = 'multimix'))
make_Z_fortran <- function(gr.out = "groups.out") {
  # Read Z from FORTRAN output. Make into R matrix.
  grout <- read.table(gr.out, header = FALSE)
  nc <- ncol(grout)
  Z <- as.matrix(grout[, -c(1, nc)])
  colnames(Z) <- NULL

  return(Z)
}