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 / lambs.gen < prev    next >
Text File  |  1995-02-16  |  759b  |  44 lines
  1. units[250]
  2. factor[lev=10] treat
  3. open 'lambs.dat';ch=2
  4. read [ch=2] treat,w[1...8],td,wg
  5. scalar mvi
  6. calc mvi=0/0
  7.  
  8. for yy=w[5...8],wg;zz=mv[5...8],mvg
  9. restrict zz;cond=yy==mvi
  10. calc zz=2
  11. restrict zz; cond=yy.ne.mvi
  12. calc zz=1
  13. restrict zz
  14. endfor
  15.  
  16. calc g[2...8]=w[2...8]-w[1]
  17. factor [lev=!(0,1)]mf[5...8],mfg
  18. calc mf[5...8],mfg=mv[5...8],mvg-1
  19.  
  20. for ff=mf[6...8],mfg;vv=g[5...8];rr=mf[5...8]
  21. restrict vv;cond=rr==0
  22. tabulate[class=treat,ff;marg=y;prin=m,c]vv
  23. model vv
  24. fit[prin=m,s,a,e] treat+ff
  25. restrict vv
  26. endfor
  27.  
  28. vari rate
  29. restrict wg,td,rate;cond=mfg==0
  30. calc rate=wg/td
  31. restrict wg,td,rate;cond=mfg.ne.0
  32. calc rate=g[7]/70
  33. restrict wg,td,rate
  34. tabulate[class=treat,mfg;marg=y;prin=m,c] rate
  35. model rate
  36. fit[prin=m,s,a,e] treat+mfg
  37. treat treat*mfg
  38. anova rate 
  39. stop
  40.  
  41.  
  42.  
  43.  
  44.