MCMC is not taking the same starting values but for different datasets
Posted: Wed Oct 14, 2015 4:07 pm
** I started a new thread instead **
I am following up on this thread, but the question is merely related.
I am using starting values and passing both FP.b and RP.b to the model as you suggested. R2mlwing manages to run this model if I use a small random sample of the original data (~ 1000 cases), but not when I want to estimate the model of the full sample (~ 1.9 million cases).
I am getting this error code:
> m1test2 <- runMLwiN(logit(AbsolutDep, cons) ~ 1 + (1 | country) + (1 |CountryClusterHouse), D = "Binomial", estoptions = list(EstM = 1, resi.store=F,
+ debugmode=F, optimat=T,
+ mcmcMeth=list(iterations=10, burnin=10),
+ mcmcOptions=list(hcen=3),
+ startval=list(FP.b = PLOSONEestimations20000IGLS[[1]]@FP , RP.b = PLOSONEestimations20000IGLS[[1]]@RP)), data = test4a, workdir = tempdir(), MLwiNPath="C:/Program Files (x86)/MLwiN v2.35/")
MLwiN is running, please wait......
/nogui option ignored
ECHO 0
Echoing is ON
STAR
iteration 0
Convergence not achieved
JOIN -0.491145551204681 '_FP_b'
JOIN 0 '_FP_v'
JOIN 1.10714721679688 1.04919278621674 1 '_RP_b'
JOIN 0 0 0 0 0 0 '_RP_v'
ECHO 0
Echoing is ON
MCMC 0 10 1 5.8 50 10 G30[1] G30[2] 2 2 2 1 1 2
error while obeying batch file C:/Users/adel/AppData/Local/Temp/Rtmp4GO7Gq/macrofile_1d002647c62.txt at line number 139:
MCMC 0 10 1 5.8 50 10 C2498] C2499] 2 2 2 1 1 2
wrong length random constraint matrix
error while obeying batch file C:/Users/adel/AppData/Local/Temp/Rtmp4GO7Gq/macrofile_1d002647c62.txt at line number 139:
MCMC 0 10 1 5.8 50 10 C2498] C2499] 2 2 2 1 1 2
wrong length random constraint matrix
.
Execution completed
Error in read.dta(MCMCfile) :
unable to open file: 'No such file or directory'
Do you have any idea what is going on? I have tried to pass only the relevant variables, in case we are having a ram problem (my PC has 64 GB of ram), but that does not work either. I am clueless about what to try next.
I am hoping you might have any ideas about what to do next.
Many thanks in advance
PS I am not sure if it is related. But I am also observing a wobbling RAM consumption behaviour when I run this model (see picture please). I have not seen anything like it before. Can these two events be related?
AN UPDATE:
1. I ran a new IGLS with the a dataset to obtain a new set of starting values. I then used those to initiate an MCMC model. But I am still getting the same error. Namely:
MLwiN is running, please wait......
/nogui option ignored
ECHO 0
Echoing is ON
STAR
iteration 0
Convergence not achieved
JOIN -0.49 '_FP_b'
JOIN 0 '_FP_v'
JOIN 1.11 1.05 1 '_RP_b'
JOIN 0 0 0 0 0 0 '_RP_v'
ECHO 0
Echoing is ON
MCMC 0 10 1 5.8 50 10 G30[1] G30[2] 2 2 2 1 1 2
error while obeying batch file C:/Users/adel/AppData/Local/Temp/Rtmp4GO7Gq/macrofile_1d00c492552.txt at line number 139:
MCMC 0 10 1 5.8 50 10 C2498] C2499] 2 2 2 1 1 2
wrong length random constraint matrix
error while obeying batch file C:/Users/adel/AppData/Local/Temp/Rtmp4GO7Gq/macrofile_1d00c492552.txt at line number 139:
MCMC 0 10 1 5.8 50 10 C2498] C2499] 2 2 2 1 1 2
wrong length random constraint matrix
.
Execution completed
Error in read.dta(MCMCfile) :
unable to open file: 'No such file or directory'
2. I rounded the starting values before passing them to runMlwin – incase that was the issue. But the problem persists.
I am following up on this thread, but the question is merely related.
I am using starting values and passing both FP.b and RP.b to the model as you suggested. R2mlwing manages to run this model if I use a small random sample of the original data (~ 1000 cases), but not when I want to estimate the model of the full sample (~ 1.9 million cases).
I am getting this error code:
> m1test2 <- runMLwiN(logit(AbsolutDep, cons) ~ 1 + (1 | country) + (1 |CountryClusterHouse), D = "Binomial", estoptions = list(EstM = 1, resi.store=F,
+ debugmode=F, optimat=T,
+ mcmcMeth=list(iterations=10, burnin=10),
+ mcmcOptions=list(hcen=3),
+ startval=list(FP.b = PLOSONEestimations20000IGLS[[1]]@FP , RP.b = PLOSONEestimations20000IGLS[[1]]@RP)), data = test4a, workdir = tempdir(), MLwiNPath="C:/Program Files (x86)/MLwiN v2.35/")
MLwiN is running, please wait......
/nogui option ignored
ECHO 0
Echoing is ON
STAR
iteration 0
Convergence not achieved
JOIN -0.491145551204681 '_FP_b'
JOIN 0 '_FP_v'
JOIN 1.10714721679688 1.04919278621674 1 '_RP_b'
JOIN 0 0 0 0 0 0 '_RP_v'
ECHO 0
Echoing is ON
MCMC 0 10 1 5.8 50 10 G30[1] G30[2] 2 2 2 1 1 2
error while obeying batch file C:/Users/adel/AppData/Local/Temp/Rtmp4GO7Gq/macrofile_1d002647c62.txt at line number 139:
MCMC 0 10 1 5.8 50 10 C2498] C2499] 2 2 2 1 1 2
wrong length random constraint matrix
error while obeying batch file C:/Users/adel/AppData/Local/Temp/Rtmp4GO7Gq/macrofile_1d002647c62.txt at line number 139:
MCMC 0 10 1 5.8 50 10 C2498] C2499] 2 2 2 1 1 2
wrong length random constraint matrix
.
Execution completed
Error in read.dta(MCMCfile) :
unable to open file: 'No such file or directory'
Do you have any idea what is going on? I have tried to pass only the relevant variables, in case we are having a ram problem (my PC has 64 GB of ram), but that does not work either. I am clueless about what to try next.
I am hoping you might have any ideas about what to do next.
Many thanks in advance
PS I am not sure if it is related. But I am also observing a wobbling RAM consumption behaviour when I run this model (see picture please). I have not seen anything like it before. Can these two events be related?
AN UPDATE:
1. I ran a new IGLS with the a dataset to obtain a new set of starting values. I then used those to initiate an MCMC model. But I am still getting the same error. Namely:
MLwiN is running, please wait......
/nogui option ignored
ECHO 0
Echoing is ON
STAR
iteration 0
Convergence not achieved
JOIN -0.49 '_FP_b'
JOIN 0 '_FP_v'
JOIN 1.11 1.05 1 '_RP_b'
JOIN 0 0 0 0 0 0 '_RP_v'
ECHO 0
Echoing is ON
MCMC 0 10 1 5.8 50 10 G30[1] G30[2] 2 2 2 1 1 2
error while obeying batch file C:/Users/adel/AppData/Local/Temp/Rtmp4GO7Gq/macrofile_1d00c492552.txt at line number 139:
MCMC 0 10 1 5.8 50 10 C2498] C2499] 2 2 2 1 1 2
wrong length random constraint matrix
error while obeying batch file C:/Users/adel/AppData/Local/Temp/Rtmp4GO7Gq/macrofile_1d00c492552.txt at line number 139:
MCMC 0 10 1 5.8 50 10 C2498] C2499] 2 2 2 1 1 2
wrong length random constraint matrix
.
Execution completed
Error in read.dta(MCMCfile) :
unable to open file: 'No such file or directory'
2. I rounded the starting values before passing them to runMlwin – incase that was the issue. But the problem persists.