# ========================================================================================================
rm(list=ls(all=TRUE))
#setwd("...")
source("mcClustExtNew.R")
load("Njki.RData")
load("LMEntryDataMatrix.RData") 
attach(LMEntryDataMatrix)
load("mccXiPrior.RData")
load("Initial_values_for_classification_H4.RData")
load("Initial_values_for_betas_H4.RData")
# ========================================================================================================
groupNr <- 4
dirToStore <- "myDir"
# ========================================================================================================
try(mcClustExtNew(      # parameter lists (all four) must be complete!!!
     Data=list(dataFile=Njk.i, 
               storeDir=dirToStore,
               priorFile= mccXiPrior,
               X = cbind( intercept=1, unEmplRDist, unskilled, skilled, whiteColl, 
                                    wageCat1Dummy, wageCat2Dummy, wageCat3Dummy, wageCat4Dummy, wageCat5Dummy,
                                    entryYear76, entryYear77, entryYear78, entryYear79, entryYear80, 
                                    entryYear81, entryYear82, entryYear83, entryYear84, entryYear85,
                                    ia.ueRD.wc1D, ia.ueRD.wc2D, ia.ueRD.wc3D, ia.ueRD.wc4D, ia.ueRD.wc5D
                        ) ),
     Prior=list(H=groupNr, 
                c=1,
                cOff=1,
                usePriorFile=TRUE,
                xiPooled=TRUE,
                N0=10,
                betaPrior = "informative", # N(0,1)
                betaPriorMean = 0,
                betaPriorVar = 1),
     Initial=list(xi.start.ind=3, 
                  pers=0.7,
                  S.i.start = InitValClass,
                  Beta.start = InitValBetas ), 
     Mcmc=list(M=15000,
               M0=10000,
               mOut=500,
               mSave=5000,
               seed=3546541) 
     ))
detach(LMEntryDataMatrix)
# ========================================================================================================
