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