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