ANOVA models for Brownian motion
|
|
- Annabelle Bernice Atkins
- 6 years ago
- Views:
Transcription
1 ANOVA models for Brownan moon Gordon Hazen Danel Apley Neehar Parkh Deparmen of Indusral Engneerng and Managemen Scences Norhwesern Unversy Transplan Oucomes Research Collaborave Norhwesern Unversy December 013 Correspondng Auhor: Gordon Hazen IEMS Deparmen, McCormck School of Engneerng Norhwesern Unversy Evanson IL Runnng Head: ANOVA for Brownan Moon 1
2 Absrac We nvesgae longudnal models havng Brownan-moon covarance srucure. We show ha any such model can be vewed as arsng from a relaed meless classcal lnear model where sample szes correspond o longudnal observaon mes. Ths relaonshp s of praccal mpac when here are closed-form ANOVA ables for he relaed classcal model. Such ables can be drecly ransformed no he analogous ables for he orgnal longudnal model. We n parcular provde complee resuls for one-way fxed and random effecs ANOVA on he drf parameer n Brownan moon, and llusrae s use n esmang heerogeney n umor growh raes. KEY WORDS: LONGITUDINAL MODELS; BROWNIAN MOTION; ANALYSIS OF VARIANCE; ANOVA; RANDOM EFFECTS; STOCHASTIC HETEROGENEOUS TUMOR GROWTH. 1. Inroducon The sascal analyss of longudnal daa has a long hsory (e.g., (Dggle, 1988; Dggle, Heagery, Lang, & Zeger, 00; Snger & Wlle, 003; Verbeke & Molenberghs, 009)). Consder he suaon n whch repeaed measuremens x 0, x1,, x on ndvdual are aken a mes J lnear model wh parameers,,, 1 : j M M j j M J. We consder here a parcular x x ( u u ) j ~ normal(0, j ). (1) where u,, 1 u are specfed consans or observed covaraes for ndvdual. We allow he ndex o be M lexcographcally nesed, so ha (1) ncludes, among many possble examples, he addvely separable model wh = (g,h): x x 0 ( ) ghj gh g h ghj ghj n whch ( 1,, M ) = (, 1,, G, 1,, H). In (1) we resrc ourselves o he specfc correlaon srucure Corr( j, j) = j j f j j () ha arses when observaons come from a Brownan moon process (e.g. (Karln & Taylor, 1975; Oksendal, 1998)) wh drf u 11 umm and volaly.
3 Random effecs versons of model (1) are mporan when measuremens are avalable from mulple ndvduals (e.g., cancer growh daa from mulple paens see secon 5). In hs case, one mgh wonder wheher he populaon of ndvduals s homogeneous all sharng he same growh rae coeffcens,..., 1 M or heerogeneous. In he laer case, one mgh suspec ha ndvdual growh rae coeffcens are sampled from a hypohecal populaon (e.g., slow- versus fas-growng umors). Such models are aracve because hey allow an analys o forecas adapvely based on an ndvdual s hsory. For nsance, n a conrbuon no nvolvng repeaed measuremens, Ayer e al (Ayer, Alagoz, & Sou, 01) examne personalzng proocols for breas cancer screenng, n whch healhy screenng oucomes for an ndvdual yeld an adapve forecas of reduced cancer ncdence. Ths allows avodance of unnecessary screenng, whch would reduce coss and false-posve morbdes whou mpacng survval. In hs paper, we show ha any model of he form (1)-() can be vewed as arsng from a relaed meless classcal lnear model n whch sample szes correspond o observaon me ncremens j =,j+1 j. Moreover, nferences from he relaed classcal model apply drecly o model (1)-(), afer modfyng he degrees of freedom n he error sum of squares. When he classcal lnear model possesses a closed-form sum-of-squares able, hs able convers drecly no a sum-of-squares able for he desred model (1)-(). In parcular, we show ha he specal case of one way Brownan ANOVA, n whch he parameers are he paen drfs hemselves, arses from a classcal one-way ANOVA ha s smple enough o solve by hand or on a spreadshee. Recasng he analyss n a famlar ANOVA forma can provde greaer nsgh han he black box resuls from sascal sofware. Any nference procedure for Brownan moon would also apply o processes such as geomerc Brownan moon obaned by ransformaon (see secon 5), snce nference could occur on ransformed daa. We beleve he praccal mpor of hs paper wll be for analyss examnng longudnal daa, bu wh lmed access o or experence wh sascal sofware or general lnear models, who wsh o quckly oban ANOVA resuls n a famlar forma whou acqurng sofware or sudyng general lnear models. We srucure he paper as follows. In Secon we ransform o a general lnear model wh consan error varances, and derve closed-form esmaes and ess. In Secon 3 we specalze hese general resuls o one-way ANOVA on a Brownan drf parameer, and also consder esmaon and esng for a random-effecs versons of hs model. In Secon 4, we examne he connecon descrbed above beween model (1)-() and relaed meless classcal lnear models. Secon 5 provdes a one-way Brownan ANOVA on umor growh daa. Secon 6 concludes. 3
4 . ANOVA for he general longudnal model In model (1)-(), le x j = x j x,j1 j = j,j1 j = 1,,J. Then we have x j = j + j j ~ ndependen normal(0, j ) (3) where u 11 umm and here are a oal of J = observaons. We rewre he equaon for he drf vecor = ( 1,..., I ) n erms of he parameer vecor = ( 1,, M ): J = U (4) Noaon We consder boh fxed effecs and random effecs models. We use he noaon x x x x x x j j J 1 j j j j J 1 j j The maeral below uses he followng marx-vecor noaon. The prme symbol denoes marx ranspose. I J wll denoe an deny marx of dmenson J. If A 1,,A I are marces, hen dag (A ) s he block dagonal marx wh blocks A. If he A have he same number of columns, hen sack (A ) s he marx A sackng A 1 ono A ono ono A I. The followng deny wll be useful: dag (A ) sack (B ) = sack (A B ) f A,B are compable under marx mulplcaon. The basc model 1 I A obaned by We rea boh fxed and random effecs, bu nally do no dsngush whch case we are consderng. Wh he ransformaon 1/ y j = x j / j he basc model (3) becomes 1/ y j = j + e j e j ~ ndependen normal(0, ) In vecor noaon, hs s 1/ y = dag ( ) + e e ~ normal(0, I J ) (5) where 4
5 y sack ( y ) y sack ( y ) j j 1/ 1/ sack j( j ) = sack ( ). The resulng model s: 1/ y = dag ( )U + e. (6) Ths s a specal case of he lnear model y = X + e n whch 1/ X = dag ( )U. (7) In he followng, we explore wha happens o ANOVA for hs lnear model. Error sum of squares The leas-squares esmaes ˆα are he unresrced values of ha mnmze he sum of squares S() subjec o = U, where 1/ 1/ 1/ S( α) yxα ydag ( ) Uα ydag ( ) μ sack y x 1/ j jyj j j j j (8) Le ˆμ =Uα ˆ The resulng error sum of squares s gven by SS E = x j yxα ˆ ˆ j j (9) j SS E has J r degrees of freedom, where r s he rank of X = dag ( because f u s he h row of U, hen X = dag ( )U = sack ( elemen 1/ 1/ 1/ )U. The laer s equal o he rank of U u ). Assumng each u s nonzero, hen each 1/ u of hs sack X s an array of rank 1 wh rows proporonal o u. If r s he rank of U, hen here are r lnearly ndependen rows u, and X mus have rank r = r. Numeraor sum of squares Suppose he null hypohess o es s H: (agans an alernave hypoheses wh no resrcons on ) and consss of q ndependen lnear resrcons on. Le α ˆ be he leas squares esmae under H, and le μ ˆ = The numeraor sum of squares for an F-es of H s SS 1/ 1/ 1/ ˆ ˆ ˆ ˆ ˆ ˆ dag ( ˆ ˆ ) dag ( ˆ ˆ Xα Xα U αα μμ) sack ( ). SS has q degrees of freedom.. (10) Uα ˆ. 5
6 For he fxed-effecs scenaro, he expecaon of he mean square MS = SS/q can be obaned from Scheffe s Rule ((Scheffe, 1959), p.39), and s gven by E[MS] = 1 + q ( ) (11) where s obaned from he formula for ˆ by replacng each erm y j by s mean replacng each erm x j by j. F Tes 1/ j. Ths s equvalen o SS E and SS are ndependen, wh SS E / havng a ch-square dsrbuon wh he approprae degrees of freedom, and SS/ beng noncenral ch-square n general, and ch-square under. So he usual F-es of he hypohess H apples. The noncenraly parameer for SS/ may be obaned from Scheffe s Rule 1 ((Scheffe, 1959), p.39): If n SS, each observaon y j s replaced by s expeced value (or x j by j ), he resul s. 1/ j Random effecs In a random effecs model, we would have ~ normal( 0,). (If only some erms of are random, hs would be refleced by zeroes n.) The dsrbuonal and mean clams above reman condonally rue gven. In parcular, gven, he quany SS E / s condonally (Jr) a dsrbuon ha does no depend on. Therefore, SS E / s uncondonally (Jr). The random-effecs counerpar o he fxed-effecs null hypohess H: s a null hypoheses H ha places ceran resrcons on. Namely, H resrcs he varances of he q ndependen lnear conrass of correspondng o o be zero. Then under H, SS E / and SS/ are condonally ndependen gven wh (Jr) and (q) dsrbuons no dependng on, so are herefore ndependen (Jr) and (q) varables under H. The F sasc herefore sll has an F-dsrbuon under H, and he same F-es apples. 3. Specal Case of One way Brownan ANOVA As we have noed, our framework encompasses one-way or mul-way ANOVA models, nesed or crossed, wh or whou covaraes, ec. In hs secon, we specalze o he smples verson of Brownan ANOVA, a one-way model, whch we defne as he general model (1)-() wh he I drfs,, 1 I he parameers of neres. In he framework of he general model, for he one-way ANOVA specal case we have = and U = I I n (4), yeldng X 6
7 = dag ( 1/ ). The general resuls from he pror secon lead drecly o many of he asserons summarzed n Table 1 below. The fxed-effecs asserons n parcular fall ou n a sraghforward manner when he leas-squares esmaes x ˆ ˆ x are subsued. Hence, n he nex subsecon, we focus aenon on dervng he random-effecs asserons n Table 1. Table 1. One-way fxed and random effecs Brownan ANOVA Fxed effecs Random effecs Model,, 1 I 1 parameers I ndependen normal(, ) Hypohess H: 1 I H: = 0 Error sum of xj squares x SS E = j j wh J I degrees of freedom j SSE MS E = J I Hypohess x x sum of SS = squares wh I1 degrees of freedom SS MS = I 1 Expeced mean squares E[MS E ] = E[MS] = 1 + ( I 1) ( ) E[MS E ] = E[MS]= +. F-es Parameer esmaes 1 Rejec H a sgnfcance level f ˆ ˆ = MSE x ˆ seˆ ˆ 1/ MS MS E > F, I 1, J I ˆ 1 1. ( I 1) MSE 1 ˆ ( MS MS E ) x ˆ ˆ ˆ ˆ ˆ seˆ ˆ In case hs quany s negave, s common pracce o replace he esmae wh zero. ˆ 1/ 7
8 A key ake-away s ha Table 1 s nearly dencal o he analogous able for a meless classcal one-way ANOVA f one subsues observaon me ncremens j n he Brownan model for whn-cell sample szes n he classcal model. Ths neresng connecon connues o hold n he general model (1)-(). We elaborae n secon 4. Dervaon of Random effecs Resuls n Table 1 In he random effecs verson, we assume he effecs are dsrbued as { : = 1,,..., I} ~ ndependen normal(, ). Regardng he fxed-effecs expecaon of he mean square MS, we have n equaon (11) from he pror secon E x 1 1 1, j j whch we denoe by he quany n Table 1, as does no depend on. In he random effecs model, we have from (11) ha E[MS] = 1 +. ( I 1) E ( ) Usng Lemma A, hs becomes E[MS]= +, where ( I 1) 1 1, as shown n he able. Here, n parallel o classcal random effecs ANOVA, when he oal observaon mes are dencal across paens, he quany wll be her common value. Combnng he las equaon wh he fac ha MS E s an unbased esmae of, we oban n he convenonal ANOVA momen esmaor for : 1 ˆ ( MS MS E ). As n classcal random effecs ANOVA, s possble ha hs esmae can be negave, n whch case s common pracce o repor a zero esmae. A leas-squares esmae of he random effecs populaon mean can be obaned as follows when we rea and as known. Subsung ~ normal(1 I, I I ) no he random effecs model y = dag ( 1/ )+ e gves y = X1 I + g g ~ normal(0, IJ+ XX) 8
9 Defnng W = ( I J + XX) 1/ X1 I, he weghed leas-squares esmae of, reang and as known, s Wz 1X( I XX) ( I XX) y 1X( I XX) y ˆ WW 1X I XX I XX X1 1X I XX X1 1/ 1/ 1 I J J I J 1/ 1/ 1 I ( J ) ( J ) I I ( J ) I The marx o be nvered here smplfes o ( I XX) dag ( I ) 1 1/ 1/ 1 J J Use Lemma A3 o wre ( ) ( ) ( ) 1/ 1/ 1 1 1/ 1/ 1 IJ IJ I 1/ 1 1/ J IJ 1 ( IJ ) 4 1/ 1/ IJ 1 1/ 1/ IJ Subsue hs no he denomnaor of he las equaon for ˆ o ge 1 and no he numeraor o ge 1/ 1/ 1/ 1 1/ 1/ 1/ 1/ 1/ ( J ) I IJ. 9
10 10 1/ 1/ 1/ 1 1/ 1/ 1/ ( ) 1 J J y y x x x x x I I The resulng quoen s he esmae ˆ lsed n Table 1. To oban an esmaed sandard error for ˆ, calculae he numeraor varance usng ndependence of he x : 1 Var Var 1 Var E[ ] E Var[ ] 1 Var E 1. x x x x Then subsue o oban ˆ se lsed n he able. 4. Equvalence o classcal ANOVA models We noed n he pror secon he near equvalence beween one-way classcal and Brownan ANOVAs when observaon me s reaed as sample sze. Ths equvalence n fac exends o any longudnal model (1)-(), ha s, o any parameerzaon = U. Such parameerzaons nclude regresson and ANCOVA models. In hs secon, we explan how hs connecon arses. We resrc ourselves here o fxed-effecs models. For comparson, we rewre our longudnal model (3) here: x j = j + j j ~ NID(0, j ) = 1,,I; j = 1,,J, ()
11 Consder he relaed classcal ANOVA model y jk = + jk jk ~ NID(0, ) = 1,,I; j = 1,,J ; k = 1,,K j (*) havng he same parameers = ( 1,..., I ) and subjec o he same srucure = U. Le ˆα, α ˆ be he leassquares esmaes 1 of under models () and (*), respecvely; le α ˆ, α ˆ be he leas squares esmaes for hypohess H: under () and (*), respecvely; and le E SS, SS, SS, * SS E be he correspondng sums of squares. Le K K be he oal number of observaons. j j In he followng s assumed we can se K j = j for all,j. Of course hs requres he j o be negers, a condon ha can be acheved by suably rescalng me n model (). Proposon 1: The leas-squares esmaes α ˆ and α ˆ under model (*) depend on he varables y jk only hrough he sample sums y. Whenever yj j x and K j = j for all,j, he followng relaonshps hold: j ˆ α = α ˆ = * SS = ˆα ˆ α SS j k SS y y SS * E jk j E where y K y 1. j j j Proof: We begn by nroducng vecor noaon for (*): y = sack j (sack k (y = sack (y ) = sack (sack j (sack k ( ))) + = sack (sack j ( 1 K j ) + 1 We assume sde condons are mposed o nsure unqueness of he leas-squares esmaes, he same sde condons n (*) as n (). For model (), he esmaes are he ones dscussed n Secon, ha s, he weghed leas-squares esmaes obaned afer ransformng he daa o have homogeneous varances. Ths echncally may be no possble f some of he j are rraonal, bu n pracce daa would be expressed o only a fne number of decmal places, hence are effecvely raonal. 11
12 = sack ( 1K ) + = dag ( 1 )+. K Wh he parameerzaon = U, we have y = dag ( 1 )U+ = X * + K where X * = dag ( 1 )U. Under model (*), he leas-square esmaes α ˆ mnmze he sum of squares K subjec o = U, where: S * ( α) yx α sack ( y) dag ( 1K ) μ sack y 1 K yjk j k yjk yj j j k y j k yjk yj Kj yj S j k j S ( α) 1 So n effec, wha mus be mnmzed s subjec o = U. Of course he resulng α ˆ depends only on he sample sums va he sample means yj. Moreover, he quany s dencal o he sum of squares of (8) (we have added he superscrp o ndcae s for he model ()) under he assumed condons yj xj and K j = j for all,j. Therefore he resulng α ˆ mus be dencal o condons. The same argumen mples ha α ˆ = α ˆ. ˆα under hese Tha follows by nong ha, whch (by he argumens n he precedng paragraph) s he same as under he assumed condons yj xj and K j = j for all,j. Fnally, subsung α α, we oban, s ˆ. SS E under he assumed condons..,, (1) Noce n Proposon 1 ha f n addon y jk does no depend on k for all,j, hen we have no only,, ) = (,, ), bu also SS SS. Under he assumpon K j = j of he proposon, we have * E E 1
13 y x { j j, jk y does no depend on k} yjk xj j for all,j Proposon 1 herefore mples ha by performng a classcal ANOVA on model (*) whle makng he subsuons K j = j, yjk xj j, we can oban all relevan sascs for he longudnal model (). Bu does he (*)- ANOVA gve us he correc sascal nferences for model ()? I s apparen canno, because (Kr) dsrbuon, whereas whch * SS E / has a SS E / has a (Jr) dsrbuon. To see wha has happened, refer o (1) above, n * SS E s decomposed no he sum of wo ndependen quanes, and, wh K J and Jr respecve degrees of freedom. Assumng y jk does no depend on k, as n he subsuon above, resuls n, = 0 and =, ~ (Jr). In general, sascal nferences for he (*)-ANOVA depend on he jon dsrbuon of he quanes (,,,, jon because conclusons depend on probablsc ndependences ha may exs, such as ndependence of and needed for hypohess esng, and ndependence of and needed for confdence nerval consrucon. One mgh worry ha he resrcon ha y jk does no depend on k could do more han merely change he dsrbuon of / from (Kr) o (Jr) perhaps could furher aler he jon dsrbuon of (,,, n some way. Tha hs s no so s a consequence of he followng resul. Proposon : The jon dsrbuon of (,,,, s ndependen of he even A = {y jk does no depend on k for all,j}. Moreover, f K j = j for all,j, hen hs dsrbuon s dencal o he jon dsrbuon of (,,,. Proof: To derve he frs clam, observe he followng: () (,,,, ) depends on only va he whnsample sums. or, equvalenly, va he whn-sample averages.. ; (), s ndependen of. and, herefore, of (,,,, ), because, s a funcon of he whn-sample varances., whch are ndependen of he whn-sample averages for he normal samples n queson; and () he saed 13
14 even A s equvalen o, 0. Combnng hese hree resuls proves ha he jon dsrbuon of (,,,, s ndependen of he even A, as assered. To show he second clam noe frs ha by Proposon 1, when K j = j for all,j, we have (,,, ) = g({x j }) and (,,,, ) = g({ y j }) for some funcon g() he same funcon n boh cases. The second clam hen follows from he easly-verfed fac ha he jon dsrbuon of { y j } under model (*) s he same as he jon dsrbuon of {x j } under model (). Wha hen s he relevance o he longudnal model () of sascal nferences from classcal (*)-ANOVA when we make he subsuons K j = j, yjk xj j? Sep back for a momen and compare wo nferenal exercses: Sascal nference for parameers under model (*) and he condon A = {y jk does no depend on k for all,j}. The daa on whch such nferences are based may be aken, by Proposon 1, o be he sample sums { yj } and sample szes {K j }. Any such nferences depend by Proposon on he same jon dsrbuon F of (,,,, ha obans n he absence of condon A, ha s, he jon dsrbuon ha obans n he usual classcal ANOVA. The only dfference s ha under condon A, we have, Kr o Jr. 0, forcng he error sum of squares o equal,, and reducng s degrees of freedom from Sascal nference for parameers under model (). Ths s based on daa {x j } and observaon mes { j }, and depends on he same (by Proposon ) jon dsrbuon F of (,,,. Because of hese parallels, nference n hese wo suaons wll be dencal on dencal daa j = K j, xj yj. Bu as we have menoned, hese equales are equvalen o j = K j, yjk xj j. Ths leads us o he followng. Concluson: Sascal nferences on obaned by subsung yjk xj j and K j = j for all,j no he classcal (*)-ANOVA and changng he degrees of freedom for from Kr o Jr wll be vald for n he longudnal model (). 14
15 Alhough hs concluson s far-reachng, s praccal mplcaons are more lmed. If exreme scalng s needed o ge all-neger j, he resulng large sample szes K j and quanes of dencal daa y jk for he (*)-ANOVA wll be nconvenen for sofware packages, and he SS E degrees of freedom n sofware oupu would need o be correced from Kr o Jr. So hs s no n general an approach we recommend for sofware packages. However, none of hese dffcules arse f here are closed-form ables for ANOVA on model (*). Closed-form ables for unbalanced models lke (*) do exs n some cases, for example, for one-way ANOVA as we have already dscussed, as well as for nesed herarchcal models when effecs are weghed by sample sze (see (Scheffe, 1959), Ch. 5.3). We dscuss he laer n he example presened below. In general, o conver a closed-form classcal ANOVA able for model (*) o a correspondng ANOVA able for model (): 1. Replace he sample szes K j by observaon mes j everywhere.. In * SS E and * SS, replace y jk and yj by x j / j. 3. Change he degrees of freedom Kr for * SS E o Jr. I should be noed ha he replacemens n Sep also enal he followng replacemens: x y K K y x x 1 1 j 1 j j j j j j x y K K y x x 1 1 j 1 j j j j j j Noe n sep 1 ha one can om any scalng of he j o neger values. Ths s because any such scalng would cancel from numeraor and denomnaor of all F-sascs. Moreover, any scaled esmae of a scaled, hence he unscaled Example: A nesed model * MS E would be unbased for he unscaled. * MS E would be an unbased The followng llusraes how a specfc unbalanced nesed classcal ANOVA ranslaes o a nesed longudnal ANOVA usng he mehods jus dscussed. Scheffe ((Scheffe, 1959),Ch. 5.3) nroduces a procedure for consrucng classcal ANOVA ables for unbalanced nesed desgns. Hs procedure apples, for example, o he followng nesed desgn. 15
16 yhjk h hjk ~ NID(0, ) (13) hjk h h h h = 1,,H; = 1,,I h ; j = 1,,J h ; k = 1,,K hj. I h h I ; J h h J ; hj. To be concree, hs model reflecs, for example, he effec h j K K due o paens. The ndces par (h,) n (13) corresponds o he sngle ndex n (*). To elmnae parameer redundancy, Scheffe mposes he followng resrcons 3 : K h h 0 h and K h h 0 for each h. The leas-square esmaes and assocaed 1 confdence nervals are: h due o hospals h and whn hospal h he effec h ˆh y h y,/ K K 1 1 h (14) ˆh y h y h,/ K K (15) 1 1 h h where. For esng he wo dfferen null hypoheses H : h 0 for all h H : h 0 for all h, he ANOVA able s: Source SS d.f. SS Kh( yh y ) H 1 h h ( h h ) h I H SS K y y Error SS E = y h j k hjk yh K I Toal SS To = y hjk y K 1 h j k 3 I s o be noed ha hese specfc resrcons are needed n order ha he sums of squares n he ANOVA able below decompose properly. 16
17 The correspondng longudnal ANOVA resuls would be obaned by followng he nsrucons 1,,3 saed earler n hs secon. The ndex here corresponds o he par (h,) here, and he longudnal ANOVA model s he same as model (13) bu wh yhjk h hjk and NID(0, hj hjk ~ NID(0, ) replaced by xhj hhj hj and hj ~ ), respecvely. The leas-square esmaes and confdence nervals (14)-(15)ranslae o ˆh x h x,/ K K 1 1 h ˆh x h x h,/ K K 1 1 h h where and we have used JI nsead of KI as he degrees of freedom for. The classcal nesed ANOVA able above ranslaes o: Source SS d.f. SS h( xh x ) H 1 Error SS E = h h ( h h ) h I H xhj h j hj xh hj J I SS x x xhj h j hj x J 1 hj Toal SS To = These resuls arse by nvokng he gudelnes above o oban he subsuons K I J I Khj hj yhjk xhj hj yh xh Kh h Kh h h h y x y x under whch he nner sum over k n SS E smplfes snce s argumen no longer depends on k. 5. Example Analyss Ebara e al. (Ebara e al., 1986) publshed daa (see Fgure 1) on he growh of unreaed hepaocellular carcnoma umors. To llusrae he mehods presened here, we perform a one-way random effecs Brownan ANOVA (Table 1) o analyze growh-rae heerogeney n hs daa. 17
18 Fgure 1. Hepaocellular carcnoma umor growh n paens from Ebara e al. (1986). Sophscaed dfferenal equaon models of umor growh are avalable (see (Araujo & Mcelwan, 004)). Here for smplcy we focus on an exponenal growh model, n whch, f z s sze a me, and s growh rae a me, hen dz d z whch yelds exponenal growh z z0e when s a consan. A sochasc verson of hs model allows growh rae o vary n random ndependen ncremens over me: d d db where s mean growh rae, s growh rae volaly, and db s an ncremen of sandard Brownan moon. Ths assumpon forces z o be geomerc Brownan moon (e.g.,(oksendal, 1998)), ha s, x = ln(z ) s Brownan moon wh drf and volaly : dz z d(ln z ) d db. Cha e al. (Cha, P.Salzman, Plevrs, & Glynn, 004) employ hs model o explore smulaon-based parameer esmaon for breas cancer umor growh. If we allow mean growh rae o depend on paen, hen our model (3) s an mmedae resul. A random-effecs model allows growh raes for paens o be randomly sampled from a normal populaon wh mean and sandard devaon. Our one-way ANOVA procedures (Table 1) allow us o esmae and and 18
19 es he hypohess H: = 0. We performed such a one-way random-effecs ANOVA based on Table 1. The procedure s smple enough ha we carred ou on a spreadshee. A word abou uns: Because x = ln(z ), we have dx dz/ z, so ha dx (and herefore x self) s a unless quany. I follows ha he quanes,,, have uns of me 1, as do all sums of squares and mean squares. To make such quanes more meanngful n our resuls below, we conver hem o percens per 100 days. The resuls from applyng he ANOVA of Table 1 are: ˆ =1.4% /100 da. ˆ = 9.64% /100 da. ˆ = 13.% /100 da. se ˆ =.53% /100 da. Therefore, an esmaed 95% confdence nerval for he mean of he populaon of mean growh raes s 13.% ± z % = 13. ± 4.96 per 100 days, and he sandard devaon of he populaon of mean growh raes s esmaed o be 9.64% per 100 days. The resuls of applyng he Table 1 ANOVA es of he null hypohess H: = 0 are Source SS df MS F p 118.7% % E 06 Error 16.6% % from whch we conclude ha he above esmae = 9.64% /100 days s sgnfcanly dfferen from 0 (p-value ). We conclude ha growh-rae heerogeney s presen n he populaon from whch hese observaons arose. 6. Concluson The prmary heorecal resul from hs paper s ha any longudnal model (1)-() can be vewed as arsng from a correspondng meless classcal ANOVA model n whch sample szes are observaon mes; and ha sascal nferences for he former can be drecly exraced from he laer once error degrees of freedom are suably revsed. The models (1)-() nclude no only ANOVA models of all ypes, bu also regresson and ANCOVA models, for example. We concede, however, ha he praccal mpac of hs resul s lmed because: () as noed n Secon 4, 19
20 he requred sofware npus for he correspondng classcal model could be cumbersome, and he sofware oupus would need o be revsed n lgh of he alered degrees of freedom; and () s possble o ransform any longudnal model (1)-() no a general lnear model, and use avalable sascal sofware o numercally derve esmaes and hypohess es resuls whou any necessy of explong he connecon o classcal ANOVA. However, our procedure does have praccal usefulness when he desred longudnal model arses from a classcal ANOVA possessng a closed-form sum-of-squares able. In ha case, he classcal able convers drecly no a sum-of-squares able for he desred longudnal model. The eases verson of hs s lkely he one-way random effecs ANOVA (Table 1) n whch s desred o esmae he populaon of drfs from whch he daa have arsen, and o es for heerogeney n hs populaon. These esmaes and ess are smple enough o be derved on a spreadshee. We beleve he praccal mpor of hs paper wll be for analyss wh lmed access o or experence wh sascal sofware or general lnear models, who wsh o quckly oban famlar ANOVA resuls whou acqurng sofware or sudyng general lnear models. References Araujo, R. P., & Mcelwan, D. L. S. (004). A Hsory of he Sudy of Sold Tumour Growh: The Conrbuon of Mahemacal Modellng. Bullen of Mahemacal Bology, 66, Ayer, T., Alagoz, O., & Sou, N. K. (01). OR Forum A POMDP approach o personalze mammography screenng decsons. Operaons Research, 60(5), Cha, Y. L., P.Salzman, Plevrs, S. K., & Glynn, P. W. (004). Smulaon-based parameer esmaon for complex models: a breas cancer naural hsory modellng llusraon. Sascal Mehods n Medcal Research, 13, Dggle, P. J. (1988). An approach o he analyss of repeaed measuremens. Bomercs, 44(4), Dggle, P. J., Heagery, P., Lang, K., & Zeger, S. L. (00). Analyss of Longudnal Daa (second edon). New York: Oxford Unversy Press. Ebara, M., M.Oho, Shnagawa, T., Sugura, N., Kmura, K., Masuan, S., e al. (1986). Naural Hsory of Mnue Hepaocellular Carcnoma Smaller Than Three Cenmeers Complcang Crrhoss: A Sudy n Paens. Gasroenerology, 90,
21 Graybll, F. A. (1983). Marces wh Applcaons n Sascs (Second Edon). Pacfc Grove, Calforna: Wadsworh & Brooks/Cole. Karln, S., & Taylor, H. M. (1975). A frs course n sochasc processes (d ed.). New York: Academc Press. Oksendal, B. (1998). Sochasc Dfferenal Equaons: An Inroducon wh Applcaons (Ffh edon ed.). Berln, Hedelberg, New York: Sprnger-Verlag. Scheffe, H. (1959). The Analyss of Varance. New York: John Wley and Sons. Snger, J. D., & Wlle, J. B. (003). Appled Longudnal Daa Analyss: Oxford Unversy Press. Verbeke, G., & Molenberghs, G. (009). Lnear Mxed Models for Longudnal Daa. New York: Sprnger Verlag. Appendx Suppose we have observaons y 1,,y J and weghs w 1,,w J. Defne he sum of squares SS = w j( y y j ) j 1 y w w j jyj where w wj. The followng resul s sandard and a consequence of elemenary algebra. j Lemma A1: SS = wy j j. j wy Lemma A: Suppose n addon ha he varables y are ndependen wh mean and varance. Then 1 E[SS] = w w w. Proof: We have E y j y has mean and varance w w, so E y w w. Then 1
22 E[SS] = E wy j j wy wj( ) w w w j j = w w w w w w 1 1 The followng s also a sandard resul (e.g., see Theorems and n (Graybll, 1983)). Lemma A3: Le he k k marx C be gven by CDab 1 where D s a nonsngular symmerc marx, a and b are each k 1 vecors, and s a scalar such ha 1aD b 0. Then C = D 1 1 ad b D ab D.
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( t) Outline of program: BGC1: Survival and event history analysis Oslo, March-May Recapitulation. The additive regression model
BGC1: Survval and even hsory analyss Oslo, March-May 212 Monday May 7h and Tuesday May 8h The addve regresson model Ørnulf Borgan Deparmen of Mahemacs Unversy of Oslo Oulne of program: Recapulaon Counng
More informationDepartment 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 informationIn 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 informationJohn Geweke a and Gianni Amisano b a Departments of Economics and Statistics, University of Iowa, USA b European Central Bank, Frankfurt, Germany
Herarchcal Markov Normal Mxure models wh Applcaons o Fnancal Asse Reurns Appendx: Proofs of Theorems and Condonal Poseror Dsrbuons John Geweke a and Gann Amsano b a Deparmens of Economcs and Sascs, Unversy
More informationTSS = 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 informationV.Abramov - FURTHER ANALYSIS OF CONFIDENCE INTERVALS FOR LARGE CLIENT/SERVER COMPUTER NETWORKS
R&RATA # Vol.) 8, March FURTHER AALYSIS OF COFIDECE ITERVALS FOR LARGE CLIET/SERVER COMPUTER ETWORKS Vyacheslav Abramov School of Mahemacal Scences, Monash Unversy, Buldng 8, Level 4, Clayon Campus, Wellngon
More informationSolution in semi infinite diffusion couples (error function analysis)
Soluon n sem nfne dffuson couples (error funcon analyss) Le us consder now he sem nfne dffuson couple of wo blocks wh concenraon of and I means ha, n a A- bnary sysem, s bondng beween wo blocks made of
More informationAppendix H: Rarefaction and extrapolation of Hill numbers for incidence data
Anne Chao Ncholas J Goell C seh lzabeh L ander K Ma Rober K Colwell and Aaron M llson 03 Rarefacon and erapolaon wh ll numbers: a framewor for samplng and esmaon n speces dversy sudes cology Monographs
More informationCS434a/541a: Pattern Recognition Prof. Olga Veksler. Lecture 4
CS434a/54a: Paern Recognon Prof. Olga Veksler Lecure 4 Oulne Normal Random Varable Properes Dscrmnan funcons Why Normal Random Varables? Analycally racable Works well when observaon comes form a corruped
More informationGENERATING CERTAIN QUINTIC IRREDUCIBLE POLYNOMIALS OVER FINITE FIELDS. Youngwoo Ahn and Kitae Kim
Korean J. Mah. 19 (2011), No. 3, pp. 263 272 GENERATING CERTAIN QUINTIC IRREDUCIBLE POLYNOMIALS OVER FINITE FIELDS Youngwoo Ahn and Kae Km Absrac. In he paper [1], an explc correspondence beween ceran
More informationOn One Analytic Method of. Constructing Program Controls
Appled Mahemacal Scences, Vol. 9, 05, no. 8, 409-407 HIKARI Ld, www.m-hkar.com hp://dx.do.org/0.988/ams.05.54349 On One Analyc Mehod of Consrucng Program Conrols A. N. Kvko, S. V. Chsyakov and Yu. E. Balyna
More informationShould Exact Index Numbers have Standard Errors? Theory and Application to Asian Growth
Should Exac Index umbers have Sandard Errors? Theory and Applcaon o Asan Growh Rober C. Feensra Marshall B. Rensdorf ovember 003 Proof of Proposon APPEDIX () Frs, we wll derve he convenonal Sao-Vara prce
More informationCS286.2 Lecture 14: Quantum de Finetti Theorems II
CS286.2 Lecure 14: Quanum de Fne Theorems II Scrbe: Mara Okounkova 1 Saemen of he heorem Recall he las saemen of he quanum de Fne heorem from he prevous lecure. Theorem 1 Quanum de Fne). Le ρ Dens C 2
More information. The geometric multiplicity is dim[ker( λi. number of linearly independent eigenvectors associated with this eigenvalue.
Lnear Algebra Lecure # Noes We connue wh he dscusson of egenvalues, egenvecors, and dagonalzably of marces We wan o know, n parcular wha condons wll assure ha a marx can be dagonalzed and wha he obsrucons
More informationJanuary Examinations 2012
Page of 5 EC79 January Examnaons No. of Pages: 5 No. of Quesons: 8 Subjec ECONOMICS (POSTGRADUATE) Tle of Paper EC79 QUANTITATIVE METHODS FOR BUSINESS AND FINANCE Tme Allowed Two Hours ( hours) Insrucons
More informationTime-interval analysis of β decay. V. Horvat and J. C. Hardy
Tme-nerval analyss of β decay V. Horva and J. C. Hardy Work on he even analyss of β decay [1] connued and resuled n he developmen of a novel mehod of bea-decay me-nerval analyss ha produces hghly accurae
More informationEcon107 Applied Econometrics Topic 5: Specification: Choosing Independent Variables (Studenmund, Chapter 6)
Econ7 Appled Economercs Topc 5: Specfcaon: Choosng Independen Varables (Sudenmund, Chaper 6 Specfcaon errors ha we wll deal wh: wrong ndependen varable; wrong funconal form. Ths lecure deals wh wrong ndependen
More informationLinear Response Theory: The connection between QFT and experiments
Phys540.nb 39 3 Lnear Response Theory: The connecon beween QFT and expermens 3.1. Basc conceps and deas Q: ow do we measure he conducvy of a meal? A: we frs nroduce a weak elecrc feld E, and hen measure
More informationComparison of Differences between Power Means 1
In. Journal of Mah. Analyss, Vol. 7, 203, no., 5-55 Comparson of Dfferences beween Power Means Chang-An Tan, Guanghua Sh and Fe Zuo College of Mahemacs and Informaon Scence Henan Normal Unversy, 453007,
More information. The geometric multiplicity is dim[ker( λi. A )], i.e. the number of linearly independent eigenvectors associated with this eigenvalue.
Mah E-b Lecure #0 Noes We connue wh he dscusson of egenvalues, egenvecors, and dagonalzably of marces We wan o know, n parcular wha condons wll assure ha a marx can be dagonalzed and wha he obsrucons are
More informationPanel Data Regression Models
Panel Daa Regresson Models Wha s Panel Daa? () Mulple dmensoned Dmensons, e.g., cross-secon and me node-o-node (c) Pongsa Pornchawseskul, Faculy of Economcs, Chulalongkorn Unversy (c) Pongsa Pornchawseskul,
More informationRobustness 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 informationA HIERARCHICAL KALMAN FILTER
A HIERARCHICAL KALMAN FILER Greg aylor aylor Fry Consulng Acuares Level 8, 3 Clarence Sree Sydney NSW Ausrala Professoral Assocae, Cenre for Acuaral Sudes Faculy of Economcs and Commerce Unversy of Melbourne
More informationUNIVERSITAT AUTÒNOMA DE BARCELONA MARCH 2017 EXAMINATION
INTERNATIONAL TRADE T. J. KEHOE UNIVERSITAT AUTÒNOMA DE BARCELONA MARCH 27 EXAMINATION Please answer wo of he hree quesons. You can consul class noes, workng papers, and arcles whle you are workng on he
More informationStandard Error of Technical Cost Incorporating Parameter Uncertainty
Sandard rror of echncal Cos Incorporang Parameer Uncerany Chrsopher Moron Insurance Ausrala Group Presened o he Acuares Insue General Insurance Semnar 3 ovember 0 Sydney hs paper has been prepared for
More information( ) () we define the interaction representation by the unitary transformation () = ()
Hgher Order Perurbaon Theory Mchael Fowler 3/7/6 The neracon Represenaon Recall ha n he frs par of hs course sequence, we dscussed he chrödnger and Hesenberg represenaons of quanum mechancs here n he chrödnger
More information[Link to MIT-Lab 6P.1 goes here.] After completing the lab, fill in the following blanks: Numerical. Simulation s Calculations
Chaper 6: Ordnary Leas Squares Esmaon Procedure he Properes Chaper 6 Oulne Cln s Assgnmen: Assess he Effec of Sudyng on Quz Scores Revew o Regresson Model o Ordnary Leas Squares () Esmaon Procedure o he
More informationNew M-Estimator Objective Function. in Simultaneous Equations Model. (A Comparative Study)
Inernaonal Mahemacal Forum, Vol. 8, 3, no., 7 - HIKARI Ld, www.m-hkar.com hp://dx.do.org/.988/mf.3.3488 New M-Esmaor Objecve Funcon n Smulaneous Equaons Model (A Comparave Sudy) Ahmed H. Youssef Professor
More information[ ] 2. [ ]3 + (Δx i + Δx i 1 ) / 2. Δx i-1 Δx i Δx i+1. TPG4160 Reservoir Simulation 2018 Lecture note 3. page 1 of 5
TPG460 Reservor Smulaon 08 page of 5 DISCRETIZATIO OF THE FOW EQUATIOS As we already have seen, fne dfference appromaons of he paral dervaves appearng n he flow equaons may be obaned from Taylor seres
More informationRELATIONSHIP BETWEEN VOLATILITY AND TRADING VOLUME: THE CASE OF HSI STOCK RETURNS DATA
RELATIONSHIP BETWEEN VOLATILITY AND TRADING VOLUME: THE CASE OF HSI STOCK RETURNS DATA Mchaela Chocholaá Unversy of Economcs Braslava, Slovaka Inroducon (1) one of he characersc feaures of sock reurns
More informationVolatility Interpolation
Volaly Inerpolaon Prelmnary Verson March 00 Jesper Andreasen and Bran Huge Danse Mares, Copenhagen wan.daddy@danseban.com brno@danseban.com Elecronc copy avalable a: hp://ssrn.com/absrac=69497 Inro Local
More informationP R = P 0. The system is shown on the next figure:
TPG460 Reservor Smulaon 08 page of INTRODUCTION TO RESERVOIR SIMULATION Analycal and numercal soluons of smple one-dmensonal, one-phase flow equaons As an nroducon o reservor smulaon, we wll revew he smples
More informationChapter 5. The linear fixed-effects estimators: matrix creation
haper 5 he lnear fed-effecs esmaors: mar creaon In hs chaper hree basc models and he daa marces needed o creae esmaors for hem are defned. he frs s ermed he cross-secon model: alhough ncorporaes some panel
More informationVariants of Pegasos. December 11, 2009
Inroducon Varans of Pegasos SooWoong Ryu bshboy@sanford.edu December, 009 Youngsoo Cho yc344@sanford.edu Developng a new SVM algorhm s ongong research opc. Among many exng SVM algorhms, we wll focus on
More informationHow about the more general "linear" scalar functions of scalars (i.e., a 1st degree polynomial of the following form with a constant term )?
lmcd Lnear ransformaon of a vecor he deas presened here are que general hey go beyond he radonal mar-vecor ype seen n lnear algebra Furhermore, hey do no deal wh bass and are equally vald for any se of
More informationOrdinary Differential Equations in Neuroscience with Matlab examples. Aim 1- Gain understanding of how to set up and solve ODE s
Ordnary Dfferenal Equaons n Neuroscence wh Malab eamples. Am - Gan undersandng of how o se up and solve ODE s Am Undersand how o se up an solve a smple eample of he Hebb rule n D Our goal a end of class
More informationGraduate Macroeconomics 2 Problem set 5. - Solutions
Graduae Macroeconomcs 2 Problem se. - Soluons Queson 1 To answer hs queson we need he frms frs order condons and he equaon ha deermnes he number of frms n equlbrum. The frms frs order condons are: F K
More informationA 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 informationDEEP UNFOLDING FOR MULTICHANNEL SOURCE SEPARATION SUPPLEMENTARY MATERIAL
DEEP UNFOLDING FOR MULTICHANNEL SOURCE SEPARATION SUPPLEMENTARY MATERIAL Sco Wsdom, John Hershey 2, Jonahan Le Roux 2, and Shnj Waanabe 2 Deparmen o Elecrcal Engneerng, Unversy o Washngon, Seale, WA, USA
More informationPhD/MA Econometrics Examination. January, 2019
Economercs Comprehensve Exam January 2019 Toal Tme: 8 hours MA sudens are requred o answer from A and B. PhD/MA Economercs Examnaon January, 2019 PhD sudens are requred o answer from A, B, and C. The answers
More informationAnalysis And Evaluation of Econometric Time Series Models: Dynamic Transfer Function Approach
1 Appeared n Proceedng of he 62 h Annual Sesson of he SLAAS (2006) pp 96. Analyss And Evaluaon of Economerc Tme Seres Models: Dynamc Transfer Funcon Approach T.M.J.A.COORAY Deparmen of Mahemacs Unversy
More informationEpistemic Game Theory: Online Appendix
Epsemc Game Theory: Onlne Appendx Edde Dekel Lucano Pomao Marcano Snscalch July 18, 2014 Prelmnares Fx a fne ype srucure T I, S, T, β I and a probably µ S T. Le T µ I, S, T µ, βµ I be a ype srucure ha
More informationComb Filters. Comb Filters
The smple flers dscussed so far are characered eher by a sngle passband and/or a sngle sopband There are applcaons where flers wh mulple passbands and sopbands are requred Thecomb fler s an example of
More information5th International Conference on Advanced Design and Manufacturing Engineering (ICADME 2015)
5h Inernaonal onference on Advanced Desgn and Manufacurng Engneerng (IADME 5 The Falure Rae Expermenal Sudy of Specal N Machne Tool hunshan He, a, *, La Pan,b and Bng Hu 3,c,,3 ollege of Mechancal and
More informationTHEORETICAL 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 informationOnline Appendix for. Strategic safety stocks in supply chains with evolving forecasts
Onlne Appendx for Sraegc safey socs n supply chans wh evolvng forecass Tor Schoenmeyr Sephen C. Graves Opsolar, Inc. 332 Hunwood Avenue Hayward, CA 94544 A. P. Sloan School of Managemen Massachuses Insue
More informationLecture 6: Learning for Control (Generalised Linear Regression)
Lecure 6: Learnng for Conrol (Generalsed Lnear Regresson) Conens: Lnear Mehods for Regresson Leas Squares, Gauss Markov heorem Recursve Leas Squares Lecure 6: RLSC - Prof. Sehu Vjayakumar Lnear Regresson
More informationCHAPTER 5: MULTIVARIATE METHODS
CHAPER 5: MULIVARIAE MEHODS Mulvarae Daa 3 Mulple measuremens (sensors) npus/feaures/arbues: -varae N nsances/observaons/eamples Each row s an eample Each column represens a feaure X a b correspons o he
More informationMechanics Physics 151
Mechancs Physcs 5 Lecure 0 Canoncal Transformaons (Chaper 9) Wha We Dd Las Tme Hamlon s Prncple n he Hamlonan formalsm Dervaon was smple δi δ Addonal end-pon consrans pq H( q, p, ) d 0 δ q ( ) δq ( ) δ
More informationBayesian Inference of the GARCH model with Rational Errors
0 Inernaonal Conference on Economcs, Busness and Markeng Managemen IPEDR vol.9 (0) (0) IACSIT Press, Sngapore Bayesan Inference of he GARCH model wh Raonal Errors Tesuya Takash + and Tng Tng Chen Hroshma
More informationMath 128b Project. Jude Yuen
Mah 8b Proec Jude Yuen . Inroducon Le { Z } be a sequence of observed ndependen vecor varables. If he elemens of Z have a on normal dsrbuon hen { Z } has a mean vecor Z and a varancecovarance marx z. Geomercally
More informatione-journal Reliability: Theory& Applications No 2 (Vol.2) Vyacheslav Abramov
June 7 e-ournal Relably: Theory& Applcaons No (Vol. CONFIDENCE INTERVALS ASSOCIATED WITH PERFORMANCE ANALYSIS OF SYMMETRIC LARGE CLOSED CLIENT/SERVER COMPUTER NETWORKS Absrac Vyacheslav Abramov School
More informationCHAPTER FOUR REPEATED MEASURES IN TOXICITY TESTING
CHAPTER FOUR REPEATED MEASURES IN TOXICITY TESTING 4. Inroducon The repeaed measures sudy s a very commonly used expermenal desgn n oxcy esng because no only allows one o nvesgae he effecs of he oxcans,
More informationLecture VI Regression
Lecure VI Regresson (Lnear Mehods for Regresson) Conens: Lnear Mehods for Regresson Leas Squares, Gauss Markov heorem Recursve Leas Squares Lecure VI: MLSC - Dr. Sehu Vjayakumar Lnear Regresson Model M
More informationDynamic Team Decision Theory. EECS 558 Project Shrutivandana Sharma and David Shuman December 10, 2005
Dynamc Team Decson Theory EECS 558 Proec Shruvandana Sharma and Davd Shuman December 0, 005 Oulne Inroducon o Team Decson Theory Decomposon of he Dynamc Team Decson Problem Equvalence of Sac and Dynamc
More informationSOME NOISELESS CODING THEOREMS OF INACCURACY MEASURE OF ORDER α AND TYPE β
SARAJEVO JOURNAL OF MATHEMATICS Vol.3 (15) (2007), 137 143 SOME NOISELESS CODING THEOREMS OF INACCURACY MEASURE OF ORDER α AND TYPE β M. A. K. BAIG AND RAYEES AHMAD DAR Absrac. In hs paper, we propose
More informationPart II CONTINUOUS TIME STOCHASTIC PROCESSES
Par II CONTINUOUS TIME STOCHASTIC PROCESSES 4 Chaper 4 For an advanced analyss of he properes of he Wener process, see: Revus D and Yor M: Connuous marngales and Brownan Moon Karazas I and Shreve S E:
More informationLet s treat the problem of the response of a system to an applied external force. Again,
Page 33 QUANTUM LNEAR RESPONSE FUNCTON Le s rea he problem of he response of a sysem o an appled exernal force. Agan, H() H f () A H + V () Exernal agen acng on nernal varable Hamlonan for equlbrum sysem
More informationLecture 18: The Laplace Transform (See Sections and 14.7 in Boas)
Lecure 8: The Lalace Transform (See Secons 88- and 47 n Boas) Recall ha our bg-cure goal s he analyss of he dfferenal equaon, ax bx cx F, where we emloy varous exansons for he drvng funcon F deendng on
More informationData 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 informationMechanics Physics 151
Mechancs Physcs 5 Lecure 9 Hamlonan Equaons of Moon (Chaper 8) Wha We Dd Las Tme Consruced Hamlonan formalsm H ( q, p, ) = q p L( q, q, ) H p = q H q = p H = L Equvalen o Lagrangan formalsm Smpler, bu
More informationABSTRACT KEYWORDS. Bonus-malus systems, frequency component, severity component. 1. INTRODUCTION
EERAIED BU-MAU YTEM ITH A FREQUECY AD A EVERITY CMET A IDIVIDUA BAI I AUTMBIE IURACE* BY RAHIM MAHMUDVAD AD HEI HAAI ABTRACT Frangos and Vronos (2001) proposed an opmal bonus-malus sysems wh a frequency
More informationCHAPTER 10: LINEAR DISCRIMINATION
CHAPER : LINEAR DISCRIMINAION Dscrmnan-based Classfcaon 3 In classfcaon h K classes (C,C,, C k ) We defned dscrmnan funcon g j (), j=,,,k hen gven an es eample, e chose (predced) s class label as C f g
More informationTHE PREDICTION OF COMPETITIVE ENVIRONMENT IN BUSINESS
THE PREICTION OF COMPETITIVE ENVIRONMENT IN BUSINESS INTROUCTION The wo dmensonal paral dfferenal equaons of second order can be used for he smulaon of compeve envronmen n busness The arcle presens he
More informationMechanics Physics 151
Mechancs Physcs 5 Lecure 9 Hamlonan Equaons of Moon (Chaper 8) Wha We Dd Las Tme Consruced Hamlonan formalsm Hqp (,,) = qp Lqq (,,) H p = q H q = p H L = Equvalen o Lagrangan formalsm Smpler, bu wce as
More informationApproximate Analytic Solution of (2+1) - Dimensional Zakharov-Kuznetsov(Zk) Equations Using Homotopy
Arcle Inernaonal Journal of Modern Mahemacal Scences, 4, (): - Inernaonal Journal of Modern Mahemacal Scences Journal homepage: www.modernscenfcpress.com/journals/jmms.aspx ISSN: 66-86X Florda, USA Approxmae
More information12d Model. Civil and Surveying Software. Drainage Analysis Module Detention/Retention Basins. Owen Thornton BE (Mech), 12d Model Programmer
d Model Cvl and Surveyng Soware Dranage Analyss Module Deenon/Reenon Basns Owen Thornon BE (Mech), d Model Programmer owen.hornon@d.com 4 January 007 Revsed: 04 Aprl 007 9 February 008 (8Cp) Ths documen
More informationAn introduction to Support Vector Machine
An nroducon o Suppor Vecor Machne 報告者 : 黃立德 References: Smon Haykn, "Neural Neworks: a comprehensve foundaon, second edon, 999, Chaper 2,6 Nello Chrsann, John Shawe-Tayer, An Inroducon o Suppor Vecor Machnes,
More informationHEAT 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 informationNPTEL Project. Econometric Modelling. Module23: Granger Causality Test. Lecture35: Granger Causality Test. Vinod Gupta School of Management
P age NPTEL Proec Economerc Modellng Vnod Gua School of Managemen Module23: Granger Causaly Tes Lecure35: Granger Causaly Tes Rudra P. Pradhan Vnod Gua School of Managemen Indan Insue of Technology Kharagur,
More informationOnline Supplement for Dynamic Multi-Technology. Production-Inventory Problem with Emissions Trading
Onlne Supplemen for Dynamc Mul-Technology Producon-Invenory Problem wh Emssons Tradng by We Zhang Zhongsheng Hua Yu Xa and Baofeng Huo Proof of Lemma For any ( qr ) Θ s easy o verfy ha he lnear programmng
More informationLi An-Ping. Beijing , P.R.China
A New Type of Cpher: DICING_csb L An-Png Bejng 100085, P.R.Chna apl0001@sna.com Absrac: In hs paper, we wll propose a new ype of cpher named DICING_csb, whch s derved from our prevous sream cpher DICING.
More informationLecture Notes 4. Univariate Forecasting and the Time Series Properties of Dynamic Economic Models
Tme Seres Seven N. Durlauf Unversy of Wsconsn Lecure Noes 4. Unvarae Forecasng and he Tme Seres Properes of Dynamc Economc Models Ths se of noes presens does hree hngs. Frs, formulas are developed o descrbe
More information2. SPATIALLY LAGGED DEPENDENT VARIABLES
2. SPATIALLY LAGGED DEPENDENT VARIABLES In hs chaper, we descrbe a sascal model ha ncorporaes spaal dependence explcly by addng a spaally lagged dependen varable y on he rgh-hand sde of he regresson equaon.
More informationSurvival Analysis and Reliability. A Note on the Mean Residual Life Function of a Parallel System
Communcaons n Sascs Theory and Mehods, 34: 475 484, 2005 Copyrgh Taylor & Francs, Inc. ISSN: 0361-0926 prn/1532-415x onlne DOI: 10.1081/STA-200047430 Survval Analyss and Relably A Noe on he Mean Resdual
More informationFall 2009 Social Sciences 7418 University of Wisconsin-Madison. Problem Set 2 Answers (4) (6) di = D (10)
Publc Affars 974 Menze D. Chnn Fall 2009 Socal Scences 7418 Unversy of Wsconsn-Madson Problem Se 2 Answers Due n lecure on Thursday, November 12. " Box n" your answers o he algebrac quesons. 1. Consder
More informationExistence and Uniqueness Results for Random Impulsive Integro-Differential Equation
Global Journal of Pure and Appled Mahemacs. ISSN 973-768 Volume 4, Number 6 (8), pp. 89-87 Research Inda Publcaons hp://www.rpublcaon.com Exsence and Unqueness Resuls for Random Impulsve Inegro-Dfferenal
More informationON THE WEAK LIMITS OF SMOOTH MAPS FOR THE DIRICHLET ENERGY BETWEEN MANIFOLDS
ON THE WEA LIMITS OF SMOOTH MAPS FOR THE DIRICHLET ENERGY BETWEEN MANIFOLDS FENGBO HANG Absrac. We denfy all he weak sequenal lms of smooh maps n W (M N). In parcular, hs mples a necessary su cen opologcal
More informationRobust and Accurate Cancer Classification with Gene Expression Profiling
Robus and Accurae Cancer Classfcaon wh Gene Expresson Proflng (Compuaonal ysems Bology, 2005) Auhor: Hafeng L, Keshu Zhang, ao Jang Oulne Background LDA (lnear dscrmnan analyss) and small sample sze problem
More informationDensity Matrix Description of NMR BCMB/CHEM 8190
Densy Marx Descrpon of NMR BCMBCHEM 89 Operaors n Marx Noaon Alernae approach o second order specra: ask abou x magnezaon nsead of energes and ranson probables. If we say wh one bass se, properes vary
More informationNotes 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 informationCH.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 informationFTCS Solution to the Heat Equation
FTCS Soluon o he Hea Equaon ME 448/548 Noes Gerald Reckenwald Porland Sae Unversy Deparmen of Mechancal Engneerng gerry@pdxedu ME 448/548: FTCS Soluon o he Hea Equaon Overvew Use he forward fne d erence
More informationChapter 6: AC Circuits
Chaper 6: AC Crcus Chaper 6: Oulne Phasors and he AC Seady Sae AC Crcus A sable, lnear crcu operang n he seady sae wh snusodal excaon (.e., snusodal seady sae. Complee response forced response naural response.
More information2.1 Constitutive Theory
Secon.. Consuve Theory.. Consuve Equaons Governng Equaons The equaons governng he behavour of maerals are (n he spaal form) dρ v & ρ + ρdv v = + ρ = Conservaon of Mass (..a) d x σ j dv dvσ + b = ρ v& +
More informationPerformance Analysis for a Network having Standby Redundant Unit with Waiting in Repair
TECHNI Inernaonal Journal of Compung Scence Communcaon Technologes VOL.5 NO. July 22 (ISSN 974-3375 erformance nalyss for a Nework havng Sby edundan Un wh ang n epar Jendra Sngh 2 abns orwal 2 Deparmen
More informationRobustness of DEWMA versus EWMA Control Charts to Non-Normal Processes
Journal of Modern Appled Sascal Mehods Volume Issue Arcle 8 5--3 Robusness of D versus Conrol Chars o Non- Processes Saad Saeed Alkahan Performance Measuremen Cener of Governmen Agences, Insue of Publc
More informationKayode Ayinde Department of Pure and Applied Mathematics, Ladoke Akintola University of Technology P. M. B. 4000, Ogbomoso, Oyo State, Nigeria
Journal of Mahemacs and Sascs 3 (4): 96-, 7 ISSN 549-3644 7 Scence Publcaons A Comparave Sudy of he Performances of he OLS and some GLS Esmaors when Sochasc egressors are boh Collnear and Correlaed wh
More informationRelative controllability of nonlinear systems with delays in control
Relave conrollably o nonlnear sysems wh delays n conrol Jerzy Klamka Insue o Conrol Engneerng, Slesan Techncal Unversy, 44- Glwce, Poland. phone/ax : 48 32 37227, {jklamka}@a.polsl.glwce.pl Keywor: Conrollably.
More informationImprovement in Estimating Population Mean using Two Auxiliary Variables in Two-Phase Sampling
Improvemen n Esmang Populaon Mean usng Two Auxlar Varables n Two-Phase amplng Rajesh ngh Deparmen of ascs, Banaras Hndu Unvers(U.P.), Inda (rsnghsa@ahoo.com) Pankaj Chauhan and Nrmala awan chool of ascs,
More information1 Constant Real Rate C 1
Consan Real Rae. Real Rae of Inees Suppose you ae equally happy wh uns of he consumpon good oday o 5 uns of he consumpon good n peod s me. C 5 Tha means you ll be pepaed o gve up uns oday n eun fo 5 uns
More informationOutline. Probabilistic Model Learning. Probabilistic Model Learning. Probabilistic Model for Time-series Data: Hidden Markov Model
Probablsc Model for Tme-seres Daa: Hdden Markov Model Hrosh Mamsuka Bonformacs Cener Kyoo Unversy Oulne Three Problems for probablsc models n machne learnng. Compung lkelhood 2. Learnng 3. Parsng (predcon
More informationBundling with Customer Self-Selection: A Simple Approach to Bundling Low Marginal Cost Goods On-Line Appendix
Bundlng wh Cusomer Self-Selecon: A Smple Approach o Bundlng Low Margnal Cos Goods On-Lne Appendx Lorn M. H Unversy of Pennsylvana, Wharon School 57 Jon M. Hunsman Hall Phladelpha, PA 94 lh@wharon.upenn.edu
More informationImprovement in Estimating Population Mean using Two Auxiliary Variables in Two-Phase Sampling
Rajesh ngh Deparmen of ascs, Banaras Hndu Unvers(U.P.), Inda Pankaj Chauhan, Nrmala awan chool of ascs, DAVV, Indore (M.P.), Inda Florenn marandache Deparmen of Mahemacs, Unvers of New Meco, Gallup, UA
More informationModeling and Solving of Multi-Product Inventory Lot-Sizing with Supplier Selection under Quantity Discounts
nernaonal ournal of Appled Engneerng Research SSN 0973-4562 Volume 13, Number 10 (2018) pp. 8708-8713 Modelng and Solvng of Mul-Produc nvenory Lo-Szng wh Suppler Selecon under Quany Dscouns Naapa anchanaruangrong
More informationScattering at an Interface: Oblique Incidence
Course Insrucor Dr. Raymond C. Rumpf Offce: A 337 Phone: (915) 747 6958 E Mal: rcrumpf@uep.edu EE 4347 Appled Elecromagnecs Topc 3g Scaerng a an Inerface: Oblque Incdence Scaerng These Oblque noes may
More informationOn computing differential transform of nonlinear non-autonomous functions and its applications
On compung dfferenal ransform of nonlnear non-auonomous funcons and s applcaons Essam. R. El-Zahar, and Abdelhalm Ebad Deparmen of Mahemacs, Faculy of Scences and Humanes, Prnce Saam Bn Abdulazz Unversy,
More informationNormal Random Variable and its discriminant functions
Noral Rando Varable and s dscrnan funcons Oulne Noral Rando Varable Properes Dscrnan funcons Why Noral Rando Varables? Analycally racable Works well when observaon coes for a corruped snle prooype 3 The
More informationEndogeneity. Is the term given to the situation when one or more of the regressors in the model are correlated with the error term such that
s row Endogeney Is he erm gven o he suaon when one or more of he regressors n he model are correlaed wh he error erm such ha E( u 0 The 3 man causes of endogeney are: Measuremen error n he rgh hand sde
More information