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 / save42.gen < prev    next >
Text File  |  1995-02-16  |  3KB  |  93 lines
  1. units [36]
  2.  
  3. factor [levels=9] t
  4. factor [levels=3] a
  5. factor [levels=3] b
  6. factor [levels=4;val=9(1...4)] block
  7. factor [levels=3;val=3(1...3)4] subblock
  8. factor [levels=3;val=(1...3)12] plot
  9.  
  10.  
  11. PROCEDURE 'expt8'
  12. OPTION NAME='SEED'; MODE=p
  13. PARAMETER NAME='T','A','B','BLOCK','SUBBLOCK','PLOT','Y';\
  14.                       MODE=p,p,p,p,p,p
  15. scalar cairn
  16.  
  17. print !t('   Lattice design');
  18.  
  19.  
  20. factor [levels=6;val=((1,2,3)3,(4,5,6)3)2] row
  21. factor [levels=6;val=3(1,2,3)2,3(4,5,6)2] col
  22.  
  23. open 'sess4.dat2'; chan=2
  24. read [print=*; channel=2] unitno,T,A,B,te,pe
  25. close chan=2
  26.  
  27. calc A=newlevels(T;!(1,1,1,2,2,2,3,3,3))
  28. calc B=newlevels(T;!(1,2,3,1,2,3,1,2,3))
  29.  
  30. random [seed=SEED; blocks=BLOCK/SUBBLOCK/PLOT] T,A,B,te
  31. calc Y=te+pe
  32.  
  33. text lm; val=!t('  ')
  34. text lh; val=!t('--------------------------------------------------------\
  35. ------')
  36. text lb; \
  37. val=!t('           .         .         .         .         .      ')
  38. text expl; val=!t('  Plot layout : 4 blocks of 3 subblocks of 3 plots')
  39. text t1; \
  40. val=!t('   B1-sub1   B1-sub2   B1-sub3  |  B3-sub1   B3-sub2   B3-sub3')
  41. text t2; \
  42. val=!t('   B1-sub1   B1-sub2   B1-sub3  |  B3-sub1   B3-sub2   B3-sub3')
  43. text t3; \
  44. val=!t('   B1-sub1   B1-sub2   B1-sub3  |  B3-sub1   B3-sub2   B3-sub3')
  45. text t4; \
  46. val=!t('   B2-sub1   B2-sub2   B2-sub3  |  B4-sub1   B4-sub2   B4-sub3')
  47. text t5; \
  48. val=!t('   B2-sub1   B2-sub2   B2-sub3  |  B4-sub1   B4-sub2   B4-sub3')
  49. text t6; \
  50. val=!t('   B2-sub1   B2-sub2   B3-sub3  |  B4-sub1   B4-sub2   B4-sub3')
  51. text expb; \
  52. val=!t('   Block 1   Block 2   Block 3   Block 4   Block 5   Block 6')
  53. text expt; val=!t('             Experimental plan  :  9 treatments')
  54. text expr; val=!t('                       Experimental results')
  55.  
  56.  
  57. for i=1
  58. tabulate [class=row,col] !(#T),!(#A),!(#B),Y; max=tmax,amax,bmax,ymax
  59. print [ip=*; ind=4; sq=y]expl
  60. print [ip=*; ind=4; sq=y]lh
  61. print [ip=*; ind=4; sq=y]t1
  62. & t2
  63. & t3
  64. & lb
  65. & t4
  66. & t5
  67. & t6
  68. print [ip=*; ind=4; sq=y]lh
  69.  
  70. print !t(' Randomising......')
  71. print [ip=*]lm
  72. print [ip=*; ind=4; sq=y] expt
  73. print [ip=*; ind=4; sq=y]lh
  74. "
  75. print [ip=*; ind=4; sq=y]expb
  76. print [ip=*; ind=4; sq=y]lh
  77. "
  78. print [ip=*;rl=*;cl=*;ind=4;sq=y] tmax; dec=0; field=10
  79. question [preamble=!t('  ',' Performing experiment'); prompt=!t(\
  80.   ' Press return for results...'); response=_resp; mode=p; default=cairn]
  81. print [ip=*; ind=4; sq=y]lh
  82. print [ip=*; ind=4; sq=y]expr
  83. print [ip=*; ind=4; sq=y]lh
  84. print [ip=*;rl=*;cl=*;ind=4;sq=y] ymax; dec=3; field=10
  85. endfor
  86. ENDPROCEDURE
  87.  
  88. open 'courseprob.lib'; chan=3; file=back
  89. store [chan=3; proc=yes; subfile=expt8; meth=over] expt8
  90. expt8 [seed=43251] t;a;b;block;subblock;plot;v1
  91. stop
  92.  
  93.