home *** CD-ROM | disk | FTP | other *** search
/ gimli.bioss.sari.ac.uk / gimli.bioss.sari.ac.uk.zip / gimli.bioss.sari.ac.uk / pub / courses / ExptDesignAnalysis / save12.gen < prev    next >
Text File  |  1995-02-16  |  2KB  |  68 lines
  1. units [36]
  2. set [case=ignore]
  3. factor [levels=6] t
  4. factor [levels=3] a
  5. factor [levels=2] b
  6. factor [levels=6] block,plot
  7.  
  8.  
  9. PROCEDURE 'expt2'
  10. OPTION NAME='SEED'; MODE=p
  11. PARAMETER NAME='T','A','B','BLOCK','PLOT','Y'; MODE=p,p,p,p,p,p
  12. scalar cairn
  13. set [case=ignore]
  14. print !t('   Randomised Block design');
  15. print !t('   Randomising....');
  16.  
  17. factor [levels=6;val=6(1...6)] row
  18. factor [levels=6;val=(1...6)6] col
  19.  
  20. open 'sess1.dat2'; chan=2
  21. read [print=*; channel=2] unitno,t,a,b,te,pe
  22. close chan=2
  23.  
  24. calc block=!(6(1...6))
  25. calc plot=!((1...6)6)
  26.  
  27. random [seed=seed; blocks=block/plot] t,a,b,te
  28. calc y=te+pe
  29.  
  30. text lm; val=!t('  ')
  31. text lh; val=!t('--------------------------------------------------------\
  32. ------')
  33. text expt; val=!t('             Experimental plan  :  6 treatments')
  34. text expb; \
  35. val=!t('   Block 1   Block 2   Block 3   Block 4   Block 5   Block 6')
  36. text expab; val=!t('             Experimental plan  :  3 X 2 factorial')
  37. text expr; val=!t('                       Experimental results')
  38.  
  39.  
  40. for i=1
  41. tabulate [class=col,row] !(#t),!(#a),!(#b),y; max=tmax,amax,bmax,ymax
  42. print [ip=*; ind=4; sq=y]lh
  43. print [ip=*; ind=4; sq=y]expt
  44. print [ip=*; ind=4; sq=y]lh
  45. print [ip=*; ind=4; sq=y]expb
  46. print [ip=*; ind=4; sq=y]lh
  47. print [ip=*;rl=*;cl=*;ind=4;sq=y] tmax; dec=0; field=10
  48. print [ip=*]lm
  49. print [ip=*; ind=4; sq=y]lh
  50. print [ip=*; ind=4; sq=y]expab
  51. print [ip=*; ind=4; sq=y]lh
  52. print [ip=*;rl=*;cl=*;ind=4;sq=y] amax,bmax; dec=0; field=8,2
  53. question [preamble=!t('  ',' Performing experiment'); prompt=!t(\
  54.   ' Press return for results...'); response=_resp; mode=p; default=cairn]
  55. print [ip=*; ind=4; sq=y]lh
  56. print [ip=*; ind=4; sq=y]expr
  57. print [ip=*; ind=4; sq=y]lh
  58. print [ip=*;rl=*;cl=*;ind=4;sq=y] ymax; dec=3; field=10
  59. endfor
  60. ENDPROCEDURE
  61.  
  62. open 'courseprob.lib'; chan=3; file=back
  63. store [chan=3; proc=yes; subfile=expt2; meth=over] expt2
  64. expt2 [seed=43251] t;a;b;block;plot;v1
  65. print t,a,b,block,plot,v1
  66. stop
  67.  
  68.