R/GWrite.R

In fchamroukhi/HDME: LASSO Regularized Mixture of Experts Models

GWRITERES = function(MAXbetak, MAXwk, MAXS, MAXLOG, MAXBIC, Y, X, K)
{
  # if(!(MAXbetak[1,1]*MAXbetak[1,2])) return()
  # if(MAXbetak[1,2] < MAXbetak[1,1])
  # {
  #   temp = MAXbetak[,1]
  #   MAXbetak[,1] = MAXbetak[,2]
  #   MAXbetak[,2] = temp
  #   MAXwk = - MAXwk
  # }
  # source("GEstep.R")
  n = dim(X)[1]
  Para = cbind(MAXbetak, t(MAXwk))
  utils::write.table(Para, file = "GPara.txt", sep = '\t', append = FALSE, col.names = FALSE, row.names = FALSE)
  utils::write.table(sqrt(MAXS), file = "GSigma.txt", append = FALSE, col.names = FALSE, row.names = FALSE)
  utils::write.table(MAXLOG, file = "GLOG.txt", append = FALSE, col.names = FALSE, row.names = FALSE)
  utils::write.table(MAXBIC, file = "GBIC.txt", append = FALSE, col.names = FALSE, row.names = FALSE)
  MAXtau = Ge.step(MAXbetak, MAXwk, MAXS, Y, X, K)
  MAXP = Pik(n, K, X, MAXwk)
  utils::write.table(MAXP, file = "GMAXP.txt", sep = '\t', append = FALSE, col.names = FALSE, row.names = FALSE)
  MAXClass = apply(MAXtau, 1, which.max)
  Data = cbind(X,Y)
  RDATA = cbind(Data, MAXClass)
  #utils::write.table(RDATA, file = "Restore Data.txt", sep = '\t', append = TRUE, col.names = FALSE, row.names = FALSE)
  #DATA = DATA[-seq(1,300),]#=========Data size = 300
  utils::write.table(RDATA, file="GRestore Data.txt",sep = '\t', append = FALSE, col.names = FALSE, row.names = FALSE)
  #NoR <<- NoR + 1
  #print(paste("NoR: ", NoR))
}