devel/scm_gof_backward_test.pl

   1 #!/usr/bin/perl
   2
   3 use lib qw(../lib);
   4 use gof_crc(test);
   5 use method::scm;
   6 use model;
   7
   8 my $debug = defined $ARGV[0] ? $ARGV[0] : 0;
   9
  10 my $model = model -> new('filename' => 'pheno.mod',
  11                          target     => 'disk' );
  12
  13 my $scm = method::scm -> new( parameters             => ['CL', 'V'],
  14                               models                 => [$model],
  15                               continuous_covariates   => ['WGT'],
  16                               categorical_covariates => ['APGR'],
  17                               fix                    => 0,
  18                               nm_version             => 6,
  19                               picky                  => 1,
  20                               logfile                => 'scmlog.txt',
  21                               threads                => 3,
  22                               search_direction       => 'backward',
  23                               gof                    => gof_crc::test,
  24                               abort_on_fail          => 0,
  25                               base_criteria_vals     =>
  26                               [{'ofv'       => 571.183,
  27                                 'estimates' => { 'IVV'  => 0.02350,
  28                                                  'IVCL' => 0.04350 },
  29                                 'stats'     => { 'TVCL' => { 'lower_relative_range' => 0.54042,
  30                                                              'upper_relative_range' => 1.77562 },
  31                                                  'TVV'  => { 'lower_relative_range' => 0.46670,
  32                                                              'upper_relative_range' => 1.90968 }}}],
  33                               included_relations     => [{'CL' => {'WGT'=>{'state'=>2}},
  34                                                           'V'  => {'WGT'=>{'state'=>2},
  35                                                                    'APGR'=>{'state'=>2}}}],
  36                               test_relations         => { 'CL' => ['WGT'],
  37                                                           'V'  => ['WGT','APGR'] } );
  38
  39
  40 $scm -> run;
  41