MZ@ !L!This program cannot be run in DOS mode. $PELZA 7(" @p@X`.text `.data @.bssp0.idataX@@.rsrc`2@PU5$AEPh Ah A ÍvU$At{jjujjt t jл؍e[^USh@x poxDA05 A5 Aà 9S@ U1ÐUBBC BASIC for Windows 3.13aBBCWrunBBCProgvCouldn't allocate virtual memoryInsufficient memoryCOMCTL32ADVAPI32SHELL32.bbcUWVS@PhjPLhjP- @E0EE@EEEEăh @PEEj EEE@h @URE܃EPj@h hjm T uj0h @hP@.j@hhTQ4 uj0h @hq@j=1t A A Aƀ A@ A@E=$AǰIf A$Ah@ h@ h@ hSj jjjjjhS ƃ,jjjV| Lj\SjPVc  t ;LVI URQ6 "u-PPRj"Q j.Pƃj\R 9wh@QjjjjjhPƃ#jjjVsLDžDžDžƅDžDžL A+; AXPh@0AhRXL$A@0A@0AIf P$AjjVj\R)LLP5 AVhSRFVr j Dž`DždhlDžpDžtx|jhpPpx)t|)‰jMQjjhhhSh@jcËt P t&uSHSD ],$AtSP1uS)S+ jjjS$uChjTPE荥([^_U ABBCWinU WVSjPj5P AT AttEܻ 8 Axu'j P5 AMj 8 A5 ABNy18 A@tVj P5 A=L At90豨@j 35 Aj 35 A@tNy5 A5 A5P Ax A~Cannot create screen bitmapCannot create screen character mapCannot create palette&Cannot allocate clipboard memoryCannot open clipboardsequencervU|WVS] uSw]w'!w:wm w3w~ly US w<6*;wrw9t& gw  w2w x e=6EPƃj&V]$Aj#AhhVǃ$VURuj0h@h@Dhjm Auj0h@h@t&j7P5#AK]S5#AAE#AE#AW5#A0Ahbhhjj5#A jHjffG11҉ȍt 8D1tt 8D1tt 8DDF~W#AWX #Auj0h@h3@jjP5#AN 0Ahj $A1fDžBEDž$P_ i #AƒPh, AV rgS5#ApE#AEE#AjX5#A0)Ѝ 4g )‰#AjZ5#A)Ѝ )‰#Aj5#AjX5#A)Ѝ ))։5#Aj 5#AjZ5#Au)Ѝ ))։5#Ao$An$A$AjhjeURAj jjMQ:jjjj%$ APEP#A#AhP1AU$A( AM1jhS@1t&EPƃj&V]$AVUR~1PMQS ;;C;;5 $AC"A"A)P5#AjEP$"A"AP5"A="A UR ;u =]$AMQzƃj5#AV5#AVVEPN UR1 }p Ap Aw~j*p Ap Aw=p Ahjp Axp A~BF$"@#@ #@'#@0#@C#@G#@P#@T#@X#@`#@d#@h#@p#@t#@x#@#@jjjMQ}yut&lhdt&\XTt&LHDt&<$AtURVSPMQL> EPVSUR@* &uj- f|u-j f}uj f}P菜j0F߃^$D$@%@%@%@%@%@%@%@%@'@'@'@'@%@%@'@'@'@'@'@'@'@'@'@'@'@'@'@'@'@'@'@'@'@'@'@'@'@'@'@'@'@'@'@'@'@'@'@'@'@'@'@'@'@'@'@'@'@'@'@'@'@'@'@'@'@'@'@'@'@'@'@'@'@'@'@'@'@'@'@%@&@&@&@&@&@&@ &@$&@(&@0&@4&@8&@@&@D&@H&@t&vrnt&fb^t&VRNt&FB>t&62.t&&"t&t&v/j f}vj f} C~<wjm f}Cx<wjQ f}Pu=j011 uj f}uo$A@u o$A1vP舙uj0ՃhP1A^0A0A 90AtjWEP`SURƃ=n$AuhPWDt&h lQhR5#At)Pp)PQRVWH %n$ASMQ9hP1A,9E"A"AP5"A( AEP19E"A"A)P5#AjUR:"A"AP5"A&="AXMQG o$Ah9Eu"A"AP5"A!A URVSP o$A   A9r : t1;w%FJ9r: t)ʍ4;vۃFPh uj0h@hP@MQt&􋅤PƋ  A ;w?J9r: t)ʍZSQVJ F F;v􋕤RaMQ uj0h@hq@EPvRj1n!At"MQVjP蓖 o$A 1J&( A"A"AMQ"AURVSP5 o$A  o$A]j jjjMQSA u=o$AuσhP1Atj5 0A5#AN#At PQ #At PD 50A5#Aǃ#At P t W  At P 5$ AXjEPzjeURw= Au ji 5!ANy!A!A$At-jjhPjjh5$A$AL A= At&j = Au1۾vX At!PX AǃX A Nyϡl AP5l Arl A= AfDžfDžDž"VDž"VfDžfDžfDžhjh/@Pjh At QL AjjjjP Ajjh@h|jl AP{{JNyDž8 Avj jPX A PI8@L@ @@ LNyF00A8 A1۾j 45 Aj 45 ANyNDž@Ph"hjƅ$AMPjhRaƅjjh5$AC$Ajjh5$A$A;5 A;5 Awfv9yj@h A)RPhZ5 A-$AtURVSPMQEPVSURx[^_UE =t=tT&T AAt&=T A+00A'w T A T A5P Ag~~~~l|88||88888|8|<<><~~<ffffff{>88xʴʴ<~~<<~~< 0``0$ff$<~~<000000llllllll~|00` 8l8vv0`0```0`00`l88l000000`00 0`xx0p0000x 0`x 8 x8x x8`x 0```xxxx| p00000000`0``0`0 0`x000xxxxxxx00000|pxxxxxx x000000xx0x0xx000 0``0 | |8l`0x ||xx ||xx8`````|| x0p000x0p0000p00000xlxx|| |x `````8|x0lx0x| x0`0000000000`0000`v8l(C) Copyright R.T.Russell 2005 t Af 0C#Ar"A0!Ap4@r!A !A !Ao t"A"A'%!Ar"A$ 0"A1҉!A!A!A!A"A"A"A At "A!A|!A5!A!AIC5!AANDABSACSADVALASCASNATNBGET BPUT COLOURCOLORCALLCASECHAINCHR$CLEAR CLOSE CLG CLS COSCOUNT CIRCLEDATADEGDEFDIVDIMDRAWENDPROC ENDWHILE ENDCASE ENDIF END ENVELOPEELSEEVALERL ERROREOF EORERR EXPEXT ELLIPSEFORFALSE FILLFNGOTOGET$GETGOSUBGCOLHIMEM INPUTIFINKEY$INKEYINTINSTR( INSTALLLINELOMEM LOCALLEFT$(LENLETLOGLNMID$(MODEMODMOVEMOUSENEXTNOTONOFFOFORIGINOROPENINOPENOUTOPENUPOSCLIOTHERWISEPRINTPAGE PRIVATEPTR PI PLOTPOINT(PROCPOS QUIT RETURN REPEATREPORT READREMRUN RADRESTORERIGHT$(RND RECTANGLESTEPSGNSINSQRSPCSTR$STRING$(SOUNDSTOP SUMSWAP SYSTANTAB(THENTIME TINTTOTRACETRUE UNTILUSRVDUVALVPOS WHILEWHEN WAIT WIDTHHIMEMLOMEMPAGEPTRTIMENo roomJump out of rangeSize neededMultiple labelMistakeMissing ,Type mismatchNot in a functionAddress out of rangeMissing "Bad DIM statementDIM spaceNot in a FN or PROCNot in a procedureBad use of arrayBad subscriptSyntax errorEscapeDivision by zeroString too longNumber too bigNegative rootLogarithm rangeAccuracy lostExponent rangeBad MODENo such variableMissing )Bad hex or binaryNo such FN/PROCBad callIncorrect argumentsNot in a FOR loopCan't match FORBad FOR variableSTEP cannot be zeroMissing TOMissing OFNot in a subroutineON syntaxON rangeNo such lineOut of dataNot in a REPEAT loopWHEN/OTHERWISE not firstMissing #Not in a WHILE loopMissing ENDCASEOF not lastMissing ENDIFBad MOUSE variableNo such system callBad librarySize mismatchDATA not LOCALCF< us< t<,t.+<(F?<)tF@<$(r<(tފFsԊC<(uC"At~%t ~$t<#uZ<$t<%t<#t81!AډFFrC<%t<$t<#t"At N#CF<(uCCSu?u^O;T$ uV;Wu SQWSQRVWU11PSVWUQR QR [X 6[X 6]_^[XMu]B_^ZY[⦉_Y[Ku^][YXJuyzPU6jHuϨyzPU6j-uyzPU5JR$6<,uFpR5l[X]$}fEuQW ff_Y9rY y )s1)v9rωJ,$!AVyl^uP>5PU5_[u;81@utـtt ttEE닰AV=!A< ׹< 5!A< ^hPI6trà y PQ$YX$R@tցR@tӁiQ@tˁN@t2T@t2_b@tY:T@t;M@t#S@u[^5!A[A[!A[_[`a[[PQRZYXt[]$xL$T$XYZQRQT$ ;YZX5(u !A>(UR1FQW_Y xBRA>,3ʋ-!MEذ$X܉-!)<,: tIt V^ta5ˍ$)r!9r)܉:PS )V^,$ \ *r퍌$9r-!A}02<,6F2@utx2 0fϸʴڸfݴܳ۰/5!uFQR1ɀ>uF u Ұ#tQRV Y!AYT $N@tT@tށR@ut<,t ;l$uql$D$nD$ \$u ۜ $-l$ \$QR ZYD$\$xӽ,u ,<,u Fgt$B $N@tT@tR@ufYY!Aϰ!0P0Pz NVYt)rQ!A t!A^tËA,<(2,<(0F",u?' uffF F*<"t< u4T!A& STOP!!Amlltڱ谪lrҳmܱmomڱ簪p簪l6аǰlܰuܰu"z屹lnvɰפɰvyyЯ3氪ɰy$ɰxw۳ѳvɰvvvyyy,xɰy9ɰ9Ⱚ氪$u F!A׹0 (kn@C:h* 52_c/1:k ^# x-·&X2x?h)SaQYӺo%(Ny9?!I@|o|樭7EW߄v |Η@½pQv0rÊAr"At߀7À:r0F+t-tNÑڑuMULu8 ynXTLTLx7szIjq_\!AT!AL$t 9<Fy "AI1ɱv+ 8S" 8û.e)u 1ɺ t  !ǃ1SWUQRȉӽ [X ]_[Ou qxLpvȜ*{~[tt81۰ûI  KRV ^[yyS r@.seoO.v7my*{O}V~[S x'[ yۀÀr x  xÀt8HȇHr1ðJ  xtہ35ϱrS QRȉw Ht}~L~㩪*[XR[QRt*sLjQd$sx[X8 V[^~8{  QRȉ;m[XtS6  RT~QRU ZYr"yQRA ZYr  QQRȉE 57{k|Aa|{7}Z*L}\~[X[Ss [s {[ y8z 8 u "A' S 0 <,uF gY tH [SVև^[u)܉19s() t Ku9rQVWF t_^Yu  3<!A y <#t@6!Ajs< t< t tfBufJы!A) c tr16 <,t FB1 9s  9sV)ˉ׍4;QY^ 4~ 0TYx R^ h [ tKtQYVցQ9YtKu^!A) ðuˁsVʉ֋=!ANO^`!X1P H Й5 Й)q  q) tt߈ʼn1؜18s͇Ӊ(t̀ r r1xs)u!)s xu=urug y8ðv ttt݈ʼn1М19r9Ӊur (r x Z_xx8Éq tt߈ʼn1М x=ur1 t%ySȉӺ1ɽ5[ r u1QR1ɺPSȉӽ s\$D$ [Xu҃tPS,s%< s1 us1^{[X sB: v$ r1J111 ucwt WEtrV 9'u111$2F2&F1F=!AQ0Yxtu8t'@'ׇ!A) C<~tz@w Att A=!A1ty-Gq y-GiM/S[1I SW׉11{u1_[ tQ9YCMEuu  y11 ss C u uMEutEtك"Au ~1Ƀ|RMu 9˱C9|كf0. t0GCu u.GKI`u!A) IFu~.GI5uEGK1Ez@[WU_9s AQ0G YSQW1ۿ t1 u!!_Y[''''ÑCtuCՀsـ s1u sp1 t8tss9i!018r')ss̀8r s) xFuÑ u. t*"A  Ҝyʀـx91y1yچ8u91y1yڑ9u99ˉstÝÍ$);-!Ar])QD$ W|$_0陿]!A׉QYQR ]$L$'[X]QRE]Uf]E]Ux"D$\$])AsIIxs)ArA)sI8u@F)uF8tN8<=mb8uF$<>s<=s<F6@uE"A u utT A0u1ɉȉϻQ@=_YXuftKftDt ttiugAu>^"S?@1[t%u:>:uFd$uʯ1ɇ@uu $[t_YYXX_YYXX̀_ZY]Xct`P YXtUPMXtIPf @Xt

r"At;= As_^ð鰱N> :t ACt;uPVFC>:^ݰijXF< t A> @> B B/ BA BB BB BB BB BBBBBBA/AAA/A//A/BBBBBBBABABABBB/BBB/AABBBBB/B/BAB/BBB/BB BB!BB"BA#$A$$A$$A$BB%BA&BA'BA(BA)BA*BA+BA,BA-BA.BA/BA0B/1BA1B@1/>3A>3@>3BA2A<3B3B/3BA3BA4BA5BA6BA7BA8BA9BA:BA;BA<BA=BA>?/////7//7//7/7/////////02345/012345/////////////////////////////////////////////////////////////////////////////////////////=&='=(=)=+=,=-=.=/=0=4=5=6=7=8=9=:=;=<A==A=>A=?A=@A=AA=BA=CA=DA=EA=FA=G//BBBBB/BBBABABABABABABABABABABABBBBBBBBBBB/B/BBBBBBBB/BB/BBB/B/BBBBBBB/B/B/B/BBBBBBBBB//?B/B/BBBBBBBBBBBBBBB/B/BBB/B/BBBBBBBBBBBBB/B/B/ BB BB BB BB BB/BB//?B/B/B/B/B/B/B/B/BB/BBBB/B/B/B/BBB B!BB"BB#B$BB%BB&BB'7 ?$""%%## !bPHff<::==;;89f'/8XOiiiGlfmswvrfrt}~|P-H-rvwssu|~}q{yupz{zxtſĿff`С``LJ```     nfofaf_f`fWhӗїӟџ@@Ӈчӏяӧѧӿѿfӧѧӯѯf,**--++()Щ422553301 ѺѣѺѻѺѳѺѫm oa`٤٬٬٫ @ggggg•1܇އܗܟޗޟܷ޷ܿ޿݇ۯߧ߇ٯ٧܏ޏݧݷݷݗݟۿߗߟ߷ٿٿٷٷݿݿܧާܯޯA"2Rbr#3Scs$4Tdt%5Ueu(8Xhx )9Yiy 0@P`p  0 @ P ` p  0 @ P ` p  0 @ P ` p  0 @ P ` p 0@P`p 0@P`p((( (0(@(P(`(p))) )0)@)P)`)p*** *0*@*P*`*p+++ +0+@+P+`+p--- -0-@-P-`-p... .0.@.P.`.p/// /0/@/P/`/pHHH H0H@HPH`HpIII I0I@IPI`IpJJJ J0J@JPJ`JpKKK K0K@KPK`KpMMM M0M@MPM`MpNNN N0N@NPN`NpOOO O0O@OPO`Op 0@P`p 0@P`p 0@P`p 0@P`p 0@P`p 0@P`p 0@P`p$ -6?@IR[mv ()*+,-./0123456789:;<=>?@ABCDEFGHIJKLMNOPQRSTUVWXYZ[\]^_hijklmnopqrstuvwxyz{|}~~#A2"A#AÇ%$Aj0jR5$At%$AÉ%$A% Aj+-$A A A AE!A!AƉ"AƉ"AƉ "AƉ$"A"AƉ8"A3"A4"A ,"APSQRZZY[X u>o$AuNo$Auv%o$A%1ۇ"A ta:<o$Aua`%$Aj ]v%$Aa-uo$A uGo$Au%o$A="Au`jao$Auyo$Auo$A%o$A8%o$Aߠ#A:#AtGS$"AË#A#Ay$A#A:#A[t o$A 1ۉ"A"Au9M%$Ajt"A tP~q"A t#PjqP5#Ai$qiRq"A tPCqo$Ap%$Ajjh5$ACs%$A]%$Aj t%$A%$At%$A;$A$As "A1һ "AP%$At%$A1һ Z)‰"AOR<,tYEPPP)RYY[zRCYeְ!AK Bۼ>,t ְ-Rр>,۰ڈֈʈݱ݈İ鬙蔼0ۀ>,ڰֈx<uFnRRZY=u1(<uF,u>tM[ZYPZٝa1Ҵ ډ1ɴ ى1Ҵ 1ɴ pFR舻R_X[ZYPYZ0Otĉ:TRR~Rx_ZYڴ 0ɰ4 - &    ROY[1u<(sFºR&YZ1Èư 蝺 x-JjjjjR%$ASjzk%$AZ1ڃs#A*#A#A*#Ats,1F-*Y[1"A u2"Au r`%$Aj *h%$Aa{tݰSQ2"Au91ۇ"A t *"AY[*됐< u Ӑ 뽐PƁÉ׊G3PÇ8< t<|t?@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@?>=;9742/,)&#    #(,059>BGLQUZ_ejotz&.6>GOX`ir{  ,8DP\iv,<L\l| /BVj~*AXp9Tp:Yx < ^  7 ^ $ N x & T  @ q >sO8t/o4xI$oWHEM`G|Q345s8v%fgz{y,WET7I90(hi12DR6UOP&kml-$!AXFYJK on"SCGHNL.j Z VBM#`acpqrtuwx'debKERNEL32USER32GDI32COMCTL32WINMMADVAPI32COMDLG32SHELL32= AtR`%$Ajjh5 AD%$Aa A`%$Ajjh5 AnD%$Aaã A`%$AjPh5$AFD%$AD$aË"Aŀut.MID%$Aj5"Ah5$AC%$A=v uðdDevice unavailable`0"A5 S%$AjSDPPjkDY ‡%$A[)T$SP%$AjPQRSD%$AX[%$ASC%$Aa%$ARD%$AhQZ uZ uJ)GB Ѝj`%$Ah jUPrrjjQSWV5#An8(@%$AaRQSWV%$AjP5$A2BK8%$AZR/Z:,uBRZSn_%$AWw h jVWvvjjQSRP5#A?6C%$AwSQ"A YZt@P`D$ %$AjPQRSfC%$AaXtPQR%$ASB%$AZYX9A%$ARB%$A3 u#A u #A)G`%$AQS5#A,7>P5#Ar>$PPSPv>%$Aa%$AXh WV5#AQSjjP>6>$PjS>D$d$[Y ]%$A䘺r@r r1ҸÃ6`%$APjA%$AD$a rBCMCCS PSRڃ%$AUjSRQjU5#A(C jSjjjU5#A6=\$C{ rC=5v=%$AZ[Y)ЉC"C2Ë"A ŀut .BMP<,uBv<,uBj<,uB^·3Á= At*`%$Ajjh5 A?%$AaSQ0`YZX= At%`%$AP5"Ah5 AR?%$Aa9ð1`rSQ0YZ P`D$ %$AjPQRS?%$AtoaXSQ!erSQYZ t3P`D$ %$AjPQRS@%$AatX9cÇ%$A?%$At tnuq_File or path not found= t {t=u@_Bad directoryu&_Too many open filest!tt tRu^Access deniedu^Invalid channel=u^Bad deviceu^Disk faulttlu^Disk read onlyul^Device faultPt=uL^File existspt'u/^Disk full^Unknown error ^Bad commandSrX 8"AݍL"AK :K!u.0Ɂ;&C#u0[ÊK [Ã<"APP`%$A5$ AjjjjRPjQS5$ Ae>=u+==uo$Au+jd5$ A=-t%$AaYX[Ç%$A\$%$AjS=SrV 8"AݍL"AK :K!u&0Ɂ;P $XK K"K#[Ã<"APP`%$A5$ AjjjjRPjQS5$ Af=<u+<=uo$Au*jd5$ A <-t%$AaYX[Ç%$A\$%$AjS< uUQWL"A t;"At ;"At5C߿<"A tC_Y`rP 8"AݍL"A;P t#1C  z%$AP:%$AaÃ[<"A1PX"A:t":*tz:ukztezt_ztY%$A`P1BBCL$AWG G4_A"A"A$P%$APP5$Aj5$AË$#AD$D$#AD$ QRjSP ZYX][)X)`jUjj0a,;l$~L$); $} $,;l$ ~T$ );T$}T$jPSRQ5$ASXP=#AuMjt$ "APP5#Af8P5#AXD$8XvRXPjjjjjjt$jjjh jj5#A5#A t(PP5#AYX"A tP{XjBy$bhhjj5#A'j5#Az%$AXy/f f"A%"A%"Af"A%"AܸL%"Aun= A "A#A"A"A"A"A="A@1=,"A@W  @@ X X``88@@ !!!!"AfC=#AfC=#AC#A"AC#A"A"ACȢ"A|#AC<t<tA< sTt<u$Ct %"Au "A lC< t< t È%"AX%"APCc$)Al%=,"A"A%$AÀ="Atf{#Afy#Ax#AÃ%$AP5$Am%$AXX[Y #A #A1#A #A"A"A"A"A"A "AI"A"A"A "A"A "AAʉى&|#A|#AE-#AÊE-#AE-#AŠER-#AZ x]9}Y; #AQ xM9}I;#AA"A "A"A"A"A;"A|Y;"A}Q"A;"A|D;"A}StretchBlt?StretchDIBitsCTextOutA*CreateFontAAnimatePalette4CreatePalette6CreatePen:CreateRectRgnArc?CreateSolidBrushBDeleteDCEDeleteObjectKEllipseLEndDocMEndPage BitBlt`ExtFloodFillbExtTextOutAdEnterCriticalSection{ExitProcessBuildCommDCBAFindCloseFindFirstFileAFindNextFileAGetCommandLineAGetCurrentDirectoryAGetCurrentThreadClearCommError GetFileAttributesACloseHandleGetFileSizeGetFullPathNameAGetLastErrorGetLocalTime GetModuleFileNameA"GetModuleHandleA:GetProcAddressMGetStartupInfoAjGetTickCount~GlobalAllocGlobalFreeGlobalLock(CopyFileAGlobalUnlockInitializeCriticalSectionLeaveCriticalSectionLoadLibraryA-CreateDirectoryAMoveFileA1CreateEventAPurgeComm4CreateFileA ReadFileRemoveDirectoryAResetEvent,SetCommStateKSetCurrentDirectoryAOSetEndOfFileSSetEventVSetFileAttributesAXSetFilePointer_SetLastError`SetLocalTimeuSetThreadPriorityySetUnhandledExceptionFilterSleepVirtualAllocDCreateProcessAVirtualFreeWaitForSingleObjectWriteFileUDeleteCriticalSectionWDeleteFileAAdjustWindowRectwDefFrameProcA{DefWindowProcADestroyCaretDestroyWindowDispatchMessageAEmptyClipboardEndPaintFillRectGetAsyncKeyStateGetClientRectGetClipboardDataGetCursorPosGetDCGetFocusGetKeyStateGetMessageA&GetWindowLongA)GetWindowRect3HideCaret=InvalidateRect?InvertRectIIsClipboardFormatAvailableWKillTimer\LoadCursorA`LoadIconAvMessageBeepwMessageBoxAMoveWindowOpenClipboardPeekMessageAPostMessageAPostQuitMessageRedrawWindowRegisterClassExARegisterHotKeyReleaseDCScreenToClientScrollDCSendMessageASetCaretPosSetClipboardDataSetCursorSetCursorPosSetRectEmptySetTimer2ClientToScreenShowCaretShowWindow3ClipCursorSystemParametersInfoA4CloseClipboardTranslateMDISysAccelTranslateMessageUnregisterHotKeyUpdateWindow<CopyRect@CreateCaretPCreateWindowExA BeginPaint@@OLE32.dll@@@@@@@@WINMM.DLL(@(@(@COMDLG32.DLL<@<@<@<@<@<@<@<@<@<@<@<@<@<@<@<@<@crtdll.dllP@P@P@P@P@P@P@P@P@P@P@P@P@P@P@P@P@P@P@P@P@P@P@P@P@P@P@P@P@P@P@P@P@P@P@P@P@P@P@P@P@P@P@P@P@P@P@P@P@P@P@P@GDI32.dlld@d@d@d@d@d@d@d@d@d@d@d@d@d@d@d@d@d@d@d@d@d@d@d@d@d@d@d@d@d@d@d@d@d@d@d@d@d@d@d@d@d@d@d@d@d@d@d@d@d@d@d@d@d@d@KERNEL32.dllx@x@x@x@x@x@x@x@x@x@x@x@x@x@x@x@x@x@x@x@x@x@x@x@x@x@x@x@x@x@x@x@x@x@x@x@x@x@x@x@x@x@x@x@x@x@x@x@x@x@x@x@x@x@x@x@x@x@x@x@USER32.dll: P:8::h:BBCWRUN`c( @0x3  0pppppppp3pppppppp30 33033 30 33wwwwwwwwzwp3wx 3pxpxpxpxpxpxpxpxpxpxpxpxpxwwwwwwwwwwxwwwwwwwwwwp?  Sample size for theta/sigma corrected 30/8/2010 to equal all of the variables in the total d.f. not in the numerator d.f. I Corrected 12/08/2011 to estimate theta/sigma only for fixed effects ):@%=10: Compile with 610x530 window (*FONT Courier New,10 Y20:0,0)'" Program CritiF by C.P. Doncaster shows test and error degrees of freedom" P<" and critical F-ratio for effects of up to three treatments, given sample" NF" sizes. Reference models are all described in Doncaster & Davey (2007)." 6P =42 =17 : TJ@ : " at line ";: SZ4:0,5)" Give the model identifier (e.g., ";:1:"3.3";:4:"): ";:1:""MDL$ (d0:MDW$=MDL$:FL$="Factor_levels.txt" Jn'" Enter the numbers of levels for each factor in turn, or enter the" mx" letter '";:1:"c";:0:"' if the factor is a covariate. ";:1:"";:0:" to use another model."' 5 TABLE(5,40),L%(4),L$(4),EF$(20),ND$(20),DD$(20) j Table of upper 5% points of the F distribution: all of column 1 and columns 2-5 for small error d.f. row 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 df 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 40 50 60 70 80 90 100 500 1000 inf 161.45,18.51,10.13, 7.71, 6.61, 5.99, 5.59, 5.32, 5.12, 4.96, 4.84, 4.75, 4.67, 4.60, 4.54, 4.49, 4.45, 4.41, 4.38, 4.35, 4.32, 4.30, 4.28, 4.26, 4.24, 4.23, 4.21, 4.20, 4.18, 4.17, 4.08, 4.03, 4.00, 3.98, 3.96, 3.95, 3.94, 3.86, 3.85, 3.84 19.00, 9.55, 6.94, 5.79, 5.14, 4.74, 4.46, 4.26, 4.10, 3.98, 3.89, 3.81, 3.74, 3.68, 3.63, 3.59, 3.55, 3.52, 3.49, 3.47, 3.44, 3.42, 3.40, 3.39, 3.37, 3.35, 3.34, 3.33, 3.32, 3.23, 3.18, 3.15, 3.13, 3.11, 3.10, 3.09, 3.01, 3.00, 3.00 $ 9.28, 6.59, 5.41 $ 6.39, 5.19 $ 5.05  I%=1 5  I%<3 K%=40 K%=5  J%=I% K% : TABLE(I%,J%): @%=5:PRINTTAB(J%*5,39+I%)TABLE(I%,J%)   : 8" model, treatments, components, n, c, b, a, effects ., effect, numerator d.f., denominator d.f. 6: 1@ "1.1"," S'( A ) " ,2,1,0,0,2,1 /J "A" ,"(a-1)" ,"(n-1)*a" T Covariate A 1^ "1.1"," S'( A ) " ,2,1,0,0,2,1 -h "A" ,"1" ,"(n-2)" r: 1| "2.1"," S'( B'( A ) ) " ,3,1,0,2,3,2 / "A" ,"(a-1)" ,"(b-1)*a" 1 "B'(A)" ,"(b-1)*a" ,"(n-1)*b*a"  Covariate A 1 "2.1"," S'( B'( A ) ) " ,3,1,0,2,3,2 - "A" ,"1" ,"(b-2)" / "B'(A)" ,"(b-2)" ,"(n-1)*b" : 1 "2.2"," S'( C'( B'( A ) ) " ,4,1,2,3,4,3 / "A" ,"(a-1)" ,"(b-1)*a" 1 "B'(A)" ,"(b-1)*a" ,"(c-1)*b*a" 3 "C'(B'(A))","(c-1)*b*a" ,"(n-1)*c*b*a"  Covariate A 1 "2.2"," S'( C'( B'( A ) ) " ,4,1,2,3,4,3 - "A" ,"1" ,"(b-2)" / "B'(A)" ,"(b-2)" ,"(c-1)*b" 1 "C'(B'(A))","(c-1)*b" ,"(n-1)*c*b" &: 10 "3.1"," S'( B | A ) " ,3,1,0,2,3,3 1: "A" ,"(a-1)" ,"(n-1)*b*a" 1D "B" ,"(b-1)" ,"(n-1)*b*a" 1N "B*A" ,"(b-1)*(a-1)" ,"(n-1)*b*a" X Covariate A 1b "3.1"," S'( B | A ) " ,3,1,0,2,3,3 /l "A" ,"1" ,"(n-2)*b" /v "B" ,"(b-1)" ,"(n-2)*b" / "B*A" ,"(b-1)" ,"(n-2)*b"  Covariate B 1 "3.1"," S'( B | A ) " ,3,1,0,2,3,3 / "A" ,"(a-1)" ,"(n-2)*a" / "B" ,"1" ,"(n-2)*a" / "B*A" ,"(a-1)" ,"(n-2)*a"  Covariates A and B 1 "3.1"," S'( B | A ) " ,3,1,0,2,3,3 - "A" ,"1" ,"(n-4)" - "B" ,"1" ,"(n-4)" - "B*A" ,"1" ,"(n-4)" : 1 "3.2"," S'( C | B | A ) " ,4,1,2,3,4,7 3 "A" ,"(a-1)" ,"(n-1)*c*b*a" 3  "B" ,"(b-1)" ,"(n-1)*c*b*a" 3 "B*A" ,"(b-1)*(a-1)" ,"(n-1)*c*b*a" 3  "C" ,"(c-1)" ,"(n-1)*c*b*a" 3* "C*A" ,"(c-1)*(a-1)" ,"(n-1)*c*b*a" 34 "C*B" ,"(c-1)*(b-1)" ,"(n-1)*c*b*a" 3> "C*B*A" ,"(c-1)*(b-1)*(a-1)","(n-1)*c*b*a" H Covariate A 1R "3.2"," S'( C | B | A ) " ,4,1,2,3,4,7 1\ "A" ,"1" ,"(n-2)*c*b" 1f "B" ,"(b-1)" ,"(n-2)*c*b" 1p "C" ,"(c-1)" ,"(n-2)*c*b" 1z "C*B" ,"(c-1)*(b-1)" ,"(n-2)*c*b" 1 "C*A" ,"(c-1)" ,"(n-2)*c*b" 1 "B*A" ,"(b-1)" ,"(n-2)*c*b" 1 "C*B*A" ,"(c-1)*(b-1)" ,"(n-2)*c*b"  Covariate B 1 "3.2"," S'( C | B | A ) " ,4,1,2,3,4,7 1 "A" ,"(a-1)" ,"(n-2)*c*a" 1 "B" ,"1" ,"(n-2)*c*a" 1 "C" ,"(c-1)" ,"(n-2)*c*a" 1 "C*B" ,"(c-1)" ,"(n-2)*c*a" 1 "C*A" ,"(c-1)*(a-1)" ,"(n-2)*c*a" 1 "B*A" ,"(a-1)" ,"(n-2)*c*a" 1 "C*B*A" ,"(c-1)*(a-1)" ,"(n-2)*c*a"  Covariate C 1 "3.2"," S'( C | B | A ) " ,4,1,2,3,4,7 1 "A" ,"(a-1)" ,"(n-2)*b*a" 1 "B" ,"(b-1)" ,"(n-2)*b*a" 1$ "B*A" ,"(b-1)*(a-1)" ,"(n-2)*b*a" 1. "C" ,"1" ,"(n-2)*b*a" 18 "C*A" ,"(a-1)" ,"(n-2)*b*a" 1B "C*B" ,"(b-1)" ,"(n-2)*b*a" 1L "C*B*A" ,"(b-1)*(a-1)" ,"(n-2)*b*a" V Covariates A and B 1` "3.2"," S'( C | B | A ) " ,4,1,2,3,4,7 /j "A" ,"1" ,"(n-4)*c" /t "B" ,"1" ,"(n-4)*c" /~ "B*A" ,"1" ,"(n-4)*c" / "C" ,"(c-1)" ,"(n-4)*c" / "C*B" ,"(c-1)" ,"(n-4)*c" / "C*A" ,"(c-1)" ,"(n-4)*c" / "C*B*A" ,"(c-1)" ,"(n-4)*c"  Covariates A and C 1 "3.2"," S'( C | B | A ) " ,4,1,2,3,4,7 / "A" ,"1" ,"(n-4)*b" / "B" ,"(b-1)" ,"(n-4)*b" / "B*A" ,"(b-1)" ,"(n-4)*b" / "C" ,"1" ,"(n-4)*b" / "C*A" ,"1" ,"(n-4)*b" / "C*B" ,"(b-1)" ,"(n-4)*b" / "C*B*A" ,"(b-1)" ,"(n-4)*b"   Covariates B and C 1 "3.2"," S'( C | B | A ) " ,4,1,2,3,4,7 / "A" ,"(a-1)" ,"(n-4)*a" /( "B" ,"1" ,"(n-4)*a" /2 "B*A" ,"(a-1)" ,"(n-4)*a" /< "C" ,"1" ,"(n-4)*a" /F "C*A" ,"(a-1)" ,"(n-4)*a" /P "C*B" ,"1" ,"(n-4)*a" /Z "C*B*A" ,"(a-1)" ,"(n-4)*a" d Covariates A, B and C 1n "3.2"," S'( C | B | A ) " ,4,1,2,3,4,7 -x "A" ,"1" ,"(n-8)" - "B" ,"1" ,"(n-8)" - "B*A" ,"1" ,"(n-8)" - "C" ,"1" ,"(n-8)" - "C*A" ,"1" ,"(n-8)" - "C*B" ,"1" ,"(n-8)" - "C*B*A" ,"1" ,"(n-8)" : 1 "3.3"," S'( C | B'( A ) ) " ,4,1,2,3,4,5 / "A" ,"(a-1)" ,"(b-1)*a" 3 "B'(A)" ,"(b-1)*a" ,"(n-1)*c*b*a" 5 "C" ,"(c-1)" ,"(c-1)*(b-1)*a" 5 "C*A" ,"(c-1)*(a-1)" ,"(c-1)*(b-1)*a" 3 "C*B'(A)" ,"(c-1)*(b-1)*a" ,"(n-1)*c*b*a"  Covariate C 1 "3.3"," S'( C | B'( A ) ) " ,4,1,2,3,4,5 / "A" ,"(a-1)" ,"(b-1)*a" 1" "B'(A)" ,"(b-1)*a" ,"(n-2)*b*a" /, "C" ,"1" ,"(b-1)*a" /6 "C*A" ,"(a-1)" ,"(b-1)*a" 1@ "C*B'(A)" ,"(b-1)*a" ,"(n-2)*b*a" J Covariate A 1T "3.3"," S'( C | B'( A ) ) " ,4,1,2,3,4,5 -^ "A" ,"1" ,"(b-2)" 1h "B'(A)" ,"(b-2)" ,"(n-1)*c*b" 3r "C" ,"(c-1)" ,"(c-1)*(b-2)" 3| "C*A" ,"(c-1)" ,"(c-1)*(b-2)" 1 "C*B'(A)" ,"(c-1)*(b-2)" ,"(n-1)*c*b"  Covariates A and C 1 "3.3"," S'( C | B'( A ) ) " ,4,1,2,3,4,5 - "A" ,"1" ,"(b-2)" / "B'(A)" ,"(b-2)" ,"(n-2)*b" - "C" ,"1" ,"(b-2)" - "C*A" ,"1" ,"(b-2)" / "C*B'(A)" ,"(b-2)" ,"(n-2)*b" : 1 "3.4"," S'( C'( B * A ) ) " ,4,1,2,3,4,4 1 "A" ,"(a-1)" ,"(c-1)*b*a" 1 "B" ,"(b-1)" ,"(c-1)*b*a" 1 "B*A" ,"(b-1)*(a-1)" ,"(c-1)*b*a" 3 "C'(B*A)" ,"(c-1)*b*a" ,"(n-1)*c*b*a"  Covariate A 1 "3.4"," S'( C'( B * A ) ) " ,4,1,2,3,4,4 /& "A" ,"1" ,"(c-2)*b" /0 "B" ,"(b-1)" ,"(c-2)*b" /: "B*A" ,"(b-1)" ,"(c-2)*b" 1D "C'(B*A)" ,"(c-2)*b" ,"(n-1)*c*b" N Covariate B 1X "3.4"," S'( C'( B * A ) ) " ,4,1,2,3,4,4 /b "A" ,"(a-1)" ,"(c-2)*a" /l "B" ,"1" ,"(c-2)*a" /v "B*A" ,"(a-1)" ,"(c-2)*a" 1 "C'(B*A)" ,"(c-2)*a" ,"(n-1)*c*a"  Covariates A and B 1 "3.4"," S'( C'( B * A ) ) " ,4,1,2,3,4,4 - "A" ,"1" ,"(c-4)" - "B" ,"1" ,"(c-4)" - "B*A" ,"1" ,"(c-4)" / "C'(B*A)" ,"(c-4)" ,"(n-1)*c" : 1 "4.1"," S'| A " ,2,1,0,0,2,1 3 "A" ,"(a-1)" ,"(n-1)*(a-1)" : 1 "4.2"," S'| B | A [Model 1]" ,3,1,0,2,3,3 3 "A" ,"(a-1)" ,"(n-1)*(a-1)" 3 "B" ,"(b-1)" ,"(n-1)*(b-1)" 9  "B*A" ,"(b-1)*(a-1)" ,"(n-1)*(b-1)*(a-1)" : 2  "4.3"," S'| C | B | A [Model 1]",4,1,2,3,4,7 3* "A" ,"(a-1)" ,"(n-1)*(a-1)" 34 "B" ,"(b-1)" ,"(n-1)*(b-1)" 9> "B*A" ,"(b-1)*(a-1)" ,"(n-1)*(b-1)*(a-1)" 3H "C" ,"(c-1)" ,"(n-1)*(c-1)" 9R "C*A" ,"(c-1)*(a-1)" ,"(n-1)*(c-1)*(a-1)" 9\ "C*B" ,"(c-1)*(b-1)" ,"(n-1)*(c-1)*(b-1)" ?f "C*B*A" ,"(c-1)*(b-1)*(a-1)","(n-1)*(c-1)*(b-1)*(a-1)" p: 1z "5.1"," B | S'| A " ,3,2,0,1,3,3 3 "A" ,"(a-1)" ,"(n-1)*(a-1)" E "B" ,"(b-1)" ,"(b-1)*(n-1)+(b-1)*(n-1)*(a-1)" E "B*A" ,"(b-1)*(a-1)" ,"(b-1)*(n-1)+(b-1)*(n-1)*(a-1)" : 1 "5.2"," C | S'| B | A " ,4,2,1,3,4,7 Q "A" ,"(a-1)" ,"(n-1)*(a-1)+(n-1)*(b-1)+(n-1)*(b-1)*(a-1)" Q "B" ,"(b-1)" ,"(n-1)*(a-1)+(n-1)*(b-1)+(n-1)*(b-1)*(a-1)" Q "B*A" ,"(b-1)*(a-1)" ,"(n-1)*(a-1)+(n-1)*(b-1)+(n-1)*(b-1)*(a-1)" o "C" ,"(c-1)" ,"(c-1)*(n-1)+(c-1)*(n-1)*(a-1)+(c-1)*(n-1)*(b-1)+(c-1)*(n-1)*(b-1)*(a-1)" o "C*A" ,"(c-1)*(a-1)" ,"(c-1)*(n-1)+(c-1)*(n-1)*(a-1)+(c-1)*(n-1)*(b-1)+(c-1)*(n-1)*(b-1)*(a-1)" o "C*B" ,"(c-1)*(b-1)" ,"(c-1)*(n-1)+(c-1)*(n-1)*(a-1)+(c-1)*(n-1)*(b-1)+(c-1)*(n-1)*(b-1)*(a-1)" o "C*B*A" ,"(c-1)*(b-1)*(a-1)","(c-1)*(n-1)+(c-1)*(n-1)*(a-1)+(c-1)*(n-1)*(b-1)+(c-1)*(n-1)*(b-1)*(a-1)" : 1 "5.3"," C | B | S'| A " ,4,3,1,2,4,7 3 "A" ,"(a-1)" ,"(n-1)*(a-1)"  "B" ,"(b-1)" ,"(b-1)*(n-1)+(b-1)*(n-1)*(a-1)+(c-1)*(n-1)+(c-1)*(n-1)*(a-1)+(c-1)*(b-1)*(n-1)+(c-1)*(b-1)*(n-1)*(a-1)" $ "B*A" ,"(b-1)*(a-1)" ,"(b-1)*(n-1)+(b-1)*(n-1)*(a-1)+(c-1)*(n-1)+(c-1)*(n-1)*(a-1)+(c-1)*(b-1)*(n-1)+(c-1)*(b-1)*(n-1)*(a-1)" . "C" ,"(c-1)" ,"(b-1)*(n-1)+(b-1)*(n-1)*(a-1)+(c-1)*(n-1)+(c-1)*(n-1)*(a-1)+(c-1)*(b-1)*(n-1)+(c-1)*(b-1)*(n-1)*(a-1)" 8 "C*A" ,"(c-1)*(a-1)" ,"(b-1)*(n-1)+(b-1)*(n-1)*(a-1)+(c-1)*(n-1)+(c-1)*(n-1)*(a-1)+(c-1)*(b-1)*(n-1)+(c-1)*(b-1)*(n-1)*(a-1)" B "C*B" ,"(c-1)*(b-1)" ,"(b-1)*(n-1)+(b-1)*(n-1)*(a-1)+(c-1)*(n-1)+(c-1)*(n-1)*(a-1)+(c-1)*(b-1)*(n-1)+(c-1)*(b-1)*(n-1)*(a-1)" L "C*B*A" ,"(c-1)*(b-1)*(a-1)","(b-1)*(n-1)+(b-1)*(n-1)*(a-1)+(c-1)*(n-1)+(c-1)*(n-1)*(a-1)+(c-1)*(b-1)*(n-1)+(c-1)*(b-1)*(n-1)*(a-1)" V : 1` "5.4"," C | B | S'| A " ,4,3,1,2,4,7 3j "A" ,"(a-1)" ,"(n-1)*(a-1)" Et "B" ,"(b-1)" ,"(b-1)*(n-1)+(b-1)*(n-1)*(a-1)" E~ "B*A" ,"(b-1)*(a-1)" ,"(b-1)*(n-1)+(b-1)*(n-1)*(a-1)" o "C" ,"(c-1)" ,"(c-1)*(n-1)+(c-1)*(n-1)*(a-1)+(c-1)*(b-1)*(n-1)+(c-1)*(b-1)*(n-1)*(a-1)" o "C*A" ,"(c-1)*(a-1)" ,"(c-1)*(n-1)+(c-1)*(n-1)*(a-1)+(c-1)*(b-1)*(n-1)+(c-1)*(b-1)*(n-1)*(a-1)" o "C*B" ,"(c-1)*(b-1)" ,"(c-1)*(n-1)+(c-1)*(n-1)*(a-1)+(c-1)*(b-1)*(n-1)+(c-1)*(b-1)*(n-1)*(a-1)" o "C*B*A" ,"(c-1)*(b-1)*(a-1)","(c-1)*(n-1)+(c-1)*(n-1)*(a-1)+(c-1)*(b-1)*(n-1)+(c-1)*(b-1)*(n-1)*(a-1)"  : 1 "5.5"," C | B | S'( A )" ,4,3,1,2,4,7 / "A" ,"(a-1)" ,"(n-1)*a" 5 "B" ,"(b-1)" ,"(b-1)*(n-1)*a" 5 "B*A" ,"(b-1)*(a-1)" ,"(b-1)*(n-1)*a" I "C" ,"(c-1)" ,"(c-1)*(n-1)*a+(c-1)*(b-1)*(n-1)*a" I "C*A" ,"(c-1)*(a-1)" ,"(c-1)*(n-1)*a+(c-1)*(b-1)*(n-1)*a" I "C*B" ,"(c-1)*(b-1)" ,"(c-1)*(n-1)*a+(c-1)*(b-1)*(n-1)*a" I "C*B*A" ,"(c-1)*(b-1)*(a-1)","(c-1)*(n-1)*a+(c-1)*(b-1)*(n-1)*a"  :  MDL$="5.6" MDW$="6.3"  : 1( "5.7"," C | B | S'( A ) " ,4,3,1,2,4,7 /2 "A" ,"(a-1)" ,"(n-1)*a" W< "B" ,"(b-1)" ,"(b-1)*(n-1)*a+(c-1)*(n-1)*a+(c-1)*(b-1)*(n-1)*a" WF "B*A" ,"(b-1)*(a-1)" ,"(b-1)*(n-1)*a+(c-1)*(n-1)*a+(c-1)*(b-1)*(n-1)*a" WP "C" ,"(c-1)" ,"(b-1)*(n-1)*a+(c-1)*(n-1)*a+(c-1)*(b-1)*(n-1)*a" WZ "C*A" ,"(c-1)*(a-1)" ,"(b-1)*(n-1)*a+(c-1)*(n-1)*a+(c-1)*(b-1)*(n-1)*a" Wd "C*B" ,"(c-1)*(b-1)" ,"(b-1)*(n-1)*a+(c-1)*(n-1)*a+(c-1)*(b-1)*(n-1)*a" Wn "C*B*A" ,"(c-1)*(b-1)*(a-1)","(b-1)*(n-1)*a+(c-1)*(n-1)*a+(c-1)*(b-1)*(n-1)*a" x :  MDL$="5.8" MDW$="6.6"  :  MDL$="5.9" MDW$="6.7"  :  MDL$="6.1" MDW$="4.1"  :  MDL$="6.2" MDW$="4.2"  : 1 "6.3"," B | S'( A ) " ,3,2,0,1,3,3 / "A" ,"(a-1)" ,"(n-1)*a" 5 "B" ,"(b-1)" ,"(b-1)*(n-1)*a" 5 "B*A" ,"(b-1)*(a-1)" ,"(b-1)*(n-1)*a"  : 1 "6.4"," C'( B )|S'( A ) " ,4,3,1,2,4,7 ; "A" ,"(a-1)" ," Quasi d.f." 7 "S'(A)" ,"(n-1)*a" ,"(c-1)*b*(n-1)*a" ;" "B" ,"(b-1)" ," Quasi d.f." ;, "B*A" ,"(b-1)*(a-1)" ," Quasi d.f." 76 "B*S'(A)" ,"(b-1)*(n-1)*a" ,"(c-1)*b*(n-1)*a" 7@ "C'(B)" ,"(c-1)*b" ,"(c-1)*b*(n-1)*a" 7J "C'(B)*A" ,"(c-1)*b*(a-1)" ,"(c-1)*b*(n-1)*a" T : 1^ "6.5"," C | B | S'( A ) " ,4,3,1,2,4,7 /h "A" ,"(a-1)" ,"(n-1)*a" 5r "B" ,"(b-1)" ,"(b-1)*(n-1)*a" 5| "B*A" ,"(b-1)*(a-1)" ,"(b-1)*(n-1)*a" 5 "C" ,"(c-1)" ,"(c-1)*(n-1)*a" 5 "C*A" ,"(c-1)*(a-1)" ,"(c-1)*(n-1)*a" ; "C*B" ,"(c-1)*(b-1)" ,"(c-1)*(b-1)*(n-1)*a" ; "C*B*A" ,"(c-1)*(b-1)*(a-1)","(c-1)*(b-1)*(n-1)*a"  : 1 "6.6"," C | S'( B'( A ) ) " ,4,2,1,3,4,5 / "A" ,"(a-1)" ,"(b-1)*a" 1 "B'(A)" ,"(b-1)*a" ,"(n-1)*b*a" 5 "C" ,"(c-1)" ,"(c-1)*(b-1)*a" 5 "C*A" ,"(c-1)*(a-1)" ,"(c-1)*(b-1)*a" 7 "C*B'(A)" ,"(c-1)*(b-1)*a" ,"(c-1)*(n-1)*b*a"  : 1 "6.7"," C | S'( B | A ) " ,4,2,1,3,4,7 1 "A" ,"(a-1)" ,"(n-1)*b*a" 1 "B" ,"(b-1)" ,"(n-1)*b*a" 1 "B*A" ,"(b-1)*(a-1)" ,"(n-1)*b*a" 7& "C" ,"(c-1)" ,"(c-1)*(n-1)*b*a" 70 "C*A" ,"(c-1)*(a-1)" ,"(c-1)*(n-1)*b*a" 7: "C*B" ,"(c-1)*(b-1)" ,"(c-1)*(n-1)*b*a" 7D "C*B*A" ,"(c-1)*(b-1)*(a-1)","(c-1)*(n-1)*b*a" N : 1X "7.1"," B'| A " ,2,0,0,1,2,1 3b "A" ,"(a-1)" ,"(b-1)*(a-1)" l : 1v "7.2"," C | B'| A " ,3,0,1,2,3,3 3 "A" ,"(a-1)" ,"(b-1)*(a-1)" 3 "C" ,"(c-1)" ,"(c-1)*(b-1)" 9 "C*A" ,"(c-1)*(a-1)" ,"(c-1)*(b-1)*(a-1)"  : % Find data on the relevant model D@A readdata  Get sample sizes J 0,10)" Model ";MDL$;" factors: ";TR$;:54,10)"Data points" 7 4:" Levels for each (";:1:"#";:4:"): "  1:TL%=1:NOc%=  I%=1 CM%  26+I%*4,11)L$  23,0,10,0,0;0;0;  L$="C" L$="c"  L$="c" (I%=1 (I%=2 (MDL$="2.1" MDL$="2.2" MDL$="3.4")) (I%=3 (MDL$="2.2" MDL$="3.3")) (MDL$)>4) dxL % L$<>"" L$(I%)=L$:L%(I%)=(L$) * L%(I%)>1 T~M 54 L$="c" (L$="" L$(I%)="c") NOc%= dxL 2> L%(I%)<10 26+I%*4,11)L$(I%);" ": DRM EH L%(I%)<100 26+I%*4,11)L$(I%);" " 26+I%*4,11)L$(I%);" " R \ Print d.f. f 23,0,10,13,0;0;0; (p 0:'" d.f."; 0z NOc% 28," ");"Theta/sigma" 39," ") 9 " Effect numerator denominator F[0.05]"; . NOc% " for 80% power" 18," ") - n=L%(n%):c=L%(c%):b=L%(b%):a=L%(a%):OK%=  I%=1 EF% W =26 " Quasi-d.f. Quasi-F": : DPO : " at line ";:  (ND$(I%))<0 OK%=  OK%= "": DdO ' " ";EF$(I%),(ND$(I%)),(DD$(I%));  :  Find critical F ! critF((ND$(I%)),(DD$(I%)))  : : Find theta/sigma for effects that are not covariates FIXED%=  J%=1 (EF$(I%)) % EF$(I%),J%,1)=(39) FIXED%= $ *. Fcrit=0 FIXED%= NOc%= TrO r8 First identify the sample size from the product of all variables in the total d.f. not in the numerator d.f. {B Use the denominator d.f. for the last variable in the model, since it always contains all variables in the total d.f. LSampleSize$="1":NV$="" -V Compile the list of numerator variables ` J%=1 (ND$(I%)) jT$=ND$(I%),J%,1) 4t T$="n" T$="a" T$="b" T$="c" NV$=NV$+T$ ~ ~ For each variable in the last denominator d.f., check that it is not in the numerator list and not already accounted for  J%=1 (DD$(EF%)) T$=DD$(EF%),J%,1) 4 T$<>"n" T$<>"a" T$<>"b" T$<>"c" T@O Good%=  K%=1 (NV$)  T$=NV$,K%,1) Good%=   K%=1 (SampleSize$) % T$=SampleSize$,K%,1) Good%=  , Good% SampleSize$=SampleSize$+"*"+T$    Find theta/sigma ;find_theta((ND$(I%)),(DD$(I%)),(SampleSize$),Fcrit) : 7( Print critical F and (if appropriate) theta/sigma 2@%=&020208 C< Fcrit=0 ," -"24," "): DPO " ",Fcrit; DF NOc% FIXED% " ";10^(RATV/2);8," ") ,22," ") P@%=10 Z: d (n OK%= 0,14)"":readdata: tXM %x Check d.f. add up to total d.f. 6 (MDW$)>4 (MDW$)<5.5 TLC$="+"+(n-1): tVO + (MDW$)>7 TLC$="+"+(b-1) TLC$=""  I%=1 EF% TLC$=TLC$+"+"+((ND$(I%)))   I%=1 EF% BAD%=  J%=1 EF%  DD$(I%)=ND$(J%) BAD%=   =26 : TbP  I%=1 TXP  J%=1 I%-1  DD$(I%)=DD$(J%) BAD%=  * BAD%= TLC$=TLC$+"+"+((DD$(I%))) " ,TLM$="" 6 I%=1 CM% (@ L%(I%)>0 TLM$=TLM$+"*"+(L%(I%)) J 'TTLC$=TLC$,2,100):TLM$=TLM$,2,100) ,^'" Total d.f. = ";TLC$;59-(TLC$)," ") .h11," ");" = ";TLM$;"-1";57-(TLM$)," ") r11," ");" = ";(TLC$); ]| (TLC$)=(TLM$)-1 59-(((TLM$)))," ") " *** ERROR IN DEGREES OF FREEDOM! ***" $ Print treatment levels to file r'" Output file '";FL$;"' contains the ";(TLM$);24-(((TLM$)))," ")'" lines of factor-level combinations." '51,11)(TLM$):OUT=(FL$):#OUT,"" #OUT,"Model "+MDL$ 3#OUT,(10)+(9)+"S"+(9)+"A"+(9)+"B"+(9)+"C"  I%=1 a  a>0 I$=(I%): dVP  a%>0 I$="c" I$=""  J%=1 b  b>0 J$=(J%): dtP  b%>0 J$="c" J$=""  K%=1 c  c>0 K$=(K%): TRQ  c%>0 K$="c" K$=""  L%=1 n 2#OUT,(10)+(9)+(L%)+(9)+I$+(9)+J$+(9)+K$ & 0 : D N#OUT,(10)+"":#OUT X trL b l: vfind_theta(p%,q%,s%,Fc)  Value of Z at beta = 0.2 Z=-0.84162 RATV=3:J=10:B=Fc*p%/q%  J%=1 5 (J=J/10::RATV=RATV-J:find_IT: IT>0 (J=J/10::RATV=RATV+J:find_IT: IT<0   : find_IT %R1=1+s%*10^RATV:R2=1+2*s%*10^RATV 4IT=(((2*q%-1)*B)-(2*p%*R1-R2/R1))/(B+R2/R1)-Z  :  : readdata   !* MD$,TR$,CM%,n%,c%,b%,a%,EF% 4 I%=1 EF% > EF$(I%),ND$(I%),DD$(I%) H R MD$=MDW$ \ f: pcritF(x,y) Rz Uses algorithm in Haines, P.D (1988, 1991 Appl. Statist.), with fudge factor  x<1 y<1 Fcrit=0:  x<3 y<6 table: +P=0.95: Critical probability 1 - alpha 1c1=51.21114354*(1/(2*(1-P)))^0.25-54.06687755 2c2=-100.7003283*(1/(2*(1-P)))^0.25+86.13944869 %A=-(c2/2)+((c2/2)^2+(c1/3)^3)^0.5 %B=-(c2/2)-((c2/2)^2+(c1/3)^3)^0.5 )u=(A)*((A))^(1/3)+(B)*((B))^(1/3) 3a1=0.278393:a2=0.230389:a3=0.000972:a4=0.078108 y0=u+(a2/(3*a4)) "R=((a3/(2*a4))^2-a2/a4+y0)^0.5 QE=((3*(a3/(2*a4))^2-R^2-2*a2/a4-(4*a2*a3/a4^2-8*a1/a4-(a3/a4)^3)/(4*R)))^0.5 x0=-a3/(4*a4)+(E-R)/2 ,a0=1-2/(9*y):b0=1-2/(9*x):c0=1-a0:d=1-b0 ZG=(-2*a0*b0-((2*a0*b0)^2-4*(2*x0^2*c0-a0^2)*(2*x0^2*d-b0^2))^0.5)/(2*(2*x0^2*c0-a0^2)) )fudge=(y/(4+x))^(1/(4*(x+0.3)*y^0.8)) "$Fcrit=(0.5+100*fudge*G^3)/100 . 8: Btable L y>100 BIT%=500 BIT%=10 Vrow=y ` y>30 row=27+(y 10) j y>100 row=37+(y 500) t y>1000 row=39 a~ y<31 Fcrit=TABLE(x,row) Fcrit=TABLE(x,row)-(TABLE(x,row)-TABLE(x,row+1))*(y BIT%)/BIT%  y>1499 Fcrit=TABLE(x,40)  bb-6