Scripts/Landsch.run

macro=; /* IMPORTANT: REMEMBER TO ERASE THE FILE all_trees.tre EACH TIME YOU RE-RUN THE ANALYSIS */ var: best_scor curr_scor rounds repls bescor bestre arbol escor level

sil=all ; lm lowmem ; macreport - ; report - ; tsave *all_trees.tre + ; k0 ; if ( argnumber > 0 )

   set level %1 ;

if ('level' == 0)

	 set repls 1 ;

lm notermpoints; lm cell 5 nest 0 ; lm iter ; end if ('level' == 1) set repls 3 ; lm notermpoints; lm cell 5 nest 1 2 ; lm iter ; end if ('level' == 2) set repls 3 ; lm termpoints; lm cell 5 nest 3 2 ; lm iter ; end if ('level' == 3) set repls 5 ; lm termpoints; lm cell 8 nest 4 2 ; lm iter ; end if ('level' == 4) set repls 5 ; lm termpoints; lm cell 10 nest 8 2 ; lm iter ;

       end

else

set repls 1 ; /* change here the number of replicates to do */ 

end loop 0 ('repls'-1)

 rseed 0 ;
 ra 1 ;
 set best_scor lmscore [ 0 ] ;
 macfloat 4 ;
 report - ; 
 set rounds 1 ; 
tbrit 0 
   progress percswap 100 Replicate #1 of 'repls' Score 'best_scor' ; 
   progress= ; 
   set rounds ++ ; 
   set curr_scor lmscore [ 0 ] ;
   if ( 'curr_scor' < 'best_scor' )
    set best_scor 'curr_scor' ;
    resetswap ;
   end
stop 
 save ; 
end
 sil=all ; 
 k0 ; 

stop ; ts/ ; p all_trees.tre ;

set bescor  100000 ;

loop=ar 0 ('repls'-1)

  set arbol #ar ; 
  set escor lmscore ['arbol'] ;
 if ('escor' < 'bescor' )
  set bescor 'escor' ;
  set bestre 'arbol' ;   
 end
 stop 
tsave *best_tree.tre ;
tch 'bestre' ;
save ;
ts/;         

lquote= ; lquote [ ;

sil-buff ; quote *Best tree replicate 'bestre' Score 'bescor' 
            Best tree kept in memory  ; 

progress/ ; sil=all ; report-; proc/;