Hi all!
I've installed psignifit3 in a Dell xps 501 running under Ubuntu 11.04. When
I ran "sudo make test" in the command line, the output said: "73 tests
failed" (the details are below).
On the other hand, when I ran "test_psignifit" within matlab (64bit 7.10.0
R2010a), the output was a fatal error:
>> test_psignifit
testing 2afc setting
MapEstimate
??? Error: Unexpected MATLAB operator.
Error in ==> MapEstimate at 156
eval ( output );
Error in ==> test_psignifit at 9
mapest = MapEstimate ( data2afc, priors );
Any help? Thanks in advance!!!
Make test output
[ OK ] PsiGumbelL 1st derivative value: 0.00456628
[ OK ] PsiGumbelL 2nd derivative value: -0.0291742
[ OK ] PsiGumbelR 1st derivative value: 0.118205
[ OK ] PsiGumbelR 2nd derivative value: -0.102208
[ OK ] PsiLogistic 1st derivative value: 0.104994
[ OK ] PsiLogistic 2nd derivative value: -0.0799625
[ OK ] Psychometric function 1st derivative w.r.t. prm 0 value: -0.0528291
[ OK ] Psychometric function 1st derivative w.r.t. prm 1 value: 0.0704388
[ OK ] Psychometric function 1st derivative w.r.t. prm 2 value: -0.208609
[ OK ] Psychometric function 2nd derivative w.r.t. prm 0 and 0 value:
0.0205253
[ OK ] Psychometric function 2nd derivative w.r.t. prm 0 and 1 value:
0.00785238
[ OK ] Psychometric function 2nd derivative w.r.t. prm 0 and 2 value:
0.110061
[ OK ] Psychometric function 2nd derivative w.r.t. prm 1 and 0 value:
0.00785238
[ OK ] Psychometric function 2nd derivative w.r.t. prm 1 and 1 value:
-0.0574291
[ OK ] Psychometric function 2nd derivative w.r.t. prm 1 and 2 value:
-0.146748
[ OK ] Psychometric function 2nd derivative w.r.t. prm 2 and 0 value:
0.110061
[ OK ] Psychometric function 2nd derivative w.r.t. prm 2 and 1 value:
-0.146748
[ OK ] Psychometric function 2nd derivative w.r.t. prm 2 and 2 value: 0
[ OK ] psi function at x=0.5 value: -1.96351
[ OK ] psi function at x=5.5 value: 1.61109
[ OK ] psi function at x=10.5 value: 2.303
[ OK ] psi function at x=15.5 value: 2.70824
[ OK ] psi function at x=20.5 value: 2.99584
[ OK ] psi function at x=25.5 value: 3.21894
[ OK ] psi function at x=30.5 value: 3.40124
[ OK ] psi function at x=35.5 value: 3.55538
[ OK ] psi function at x=40.5 value: 3.68891
[ OK ] psi function at x=45.5 value: 3.80668
[ OK ] psi function at x=50.5 value: 3.91204
[ OK ] digamma function at x=0.5 value: 4.9348
[ OK ] digamma function at x=5.5 value: 0.199342
[ OK ] digamma function at x=10.5 value: 0.099917
[ OK ] digamma function at x=15.5 value: 0.066642
[ OK ] digamma function at x=20.5 value: 0.0499896
[ OK ] digamma function at x=25.5 value: 0.0399947
[ OK ] digamma function at x=30.5 value: 0.0333302
[ OK ] digamma function at x=35.5 value: 0.0285695
[ OK ] digamma function at x=40.5 value: 0.0249987
[ OK ] digamma function at x=45.5 value: 0.0222213
[ OK ] digamma function at x=50.5 value: 0.0199993
=============> Solutions of optimizer <============
[ OK ] OptimizerSolution 2AFC alpha value: 3.29751
[ OK ] OptimizerSolution 2AFC beta value: 0.959018
[ OK ] OptimizerSolution 2AFC lambda value: 0.0191278
[ OK ] OptimizerSolution 2AFC deviance value: 3.9847
[ OK ] OptimizerSolution 2AFC deviance sum value: 3.9847
[ OK ] OptimizerSolution 2AFC Rpd value: 0.155845
[ OK ] OptimizerSolution 2AFC Rkd value: -0.319545
[FAIL] OptimizerSolution Y/N alpha value: 3.42566 should be: 3.43942
[FAIL] OptimizerSolution Y/N beta value: 1.02624 should be: 0.988357
[ OK ] OptimizerSolution Y/N lambda value: 0.000437318
[ OK ] OptimizerSolution Y/N gamma value: 0.0240525
[FAIL] OptimizerSolution Y/N deviance value: 2.21098 should be: 2.08172
[ OK ] OptimizerSolution Y/N deviance sum value: 2.21098
[FAIL] OptimizerSolution Y/N Rpd value: 0.339664 should be: 0.217146
[FAIL] OptimizerSolution Y/N Rkd value: -0.350166 should be: -0.477967
[ OK ] Optimizer Solution gamma=lambda, alpha value: 3.32555
[ OK ] Optimizer Solution gamma=lambda, beta value: 1.06769
[FAIL] Optimizer Solution gamma=lambda, lambda value: 0.000327229 should
be: 0.000506492
=============> Bootstrap properties <============
[ OK ] Acceleration constant (threshold) value: 0.00033244 < 0.018662
[FAIL] Bias (threshold) value: -0.0677307 should be: -0.118085
[FAIL] th(.1) value: 2.65266 should be: 2.58437
[FAIL] th(.9) value: 3.89757 should be: 3.83427
[ OK ] Acceleration constant (slope) value: 0.00033244
[ OK ] Bias (slope) value: -0.0401168
[ OK ] sl(.1) value: 0.172577
[FAIL] sl(.9) value: 0.497512 should be: 0.545639
[FAIL] Deviance limits value: 9.67016 should be: 8.95713
[ OK ] Rpd( 2.5%) value: -0.532793
[FAIL] Rpd(97.5%) value: 0.632072 should be: 0.495472
[ OK ] Rkd( 2.5%) value: -0.922931
[ OK ] Rkd(97.5%) value: 0.6258
[ OK ] influential 0 good
[ OK ] outliers 0 good
[ OK ] influential 1 good
[ OK ] outliers 1 good
[ OK ] influential 2 good
[ OK ] outliers 2 good
[ OK ] influential 3 good
[ OK ] outliers 3 good
[ OK ] influential 4 good
[ OK ] outliers 4 good
[ OK ] influential 5 good
[ OK ] outliers 5 good
=============> Properties of sigmoids <============
[ OK ] Phi(0) value: 0.5
[ OK ] invPhi(0.5) value: 0
[ OK ] invPhi(Phi(0.3)) value: 0.3
[ OK ] Phi(invPhi(0.3)) value: 0.3
[ OK ] PsiLogistic->f( 0) value: 0.5
[ OK ] PsiLogistic->f(-3) value: 0.0474259 < 0.05
[ OK ] PsiLogistic->f( 3) value: 0.952574 > 0.95
[ OK ] PsiLogistic( 3)-PsiLogistic(-3) value: 0.952574
[ OK ] PsiLogistic->df(3) value: 0.0451767 < 0.05
[ OK ] PsiLogistic->df(-3) value: 0.0451767 < 0.05
[ OK ] PsiLogistic monotonically increasing value: 0.00664806 > 0
[ OK ] PsiLogistic->ddf(3) value: -0.0408916 < 0
[ OK ] PsiLogistic->ddf(-3) value: 0.0408916 > 0
[ OK ] PsiGauss->f( 0) value: 0.5
[ OK ] PsiGauss->f(-3) value: 0.0013499 < 0.01
[ OK ] PsiGauss->f( 3) value: 0.99865 > 0.99
[ OK ] PsiGauss( 3)-PsiGauss(-3) value: 0.99865
[ OK ] PsiGaussian monotonically increasing value: 1.48672e-06 > 0
[ OK ] PsiGauss->df(3) value: 0.00443185 < 0.01
[ OK ] PsiGauss->df(-3) value: 0.00443185 < 0.01
[ OK ] PsiGauss->ddf(3) value: -0.0132955 < 0
[ OK ] PsiGauss->ddf(-3) value: 0.0132955 > 0
[ OK ] PsiGumbelL->f(0) value: 0.632121
[ OK ] PsiGumbelL->f(0) value: 0.048568
[ OK ] PsiGumbelL->f(0) value: 1
[ OK ] PsiGumbelL( 3 )-PsiGumbelL (-3 ) value: 1 > 0.951432
[ OK ] PsiGumbelL monotonically increasing value: 5.20543e-63 > 0
[ OK ] PsiGumbelL->df(3) value: 3.80054e-08 < 0.01
[ OK ] PsiGumbelL->df(-3) value: 0.047369 < 0.05
[ OK ] PsiGumbelL->ddf(3) value: -7.25354e-07 < 0
[ OK ] PsiGumbelL->ddf(-3) value: 0.0450106 > 0
[ OK ] PsiGumbelR->f(0) value: 0.367879
[ OK ] PsiGumbelR->f(0) value: 1.89218e-09
[ OK ] PsiGumbelR->f(0) value: 0.951432
[ OK ] PsiGumbelR( -3 )-PsiGumbelR (3 ) value: 1 > 0.951432
[ OK ] PsiGumbelR monotonically increasing value: 5.20543e-63 > 0
[ OK ] PsiGumbelR->df(3) value: 0.047369 < 0.05
[ OK ] PsiGumbelR->df(-3) value: 3.80054e-08 < 0.01
[ OK ] PsiGumbelR->ddf(3) value: -0.0450106 < 0
[ OK ] PsiGumbelR->ddf(-3) value: 7.25354e-07 > 0
[ OK ] PsiCauchy->f(0) value: 0.5
[ OK ] PsiCauchy->f(-z(0.1)*(-.5) ) value: 0.1
[ OK ] PsiCauchy->f( z(0.9)*.5 ) value: 0.9
[ OK ] PsiCauchy monotonically increasing value: 0.0122427 > 0
[ OK ] PsiCauchy->df(3) value: 0.031831 < 0.05
[ OK ] PsiCauchy->df(-3) value: 0.031831 < 0.05
[ OK ] PsiCauchy->ddf(3) value: -0.0190986 < 0
[ OK ] PsiCauchy->ddf(-3) value: 0.0190986 > 0
[ OK ] PsiId->f(-1) value: -1
[ OK ] PsiId->f(0) value: 0
[ OK ] PsiId->f( 1) value: 1
[ OK ] PsiId->df(-1) value: 1
[ OK ] PsiId->df( 0) value: 1
[ OK ] PsiId->df( 1) value: 1
[ OK ] PsiId->ddf(-1) value: 0
[ OK ] PsiId->ddf( 0) value: 0
[ OK ] PsiId->ddf( 1) value: 0
[ OK ] PsiId->inv(PsiId->f(-1)) value: -1
[ OK ] PsiId->inv(PsiId->f( 0)) value: 0
[ OK ] PsiId->inv(PsiId->f( 1)) value: 1
=============> Tests of core objects <============
[ OK ] abCore at threshold value: 0
[ OK ] abCore derivative stimulus value: 0.5
[ OK ] abCore derivative 0 value: -0.5
[ OK ] abCore derivative 1 value: -0
[ OK ] abCore 2nd derivative 0,0 value: 0
[ OK ] abCore 2nd derivative 0,1 value: 0.25
[ OK ] abCore 2nd derivative 1,1 value: 0
[ OK ] abCore inversion g(inv(2)) value: 2
[ OK ] abCore inversion inv(g(2)) value: 2
[ OK ] abCore inversion dinv(2,0) value: 1
[ OK ] abCore inversion dinv(2,1) value: 2
[ OK ] mwCore at threshold value: 0
[ OK ] mwCore derivative stimulus value: 2.19722
[ OK ] mwCore derivative 0 value: -2.19722
[ OK ] mwCore derivative 1 value: -0
[ OK ] mwCore 2nd derivative 0,0 value: 0
[ OK ] mwCore 2nd derivative 0,1 value: 1.09861
[ OK ] mwCore 2nd derivative 1,1 value: 0
[ OK ] mwCore inversion g(inv(2)) value: 2
[ OK ] mwCore inversion inv(g(2)) value: 2
[ OK ] mwCore inversion dinv(2,0) value: 1
[ OK ] mwCore inversion dinv(2,1) value: 0.45512
[ OK ] mwCore (m) for PsiLogistic value: 0.5
[ OK ] mwCore (w) for PsiLogistic value: 2
[ OK ] mwCore (m) for PsiGauss value: 0.5
[ OK ] mwCore (w) for PsiGauss value: 2
[ OK ] mwCore (m) for PsiGumbelL value: 0.5
[ OK ] mwCore (w) for PsiGumbelL value: 2
[ OK ] mwCore (m) for PsiCauchy value: 0.5
[ OK ] mwCore (w) for PsiCauchy value: 2
[ OK ] mwCore (m) for PsiExponential value: 0.5
[ OK ] mwCore (w) for PsiExponential value: 2
[ OK ] mwCore (m) for PsiGumbelR value: 0.5
[ OK ] mwCore (w) for PsiGumbelR value: 2
[ OK ] linearCore at threshold value: 0
[ OK ] linearCore derivative stimulus value: 3
[ OK ] linearCore derivative(0) at threshold value: -0.666667
[ OK ] linearCore derivative(1) at threshold value: 1
[ OK ] linearCore derivative(0,0) at threshold value: 0
[ OK ] linearCore derivative(0,0) at threshold value: 0
[ OK ] linearCore derivative(0,0) at threshold value: 0
[ OK ] linearCore 2nd derivative symmetry at threshold value: 0
[ OK ] linearCore inverse threshold value: -0.666667
[ OK ] linearCore inverse g(inv(2)) value: 2
[ OK ] linearCore inverse inv(g(2)) value: 2
[ OK ] linearCore inverse derivative(0) value: 0
[ OK ] linearCore inverse derivative(1) value: -0.333333
[ OK ] logCore at threshold value: 0
[ OK ] linearCore derivative stimulus value: 1
[ OK ] logCore derivative(0) at threshold value: -0.666667
[ OK ] logCore derivative(1) at threshold value: 1
[ OK ] logCore derivative(0,0) at threshold value: 0
[ OK ] logCore derivative(1,0) at threshold value: 0
[ OK ] logCore derivative(1,1) at threshold value: 0
[ OK ] logCore 2nd derivative symmetry at threshold value: 0
[ OK ] logCore inverse value: 0.513417
[ OK ] logCore inversion g(inv(2)) value: 2
[ OK ] logCore inversion inv(g(2)) value: 2
[ OK ] logCore inversion dinv(2,0) value: 0
[ OK ] logCore inversion dinv(2,1) value: -0.333333
[ OK ] NakaRushton at threshold value: 0.5
[ OK ] NakaRushton derivative stimulus value: 0.24
[ OK ] NakaRushton inverse value: 4
[ OK ] NakaRushton inversion inv(g(2)) value: 2
[ OK ] NakaRushton inversion g(inv(0.5)) value: 0.5
[ OK ] NakaRushton inversion dinv(2,0) value: 0.5
[ OK ] NakaRushton inversion dinv(2,1) value: 0.693147
=============> MCMC <============
[ OK ] Hybrid MCMC alpha value: 3.5093
[ OK ] Hybrid MCMC beta value: 0.941225
[ OK ] Hybrid MCMC lambda value: 0.0242945
[ OK ] Metropolis Hastings alpha value: 3.13955
[ OK ] Metropolis Hastings beta value: 1.03687
[ OK ] Metropolis Hastings lambda value: 0.0199283
[ OK ] Generic Metropolis MCMC alpha value: 3.2443
[ OK ] Generic Metropolis MCMC beta value: 0.980343
[ OK ] Generic Metropolis MCMC lambda value: 0.0276917
=============> Priors <============
[ OK ] Flat prior at 0 value: 1
[ OK ] Flat prior derivative at 0 value: 0
[ OK ] Uniform prior at -0.5 value: 0
[ OK ] Uniform prior at 0.5 value: 1
[ OK ] Uniform prior at 1.5 value: 0
[ OK ] Uniform prior derivative at -0.5 value: 0
[ OK ] Uniform prior derivative at 0.5 value: 0
[ OK ] Uniform prior derivative at 1.5 value: 0
[ OK ] Gaussian prior at -1 value: 0.241971
[ OK ] Gaussian prior at 0 value: 0.398942
[ OK ] Gaussian prior at 1 value: 0.241971
[ OK ] Gaussian prior derivative at -1 value: 0.241971
[ OK ] Gaussian prior derivative at 0 value: -0
[ OK ] Gaussian prior derivative at 1 value: -0.241971
[ OK ] BetaPrior at 0 value: 0
[ OK ] BetaPrior at 0.1 value: 1.68095
[ OK ] BetaPrior at 0.5 value: 1.1601
[ OK ] BetaPrior at 1.1 value: 0
[ OK ] BetaPrior derivative at 0 value: 0
[ OK ] BetaPrior derivative at 0.1 value: 12.1402
[ OK ] BetaPrior derivative at 0.5 value: 5.80049
[ OK ] BetaPrior derivative at 1.1 value: 0
[ OK ] GammaPrior at -0.5 value: 0
[ OK ] GammaPrior at 0.5 value: 0.12998
[ OK ] GammaPrior at 1.0 value: 0.1556
[ OK ] GammaPrior at 1.5 value: 0.161314
[ OK ] GammaPrior derivative at -0.5 value: 0
[ OK ] GammaPrior derivative at 0.5 value: 0.0866532
[ OK ] GammaPrior derivative at 1.0 value: 0.0259333
[ OK ] GammaPrior derivative at 1.5 value: 0
[ OK ] nGammaPrior at 0.5 value: 0
[ OK ] nGammaPrior at -0.5 value: 0.12998
[ OK ] nGammaPrior at -1.0 value: 0.1556
[ OK ] nGammaPrior at -1.5 value: 0.161314
[ OK ] nGammaPrior derivative at 0.5 value: -0
[FAIL] nGammaPrior derivative at -0.5 value: -0.0866532 should be:
0.0866532
[FAIL] nGammaPrior derivative at -1.0 value: -0.0259333 should be:
0.0259333
[ OK ] nGammaPrior derivative at -1.5 value: -0
=============> Linear algebra routines <============
[ OK ] Inverse (0,0) value: 1.80489
[ OK ] Inverse (1,0) value: -0.677728
[ OK ] Inverse (2,0) value: 0.00780493
[ OK ] Inverse (0,1) value: -0.677728
[ OK ] Inverse (1,1) value: 1.20944
[ OK ] Inverse (2,1) value: 0.753826
[ OK ] Inverse (0,2) value: 0.00780493
[ OK ] Inverse (1,2) value: 0.753826
[ OK ] Inverse (2,2) value: 2.48652
[ OK ] solving Ax=b, x[0] value: 1.46603
[ OK ] solving Ax=b, x[1] value: -0.0730086
[ OK ] solving Ax=b, x[2] value: 0.384718
[ OK ] Cholesky (0,0) value: 0.866025
[ OK ] Cholesky (1,0) value: 0.600444
[ OK ] Cholesky (2,0) value: -0.184752
[ OK ] Cholesky (0,1) value: 0
[ OK ] Cholesky (1,1) value: 1.00969
[ OK ] Cholesky (2,1) value: -0.306102
[ OK ] Cholesky (0,2) value: 0
[ OK ] Cholesky (1,2) value: 0
[ OK ] Cholesky (2,2) value: 0.634168
[ OK ] LU (0,0) value: 0.75
[ OK ] LU (1,0) value: 0.693333
[ OK ] LU (2,0) value: -0.213333
[ OK ] LU (0,1) value: 0.52
[ OK ] LU (1,1) value: 1.01947
[ OK ] LU (2,1) value: -0.303165
[ OK ] LU (0,2) value: -0.16
[ OK ] LU (1,2) value: -0.309067
[ OK ] LU (2,2) value: 0.402168
[ OK ] Ax=b, b[0] value: 1
[ OK ] Ax=b, b[1] value: 0.5
[ OK ] Ax=b, b[2] value: 2.77556e-17
[ OK ] M should be symmetric value: 1
[ OK ] matrix scaling value: 1.5
[ OK ] matrix scaling value: 1.04
[ OK ] matrix scaling value: -0.32
[ OK ] matrix scaling value: 1.04
[ OK ] matrix scaling value: 2.76
[ OK ] matrix scaling value: -0.84
[ OK ] matrix scaling value: -0.32
[ OK ] matrix scaling value: -0.84
[ OK ] matrix scaling value: 1.06
[ OK ] pivot Ax=b, x[0] value: -5.24051
[ OK ] pivot Ax=b, x[1] value: 1.32526
[ OK ] pivot Ax=b, x[2] value: 0.334347
=============> Testing return bug in jackknifedata <============
[ OK ] optimizer hung value: 0
=============> Initial parameter heuristics <============
[ OK ] PsiPsychometric->getStart() for increasing data value: 0.608 > 0
[ OK ] PsiPsychometric->getStart() for decreasing data value: -0.202335 <
0
=============> Finding good starting values <============
[FAIL] yes-no: Starting value for alpha value: 3.33333 should be: 3.49558
[FAIL] yes-no: Starting value for beta value: 1.03704 should be: 0.898865
[FAIL] yes-no: Starting value for lambda value: 0.00185185 should be:
0.00555556
[FAIL] yes-no: Starting value for gamma value: 0.0166667 should be:
0.0444444
[ OK ] yes-no: minimum of alpha range value: 0
[ OK ] yes-no: maximum of alpha range value: 10
[ OK ] yes-no: minimum of alpha range value: 0
[ OK ] yes-no: maximum of alpha range value: 10
[FAIL] yes-no: minimum of beta range value: 2 should be: 0.45512
[FAIL] yes-no: maximum of beta range value: 10 should be: 2.2756
[FAIL] yes-no: minimum of beta range value: 2 should be: 0.45512
[FAIL] yes-no: maximum of beta range value: 10 should be: 2.2756
[FAIL] yes-no: minimum of lambda range value: 2 should be: 0
[ OK ] yes-no: maximum of lambda range value: 0.1
[FAIL] yes-no: minimum of lambda range value: 2 should be: 0
[ OK ] yes-no: maximum of lambda range value: 0.1
[FAIL] yes-no: minimum of gamma range value: 2 should be: 0
[ OK ] yes-no: maximum of gamma range value: 0.1
[FAIL] yes-no: minimum of gamma range value: 2 should be: 0
[ OK ] yes-no: maximum of gamma range value: 0.1
[FAIL] 2afc: Starting value for alpha value: 3.33333 should be: 3.29584
[FAIL] 2afc: Starting value for beta value: 0.888889 should be: 0.988751
[FAIL] 2afc: Starting value for lambda value: 0.0203704 should be:
0.0185185
[ OK ] 2afc: minimum of alpha range value: 0
[ OK ] 2afc: maximum of alpha range value: 10
[ OK ] 2afc: minimum of alpha range value: 0
[ OK ] 2afc: maximum of alpha range value: 10
[FAIL] 2afc: minimum of beta range value: 2 should be: 0.45512
[FAIL] 2afc: maximum of beta range value: 10 should be: 2.2756
[FAIL] 2afc: minimum of beta range value: 2 should be: 0.45512
[FAIL] 2afc: maximum of beta range value: 10 should be: 2.2756
[FAIL] 2afc: minimum of lambda range value: 2 should be: 0
[ OK ] 2afc: maximum of lambda range value: 0.1
[FAIL] 2afc: minimum of lambda range value: 2 should be: 0
[ OK ] 2afc: maximum of lambda range value: 0.1
[ OK ] grid parameter 0 lower limit value: 0
[ OK ] grid parameter 0 upper limit value: 0.5
[ OK ] grid parameter 1 lower limit value: 0
[ OK ] grid parameter 1 upper limit value: 0.5
[ OK ] grid parameter 2 lower limit value: 0
[ OK ] grid parameter 2 upper limit value: 0.05
[ OK ] grid.empty() on nonempty grid value: 0
[ OK ] grid.empty() on empty grid value: 1
[ OK ] grid.dimension() on 3d grid value: 3
[ OK ] grid.get_gridsize() on small grid value: 5
[ OK ] grid.front()[0] value: 0
[ OK ] grid.front()[1] value: 0.125
[ OK ] grid.front()[2] value: 0.25
[ OK ] grid.front()[3] value: 0.375
[ OK ] grid.front()[4] value: 0.5
[ OK ] shifted grid.front()[0] value: -0.125
[ OK ] shifted grid.front()[1] value: 0
[ OK ] shifted grid.front()[2] value: 0.125
[ OK ] shifted grid.front()[3] value: 0.25
[ OK ] shifted grid.front()[4] value: 0.375
[ OK ] shrunken grid.front()[0] value: 0
[ OK ] shrunken grid.front()[1] value: 0.0625
[ OK ] shrunken grid.front()[2] value: 0.125
[ OK ] shrunken grid.front()[3] value: 0.1875
[ OK ] shrunken grid.front()[4] value: 0.25
[ OK ] subgrid dimension value: 2
[ OK ] linspace 0 value: 0
[ OK ] linspace 1 value: 0.5
[ OK ] linspace 2 value: 1
[ OK ] Number of generated gridpoints from 2x2x2 grid value: 8
[ OK ] gridpoint 0 first param value: 0
[ OK ] gridpoint 0 second param value: 0
[ OK ] gridpoint 0 third param value: 0
[ OK ] gridpoint 1 first param value: 0
[ OK ] gridpoint 1 second param value: 0
[ OK ] gridpoint 1 third param value: 0.05
[ OK ] gridpoint 2 first param value: 0
[ OK ] gridpoint 2 second param value: 0.5
[ OK ] gridpoint 2 third param value: 0
[ OK ] gridpoint 3 first param value: 0
[ OK ] gridpoint 3 second param value: 0.5
[ OK ] gridpoint 3 third param value: 0.05
[ OK ] gridpoint 4 first param value: 0.5
[ OK ] gridpoint 4 second param value: 0
[ OK ] gridpoint 4 third param value: 0
[ OK ] gridpoint 5 first param value: 0.5
[ OK ] gridpoint 5 second param value: 0
[ OK ] gridpoint 5 third param value: 0.05
[ OK ] gridpoint 6 first param value: 0.5
[ OK ] gridpoint 6 second param value: 0.5
[ OK ] gridpoint 6 third param value: 0
[ OK ] gridpoint 7 first param value: 0.5
[ OK ] gridpoint 7 second param value: 0.5
[ OK ] gridpoint 7 third param value: 0.05
[FAIL] Best fit on grid value: 35.4381 should be: 23.8999
[FAIL] Second best fit on grid value: 43.346 should be: 32.0023
[FAIL] Best fitting alpha on grid value: 0.5 should be: 0
[ OK ] Best fitting beta on grid value: 0.5
[ OK ] Best fitting lambda on grid value: 0.05
[FAIL] Second best fitting alpha on grid value: 0 should be: 0.5
[ OK ] Second best fitting beta on grid value: 0.5
[ OK ] Second best fitting lambda on grid value: 0.05
[FAIL] gridpoint 0 first param value: 0.25 should be: -0.25
[ OK ] gridpoint 0 second param value: 0.25
[ OK ] gridpoint 0 third param value: 0.025
[FAIL] gridpoint 1 first param value: 0.25 should be: -0.25
[ OK ] gridpoint 1 second param value: 0.25
[ OK ] gridpoint 1 third param value: 0.075
[FAIL] gridpoint 2 first param value: 0.25 should be: -0.25
[ OK ] gridpoint 2 second param value: 0.75
[ OK ] gridpoint 2 third param value: 0.025
[FAIL] gridpoint 3 first param value: 0.25 should be: -0.25
[ OK ] gridpoint 3 second param value: 0.75
[ OK ] gridpoint 3 third param value: 0.075
[FAIL] gridpoint 4 first param value: 0.75 should be: 0.25
[ OK ] gridpoint 4 second param value: 0.25
[ OK ] gridpoint 4 third param value: 0.025
[FAIL] gridpoint 5 first param value: 0.75 should be: 0.25
[ OK ] gridpoint 5 second param value: 0.25
[ OK ] gridpoint 5 third param value: 0.075
[FAIL] gridpoint 6 first param value: 0.75 should be: 0.25
[ OK ] gridpoint 6 second param value: 0.75
[ OK ] gridpoint 6 third param value: 0.025
[FAIL] gridpoint 7 first param value: 0.75 should be: 0.25
[ OK ] gridpoint 7 second param value: 0.75
[ OK ] gridpoint 7 third param value: 0.075
=============> Approximate numerical integration <============
[FAIL] Lower bound of m value: -2.22222 should be: 0
[FAIL] Upper bound of m value: 8.88889 should be: 10
[ OK ] Length of m grid value: 7
[FAIL] Lower bound of w value: -0.538267 should be: 2
[FAIL] Upper bound of w value: 8.35062 should be: 10
[ OK ] Length of w grid value: 7
[FAIL] Lower bound of lm value: -0.0351852 should be: 0
[FAIL] Upper bound of lm value: 0.0759259 should be: 0.1
[ OK ] Length of lm grid value: 7
[ OK ] cdf_grid Gauss value: -1.28155
[ OK ] cdf_grid Gauss value: -0.727913
[ OK ] cdf_grid Gauss value: -0.340695
[ OK ] cdf_grid Gauss value: 4.27646e-18
[ OK ] cdf_grid Gauss value: 0.340695
[ OK ] cdf_grid Gauss value: 0.727913
[ OK ] cdf_grid Gauss value: 1.28155
[FAIL] fit_posterior Gauss mean value: -0.00547325 should be: 0
[FAIL] fit_posterior Gauss std value: 0.00115231 should be: 1
[FAIL] cdf_grid Gamma value: inf should be: 1.74477
[FAIL] cdf_grid Gamma value: -nan should be: 2.45729
[FAIL] cdf_grid Gamma value: -nan should be: 3.06205
[FAIL] cdf_grid Gamma value: -nan should be: 3.67206
[FAIL] cdf_grid Gamma value: -nan should be: 4.35885
[FAIL] cdf_grid Gamma value: -nan should be: 5.23689
[FAIL] cdf_grid Gamma value: -nan should be: 6.68078
[FAIL] fit_posterior Gamma shape value: 2 should be: 4
[FAIL] fit_posterior Gamma scale value: 1.815 should be: 1
[ OK ] cdf_grid Beta value: 0.0256173
[ OK ] cdf_grid Beta value: 0.0437104
[ OK ] cdf_grid Beta value: 0.0606308
[ OK ] cdf_grid Beta value: 0.0786438
[ OK ] cdf_grid Beta value: 0.0996593
[ OK ] cdf_grid Beta value: 0.127175
[ OK ] cdf_grid Beta value: 0.172935
[FAIL] fit_posterior Beta alpha value: 1.99882 should be: 2
[FAIL] fit_posterior Beta beta value: 19.9891 should be: 20
[FAIL] Posterior for m -- mu value: 3.39015 should be: 3.28
[FAIL] Posterior for m -- sg value: 0.451554 should be: 0.436
[FAIL] Posterior for w -- k value: 10.4441 should be: 12.732
[FAIL] Posterior for w -- th value: 0.440721 should be: 0.362
[FAIL] Posterior for lm -- al value: 2.54935 should be: 3.17
[FAIL] Posterior for lm -- bt value: 106.59 should be: 111.26
H = 0.963892
[ OK ] Sampled and fitted posterior mean for m value: 3.21063
[FAIL] Sampled and fitted posterior mean for w value: 4.58142 should be:
4.60294
[ OK ] Sampled and fitted posterior mean for lm value: 0.0270066
[ OK ] number of duplicates in SIR value: 0.04 < 0.1
[ OK ] Rkd is set value: 1
[ OK ] Rpd is set value: 1
[ OK ] logratio is set value: 1
[ OK ] ppRkd is set value: 1
[ OK ] ppRpd is set value: 1
[ OK ] ppData is set value: 6
[ OK ] Deviance is set value: 1
[ OK ] ppDeviance is set value: 1
==> 73 tests failed
--
Andres
|