Experimental Design and the Analysis of Variance

Size: px
Start display at page:

Download "Experimental Design and the Analysis of Variance"

Transcription

1 Expermenl Desgn nd he nlyss of Vrnce

2 Comprng > Groups - Numerc Responses Exenson of Mehods used o Compre Groups Independen Smples nd Pred D Desgns Norml nd non-norml d dsruons D Desgn Independen Smples (CRD) Pred D (RBD) Norml -Tes 1-Wy NOV -Tes -Wy NOV Nonnorml Kruskl- Wlls Tes redmn s Tes

3 Compleely Rndomzed Desgn (CRD) Conrolled Expermens - Sujecs ssgned rndom o one of he remens o e compred Oservonl Sudes - Sujecs re smpled from exsng groups Sscl model y j s mesuremen from he j h sujec from group : yj j j where s he overll men, s he effec of remen, j s rndom error, nd s he populon men for group

4 1-Wy NOV for Norml D (CRD) or ech group on he men, sndrd devon, nd smple sze: 1 ) (.. j j j j n y y s n y y On he overll men nd smple sze N y N y n y n y n n N j j

5 nlyss of Vrnce - Sums of Squres Tol Vron n ( ) 1 1 j1 j.. Tol TSS y y df N Beween Group (Smple) Vron SST n ( y j y n y y dft ) ( 1...) 1 Whn Group (Smple) Vron SSE n ( y j j y n 1 1.) ( 1) 1 s df E N TSS SST SSE df Tol df T df E

6 nlyss of Vrnce Tle nd -Tes Source of Vron Sum of Squres Degrres of reedom Men Squre Tremens SST -1 MST=SST/(-1) =MST/MSE Error SSE N- MSE=SSE/(N-) Tol TSS N-1 ssumpon: ll dsruons norml wh common vrnce H 0 : No dfferences mong Group Mens ( 1 =0) H : Group mens re no ll equl (No ll re 0) T. S.: R. R.: os os MST MSE P vl : P(, 1, N os ) ( Tle 9)

7 Expeced Men Squres Model: y j = + + j wh j ~ N(0,s ), Sn = 0: 1 ) ( ) ( rue), s ( oherwse 1 ) ( ) ( rue, s 0 : When 1) ( 1 1 ) ( ) ( 1 ) ( ) ( MSE E MST E H MSE E MST E H n n MSE E MST E n MST E MSE E s s s s s

8 Expeced Men Squres 3 cors effec mgnude of -ssc (for fxed ) True group effecs ( 1,, ) Group smple szes (n 1,,n ) Whn group vrnce (s ) os = MST/MSE When H 0 s rue ( 1 = = =0), E(MST)/E(MSE)=1 Mrgnl Effecs of ech fcor (ll oher fcors fxed) s spred n ( 1,, ) E(MST)/E(MSE) s (n 1,,n ) E(MST)/E(MSE) (when H 0 flse) s s E(MST)/E(MSE) (when H 0 flse)

9 ) =100, 1 =-0, =0, 3 =0, s = 0 E( MST ) E( MSE) B) =100, 1 =-0, =0, 3 =0, s = n B C D C) =100, 1 =-5, =0, 3 =5, s = 0 D) =100, 1 =-5, =0, 3 =5, s = 5

10 Exmple - Sesonl De Perns n Rvens Tremens - = 4 sesons of yer (3 replces ech) Wner: Novemer, Decemer, Jnury Sprng: erury, Mrch, prl Summer: My, June, July ll: ugus, Sepemer, Ocoer Response (Y) - Vegeon (percen of ol pelle wegh) Trnsformon (or pproxme normly): Y ' rcsn Y 100 Source: K.. Engel nd L.S. Young (1989). Spl nd Temporl Perns n he De of Common Rvens n Souhwesern Idho, The Condor, 91:37-378

11 Sesonl De Perns n Rvens - D/Mens y y y y y Y Wner(=1) ll(=) Summer(=3) ll (=4) j= j= j= Y' Wner(=1) ll(=) Summer(=3) ll (=4) j= j= j=

12 Sesonl De Perns n Rvens - D/Mens Plo of Trnsformed D y Seson Trnsformed % Vegeon Seson

13 Sesonl De Perns n Rvens - NOV Tol Vron TSS Beween Group SST :(df Tol ( ) Vron 3 ( ) Whn Group Vron : ) :(df (df SSE ( ) E... ( ) T 4-1 3)... ( ) 1-4 8)... ( ) NOV Source of Vron SS df MS P-vlue cr Beween Groups Whn Groups Tol Do no conclude h sesons dffer wh respec o vegeon nke

14 Sesonl De Perns n Rvens - Spredshee Monh Seson Y' Seson MenOverll Men TSS SST SSE NOV DEC JN EB MR PR MY JUN JUL UG SEP OCT Sum Tol SS Beween Seson SS Whn Seson SS (Y -Overll Men) (Group Men-Overll Men) (Y -Group Men)

15 CRD wh Non-Norml D Kruskl-Wlls Tes Exenson of Wlcoxon Rnk-Sum Tes o > Groups Procedure: Rnk he oservons cross groups from smlles (1) o lrges ( N = n n ), djusng for es Compue he rnk sums for ech group: T 1,...,T. Noe h T T = N(N+1)/

16 Kruskl-Wlls Tes H 0 : The populon dsruons re dencl (M 1 =...=M ) H : No ll dsruons re dencl (No ll M re equl) 1 T N( N 1) 1 n T. S.: H 3( N 1) R. R.: H, 1 P vl : P( H ) n djusmen o H s suggesed when here re mny es n he d. ormul s gven on pge 344 of O&L.

17 Exmple - Sesonl De Perns n Rvens Monh Seson Y' Rnk NOV DEC JN EB MR PR MY JUN JUL UG SEP OCT T 1 = = 6 T = = 4.5 T 3 = = 8 T 4 = = 19.5 H 0 : Nosesonl dfference H :Sesonl Dfference s T. S.: H 1 (6) 1(1 1) 3 R. R.( 0.05) : H.05,41 (4.5) (8) 3 (19.5) 3 3(1 1) P vlue : P( H 5.1).163

18 s Trnsformons for Consn Vrnce k y y or y y Used for Posson Dsruon k 1 T T s k y ln y or y ln y 1 T ^ 1 s k 1 y sn y Used for Bnoml Dsruon k, y T T 1 n Box-Cox Trnsformon: Power Trnsformon used o on normly nd ofen consn vrnce y T y 0 ln y 0

19 Welch s Tes Unequl Vrnces 1 1 * pprox * 1, ~ 1 W W W W W W W W n w w w s wy w y w w C m C n w C m

20 C m * W Exmple Sesonl De Perns n Rvens W W Seson n_ s_ yr_ s_^ w_ w_y_ w_y_^c_w Sum wy w y 1 1 w w n 1 w 4 1 C W * C W W mw 0.760(1.3390) qf (.95,3, 4.3) n R 0.05,3,4.3

21 Lner Conrss Lner funcons of he remen mens (populon nd smple) such h he coeffcens sum o 0. Used o compre groups or prs of remen mens, sed on reserch queson(s) of neres Populon Conrs: l... where 0 ^ 1 1 Esmed Conrs: l y... y y s s... s V l V l MSE SE l MSE 1 ^ ^ ^ ^ ^ 1 n1 n 1 n 1 n 1 n ^ ^ MSE ^ ^ MSE n... n n V l SE l n 1 n 1

22 Orhogonl Conrss & Sums of Squres Two Conrss: l l l1, l re orhogonl f: 0 for lnced d, 0 n 1 1 Conrs Sum of Squres: SSC ^ y l 1 n n 1 1 for lnced d, SSC ^ nl 1 mong remens, we cn on 1 prwse orhogonl Conrss: l,..., l 1 1 Then, we cn decompose he Beween Tremen Sum of Squres no he Conrss: SST SSC... SSC where ech of he Conrs Sums of Squres hs 1 degree of freedom 1 1 or ny conrs: Tesng H : l 0 H : l 0 0k k k k SSCk l -Tes: TS : k RR : k,1, N -Tes: k ^ MSE SE ^ k ^ l k RR : k /; N

23 Smulneous Tess of Mulple Conrss Usng m conrss for comprsons mong remens Ech conrs o e esed sgnfcnce level, whch we lel s I for ndvdul comprson Type I error re Proly of mkng les one flse rejecon of one of he m null hypoheses s he expermenwse Type I error re, whch we lel s E Tess re no ndependen unless he error (Whn Group) degrees re nfne, however Bonferron nequly mples h E m I Choose I = E / m

24 Scheffe s Mehod for ll Conrss Cn e used for ny numer of conrss, even hose suggesed y d. Conservve (Wde CI s, Low Power) l s ^ ^ ^ l y SE l MSE 1 1 n H : l 0 H : l ^ ^ ^ Rejec H f l SE l ( 1) df 1, df df N 0, df, df 1 E 1 Error Smulneous 1 ^ ^ ^ 100% Confdence Inervls: lse l ( 1) E, df, df 1

25 Pos-hoc Comprsons of Tremens If dfferences n group mens re deermned from he - es, reserchers wn o compre prs of groups. Three populr mehods nclude: sher s LSD - Upon rejecng he null hypohess of no dfferences n group mens, LSD mehod s equvlen o dong prwse comprsons mong ll prs of groups s n Chper 6. Tukey s Mehod - Specfclly compres ll (-1)/ prs of groups. Ulzes specl le (Tle 11, p. 701). Bonferron s Mehod - djuss ndvdul comprson error res so h ll conclusons wll e correc desred confdence/sgnfcnce level. ny numer of comprsons cn e mde. Very generl pproch cn e ppled o ny nferenl prolem

26 sher s Les Sgnfcn Dfference Procedure Proeced Verson s o only pply mehod fer sgnfcn resul n overll -es or ech pr of groups, compue he les sgnfcn dfference (LSD) h he smple mens need o dffer y o conclude he populon mens re no equl 1 1 LSDj / MSE wh df N n n j Conclude f y y LSD j. j. j

27 Tukey s W Procedure More conservve hn sher s LSD (mnmum sgnfcn dfference nd confdence nervl wdh re hgher). Derved so h he proly h les one flse dfference s deeced s (expermenwse error re) MSE Wj q (, ) q gven n Tle 11, p. 701 wh N - n Conclude f y y W j. j. j Tukey's Confdence Inervl:. j. y y W j When he smple szes re unequl, use q (, ) 1 1 Wj MSE n n j

28 Bonferron s Mehod (Mos Generl) Wsh o mke C comprsons of prs of groups wh smulneous confdence nervls or -sded ess When ll pr of remens re o e compred, C = (-1)/ Wn he overll confdence level for ll nervls o e correc o e 95% or he overll ype I error re for ll ess o e 0.05 or confdence nervls, consruc (1-(0.05/C))100% CIs for he dfference n ech pr of group mens (wder hn 95% CIs) Conduc ech es =0.05/C sgnfcnce level (rejecon regon cu-offs more exreme hn when =0.05) Crcl -vlues re gven n le on clss wese, we wll use noon: /,C, where C=#Comprsons, = df

29 Bonferron s Mehod (Mos Generl) B j /, C, v ( gven on clss wese wh Conclude MSE j f 1 n y. 1 n j y j. v N-) B j Bonferron 's Confdence Inervl : y y Bj. j.

30 Exmple - Sesonl De Perns n Rvens Noe: No dfferences were found, hese clculons re only for demonsron purposes MSE n 3.05,8.306 q.05, 4, df , C6, df 8 LSD W B j j j ( ) ( ) 1 3 ( ) E E Comprson( vs j) Group Men Group j MenDfference 1 vs vs vs vs vs vs

