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