Changeset b41310e5 in tol


Ignore:
Timestamp:
Aug 24, 2009, 8:30:00 AM (17 years ago)
Author:
Víctor de Buen <vdebuen@…>
Branches:
1941-docker-image-to-test, devel, master
Children:
bd637b5
Parents:
a0c0550a
Message:

Reordering code

git-svn-id: ​https://www.tol-project.org/svn/tolp/trunk@1239 be4846ae-f051-483a-831f-b64e7d08527f

File:
1 edited

Legend:

Unmodified
Added
Removed
  • tol_tests/tol/BSR/test_13/test.tol

    ra0c0550a rb41310e5  
    1111
    1212Real numErr0 = Copy(NError);
     13
     14BysMcmc::Config bsr.config =
     15[[
     16//MCMC dimensions
     17  Real mcmc.burnin = 100;
     18  Real mcmc.sampleLength = 2000;
     19  Real mcmc.cacheLength = 100;
     20
     21//Partial sampling
     22  Real partialSampling.Sequential = False;
     23
     24//Basic master configuration
     25  Real basic.cholesky.epsilon = 1.E-13;
     26  Real basic.cholesky.warningFreq = 100;
     27  Real basic.truncMNormal.gibbsNumIter = 5;
     28  Real bsr.iniVal.try.givenByUser = 1;
     29  Real bsr.iniVal.try.zero = 2;
     30  Real bsr.iniVal.try.cnstrMinNorm = 3;
     31  Real bsr.iniVal.try.cnstrMLE = 4;
     32
     33  Text bsr.arimaFilter = BysMcmc::Options::Arima.Filter::FastChol;
     34  Real bsr.arimaSkipIter = 1;
     35
     36  Real bsr.linBlk.MH.useAfterIter   = 10;
     37  Real bsr.linBlk.MH.maxRejected    = 3;
     38  Real bsr.linBlk.MH.maxConsecutive = 100;
     39
     40  Real bsr.linBlk.X.drop = 1.E-15;
     41  Real bsr.linBlk.usePreconditioning = 0;
     42
     43//Report configuration
     44  Real report.raftery.diag.q   = 0.25;
     45  Real report.raftery.diag.r   = 0.07;   
     46  Real report.raftery.diag.s   = 0.90;   
     47  Real report.raftery.diag.eps = 0.01;
     48  Real report.acf.lag          = 20,
     49  Real report.histogram.parts  = 100;
     50  Real report.kerDens.points   = 0;
     51//Generic flags
     52  Real do.resume          = 0;
     53  Real do.report          = True;
     54  Real do.eval            = True;
     55  Real do.linear.effects  = False;
     56  Real do.mcmc.profiler   = True
     57]];
     58
     59
     60BysMcmc::McmcRecover recover = BysMcmc::Get.Recover(
     61  "../test_11/", "Test11", "Test11.01",
     62  bsr.config);
     63
    1364
    1465Set parsed = BSR.Parse("model.bsr");
    … …  
    117168
    118169
    119 BysMcmc::Config bsr.config =
    120 [[
    121 //MCMC dimensions
    122   Real mcmc.burnin = 100;
    123   Real mcmc.sampleLength = 2000;
    124   Real mcmc.cacheLength = 100;
    125 
    126 //Partial sampling
    127   Real partialSampling.Sequential = False;
    128 
    129 //Basic master configuration
    130   Real basic.cholesky.epsilon = 1.E-13;
    131   Real basic.cholesky.warningFreq = 100;
    132   Real basic.truncMNormal.gibbsNumIter = 5;
    133   Real bsr.iniVal.try.givenByUser = 1;
    134   Real bsr.iniVal.try.zero = 2;
    135   Real bsr.iniVal.try.cnstrMinNorm = 3;
    136   Real bsr.iniVal.try.cnstrMLE = 4;
    137 
    138   Text bsr.arimaFilter = BysMcmc::Options::Arima.Filter::FastChol;
    139   Real bsr.arimaSkipIter = 1;
    140 
    141   Real bsr.linBlk.MH.useAfterIter   = 10;
    142   Real bsr.linBlk.MH.maxRejected    = 3;
    143   Real bsr.linBlk.MH.maxConsecutive = 100;
    144 
    145   Real bsr.linBlk.X.drop = 1.E-15;
    146   Real bsr.linBlk.usePreconditioning = 0;
    147 
    148 //Report configuration
    149   Real report.raftery.diag.q   = 0.25;
    150   Real report.raftery.diag.r   = 0.07;   
    151   Real report.raftery.diag.s   = 0.90;   
    152   Real report.raftery.diag.eps = 0.01;
    153   Real report.acf.lag          = 20,
    154   Real report.histogram.parts  = 100;
    155   Real report.kerDens.points   = 0;
    156 //Generic flags
    157   Real do.resume          = 0;
    158   Real do.report          = True;
    159   Real do.eval            = True;
    160   Real do.linear.effects  = False;
    161   Real do.mcmc.profiler   = True
    162 ]];
    163 
    164 
    165 BysMcmc::McmcRecover recover = BysMcmc::Get.Recover(
    166   "../test_11/", "Test11", "Test11.01",
    167   bsr.config);
    168 
    169170Set fixedParam = Copy(Empty);
    170171VMatrix fixedValues = Rand(0,1,0,0);
    … …  
    188189Real cycler::_.sampler::_.arm.blk::refreshFrequency := 0;
    189190
    190 //Real cycler::initialize(0);
    191 
     191Real cycler::initialize(0);
     192
     193/*
    192194BysMcmc::Estimation estim = BysMcmc::Estim(cycler);
    193195
Note: See TracChangeset for help on using the changeset viewer.