31 Rndomzed Block Desgn (RBD) > Tremens (groups) o e compred Blocks of homogeneous uns re smpled. Blocks cn e ndvdul sujecs. Blocks re mde up of suuns Suuns whn lock receve one remen. When sujecs re locks, receve remens n rndom order. Oucome when Tremen s ssgned o Block j s leled Y j Effec of Tr s leled Effec of Block j s leled j Rndom error erm s leled j Effcency gn from removng lock-o-lock vrly from expermenl error

32 Y 1 Rndomzed Complee Block Desgns j j j j j 0 E( ) 0 V( ) s j j Noe: 1 1 Y 1 Y11... Y Y Y Y Tes for dfferences mong remen effecs: H 0 : ( 1... ) H : No ll = 0 (No ll re equl)

33 RBD - NOV -Tes (Norml D) D Srucure: ( Tremens, Blocks) Men for Tremen : Men for Sujec (Block) j: Overll Men: y.. y. Overll smple sze: N = y. j NOV:Tremen, Block, nd Error Sums of Squres 1 j1 1 j1 1 j1 TSS y y df j.. SST y y df 1 SSB y y df 1 s EMST j.... j.. 1 Tol T B SSE y y y y TSS SST SSB df ( 1)( 1) E MSE j s 1 E

34 RBD - NOV -Tes (Norml D) NOV Tle: Source SS df MS Tremens SST -1 MST = SST/(-1) = MST/MSE Blocks SSB -1 MSB = SSB/(-1) Error SSE (-1)(-1) MSE = SSE/[(-1)(-1)] Tol TSS -1 H 0 : ( 1... ) H : No ll = 0 (No ll re equl) T. S.: R. R.: P vl os os : P( MST MSE, 1,( 1)( 1) os )

35 Prwse Comprson of Tremen Mens Tukey s Mehod- q n Sudenzed Rnge Tle wh = (-1)(-1) W j q (, v) MSE Conclude j f y. y j. W j Tukey' s Confdence Inervl : y y Wj. j. Bonferron s Mehod - -vlues from le on clss wese wh = (-1)(-1) nd C=(-1)/ B j /, C, v Conclude MSE j f y. y j. B j Bonferron 's Confdence Inervl : y y Bj. j.

36 Expeced Men Squres / Relve Effcency Expeced Men Squres: s wh CRD, he Expeced Men Squres for Tremen nd Error re funcons of he smple szes (, he numer of locks), he rue remen effecs ( 1,, ) nd he vrnce of he rndom error erms (s ) By ssgnng ll remens o uns whn locks, error vrnce s (much) smller for RBD hn CRD (whch comnes lock vron&rndom error no error erm) Relve Effcency of RBD o CRD (how mny mes s mny replces would e needed for CRD o hve s precse of esmes of remen mens s RBD does): RE ( RCB, CR) MSE MSE CR RCB ( 1) MSB ( 1) MSE ( 1) MSE

37 Exmple - Cffene nd Endurnce Tremens: =4 Doses of Cffene: 0, 5, 9, 13 mg Blocks: =9 Well-condoned cyclss Response: y j =Mnues o exhuson for cycls dose D: Dose \ Sujec

38 Plo of Y versus Sujec y Dose Tme o Exhuson mg 5 mg 9mg 13 mg Cycls

39 Exmple - Cffene nd Endurnce Sujec\Dose 0mg 5mg 9mg 13mg Suj MenSuj Dev Sqr Dev Dose Men Dose Dev Squred Dev TSS TSS ( ) SST TSS 9 ( ) SSB 4 ( ) SST SSB ( ) ( ) ( ) SSE ( ) df Tol 9(103.68) ( ) (9) 1 35 ( ) df E df T df (4 1)(9 1) 4 B

40 Exmple - Cffene nd Endurnce Source df SS MS Dose Cycls Error Tol H H : No Cffene : Dfference R. R.( P 0 T. S.: os vlue Conclude MST MSE 0.05) : : P( s h rue Dose Exs os mens Effec mong.05,3,4 5.9).0036 ( 1 Doses re no ll (rom EXCEL) equl 4 0)

41 Exmple - Cffene nd Endurnce Tukey' sw : q.05,4, W Bonferron 's B :.05/,6,4.875 B Doses Hgh Men Low Men Dfference Conclude 5mg vs 0mg mg vs 0mg mg vs 0mg mg vs 5mg NSD 13mg vs 5mg NSD 13mg vs 9mg NSD

42 Exmple - Cffene nd Endurnce Relve Effcency of Rndomzed Block o 4 9 MSB MSE 5.57 ( 1) MSB ( 1) MSE RE ( RCB, CR) ( 1) MSE Compleely Rndomzed 8(694.75) 9(3)(5.57) (9(4) 1)(5.57) Desgn : Would hve needed 3.79 mes s mny cyclss per dose o hve he sme precson on he esmes of men endurnce me. 9(3.79) 35 cyclss per dose 4(35) = 140 ol cyclss

43 RBD -- Non-Norml D redmn s Tes When d re non-norml, es s sed on rnks Procedure o on es ssc: Rnk he remens whn ech lock (1=smlles, =lrges) djusng for es Compue rnk sums for remens (T ) cross locks H 0 : The populons re dencl (M 1 =...=M ) H : Dfferences exs mong he group medns 1 T S T ( 1)..: r 3 ( 1) 1 R. R.: r, 1 P vl P : ( r )

44 Exmple - Cffene nd Endurnce H H Sujec\Dose 0mg 5mg 9mg 13mg Rnks 0mg 5mg 9mg 13mg Tol T S : No Dose Dfferences : Dose Dfferences Exs (4)(4 1) 180..: r (10) (8) 3(9)(4 1) R. R.( 0.05) : P r.05,41 -vlue: P( 14.).006 (rom EXCEL) Conclude Medns re no ll equl

45 Ln Squre Desgn Desgn used o compre remens when here re wo sources of exrneous vron (ypes of locks), ech oserved levels Bes sued for nlyses when 10 Clssc Exmple: Cr Tre Comprson Tremens: 4 Brnds of res (,B,C,D) Exrneous Source 1: Cr (1,,3,4) Exrneous Source : Poson (Drver ron, Pssenger ron, Drver Rer, Pssenger Rer) Cr\Poson D P DR PR 1 B C D B C D 3 C D B 4 D B C

46 Ln Squre Desgn - Model Model ( remens, rows, columns, N= ) : y jk k j jk Overll Men Effec of Tremen k y y k k.. k... j jk Effec due o row y y ^ ^ ^... ^..... Effec due o Column j y y Rndom Error Term y j. j....

47 Ln Squre Desgn - NOV & -Tes Tol Sum of Squres : TSS Tremen Sum of Squres SST Row Sum of Squres SSR 1 j1 1 Column Sum of Squres SSC y.. j1 y k 1 jk y. j. y y.. k... y... y... y... df R df 1 df df C T Error Sum of Squres SSE TSS SST SSR SSC df E ( 1) 3( 1) ( 1)( ) H 0 : 1 = = = 0 H : No ll k = 0 TS: os = MST/MSE = (SST/(-1))/(SSE/((-1)(-))) RR: os, -1, (-1)(-)

48 Prwse Comprson of Tremen Mens Tukey s Mehod- q n Sudenzed Rnge Tle wh = (-1)(-) W j q (, v) MSE Conclude j f y. y j. W j Tukey' s Confdence Inervl : y y Wj. j. Bonferron s Mehod - -vlues from le on clss wese wh = (-1)(-) nd C=(-1)/ B j /, C, v Conclude MSE j f y. y j. B j Bonferron 's Confdence Inervl : y y Bj. j.

49 Expeced Men Squres / Relve Effcency Expeced Men Squres: s wh CRD, he Expeced Men Squres for Tremen nd Error re funcons of he smple szes (, he numer of locks), he rue remen effecs ( 1,, ) nd he vrnce of he rndom error erms (s ) By ssgnng ll remens o uns whn locks, error vrnce s (much) smller for LS hn CRD (whch comnes lock vron&rndom error no error erm) Relve Effcency of LS o CRD (how mny mes s mny replces would e needed for CRD o hve s precse of esmes of remen mens s LS does): RE ( LS, CR) MSE MSE CR LS MSR MSC ( 1) MSE ( 1) MSE

50 Power pproch o Smple Sze Choce R Code When he mens re no ll equl, he -ssc s non-cenrl : * 1 1 ~ 1,, where where N n 1 1 When ll smple szes re equl: where s The power of he es, when conduced he sgnfcnce level of : n * * Pr 1 ; 1, N ~ 1, N, In R: 1 ; 1, N qf (1, 1, N ) n Power = 1 1 pf qf (1, 1, N ), 1, N, s N

51 -Wy NOV nomnl or ordnl fcors re eleved o e reled o qunve response ddve Effecs: The effecs of he levels of ech fcor do no depend on he levels of he oher fcor. Inercon: The effecs of levels of ech fcor depend on he levels of he oher fcor Noon: j s he men response when fcor s level nd cor B j

52 j -Wy NOV - Model y 1,..., j 1,..., k 1,..., n jk j j jk h y Mesuremen on k un recevng cors level, B level j jk Overll Men jk Effec of Effec of j h h level of fcor level of fcor B h h j Inercon effec when level of nd j level of B re comned Rndom Error Terms Model depends on wheher ll levels of neres for fcor re ncluded n expermen: xed Effecs: ll levels of fcors nd B ncluded Rndom Effecs: Suse of levels ncluded for fcors nd B Mxed Effecs: One fcor hs ll levels, oher fcor suse

53 xed Effecs Model cor : Effecs re fxed consns nd sum o 0 cor B: Effecs re fxed consns nd sum o 0 Inercon: Effecs re fxed consns nd sum o 0 over ll levels of fcor B, for ech level of fcor, nd vce vers Error Terms: Rndom Vrles h re ssumed o e ndependen nd normlly dsrued wh men 0, vrnce s 1 0, j 0 j 0 j j 0 jk ~ N 0, s j1 1 j1

54 Exmple - Thldomde for IDS Response: 8-dy wegh gn n IDS pens cor : Drug: Thldomde/Plceo cor B: TB Sus of Pen: TB + /TB - Sujecs: 3 pens (16 TB + nd 16 TB - ). Rndom ssgnmen of 8 from ech group o ech drug). D: Thldomde/TB + : 9,6,4.5,,.5,3,1,1.5 Thldomde/TB - :.5,3.5,4,1,0.5,4,1.5, Plceo/TB + : 0,1,-1,-,-3,-3,0.5,-.5 Plceo/TB - : -0.5,0,.5,0.5,-1.5,0,1,3.5

55 NOV pproch Tol Vron (TSS) s proned no 4 componens: cor : Vron n mens mong levels of cor B: Vron n mens mong levels of B Inercon: Vron n mens mong comnons of levels of nd B h re no due o or B lone Error: Vron mong sujecs whn he sme comnons of levels of nd B (Whn SS)

56 nlyss of Vrnce Tol Vron: TSS y y df n 1 1 j1 k 1 cor Sum of Squres: SS n y y df 1 1 j cor B Sum of Squres: SSB n y y df 1 Inercon Sum of Squres: n jk. j.... Tol B SSB n y y y y df ( 1)( 1) n 1 j1 k 1 1 j1 j.... j.... Error Sum of Squres: SSE y y df ( n 1) TSS = SS + SSB + SSB + SSE df Tol = df + df B + df B + df E jk j. E B

57 NOV pproch - xed Effecs Source df SS MS cor -1 SS MS=SS/(-1) =MS/MSE cor B -1 SSB MSB=SSB/(-1) B=MSB/MSE Inercon (-1)(-1) SSB MSB=SSB/[(-1)(-1)] B=MSB/MSE Error (n-1) SSE MSE=SSE/[(n-1)] Tol n-1 TSS Procedure: rs es for nercon effecs If nercon es no sgnfcn, es for cor nd B effecs Tes for Inercon: Tes for cor Tes for cor B H :... 0 H :... 0 H : H : No ll 0 H : No ll 0 H : No ll 0 j MSB MS MSB TS : B TS : TS : B MSE MSE MSE RR : RR : RR : B,( 1)( 1), ( n1),( 1), ( n1) B ( 1), ( n1) j

