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 >
Wrap
Text File
|
1995-02-16
|
2KB
|
75 lines
units [36]
set [case=ignore]
factor [levels=6] t
factor [levels=3] a
factor [levels=2] b
"factor [levels=6] row,column"
PROCEDURE 'expt3'
OPTION NAME='SEED'; MODE=p
PARAMETER NAME='T','A','B','Y'; MODE=p,p,p,p
scalar cairn
set [case=ignore]
print !t(' Latin square design');
print !t(' Randomising....');
factor [levels=6;val=6(1...6)] row
factor [levels=6;val=(1...6)6] col
"
factor [levels=6;val=6(1...6)] t
factor [levels=3] a
factor [levels=2] b
factor [levels=6] block,plot
"
open 'sess2.dat1'; chan=2
read [print=*; channel=2] unitno,t,a,b,te,pe
close chan=2
"calc t=!(6(1...6))
calc a=newlevels(t;!(1,1,2,2,3,3))
calc b=newlevels(t;!(1,2,1,2,1,2))
calc row=!((1...6)6)
calc column=!(6(1...6))
"
random [seed=seed; blocks=row*col] t,a,b,te
calc y=te+pe
text lm; val=!t(' ')
text lh; val=!t('--------------------------------------------------------\
------')
text expt; val=!t(' Experimental plan : 6 treatments')
text expb; \
val=!t(' Col.1 Col.2 Col.3 Col.4 Col.5 Col.6')
text expab; val=!t(' Experimental plan : 3 X 2 factorial')
text expr; val=!t(' Experimental results')
for i=1
tabulate [class=row,col] !(#t),!(#a),!(#b),y; max=tmax,amax,bmax,ymax
print [ip=*; ind=4; sq=y]lh
print [ip=*; ind=4; sq=y]expt
print [ip=*; ind=4; sq=y]lh
print [ip=*; ind=4; sq=y]expb
print [ip=*; ind=4; sq=y]lh
print [ip=*;rl=*;cl=*;ind=4;sq=y] tmax; dec=0; field=10
print [ip=*]lm
print [ip=*; ind=4; sq=y]lh
print [ip=*; ind=4; sq=y]expab
print [ip=*; ind=4; sq=y]lh
print [ip=*;rl=*;cl=*;ind=4;sq=y] amax,bmax; dec=0; field=8,2
question [preamble=!t(' ',' Performing experiment'); prompt=!t(\
' Press return for results...'); response=_resp; mode=p; default=cairn]
print [ip=*; ind=4; sq=y]lh
print [ip=*; ind=4; sq=y]expr
print [ip=*; ind=4; sq=y]lh
print [ip=*;rl=*;cl=*;ind=4;sq=y] ymax; dec=3; field=10
endfor
ENDPROCEDURE
open 'courseprob.lib'; chan=3; file=back
store [chan=3; proc=yes; subfile=expt3; meth=over] expt3
expt3 [seed=43251] t;a;b;v1
stop