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 / save24.gen < prev    next >
Text File  |  1995-02-16  |  3KB  |  96 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. factor [levels=3; values=12(1,2,3)] block
  7. factor [levels=6; values=2(1...6)3] mainplot
  8. factor [levels=2; values=(1,2)18] subplot
  9.  
  10. PROCEDURE 'expt6'
  11. OPTION NAME='SEED'; MODE=p
  12. PARAMETER NAME='T','A','B','C','BLOCK','MAINPLOT','SUBPLOT','Y'; \
  13.                                          MODE=p,p,p,p,p,p
  14. scalar cairn
  15. set [case=ignore]
  16. print !t('   Split plot design');
  17.  
  18. factor [levels=6;val=(1...6)6] row
  19. factor [levels=6;val=6(1...6)] col
  20.  
  21. open 'sess2.dat4'; chan=2
  22. read [print=*; channel=2] unitno,t,a,b,c,te,pe
  23. close chan=2
  24.  
  25.  
  26.  
  27. random [seed=seed; blocks=block/mainplot/subplot] t,a,b,c,te
  28. calc y=te+pe
  29.  
  30. text lm; val=!t('  ')
  31. text lh; val=!t('--------------------------------------------------------\
  32. ------')
  33. text lb; \
  34. val=!t('           .         .         .         .         .      ')
  35. text expt; val=!t('  Plot layout : 3 blocks of 6 main plots of 2 subplots')
  36. text t1; \
  37. val=!t('   M1-sub1   M4-sub1   M1-sub1   M4-sub1   M1-sub1   M4-sub1')
  38. text t2; \
  39. val=!t('   M1-sub2   M4-sub2   M1-sub2   M4-sub2   M1-sub2   M4-sub2')
  40. text t3; \
  41. val=!t('   M2-sub1   M5-sub1   M2-sub1   M5-sub1   M2-sub1   M5-sub1')
  42. text t4; \
  43. val=!t('   M2-sub2   M5-sub2   M2-sub2   M5-sub2   M2-sub2   M5-sub2')
  44. text t5; \
  45. val=!t('   M3-sub1   M6-sub1   M3-sub1   M6-sub1   M3-sub1   M6-sub1')
  46. text t6; \
  47. val=!t('   M3-sub2   M6-sub2   M3-sub2   M6-sub2   M3-sub2   M6-sub2')
  48. text expb; \
  49. val=!t('          Block 1      |      Block 2      |      Block 3     ')
  50.  
  51. text expab; val=!t('          Experimental plan  :  3 X 2 X 2 factorial')
  52. text expr; val=!t('                       Experimental results')
  53.  
  54.  
  55. for i=1
  56. tabulate [class=row,col] !(#t),!(#a),!(#b),!(#c),y; \
  57.                                     max=tmax,amax,bmax,cmax,ymax
  58.  
  59. print [ip=*; ind=4; sq=y]expt
  60. print [ip=*; ind=4; sq=y]lh
  61. & expb
  62. & lh
  63. print [ip=*; ind=4; sq=y]t1
  64. & t2
  65. print [ip=*; ind=4; sq=y]lb
  66. & t3
  67. & t4
  68. & lb
  69. & t5
  70. & t6
  71. print [ip=*; ind=4; sq=y]lh
  72.  
  73. print !t(' Randomising......')
  74.  
  75. print [ip=*; ind=4; sq=y]lh
  76. print [ip=*; ind=4; sq=y]expab
  77. print [ip=*; ind=4; sq=y]lh
  78. print [ip=*;rl=*;cl=*;ind=4;sq=y] amax,bmax,cmax; dec=0; field=6,2,2
  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=expt6; meth=over] expt6
  90. expt6 [seed=43251] t;a;b;c;block;mainplot;subplot;v1
  91. block block/mainplot/subplot
  92. treat a*b*c
  93. anova v1
  94. stop
  95.  
  96.