58 Exmple - Thldomde for IDS 7.5 Indvdul Pens Negve Posve Group Mens wgn menwg Plceo T hldomde Plceo T hldomde drug drug Repor WTGIN GROUP TB+/Thldomde TB-/Thldomde TB+/Plceo TB-/Plceo Tol Men N Sd. Devon

59 Exmple - Thldomde for IDS Dependen Vrle: WTGIN Source Correced Model Inercep DRUG TB DRUG * TB Error Tol Correced Tol Tess of Beween-Sujecs Effecs Type III Sum of Squres df Men Squre Sg R Squred =.5 (djused R Squred =.471) There s sgnfcn Drug*TB nercon ( DT =5.897, P=.0) The Drug effec depends on TB sus (nd vce vers)

60 NOV ddve Model If he Inercon s no sgnfcn, he Inercon erm cn e removed, nd n ddve model cn e f. n j Y SSE Y Y Y Y df n 1 jk j jk jk E 1 j1 k 1 Source df SS MS cor -1 SS MS=SS/(-1) =MS/MSE cor B -1 SSB MSB=SSB/(-1) B =MSB/MSE Error n--+1 SSE MSE =SSE /[n--+1] Tol n-1 TSS

61 Comprng Mn Effecs (No Inercon) Tukey s Mehod- q n Sudenzed Rnge Tle wh = n--+1 MSE B W q (, v) W q (, v) j j n Conclude: f y y W f y y W j.. j.. j... j... j..... j. Tukey's CI: ( ) : y y W ( ) : y y W j j j j MSE n Bonferron s Mehod - -vlues n Bonferron le wh =n--+1 MSE B MSE B j /, ( 1)/, v B j /, ( 1)/, v n n Conclude: f y y B f y y B j.. j.. j... j... j..... j. Bonferron's CI:( α - α ): y y B ( - ): y y B j j j j B j j B j B B j

62 Comprng Mn Effecs (Inercon) Tukey s Mehod- q n Sudenzed Rnge Tle wh = (n-1) W q (, v) j MSE n h Whn k level of cor B, Conclude: f y y W j k. jk. j j k. jk. Tukey's CI: ( ) : y y W Smlr for cor B n Bonferron s Mehod - -vlues n Bonferron le wh = (n-1) B j /, ( 1)/, v MSE n h Whn k level of B, Conclude: Bonferron's CI:( α - α ): y y B j k. jk. j f y y B j k. jk. j j

63 Mscellneous Topcs -cor NOV cn e conduced n Rndomzed Block Desgn, where ech lock s mde up of expermenl uns. nlyss s drec exenson of RBD wh 1-fcor NOV corl Expermens cn e conduced wh ny numer of fcors. Hgher order nercons cn e formed (for nsnce, he B nercon effecs my dffer for vrous levels of fcor C). When expermens re no lnced, clculons re mmensely messer nd you mus use sscl sofwre pckges for clculons

64 Unequl Smple Szes When smple szes re unequl, clculons nd prmeer nerpreons (especlly mrgnl ones) ecome messer Oservonl sudes ofen hve unequl smple szes due o vlly of smplng uns for cern comnons of fcor levels (vllgers of cern ypes n rurl sudy for nsnce) Expermenl sudes, even when plnned wh equl smple szes cn end up unlnced hrough echncl prolems or drop ous Some condons my e cheper o mesure hn ohers, nd wll hve lrger smple szes Some suons hve prculr conrss of hgher mpornce

65 Regresson pproch - I Smple Szes: # of Cses when cor s level, j: n nj j j j j T j j jk j j1 1 1 j1 k 1 nj n n n n n n Y Y Y j Y N Model: Y ~ 0, s (ndependen) jk j j jk jk Resrcons on Effecs: j j j 1 j1 1 j , 1 j 1 j j 1 j

66 Regresson Model: Regresson pproch - II Y X... X X... X jk 1 jk1 1 jk, 1 1 jk 1 jk,... X X X X 11 jk1 jk 1, 1 jk, 1 jk, jk 1 f cse from level 1 of fcor where: X jk1 1 f cse from level of fcor 0 oherwse 1 f cse from level -1 of fcor X jk, 1 1 f cse from level of fcor 0 oherwse 1 f cse from level 1 of fcor B where: X jk 1 f cse from level of fcor B 0 oherwse 1 f cse from level -1 of fcor B X jk, 1 f cse from level of fcor B 0 oherwse

67 Regresson pproch Exmple I Wrer Type (B) Syle (cor ) Poes (B1) Yer Pek Novelss (B) Yer Pek Concepulss (1) Elo 3 zgerld 9 (nders) Cummngs 6 Hemngwy 30 Plh 30 Melvlle 3 Pound 30 Lwrence 35 Wlur 34 Joyce 40 n & Men n11= n1= Expermenlss () Bshop 9 Jmes 38 (Seekers) Moore 3 ulkner 39 Wllms 40 Dckens 41 Lowell 41 Woolf 45 Sevens 4 Conrd 47 ros 48 Twn 50 Hrdy 51 n & Men n1= n= Y X X X X jk 1 jk1 1 jk 11 jk1 jk jk X_jk1 X_jk X1X

68 Tesng Sreges Models Model 1: ll cor, cor B, nd Inercon B Effecs Model :ll cor, cor B Effecs (Remove Inercon) Model 3: ll cor B,Inercon B Effecs (Remove ) Model 4:ll cor,inercon B Effecs (Remove B) To es for Inercon Effecs, Model 1 s ull Model, Model s Reduced df Numeror =(-1)(-1) df den =n T - Tesng for cor Effecs, ull=model 1, Reduced=Model 3 df Numeror =(-1) df den =n T - Tesng for cor B Effecs, ull=model 1, Reduced=Model 4 df Numeror =(-1) df den =n T -

69 ^ Regresson pproch Exmple - Connued Y X X X X jk 1 jk1 1 jk 11 jk1 jk jk Model 1: jk 1 jk1 1 jk 11 jk1 jk Y X.59X 0.9X X 1 1 Model : E Y X X Y X.63X Model 3: E Y X X X X E Y jk 1 jk1 1 jk 1 jk ^ X X X Y X 0.47X X ^ 1 jk 11 jk1 jk 1 Model 4: E Y X X X Y X 0.6X X jk 1 jk1 11 jk1 jk 1 1 NOV Model1 Model Model3 Model4 df SS df SS df SS df SS Regresson Resdul Tol ^

70 H Regresson pproch Exmple - Connued : 0 H : Inercon Exss SSE R df R 0 SSE df 19 E SSE R SSE df * E R dfe : B SSE dfe * RR B TS :.95,1, E NOV Model1 Model Model3 Model4 df SS df SS df SS df SS Regresson Resdul Tol

71 H Regresson pproch Exmple - Connued E * * * B : 0 H : cor Effecs Exs: SSE R df R 1.80 RR :.95,1, H : 0 H : cor B Effecs Exs: 0 1 E SSE R df R * 5.16 RR :.95,1, NOV Model1 Model Model3 Model4 df SS df SS df SS df SS Regresson Resdul Tol B

72 Esmng Tremen nd cor Level Mens/Conrss Tremen Mens: ^ jk k 1 Prmeer: j Esmor: j = Yj Esmed Sndrd Error: nj cor Mens: j ^ j1 j1 n j Y ^ s Prmeer: = Esmor: Esmed Sndrd Error: cor B Mens: Y Prmeer: = Esmor: Esmed Sndrd Error: j Y j j ^ MSE n j j ^ ^ 1 1 MSE 1 j s j 1 nj Conrs or Lner uncon of cor Mens: s j MSE MSE Prmeer: L c Esmor: L c Esmed Sndrd Error: s L c ^ ^ ^ j1 nj Conrs or Lner uncon of cor B Mens: ^ ^ ^ MSE Prmeer: L c Esmor: LB c j j Esmed Sndrd Error: s B L c B j j j1 Conrs or Lner uncon of Tremen Mens: j j1 j1 1 nj ^ ^ B jj B j j d Sndrd Error: sl B 1 j1 1 j1 Prmeer: L c Esmor: L c Y Esme j1 1 n j MSE 1 1 c n j 1 j1 j

73 Sndrd Error Mulplers Sngle Comprsons: / ; N Generl Mulple Comprsons of Tremen (Cell) Mens : Scheffe: S 1 ; 1, N Bonferron: B g, N g # of comprsons Tukey (ll prs of remen mens): T 1 q ;, N Generl Mulple Comprsons of cor Level Mens : B B g N g Scheffe: cor : S 1 ; 1, N cor B: S 1 1 ; 1, N Bonferron: cor or cor B :, # of comprsons 1 1 Tukey: cor : T q ;, N cor B: TB q ;, N

74 Creve Lfe Cycles Comprng Tremen Mens Comprng ll 4 Tremen Mens(hough no nercon ws presen): MSE 9.3 MSE 9.3 Y 8.60 n 5 s Y.4 Y 33.0 n 5 s Y n11 n1 MSE sy Y 1 n s Y 1 Y n n T q 0.95, 4, sy j Y ' j ' MSE n j n Y 11 Y 1 s Y 11 Y 1 HSD Y 11 Y s (3.43) Y Y HSD.81(3.8) Y 11 Y s Y 11 Y HSD Y 1 Y 1 s Y 1 Y 1 HSD ' j ' MSE n (3.17) ( 3.8) Y 1 Y s Y 1 Y HSD Y 1 Y s Y 1 Y HSD (3.17) (3.01) Concepulss/Poes Concepulss/Novelss Expermenlss/Poes Expermenlss/Novelss

75 Creve Lfe Cycles Comprng cor Level Mens cor (Syle): 1 j ^ j1 Y11 Y j ^ j1 Y 1 Y ^ ^ 1 Y Y s ^ ^ ( 1) , % CI for (Concepulss - Expermenlss): (.093)(3.3) , cor B (Wrer Type): 1 ^ 1 Y11 Y ^ ^ j1 Y Y Y ^ ^ 1 ^ s Y ( 1) , % CI for (Poes - Novelss): (.093)(3.3) ,1.57

76 1-Wy Rndom Effecs Model Tremen Levels n Expermen re Smple from Populon of Levels Effecs re Rndom Vrles (No xed Consns) Y 1,..., ; j 1,..., n j j NID ~ NID 0, s ~ 0, s, ndependen j j NOV oned s n xed Effecs Model: T 1 SST SST n Y Y df E MST E s ns 1 1 n SSE SSE Yj Y dfe n 1 EMSE E s 1 j1 n 1 ^ ^ MST Esmng Populon Men: Y SE 1 100% CI: n ^ ^ MST MSE Pon Esmes of Vrnce Componens: s MSE s n MST Tesng for Tr Effecs: H : s 0 H : s 0 TS : RR : MSE ^ /; 1 0 ; l, n1 MST n

77 Rndomzed Complee Block Desgns Sujecs s Blocks Y 1,..., ; j 1,..., j j j j j 1 0 E( ) 0 V ( ) s Rndom Blocks: ~ N 0, s j j j Noe: 1 1 Y 1 Y11... Y Y 1 1 Y Y Tes for dfferences mong remen effecs: H 0 : H : No ll = 0 TS: = MST/MSE RR: ;-1,(-1)(-1)

