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 / save11.gen < prev    next >
Text File  |  1995-02-16  |  2KB  |  60 lines
  1. units [36]
  2. set [case=ignore]
  3. factor [levels=6] t
  4. factor [levels=3] a
  5. factor [levels=2] b
  6.  
  7.  
  8.  
  9. PROCEDURE 'expt1'
  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('   Completely randomised 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. open 'sess1.dat1'; chan=2
  21. read [print=*; channel=2] unitno,t,a,b,te,pe
  22. close chan=2
  23. random [seed=seed] t,a,b,te
  24. calc y=te+pe
  25.  
  26. text lm; val=!t('  ')
  27. text lh; val=!t('--------------------------------------------------------\
  28. ------')
  29. text expt; val=!t('             Experimental plan  :  6 treatments')
  30. text expab; val=!t('             Experimental plan  :  3 X 2 factorial')
  31. text expr; val=!t('                       Experimental results')
  32.  
  33.  
  34. for i=1
  35. tabulate [class=col,row] !(#t),!(#a),!(#b),y; max=tmax,amax,bmax,ymax
  36. print [ip=*; ind=4; sq=y]lh
  37. print [ip=*; ind=4; sq=y]expt
  38. print [ip=*; ind=4; sq=y]lh
  39. print [ip=*;rl=*;cl=*;ind=4;sq=y] tmax; dec=0; field=10
  40. print [ip=*]lm
  41. print [ip=*; ind=4; sq=y]lh
  42. print [ip=*; ind=4; sq=y]expab
  43. print [ip=*; ind=4; sq=y]lh
  44. print [ip=*;rl=*;cl=*;ind=4;sq=y] amax,bmax; dec=0; field=8,2
  45. question [preamble=!t('  ',' Performing experiment'); prompt=!t(\
  46.   ' Press return for results...'); response=_resp; mode=p; default=cairn]
  47. print [ip=*; ind=4; sq=y]lh
  48. print [ip=*; ind=4; sq=y]expr
  49. print [ip=*; ind=4; sq=y]lh
  50. print [ip=*;rl=*;cl=*;ind=4;sq=y] ymax; dec=3; field=10
  51. endfor
  52. ENDPROCEDURE
  53.  
  54. open 'courseprob.lib'; chan=3; file=back
  55. store [chan=3; proc=yes; subfile=expt1; meth=over] expt1
  56. expt1 [seed=43251] t;a;b;v1
  57. print t,a,b,v1
  58. stop
  59.  
  60.