78 ssume: Mxed Effecs Models cor xed (ll levels of neres n sudy) 1 0 cor B Rndom (Smple of levels used n sudy) j ~ N(0,s ) (Independen) B Inercon erms Rndom ) j ~ N(0,s (Independen) nlyss of Vrnce s compued excly s n xed Effecs cse (Sums of Squres, df s, MS s) Error erms for ess chnge (See nex slde).

79 Expeced Men Squres for -Wy NOV Men Squre df xed Model Rndom Model Mxed Model ( xed, B xed) ( Rndom, B Rndom) ( xed, B Rndom) n n 1 1 MS -1 s + s ns ns s ns 1 1 MSB -1 s + n j j1 s ns ns s ns ns 1 n j 1 j1 MSB 1 1 s + s ns s ns 11 MSE n 1 s s s

80 NOV pproch Mxed Effecs Source df SS MS cor -1 SS MS=SS/(-1) =MS/MSB cor B -1 SSB MSB=SSB/(-1) B=MSB/MSB Inercon (-1)(-1) SSB MSB=SSB/[(-1)(-1)] B=MSB/MSE Error (n-1) SSE MSE=SSE/[(n-1)] Tol n-1 TSS Procedure: rs es for nercon effecs If nercon es no sgnfcn, es for cor nd B effecs Tes for Inercon: Tes for cor Tes for cor B H : s 0 H :... 0 H : s 0 H : s 0 H : No ll 0 H : s 0 MSB MS MSB TS : B TS : TS : B MSE MSB MSB RR : RR : RR : B,( 1)( 1), ( n1),( 1),( 1)( 1) B,( 1),( 1)( 1)

81 Comprng Mn Effecs for (No Inercon) Tukey s Mehod- q n Sudenzed Rnge Tle wh = (-1)(-1) W q (, v) j MSB n Conclude: f y y W j.. j.. Tukey's CI: ( ) : y y W j.. j.. Bonferron s Mehod - -vlues n Bonferron le wh = (-1)(-1) B j Conclude: /, ( 1)/, v MSB n f y y B j.. j.. Bonferron's CI:( α - α ): y y B j.. j.. j j j j

82 ssume: Rndom Effecs Models cor Rndom (Smple of levels used n sudy) ~ N(0,s ) (Independen) cor B Rndom (Smple of levels used n sudy) j ~ N(0,s ) (Independen) B Inercon erms Rndom ) j ~ N(0,s (Independen) nlyss of Vrnce s compued excly s n xed Effecs cse (Sums of Squres, df s, MS s) Error erms for ess chnge (See nex slde).

83 NOV pproch Rndom Effecs Source df SS MS cor -1 SS MS=SS/(-1) =MS/MSB cor B -1 SSB MSB=SSB/(-1) B =MSB/MSB Inercon (-1)(-1) SSB MSB=SSB/[(-1)(-1)] B =MSB/MSE Error (n-1) SSE MSE=SSE/[(n-1)] Tol n-1 TSS Procedure: rs es for nercon effecs If nercon es no sgnfcn, es for cor nd B effecs Tes for Inercon: Tes for cor Tes for cor B H : s 0 H : s 0 H : s 0 H : s 0 H : s 0 H : s 0 MSB MS MSB TS : B TS : TS : B MSE MSB MSB RR : RR : RR : B,( 1)( 1), ( n1),( 1),( 1)( 1) B ( 1),( 1)( 1)

84 Nesed Desgns Desgns where levels of one fcor re nesed (s opposed o crossed) wr oher fcor Exmples Include: Clssrooms nesed whn schools Lers nesed whn eed Vrees Hr swches nesed whn shmpoo ypes Swmps of vryng szes (e.g. lrge, medum, smll) Resurns nesed whn nonl chns

85 Nesed Desgn - Model Y 1,... j 1,..., k 1,..., n jk j() jk where: Y jk jk j () h h h Response for k rep of cor level, B j level whn Overll Men h Effec of level of (xed or Rndom) h h Effec of j level of B whn level of (xed or Rndom) h Rndom error erm for k rep when s, B s j()

86 Nesed Desgn - NOV Tol Vron: n TSS Y Y df n 1 cor : jk... Tol 1 j1 k 1 1 SS n Y.. Y... df 1 1 cor B Nesed Whn SSB( ) n Y Y df Error: SSE Y Y j... B( ) 1 j1 1 n jk j. E 1 j1 k1 1 df ( n 1)

87 Esmors, nlyss of Vrnce, -ess ^ ^ ^ Y Y Y Y Y j () j ^ ^ ^ ^ ed Vlues: Y = Y Y Y Y Y Y jk j () j j ^ n jk jk jk jk j jk j E 1 j1 k 1 Resduls: e Y Y Y Y SSE Y Y df n 1 cor Sum of Squres: SS n Y Y df j1 j B( ) cor B Whn Sum of Squres: SSB( ) n Y Y df 1 xed, B xed xed, B Rndom Rndom, B Rndom n n ( ) ( ) 1 1 E MS s E MS s ns E MS s ns ns 1 1 E MSB( ) n j () 1 j1 EMSB( ) s ns ( ) EMSB( ) s ns ( ) s 1 s EMSE E MSE s s E MSE MS MS MS H0 : TS : H0 : TS : H0 : s 0 TS : MSE MSB( ) MSB( ) MSB( ) MSB( ) MSB( ) H0 : 1(1)... ( ) 0 TS : B ( ) H0 : s ( ) 0 TS : B ( ) H0 : s ( ) 0 TS : B ( ) MSE MSE MSE

88 cors nd B xed 0 0 1,..., ~ N 0, s j() jk 1 j1 Tess for Dfferences mong cor Effecs H :... 0 H : No ll MS Tes Ssc: P-vlue: P MSE Rejecon Regon:, 1,( n1) Tess for Dfferences mong cor B Effecs H : 0, j H : No ll 0 0 j( ) j( ) MSB( ) Tes Ssc: P-vlue: P MSE Rejecon Regon: B( ),,( n1) B( ) B( )

89 xed Effecs Model ( nd B xed) MS cor : H0 : TS : RR : 1 ; 1, n 1 MSE s MSE EY s Y s Y n n Conrss mong Mens of : c 0 1 ^ MSE ^ L c c L c Y s L c n 1 ^ MSE 1 100% CI for L: L 1 / ; n 1 c n ll Possle Comprsons mong prs of mens of : c 1 1 C 1 MSE Tukey: Y Y ' q 1 ;, n 1 Bonferron: Y Y ' 1 ; n 1 n C MSE n MSB( ) cor B(): H0 : 1(1)... ( ) 0 TS : B ( ) RR : B ( ) 1 ; 1, n 1 MSE s MSE EY j j j() s Y j s Y j n n 1 ll Possle Comprsons mong prs of mens of B whn gven level ( ) of : CB( ) MSE MSE Tukey: Y j Y j ' q 1 ;, n 1 Bonferron: Y j Y j ' 1 ; n 1 n C B( ) n

90 Comprng Mn Effecs for Tukey s Mehod- q n Sudenzed Rnge Tle wh = (n-1)s W q (, v) j MSE 1 1 n n Conclude: f y y W j.. j.. Tukey's CI: ( ) : y y W j.. j.. Bonferron s Mehod - -vlues n Bonferron le wh = (n-1)s 1 1 B j /, ( 1)/, v MSE n n j Conclude: j.. j.. j f y y B Bonferron's CI:( α - α ): y y B j.. j.. j j j j

91 Comprng Effecs for cor B Whn Tukey s Mehod- q n Sudenzed Rnge Tle wh = (n-1)s B W q (, v) j ( k ) k Conclude: MSE n f y y W ( k ) j( k ) k. kj. Tukey's CI: ( ) : y y W ( k ) j( k ) k. kj. Bonferron s Mehod - -vlues n Bonferron le wh = (n-1)s B B j ( k ) /, k ( k 1)/, v MSE n Conclude: ( k ) j( k ) k. kj. f y y B ( k ) j( k ) k. kj. j ( k ) B j ( k ) Bonferron's CI:( - ): y y B B B j ( k ) B j ( k )

92 cor xed nd B Rndom r B B B jk j P MSE MSB H j H P MSB MS H H N N 1),(, ) ( ) ( ) ( 0 1,, 1 0 ) ( 1 : Regon Rejecon : P - vlue ) ( : Tes Ssc 0 :, 0 : cor B Effecs mong s for Dfference Tess : Regon Rejecon : P - vlue ) ( : Tes Ssc 0 No ll : 0... : cor Effecs mong s for Dfference Tess 0, ~ 0, ~ 0 s s s s

93 Mxed Effecs Model ( xed nd B Rndom) MS cor : H0 : TS : RR : 1 ; 1, 1 MSB( ) s ns ( ) MSB( ) E Y s Y s Y n n 1 C ll Possle Comprsons mong prs of mens of : c 1 Tukey: ( ) MSB( ) q 1 ;, 1 Bonferron: Y Y 1 ; 1 n C n MSB Y Y' ' MSB( ) cor B(): H0 : s( ) 0 TS : B ( ) RR : B ( ) 1 ; 1, n 1 MSE 1 1 EMSB( ) s ns ( ) EMSE s EMSB( ) EMSE s n n 1 1 s( ) ( ) ( ) pproxme df (Serhwe): ( ) s MSB MSE df n n MSB( ) MSE n n 1 n 1 pproxme 1 100% CI or s : ( ) df s 1 ; df ( ) ( ) ( ) df ( ) s ( ), ; df ( ) ( )

94 Comprng Mn Effecs for (B Rndom) Tukey s Mehod- q n Sudenzed Rnge Tle wh = S - W q (, v) j MSB( ) 1 1 n n Conclude: f y y W j.. j.. Tukey's CI: ( ) : y y W j.. j.. Bonferron s Mehod - -vlues n Bonferron le wh = S - B /, ( 1)/, ( ) j v MSB n n j Conclude: j.. j.. j.. j f y y B Bonferron's CI:( α - α ): y y B j j j j j

95 cors nd B Rndom r B B B jk j P MSE MSB H j H P MSB MS H H N N N 1),(, ) ( ) ( ) ( 0 1,, 0 ) ( : Regon Rejecon : vlue P - ) ( : Tes Ssc 0 :, 0 : cor B Effecs mong s for Dfference Tess : Regon Rejecon : vlue P - ) ( : Tes Ssc 0 : 0 : cor Effecs mong s for Dfference Tess 0, ~ 0, ~ 0, ~ s s s s s s s

96 Rndom Effecs Model ( nd B Rndom) MS MSB( ) cor : H0 : s 0 TS : RR : 1 ; 1, 1 ( ) E MS s ns ns E MSB s ns ( ) ( ) 1 1 EMS) EMSB( ) s n n 1 1 s MS MSB( ) pproxme df (Serhwe): df n n df s df s pproxme 1100% CI or s :, 1 ; df ; df cor B Whn s sme s n Mxed Effecs Model s MS MSB( ) n n 1 1

97 Elemens of Spl-Plo Desgns Spl-Plo Expermen: corl desgn wh les fcors, where expermenl uns wr fcors dffer n sze or oservonl pons. Whole plo: Lrges expermenl un Whole Plo cor: cor h hs levels ssgned o whole plos. Cn e exended o or more fcors Suplo: Expermenl uns h he whole plo s spl no (where oservons re mde) Suplo cor: cor h hs levels ssgned o suplos Blocks: ggreges of whole plos h receve ll levels of whole plo fcor

98 Spl Plo Desgn Block 1 Block Block 3 Block 4 =1, B=1 =1, B=1 =1, B=1 =1, B=1 =1, B= =1, B= =1, B= =1, B= =1, B=3 =1, B=3 =1, B=3 =1, B=3 =1, B=4 =1, B=4 =1, B=4 =1, B=4 =, B=1 =, B=1 =, B=1 =, B=1 =, B= =, B= =, B= =, B= =, B=3 =, B=3 =, B=3 =, B=3 =, B=4 =, B=4 =, B=4 =, B=4 =3, B=1 =3, B=1 =3, B=1 =3, B=1 =3, B= =3, B= =3, B= =3, B= =3, B=3 =3, B=3 =3, B=3 =3, B=3 =3, B=4 =3, B=4 =3, B=4 =3, B=4 Noe: Whn ech lock we would ssgn rndom he 3 levels of o he whole plos nd he 4 levels of B o he suplos whn whole plos

99 Exmples grculure: Vrees of crop or gs my need o e grown n lrge res, whle vrees of ferlzer or vryng growh perods my e oserved n suses of he re. Engneerng: My need long heng perods for process nd my e le o compre severl formulons of y-produc whn ech level of he heng fcor. Behvorl Scences: Mny sudes nvolve repeed mesuremens on he sme sujecs nd re nlyzed s spl-plo (See Repeed Mesures lecure)

100 Desgn Srucure Blocks: groups of expermenl uns o e exposed o ll comnons of whole plo nd suplo fcors Whole plos: expermenl uns o whch he whole plo fcor levels wll e ssgned o rndom whn locks Suplos: c suuns whn whole plos o whch he suplo fcor levels wll e ssgned o rndom. ully lnced expermen wll hve n=c oservons

101 D Elemens (xed cors, Rndom Blocks) Y jk : Oservon from wp, lock j, nd sp k : Overll men level : Effec of h level of whole plo fcor (xed) j : Effec of j h lock (Rndom) ( ) j : Rndom error correspondng o whole plo elemens n lock j where wp s ppled k : Effec of k h level of suplo fcor (xed) ) k : Inercon wn wp nd sp k (c ) jk : Inercon wn lock j nd sp k (ofen se o 0) jk : Rndom Error= (c ) jk + (c ) jk Noe h f lock/sp nercon s ssumed o e 0, represens he lock/sp whn wp nercon

102 Model nd Common ssumpons Y jk = + + j + ( ) j + k + ( ) k + jk 1 j ~ 0 NID(0, s ) ( ) j ~ NID(0, s ) c k 1 k 0 1 ( ) k c k 1 ( ) k 0 jk ~ COV NID(0. ( j s,( ) j ) ) COV ( j, jk ) COV (( ) j, jk ) 0

103 Tess for xed Effecs ) ~ ( Ssc: Tes, 0 ) ( : SP Inercon : WP ) ~ ( Ssc: Tes 0 : Suplo Tr Effecs : ) ~ ( Ssc: Tes 0 : Plo Tr Effecs : Whole 1) 1)( ( 1). 1)( ( 0 1) 1)( ( ) 1)( 1.( * 1 0 c c SP WP SP WP ERROR SP WP SP WP k c c SP SP ERROR SP SP c WP WP WP BLOCK WP WP P P MS MS k H P P MS MS H P P MS MS H

104 Comprng cor Levels Whole Plo cor Levels : 95% CI for 95% CI for : Y.. Y '.. Su Plo cor Levels : ( ) ( ) : Y. k Y. k ' MSBLOCK WP ( c 1) MSERROR Y. k Y '. k ( c 1) MSERROR MSBLOCK WP ( c 1) MS MS ( 1)( c 1) BLOCK WP ( 1)( 1) MS c BLOCK WP MSERROR 95% CI for ( k k ') : Y.. k Y.. k ' Su Plo Effecs Whn sme whole plo (Inercon) : Whole Plo Effecs whn sme su plo (Inercon) : ^ k ' k ' ERROR k c k ' (df MS ERROR gven elow)

105 Repeed Mesures Desgns Tremens/Condons o compre N sujecs o e ncluded n sudy (ech sujec wll receve only one remen) n sujecs receve r : n = N me perods of d wll e oned Effecs of r, me nd rxme nercon of prmry neres. Beween Sujec cor: Tremen Whn Sujec cors: Tme, TrxTme

106 Model Y jk j( ) overll men k ( ) k jk effec of r j( ) effec of j h 1 0 sujec n r j( ) ~ NID 0, s k effec of k h me perod k 1 k 0 ( ) jk k nercon eween r rndom error erm jk ~ nd me k NID 0, s 1 ( ) k k 1 ( ) k 0 Noe he rndom error erm s cully he nercon eween sujecs (whn remens) nd me

107 Tess for xed Effecs ) ~ ( Ssc: Tes, 0 ) ( : Tremen/Tme Inercon : ) ~ ( Ssc: Tes 0 : Tme Effecs : ) ~ ( Ssc: Tes 0 : Tremen Effecs : 1) 1)( ( 1), 1)( ( 0 1) 1)( ( 1, 1 0 1) ( 1, ) ( 1 0 n TIME TRT TIME TRT ERROR TIME TRT TIME TRT k n TIME TIME ERROR TIME TIME n TRTS TRTS TRTS SUBJECTS TRTS TRTS P P MS MS k H P P MS MS H P P MS MS H

108 Y Comprng cor Levels Comprng Tremen Levels : 95% CI for : Comprng Tme Levels : wh pproxme df : ( 1) MSERROR MSSUBJECT ( TRT ) ( 1) MS MSSUBJECT ( TRT ) ( n 1)( 1) Y MSERROR 95% CI for k k ' : Y.. k Y.. k ' n Comprng Tremen Levels Whn Tme Levels : Y ^. k '. k ' MS ERROR Y.. '.. SUBJECTS ( TRTS ) n ( ( n 1) MS 1) MS SUBJECTS ( TRTS ) n ERROR

TSS = SST + SSE An orthogonal partition of the total SS

TSS = SST + SSE An orthogonal partition of the total SS ANOVA: Topc 4. Orhogonal conrass [ST&D p. 183] H 0 : µ 1 = µ =... = µ H 1 : The mean of a leas one reamen group s dfferen To es hs hypohess, a basc ANOVA allocaes he varaon among reamen means (SST) equally

More information

Simplified Variance Estimation for Three-Stage Random Sampling

Simplified Variance Estimation for Three-Stage Random Sampling Deprmen of ppled Sscs Johnnes Kepler Unversy Lnz IFS Reserch Pper Seres 04-67 Smplfed rnce Esmon for Three-Sge Rndom Smplng ndres Quember Ocober 04 Smplfed rnce Esmon for Three-Sge Rndom Smplng ndres Quember

More information

Interval Estimation. Consider a random variable X with a mean of X. Let X be distributed as X X

Interval Estimation. Consider a random variable X with a mean of X. Let X be distributed as X X ECON 37: Ecoomercs Hypohess Tesg Iervl Esmo Wh we hve doe so fr s o udersd how we c ob esmors of ecoomcs reloshp we wsh o sudy. The queso s how comforble re we wh our esmors? We frs exme how o produce

More information

In the complete model, these slopes are ANALYSIS OF VARIANCE FOR THE COMPLETE TWO-WAY MODEL. (! i+1 -! i ) + [(!") i+1,q - [(!

In the complete model, these slopes are ANALYSIS OF VARIANCE FOR THE COMPLETE TWO-WAY MODEL. (! i+1 -! i ) + [(!) i+1,q - [(! ANALYSIS OF VARIANCE FOR THE COMPLETE TWO-WAY MODEL The frs hng o es n wo-way ANOVA: Is here neracon? "No neracon" means: The man effecs model would f. Ths n urn means: In he neracon plo (wh A on he horzonal

More information

Minimum Squared Error

Minimum Squared Error Minimum Squred Error LDF: Minimum Squred-Error Procedures Ide: conver o esier nd eer undersood prolem Percepron y i > for ll smples y i solve sysem of liner inequliies MSE procedure y i = i for ll smples

More information

Minimum Squared Error

Minimum Squared Error Minimum Squred Error LDF: Minimum Squred-Error Procedures Ide: conver o esier nd eer undersood prolem Percepron y i > 0 for ll smples y i solve sysem of liner inequliies MSE procedure y i i for ll smples

More information

4.8 Improper Integrals

4.8 Improper Integrals 4.8 Improper Inegrls Well you ve mde i hrough ll he inegrion echniques. Congrs! Unforunely for us, we sill need o cover one more inegrl. They re clled Improper Inegrls. A his poin, we ve only del wih inegrls

More information

EEM 486: Computer Architecture

EEM 486: Computer Architecture EEM 486: Compuer Archecure Lecure 4 ALU EEM 486 MIPS Arhmec Insrucons R-ype I-ype Insrucon Exmpe Menng Commen dd dd $,$2,$3 $ = $2 + $3 sub sub $,$2,$3 $ = $2 - $3 3 opernds; overfow deeced 3 opernds;

More information

ENGI 4421 Probability and Statistics Faculty of Engineering and Applied Science Problem Set 10 Solutions Chi-Square Tests; Simple Linear Regression

ENGI 4421 Probability and Statistics Faculty of Engineering and Applied Science Problem Set 10 Solutions Chi-Square Tests; Simple Linear Regression ENGI 441 Prbbly nd Sscs Fculy f Engneerng nd Appled Scence Prblem Se 10 Sluns Ch-Squre Tess; Smple Lner Regressn 1. Is he fllwng se f bservns f bjecs n egh dfferen drecns cnssen wh unfrm dsrbun? Drecn

More information

Applied Statistics Qualifier Examination

Applied Statistics Qualifier Examination Appled Sttstcs Qulfer Exmnton Qul_june_8 Fll 8 Instructons: () The exmnton contns 4 Questons. You re to nswer 3 out of 4 of them. () You my use ny books nd clss notes tht you mght fnd helpful n solvng

More information

September 20 Homework Solutions

September 20 Homework Solutions College of Engineering nd Compuer Science Mechnicl Engineering Deprmen Mechnicl Engineering A Seminr in Engineering Anlysis Fll 7 Number 66 Insrucor: Lrry Creo Sepember Homework Soluions Find he specrum

More information

Physics 201 Lecture 2

Physics 201 Lecture 2 Physcs 1 Lecure Lecure Chper.1-. Dene Poson, Dsplcemen & Dsnce Dsngush Tme nd Tme Inerl Dene Velocy (Aerge nd Insnneous), Speed Dene Acceleron Undersnd lgebrclly, hrough ecors, nd grphclly he relonshps

More information

Supporting information How to concatenate the local attractors of subnetworks in the HPFP

Supporting information How to concatenate the local attractors of subnetworks in the HPFP n Effcen lgorh for Idenfyng Prry Phenoype rcors of Lrge-Scle Boolen Newor Sng-Mo Choo nd Kwng-Hyun Cho Depren of Mhecs Unversy of Ulsn Ulsn 446 Republc of Kore Depren of Bo nd Brn Engneerng Kore dvnced

More information

We are estimating the density of long distant migrant (LDM) birds in wetlands along Lake Michigan.

We are estimating the density of long distant migrant (LDM) birds in wetlands along Lake Michigan. Ch 17 Random ffecs and Mxed Models 17. Random ffecs Models We are esmang he densy of long dsan mgran (LDM) brds n welands along Lake Mchgan. μ + = LDM per hecaren h weland ~ N(0, ) The varably of expeced

More information

rank Additionally system of equation only independent atfect Gawp (A) possible ( Alb ) easily process form rang A. Proposition with Definition

rank Additionally system of equation only independent atfect Gawp (A) possible ( Alb ) easily process form rang A. Proposition with Definition Defiion nexivnol numer ler dependen rows mrix sid row Gwp elimion mehod does no fec h numer end process i possile esily red rng fc for mrix form der zz rn rnk wih m dcussion i holds rr o Proposiion ler

More information

AN INTRODUCTORY GUIDELINE FOR THE USE OF BAYESIAN STATISTICAL METHODS IN THE ANALYSIS OF ROAD TRAFFIC ACCIDENT DATA

AN INTRODUCTORY GUIDELINE FOR THE USE OF BAYESIAN STATISTICAL METHODS IN THE ANALYSIS OF ROAD TRAFFIC ACCIDENT DATA AN INTRODUCTORY GUIDELINE FOR THE USE OF BAYESIAN STATISTICAL METHODS IN THE ANALYSIS OF ROAD TRAFFIC ACCIDENT DATA CJ Molle Pr Eng Chef Engneer : Trffc Engneerng Deprmen of Economc Affrs, Agrculure nd

More information

e t dt e t dt = lim e t dt T (1 e T ) = 1

e t dt e t dt = lim e t dt T (1 e T ) = 1 Improper Inegrls There re wo ypes of improper inegrls - hose wih infinie limis of inegrion, nd hose wih inegrnds h pproch some poin wihin he limis of inegrion. Firs we will consider inegrls wih infinie

More information

Data Collection Definitions of Variables - Conceptualize vs Operationalize Sample Selection Criteria Source of Data Consistency of Data

Data Collection Definitions of Variables - Conceptualize vs Operationalize Sample Selection Criteria Source of Data Consistency of Data Apply Sascs and Economercs n Fnancal Research Obj. of Sudy & Hypoheses Tesng From framework objecves of sudy are needed o clarfy, hen, n research mehodology he hypoheses esng are saed, ncludng esng mehods.

More information

Student Activity 3: Single Factor ANOVA

Student Activity 3: Single Factor ANOVA MATH 40 Student Activity 3: Single Fctor ANOVA Some Bsic Concepts In designed experiment, two or more tretments, or combintions of tretments, is pplied to experimentl units The number of tretments, whether

More information

Decompression diagram sampler_src (source files and makefiles) bin (binary files) --- sh (sample shells) --- input (sample input files)

Decompression diagram sampler_src (source files and makefiles) bin (binary files) --- sh (sample shells) --- input (sample input files) . Iroduco Probblsc oe-moh forecs gudce s mde b 50 esemble members mproved b Model Oupu scs (MO). scl equo s mde b usg hdcs d d observo d. We selec some prmeers for modfg forecs o use mulple regresso formul.

More information

FINANCIAL ECONOMETRICS

FINANCIAL ECONOMETRICS FINANCIAL ECONOMETRICS SPRING 07 WEEK IV NONLINEAR MODELS Prof. Dr. Burç ÜLENGİN Nonlner NONLINEARITY EXISTS IN FINANCIAL TIME SERIES ESPECIALLY IN VOLATILITY AND HIGH FREQUENCY DATA LINEAR MODEL IS DEFINED

More information

Advanced Electromechanical Systems (ELE 847)

Advanced Electromechanical Systems (ELE 847) (ELE 847) Dr. Smr ouro-rener Topc 1.4: DC moor speed conrol Torono, 2009 Moor Speed Conrol (open loop conrol) Consder he followng crcu dgrm n V n V bn T1 T 5 T3 V dc r L AA e r f L FF f o V f V cn T 4

More information

Hidden Markov Model. a ij. Observation : O1,O2,... States in time : q1, q2,... All states : s1, s2,..., sn

Hidden Markov Model. a ij. Observation : O1,O2,... States in time : q1, q2,... All states : s1, s2,..., sn Hdden Mrkov Model S S servon : 2... Ses n me : 2... All ses : s s2... s 2 3 2 3 2 Hdden Mrkov Model Con d Dscree Mrkov Model 2 z k s s s s s s Degree Mrkov Model Hdden Mrkov Model Con d : rnson roly from

More information

3. Renewal Limit Theorems

3. Renewal Limit Theorems Virul Lborories > 14. Renewl Processes > 1 2 3 3. Renewl Limi Theorems In he inroducion o renewl processes, we noed h he rrivl ime process nd he couning process re inverses, in sens The rrivl ime process

More information

Properties of Logarithms. Solving Exponential and Logarithmic Equations. Properties of Logarithms. Properties of Logarithms. ( x)

Properties of Logarithms. Solving Exponential and Logarithmic Equations. Properties of Logarithms. Properties of Logarithms. ( x) Properies of Logrihms Solving Eponenil nd Logrihmic Equions Properies of Logrihms Produc Rule ( ) log mn = log m + log n ( ) log = log + log Properies of Logrihms Quoien Rule log m = logm logn n log7 =

More information

ENGR 1990 Engineering Mathematics The Integral of a Function as a Function

ENGR 1990 Engineering Mathematics The Integral of a Function as a Function ENGR 1990 Engineering Mhemics The Inegrl of Funcion s Funcion Previously, we lerned how o esime he inegrl of funcion f( ) over some inervl y dding he res of finie se of rpezoids h represen he re under

More information

Introduction to Numerical Integration Part II

Introduction to Numerical Integration Part II Introducton to umercl Integrton Prt II CS 75/Mth 75 Brn T. Smth, UM, CS Dept. Sprng, 998 4/9/998 qud_ Intro to Gussn Qudrture s eore, the generl tretment chnges the ntegrton prolem to ndng the ntegrl w

More information

A Kalman filtering simulation

A Kalman filtering simulation A Klmn filering simulion The performnce of Klmn filering hs been esed on he bsis of wo differen dynmicl models, ssuming eiher moion wih consn elociy or wih consn ccelerion. The former is epeced o beer

More information

FM Applications of Integration 1.Centroid of Area

FM Applications of Integration 1.Centroid of Area FM Applicions of Inegrion.Cenroid of Are The cenroid of ody is is geomeric cenre. For n ojec mde of uniform meril, he cenroid coincides wih he poin which he ody cn e suppored in perfecly lnced se ie, is

More information

1.0 Electrical Systems

1.0 Electrical Systems . Elecricl Sysems The ypes of dynmicl sysems we will e sudying cn e modeled in erms of lgeric equions, differenil equions, or inegrl equions. We will egin y looking fmilir mhemicl models of idel resisors,

More information

Department of Economics University of Toronto

Department of Economics University of Toronto Deparmen of Economcs Unversy of Torono ECO408F M.A. Economercs Lecure Noes on Heeroskedascy Heeroskedascy o Ths lecure nvolves lookng a modfcaons we need o make o deal wh he regresson model when some of

More information

F-Tests and Analysis of Variance (ANOVA) in the Simple Linear Regression Model. 1. Introduction

F-Tests and Analysis of Variance (ANOVA) in the Simple Linear Regression Model. 1. Introduction ECOOMICS 35* -- OTE 9 ECO 35* -- OTE 9 F-Tess and Analyss of Varance (AOVA n he Smple Lnear Regresson Model Inroducon The smple lnear regresson model s gven by he followng populaon regresson equaon, or

More information

5.1-The Initial-Value Problems For Ordinary Differential Equations

5.1-The Initial-Value Problems For Ordinary Differential Equations 5.-The Iniil-Vlue Problems For Ordinry Differenil Equions Consider solving iniil-vlue problems for ordinry differenil equions: (*) y f, y, b, y. If we know he generl soluion y of he ordinry differenil

More information

Motion. Part 2: Constant Acceleration. Acceleration. October Lab Physics. Ms. Levine 1. Acceleration. Acceleration. Units for Acceleration.

Motion. Part 2: Constant Acceleration. Acceleration. October Lab Physics. Ms. Levine 1. Acceleration. Acceleration. Units for Acceleration. Moion Accelerion Pr : Consn Accelerion Accelerion Accelerion Accelerion is he re of chnge of velociy. = v - vo = Δv Δ ccelerion = = v - vo chnge of velociy elpsed ime Accelerion is vecor, lhough in one-dimensionl

More information

THEORETICAL AUTOCORRELATIONS. ) if often denoted by γ. Note that

THEORETICAL AUTOCORRELATIONS. ) if often denoted by γ. Note that THEORETICAL AUTOCORRELATIONS Cov( y, y ) E( y E( y))( y E( y)) ρ = = Var( y) E( y E( y)) =,, L ρ = and Cov( y, y ) s ofen denoed by whle Var( y ) f ofen denoed by γ. Noe ha γ = γ and ρ = ρ and because

More information

14.3 comparing two populations: based on independent samples

14.3 comparing two populations: based on independent samples Chpter4 Nonprmetric Sttistics Introduction: : methods for mking inferences bout popultion prmeters (confidence intervl nd hypothesis testing) rely on the ssumptions bout probbility distribution of smpled

More information

Chapter 9: Inferences based on Two samples: Confidence intervals and tests of hypotheses

Chapter 9: Inferences based on Two samples: Confidence intervals and tests of hypotheses Chpter 9: Inferences bsed on Two smples: Confidence intervls nd tests of hypotheses 9.1 The trget prmeter : difference between two popultion mens : difference between two popultion proportions : rtio of

More information

Origin Destination Transportation Models: Methods

Origin Destination Transportation Models: Methods In Jr. of Mhemcl Scences & Applcons Vol. 2, No. 2, My 2012 Copyrgh Mnd Reder Publcons ISSN No: 2230-9888 www.journlshub.com Orgn Desnon rnsporon Models: Mehods Jyo Gup nd 1 N H. Shh Deprmen of Mhemcs,

More information

CISE 301: Numerical Methods Lecture 5, Topic 4 Least Squares, Curve Fitting

CISE 301: Numerical Methods Lecture 5, Topic 4 Least Squares, Curve Fitting CISE 3: umercl Methods Lecture 5 Topc 4 Lest Squres Curve Fttng Dr. Amr Khouh Term Red Chpter 7 of the tetoo c Khouh CISE3_Topc4_Lest Squre Motvton Gven set of epermentl dt 3 5. 5.9 6.3 The reltonshp etween

More information

Robustness Experiments with Two Variance Components

Robustness Experiments with Two Variance Components Naonal Insue of Sandards and Technology (NIST) Informaon Technology Laboraory (ITL) Sascal Engneerng Dvson (SED) Robusness Expermens wh Two Varance Componens by Ana Ivelsse Avlés avles@ns.gov Conference

More information

Solutions to Problems from Chapter 2

Solutions to Problems from Chapter 2 Soluions o Problems rom Chper Problem. The signls u() :5sgn(), u () :5sgn(), nd u h () :5sgn() re ploed respecively in Figures.,b,c. Noe h u h () :5sgn() :5; 8 including, bu u () :5sgn() is undeined..5

More information

Farey Fractions. Rickard Fernström. U.U.D.M. Project Report 2017:24. Department of Mathematics Uppsala University

Farey Fractions. Rickard Fernström. U.U.D.M. Project Report 2017:24. Department of Mathematics Uppsala University U.U.D.M. Project Report 07:4 Frey Frctions Rickrd Fernström Exmensrete i mtemtik, 5 hp Hledre: Andres Strömergsson Exmintor: Jörgen Östensson Juni 07 Deprtment of Mthemtics Uppsl University Frey Frctions

More information

Chapter 5 Supplemental Text Material R S T. ij i j ij ijk

Chapter 5 Supplemental Text Material R S T. ij i j ij ijk Chpter 5 Supplementl Text Mterl 5-. Expected Men Squres n the Two-fctor Fctorl Consder the two-fctor fxed effects model y = µ + τ + β + ( τβ) + ε k R S T =,,, =,,, k =,,, n gven s Equton (5-) n the textook.

More information

II The Z Transform. Topics to be covered. 1. Introduction. 2. The Z transform. 3. Z transforms of elementary functions

II The Z Transform. Topics to be covered. 1. Introduction. 2. The Z transform. 3. Z transforms of elementary functions II The Z Trnsfor Tocs o e covered. Inroducon. The Z rnsfor 3. Z rnsfors of eleenry funcons 4. Proeres nd Theory of rnsfor 5. The nverse rnsfor 6. Z rnsfor for solvng dfference equons II. Inroducon The

More information

Section 11.5 Estimation of difference of two proportions

Section 11.5 Estimation of difference of two proportions ection.5 Estimtion of difference of two proportions As seen in estimtion of difference of two mens for nonnorml popultion bsed on lrge smple sizes, one cn use CLT in the pproximtion of the distribution

More information

A First Guide to Hypothesis Testing in Linear Regression Models. A Generic Linear Regression Model: Scalar Formulation

A First Guide to Hypothesis Testing in Linear Regression Models. A Generic Linear Regression Model: Scalar Formulation ECON 5* -- A rs Gude o Hypoess Tesng MG Abbo A rs Gude o Hypoess Tesng n Lnear Regresson Models A Generc Lnear Regresson Model: Scalar mulaon e - populaon ( sample observaon, e scalar fmulaon of e PRE

More information

Continuous Random Variable X:

Continuous Random Variable X: Continuous Rndom Vrile : The continuous rndom vrile hs its vlues in n intervl, nd it hs proility distriution unction or proility density unction p.d. stisies:, 0 & d Which does men tht the totl re under

More information

Single-loop System Reliability-Based Design & Topology Optimization (SRBDO/SRBTO): A Matrix-based System Reliability (MSR) Method

Single-loop System Reliability-Based Design & Topology Optimization (SRBDO/SRBTO): A Matrix-based System Reliability (MSR) Method 10 h US Naonal Congress on Compuaonal Mechancs Columbus, Oho 16-19, 2009 Sngle-loop Sysem Relably-Based Desgn & Topology Opmzaon (SRBDO/SRBTO): A Marx-based Sysem Relably (MSR) Mehod Tam Nguyen, Junho

More information

( ) ( ) ( ) ( ) ( ) ( y )

( ) ( ) ( ) ( ) ( ) ( y ) 8. Lengh of Plne Curve The mos fmous heorem in ll of mhemics is he Pyhgoren Theorem. I s formulion s he disnce formul is used o find he lenghs of line segmens in he coordine plne. In his secion you ll

More information

Contraction Mapping Principle Approach to Differential Equations

Contraction Mapping Principle Approach to Differential Equations epl Journl of Science echnology 0 (009) 49-53 Conrcion pping Principle pproch o Differenil Equions Bishnu P. Dhungn Deprmen of hemics, hendr Rn Cmpus ribhuvn Universiy, Khmu epl bsrc Using n eension of

More information

Background and Motivation: Importance of Pressure Measurements

Background and Motivation: Importance of Pressure Measurements Imornce of Pressre Mesremens: Pressre s rmry concern for mny engneerng lcons e.g. lf nd form drg. Cvon : Pressre s of fndmenl mornce n ndersndng nd modelng cvon. Trblence: Velocy-Pressre-Grden ensor whch

More information

f t f a f x dx By Lin McMullin f x dx= f b f a. 2

f t f a f x dx By Lin McMullin f x dx= f b f a. 2 Accumulion: Thoughs On () By Lin McMullin f f f d = + The gols of he AP* Clculus progrm include he semen, Sudens should undersnd he definie inegrl s he ne ccumulion of chnge. 1 The Topicl Ouline includes

More information

Chapter Simpson s 1/3 Rule of Integration. ( x)

Chapter Simpson s 1/3 Rule of Integration. ( x) Cper 7. Smpso s / Rule o Iegro Aer redg s per, you sould e le o. derve e ormul or Smpso s / rule o egro,. use Smpso s / rule o solve egrls,. develop e ormul or mulple-segme Smpso s / rule o egro,. use

More information

INVESTIGATION OF HABITABILITY INDICES OF YTU GULET SERIES IN VARIOUS SEA STATES

INVESTIGATION OF HABITABILITY INDICES OF YTU GULET SERIES IN VARIOUS SEA STATES Brodogrdnj/Shpuldng Volume 65 Numer 3, 214 Ferd Ckc Muhsn Aydn ISSN 7-215X eissn 1845-5859 INVESTIGATION OF HABITABILITY INDICES OF YTU GULET SERIES IN VARIOUS SEA STATES UDC 629.5(5) Professonl pper Summry

More information

Average & instantaneous velocity and acceleration Motion with constant acceleration

Average & instantaneous velocity and acceleration Motion with constant acceleration Physics 7: Lecure Reminders Discussion nd Lb secions sr meeing ne week Fill ou Pink dd/drop form if you need o swich o differen secion h is FULL. Do i TODAY. Homework Ch. : 5, 7,, 3,, nd 6 Ch.: 6,, 3 Submission

More information

Privacy-Preserving Bayesian Network Parameter Learning

Privacy-Preserving Bayesian Network Parameter Learning 4h WSEAS In. Conf. on COMUTATIONAL INTELLIGENCE, MAN-MACHINE SYSTEMS nd CYBERNETICS Mm, Flord, USA, November 7-9, 005 pp46-5) rvcy-reservng Byesn Nework rmeer Lernng JIANJIE MA. SIVAUMAR School of EECS,

More information

THE EXISTENCE OF SOLUTIONS FOR A CLASS OF IMPULSIVE FRACTIONAL Q-DIFFERENCE EQUATIONS

THE EXISTENCE OF SOLUTIONS FOR A CLASS OF IMPULSIVE FRACTIONAL Q-DIFFERENCE EQUATIONS Europen Journl of Mhemcs nd Compuer Scence Vol 4 No, 7 SSN 59-995 THE EXSTENCE OF SOLUTONS FOR A CLASS OF MPULSVE FRACTONAL Q-DFFERENCE EQUATONS Shuyun Wn, Yu Tng, Q GE Deprmen of Mhemcs, Ynbn Unversy,

More information

MIXED MODELS (Sections ) I) In the unrestricted model, interactions are treated as in the random effects model:

MIXED MODELS (Sections ) I) In the unrestricted model, interactions are treated as in the random effects model: 1 2 MIXED MODELS (Sections 17.7 17.8) Exmple: Suppose tht in the fiber breking strength exmple, the four mchines used were the only ones of interest, but the interest ws over wide rnge of opertors, nd

More information

378 Relations Solutions for Chapter 16. Section 16.1 Exercises. 3. Let A = {0,1,2,3,4,5}. Write out the relation R that expresses on A.

378 Relations Solutions for Chapter 16. Section 16.1 Exercises. 3. Let A = {0,1,2,3,4,5}. Write out the relation R that expresses on A. 378 Reltions 16.7 Solutions for Chpter 16 Section 16.1 Exercises 1. Let A = {0,1,2,3,4,5}. Write out the reltion R tht expresses > on A. Then illustrte it with digrm. 2 1 R = { (5,4),(5,3),(5,2),(5,1),(5,0),(4,3),(4,2),(4,1),

More information

1. Find a basis for the row space of each of the following matrices. Your basis should consist of rows of the original matrix.

1. Find a basis for the row space of each of the following matrices. Your basis should consist of rows of the original matrix. Mh 7 Exm - Prcice Prolem Solions. Find sis for he row spce of ech of he following mrices. Yor sis shold consis of rows of he originl mrix. 4 () 7 7 8 () Since we wn sis for he row spce consising of rows

More information

The steps of the hypothesis test

The steps of the hypothesis test ttisticl Methods I (EXT 7005) Pge 78 Mosquito species Time of dy A B C Mid morning 0.0088 5.4900 5.5000 Mid Afternoon.3400 0.0300 0.8700 Dusk 0.600 5.400 3.000 The Chi squre test sttistic is the sum of

More information

Chapter 2: Evaluative Feedback

Chapter 2: Evaluative Feedback Chper 2: Evluive Feedbck Evluing cions vs. insrucing by giving correc cions Pure evluive feedbck depends olly on he cion ken. Pure insrucive feedbck depends no ll on he cion ken. Supervised lerning is

More information

Math 8 Winter 2015 Applications of Integration

Math 8 Winter 2015 Applications of Integration Mth 8 Winter 205 Applictions of Integrtion Here re few importnt pplictions of integrtion. The pplictions you my see on n exm in this course include only the Net Chnge Theorem (which is relly just the Fundmentl

More information

DCDM BUSINESS SCHOOL NUMERICAL METHODS (COS 233-8) Solutions to Assignment 3. x f(x)

DCDM BUSINESS SCHOOL NUMERICAL METHODS (COS 233-8) Solutions to Assignment 3. x f(x) DCDM BUSINESS SCHOOL NUMEICAL METHODS (COS -8) Solutons to Assgnment Queston Consder the followng dt: 5 f() 8 7 5 () Set up dfference tble through fourth dfferences. (b) Wht s the mnmum degree tht n nterpoltng

More information

Electromagnetic Transient Simulation of Large Power Transformer Internal Fault

Electromagnetic Transient Simulation of Large Power Transformer Internal Fault Inernonl Conference on Advnces n Energy nd Envronmenl Scence (ICAEES 5) Elecromgnec Trnsen Smulon of rge Power Trnsformer Inernl Ful Jun u,, Shwu Xo,, Qngsen Sun,c, Huxng Wng,d nd e Yng,e School of Elecrcl

More information

Section 7.2 Two-way ANOVA with random effect(s)

Section 7.2 Two-way ANOVA with random effect(s) Secto 7. Two-wy ANOVA wth rdom effect(s) 1 1. Model wth Two Rdom ffects The fctors hgher-wy ANOVAs c g e cosdered fxed or rdom depedg o the cotext of the study. or ech fctor: Are the levels of tht fctor

More information

Rank One Update And the Google Matrix by Al Bernstein Signal Science, LLC

Rank One Update And the Google Matrix by Al Bernstein Signal Science, LLC Introducton Rnk One Updte And the Google Mtrx y Al Bernsten Sgnl Scence, LLC www.sgnlscence.net here re two dfferent wys to perform mtrx multplctons. he frst uses dot product formulton nd the second uses

More information

CH.3. COMPATIBILITY EQUATIONS. Continuum Mechanics Course (MMC) - ETSECCPB - UPC

CH.3. COMPATIBILITY EQUATIONS. Continuum Mechanics Course (MMC) - ETSECCPB - UPC CH.3. COMPATIBILITY EQUATIONS Connuum Mechancs Course (MMC) - ETSECCPB - UPC Overvew Compably Condons Compably Equaons of a Poenal Vecor Feld Compably Condons for Infnesmal Srans Inegraon of he Infnesmal

More information

Magnetostatics Bar Magnet. Magnetostatics Oersted s Experiment

Magnetostatics Bar Magnet. Magnetostatics Oersted s Experiment Mgneosics Br Mgne As fr bck s 4500 yers go, he Chinese discovered h cerin ypes of iron ore could rc ech oher nd cerin mels. Iron filings "mp" of br mgne s field Crefully suspended slivers of his mel were

More information

MTH 146 Class 11 Notes

MTH 146 Class 11 Notes 8.- Are of Surfce of Revoluion MTH 6 Clss Noes Suppose we wish o revolve curve C round n is nd find he surfce re of he resuling solid. Suppose f( ) is nonnegive funcion wih coninuous firs derivive on he

More information

The Characterization of Jones Polynomial. for Some Knots

The Characterization of Jones Polynomial. for Some Knots Inernon Mhemc Forum,, 8, no, 9 - The Chrceron of Jones Poynom for Some Knos Mur Cncn Yuuncu Y Ünversy, Fcuy of rs nd Scences Mhemcs Deprmen, 8, n, Turkey m_cencen@yhoocom İsm Yr Non Educon Mnsry, 8, n,

More information

P 3 (x) = f(0) + f (0)x + f (0) 2. x 2 + f (0) . In the problem set, you are asked to show, in general, the n th order term is a n = f (n) (0)

P 3 (x) = f(0) + f (0)x + f (0) 2. x 2 + f (0) . In the problem set, you are asked to show, in general, the n th order term is a n = f (n) (0) 1 Tylor polynomils In Section 3.5, we discussed how to pproximte function f(x) round point in terms of its first derivtive f (x) evluted t, tht is using the liner pproximtion f() + f ()(x ). We clled this

More information

For the percentage of full time students at RCC the symbols would be:

For the percentage of full time students at RCC the symbols would be: Mth 17/171 Chpter 7- ypothesis Testing with One Smple This chpter is s simple s the previous one, except it is more interesting In this chpter we will test clims concerning the sme prmeters tht we worked

More information

Comparing Means: t-tests for One Sample & Two Related Samples

Comparing Means: t-tests for One Sample & Two Related Samples Comparing Means: -Tess for One Sample & Two Relaed Samples Using he z-tes: Assumpions -Tess for One Sample & Two Relaed Samples The z-es (of a sample mean agains a populaion mean) is based on he assumpion

More information

How do we solve these things, especially when they get complicated? How do we know when a system has a solution, and when is it unique?

How do we solve these things, especially when they get complicated? How do we know when a system has a solution, and when is it unique? XII. LINEAR ALGEBRA: SOLVING SYSTEMS OF EQUATIONS Tody we re going to tlk out solving systems of liner equtions. These re prolems tht give couple of equtions with couple of unknowns, like: 6= x + x 7=

More information

Principle Component Analysis

Principle Component Analysis Prncple Component Anlyss Jng Go SUNY Bufflo Why Dmensonlty Reducton? We hve too mny dmensons o reson bout or obtn nsghts from o vsulze oo much nose n the dt Need to reduce them to smller set of fctors

More information

T-Match: Matching Techniques For Driving Yagi-Uda Antennas: T-Match. 2a s. Z in. (Sections 9.5 & 9.7 of Balanis)

T-Match: Matching Techniques For Driving Yagi-Uda Antennas: T-Match. 2a s. Z in. (Sections 9.5 & 9.7 of Balanis) 3/0/018 _mch.doc Pge 1 of 6 T-Mch: Mching Techniques For Driving Ygi-Ud Anenns: T-Mch (Secions 9.5 & 9.7 of Blnis) l s l / l / in The T-Mch is shun-mching echnique h cn be used o feed he driven elemen

More information

MODELLING AND EXPERIMENTAL ANALYSIS OF MOTORCYCLE DYNAMICS USING MATLAB

MODELLING AND EXPERIMENTAL ANALYSIS OF MOTORCYCLE DYNAMICS USING MATLAB MODELLING AND EXPERIMENTAL ANALYSIS OF MOTORCYCLE DYNAMICS USING MATLAB P. Florn, P. Vrání, R. Čermá Fculy of Mechncl Engneerng, Unversy of Wes Bohem Asrc The frs pr of hs pper s devoed o mhemcl modellng

More information

1B40 Practical Skills

1B40 Practical Skills B40 Prcticl Skills Comining uncertinties from severl quntities error propgtion We usully encounter situtions where the result of n experiment is given in terms of two (or more) quntities. We then need

More information

SOME USEFUL MATHEMATICS

SOME USEFUL MATHEMATICS SOME USEFU MAHEMAICS SOME USEFU MAHEMAICS I is esy o mesure n preic he behvior of n elecricl circui h conins only c volges n currens. However, mos useful elecricl signls h crry informion vry wih ime. Since

More information

Chapter 2 Linear Mo on

Chapter 2 Linear Mo on Chper Lner M n .1 Aerge Velcy The erge elcy prcle s dened s The erge elcy depends nly n he nl nd he nl psns he prcle. Ths mens h prcle srs rm pn nd reurn bck he sme pn, s dsplcemen, nd s s erge elcy s

More information

Comparison Procedures

Comparison Procedures Comprison Procedures Single Fctor, Between-Subects Cse /8/ Comprison Procedures, One-Fctor ANOVA, Between Subects Two Comprison Strtegies post hoc (fter-the-fct) pproch You re interested in discovering

More information

An integral having either an infinite limit of integration or an unbounded integrand is called improper. Here are two examples.

An integral having either an infinite limit of integration or an unbounded integrand is called improper. Here are two examples. Improper Inegrls To his poin we hve only considered inegrls f(x) wih he is of inegrion nd b finie nd he inegrnd f(x) bounded (nd in fc coninuous excep possibly for finiely mny jump disconinuiies) An inegrl

More information

The solution is often represented as a vector: 2xI + 4X2 + 2X3 + 4X4 + 2X5 = 4 2xI + 4X2 + 3X3 + 3X4 + 3X5 = 4. 3xI + 6X2 + 6X3 + 3X4 + 6X5 = 6.

The solution is often represented as a vector: 2xI + 4X2 + 2X3 + 4X4 + 2X5 = 4 2xI + 4X2 + 3X3 + 3X4 + 3X5 = 4. 3xI + 6X2 + 6X3 + 3X4 + 6X5 = 6. [~ o o :- o o ill] i 1. Mrices, Vecors, nd Guss-Jordn Eliminion 1 x y = = - z= The soluion is ofen represened s vecor: n his exmple, he process of eliminion works very smoohly. We cn elimine ll enries

More information

Section 6: Area, Volume, and Average Value

Section 6: Area, Volume, and Average Value Chpter The Integrl Applied Clculus Section 6: Are, Volume, nd Averge Vlue Are We hve lredy used integrls to find the re etween the grph of function nd the horizontl xis. Integrls cn lso e used to find

More information

Continuous Random Variables

Continuous Random Variables CPSC 53 Systems Modeling nd Simultion Continuous Rndom Vriles Dr. Anirn Mhnti Deprtment of Computer Science University of Clgry mhnti@cpsc.uclgry.c Definitions A rndom vrile is sid to e continuous if there

More information

0 for t < 0 1 for t > 0

0 for t < 0 1 for t > 0 8.0 Sep nd del funcions Auhor: Jeremy Orloff The uni Sep Funcion We define he uni sep funcion by u() = 0 for < 0 for > 0 I is clled he uni sep funcion becuse i kes uni sep = 0. I is someimes clled he Heviside

More information

1.B Appendix to Chapter 1

1.B Appendix to Chapter 1 Secon.B.B Append o Chper.B. The Ordnr Clcl Here re led ome mporn concep rom he ordnr clcl. The Dervve Conder ncon o one ndependen vrble. The dervve o dened b d d lm lm.b. where he ncremen n de o n ncremen

More information

Review of Gaussian Quadrature method

Review of Gaussian Quadrature method Review of Gussin Qudrture method Nsser M. Asi Spring 006 compiled on Sundy Decemer 1, 017 t 09:1 PM 1 The prolem To find numericl vlue for the integrl of rel vlued function of rel vrile over specific rnge

More information

Notes on the stability of dynamic systems and the use of Eigen Values.

Notes on the stability of dynamic systems and the use of Eigen Values. Noes on he sabl of dnamc ssems and he use of Egen Values. Source: Macro II course noes, Dr. Davd Bessler s Tme Seres course noes, zarads (999) Ineremporal Macroeconomcs chaper 4 & Techncal ppend, and Hamlon

More information

1 jordan.mcd Eigenvalue-eigenvector approach to solving first order ODEs. -- Jordan normal (canonical) form. Instructor: Nam Sun Wang

1 jordan.mcd Eigenvalue-eigenvector approach to solving first order ODEs. -- Jordan normal (canonical) form. Instructor: Nam Sun Wang jordnmcd Eigenvlue-eigenvecor pproch o solving firs order ODEs -- ordn norml (cnonicl) form Insrucor: Nm Sun Wng Consider he following se of coupled firs order ODEs d d x x 5 x x d d x d d x x x 5 x x

More information

QUADRATIC EQUATIONS OBJECTIVE PROBLEMS

QUADRATIC EQUATIONS OBJECTIVE PROBLEMS QUADRATIC EQUATIONS OBJECTIVE PROBLEMS +. The solution of the eqution will e (), () 0,, 5, 5. The roots of the given eqution ( p q) ( q r) ( r p) 0 + + re p q r p (), r p p q, q r p q (), (d), q r p q.

More information

Trigonometry. Trigonometry. Solutions. Curriculum Ready ACMMG: 223, 224, 245.

Trigonometry. Trigonometry. Solutions. Curriculum Ready ACMMG: 223, 224, 245. Trgonometry Trgonometry Solutons Currulum Redy CMMG:, 4, 4 www.mthlets.om Trgonometry Solutons Bss Pge questons. Identfy f the followng trngles re rght ngled or not. Trngles,, d, e re rght ngled ndted

More information

J i-1 i. J i i+1. Numerical integration of the diffusion equation (I) Finite difference method. Spatial Discretization. Internal nodes.

J i-1 i. J i i+1. Numerical integration of the diffusion equation (I) Finite difference method. Spatial Discretization. Internal nodes. umercal negraon of he dffuson equaon (I) Fne dfference mehod. Spaal screaon. Inernal nodes. R L V For hermal conducon le s dscree he spaal doman no small fne spans, =,,: Balance of parcles for an nernal

More information

HEAT CONDUCTION PROBLEM IN A TWO-LAYERED HOLLOW CYLINDER BY USING THE GREEN S FUNCTION METHOD

HEAT CONDUCTION PROBLEM IN A TWO-LAYERED HOLLOW CYLINDER BY USING THE GREEN S FUNCTION METHOD Journal of Appled Mahemacs and Compuaonal Mechancs 3, (), 45-5 HEAT CONDUCTION PROBLEM IN A TWO-LAYERED HOLLOW CYLINDER BY USING THE GREEN S FUNCTION METHOD Sansław Kukla, Urszula Sedlecka Insue of Mahemacs,

More information

S Radio transmission and network access Exercise 1-2

S Radio transmission and network access Exercise 1-2 S-7.330 Rdio rnsmission nd nework ccess Exercise 1 - P1 In four-symbol digil sysem wih eqully probble symbols he pulses in he figure re used in rnsmission over AWGN-chnnel. s () s () s () s () 1 3 4 )

More information

Modeling and Predicting Sequences: HMM and (may be) CRF. Amr Ahmed Feb 25

Modeling and Predicting Sequences: HMM and (may be) CRF. Amr Ahmed Feb 25 Modelg d redcg Sequeces: HMM d m be CRF Amr Ahmed 070 Feb 25 Bg cure redcg Sgle Lbel Ipu : A se of feures: - Bg of words docume - Oupu : Clss lbel - Topc of he docume - redcg Sequece of Lbels Noo Noe:

More information

Research Article Oscillatory Criteria for Higher Order Functional Differential Equations with Damping

Research Article Oscillatory Criteria for Higher Order Functional Differential Equations with Damping Journl of Funcon Spces nd Applcons Volume 2013, Arcle ID 968356, 5 pges hp://dx.do.org/10.1155/2013/968356 Reserch Arcle Oscllory Crer for Hgher Order Funconl Dfferenl Equons wh Dmpng Pegung Wng 1 nd H

More information

1 PYTHAGORAS THEOREM 1. Given a right angled triangle, the square of the hypotenuse is equal to the sum of the squares of the other two sides.

1 PYTHAGORAS THEOREM 1. Given a right angled triangle, the square of the hypotenuse is equal to the sum of the squares of the other two sides. 1 PYTHAGORAS THEOREM 1 1 Pythgors Theorem In this setion we will present geometri proof of the fmous theorem of Pythgors. Given right ngled tringle, the squre of the hypotenuse is equl to the sum of the

More information

Matrix Eigenvalues and Eigenvectors September 13, 2017

Matrix Eigenvalues and Eigenvectors September 13, 2017 Mtri Eigenvlues nd Eigenvectors September, 7 Mtri Eigenvlues nd Eigenvectors Lrry Cretto Mechnicl Engineering 5A Seminr in Engineering Anlysis September, 7 Outline Review lst lecture Definition of eigenvlues

More information