Approximations to the Loglikelihood Function in the Nonlinear Mixed Effects Model
|
|
- Christiana George
- 6 years ago
- Views:
Transcription
1 Approxmatons to the Lolkelhood Functon n the Nonlnear Mxed Effects Model José C. Pnhero and Doulas M. Bates Department of Statstcs Unversty of Wsconsn, Madson Nonlnear mxed effects models have receved a reat deal of attenton n the statstcal lterature n recent years because of the flexblty they offer n handln unbalanced repeated measures data that arse n dfferent areas of nvestaton, such as pharmacoknetcs and economcs. Several dfferent methods for estmatn the parameters n nonlnear mxed effects model have been proposed. We concentrate here on two of them: maxmum lkelhood and restrcted maxmum lkelhood. A rather complex numercal ssue for (restrcted) maxmum lkelhood estmaton s the evaluaton of the lolkelhood functon of the data, snce t nvolves the evaluaton of a multple nteral that n most cases does not have a closed form expresson. We consder here four dfferent approxmatons to the lolkelhood, comparn ther computatonal and statstcal propertes. We conclude that the alternatn approxmaton suested by Lndstrom and Bates (1990), the Laplacan approxmaton, and Gaussan quadrature centered at the condtonal modes of the random effects are qute accurate and computatonally effcent. Gaussan quadrature centered at the expected value of the random effects s qute naccurate for a smaller number of abscssas and computatonally neffcent for a larer number of abscssas. Importance sampln s accurate but qute neffcent computatonally. Keywords: Nonlnear mxed effects models, maxmum lkelhood estmaton, Laplacan approxmaton, Gaussan quadrature, mportance sampln. 1 Introducton. Several dfferent nonlnear mxed effects models and estmaton methods for ther parameters have been proposed n recent years (Shener and Beal, 1980; Mallet, Mentre, Stemer and Lokek, 1988; Lndstrom and Bates, 1990; Vonesh and Carter, 199; Davdan and Gallant, 199; Wakefeld, Smth, Racne-Poon and Gelfand, 1994). We consder here a slhtly modfed verson of the model proposed n Lndstrom 1
2 and Bates (1990). Ths model can be vewed as a herarchcal model that n some ways eneralzes both the lnear mxed effects model of Lard and Ware (198) and the usual nonlnear model for ndependent data (Bates and Watts, 1988). In the frst stae the jth observaton on the th cluster s modeled as y j = f ( j ; x j ) + j ; = 1; : : :; M; j = 1; : : :; n (1) where f s a nonlnear functon of a cluster-specfc parameter vector j and the predctor vector x j, j s a normally dstrbuted nose term, M s the total number of clusters, and n s the number of observatons n the th cluster. In the second stae the cluster-specfc parameter vector s modeled as j = A j + B j b ; b N (0; D); where s a p-dmensonal vector of fxed populaton parameters, b s a q-dmensonal random effects vector assocated wth the th cluster (not varyn wth j), A j and B j are desn matrces for the fxed and random effects respectvely, and D s a (eneral) varance-covarance matrx. It s further assumed that observatons made on dfferent clusters are ndependent and that the j are..d. N (0; ) and ndependent of the b. We consder estmaton of the model s parameters by ether maxmum lkelhood, or restrcted maxmum lkelhood, based on the marnal densty of y Z p(y j ; D; ) = p(y j b; ; D; ) p(b) db () In eneral ths nteral does not have a closed-form expresson when the model functon f s nonlnear n b so dfferent approxmatons have been proposed for estmatn t. Some of these methods consst of takn a frst order Taylor expanson of the model functon f around the expected value of the random effects (Shener and Beal, 1980; Vonesh and Carter, 199), or around the condtonal (on D) modes of the random effects (Lndstrom and Bates, 1990). Others have proposed the use of Gaussan quadrature rules (Davdan and Gallant, 199). We consder here four dfferent approxmatons to the lolkelhood (): Lndstrom and Bates (1990) s alternatn method, a modfed Laplacan approxmaton (Terney and Kadane, 1986), mportance sampln (Geweke, 1989), and Gaussan quadrature (Davdan and Gallant, 199). We compare
3 them based on ther computatonal and statstcal propertes, usn both real data examples and smulaton results. Secton contans a descrpton of the dfferent approxmatons to the lolkelhood as appled to the nonlnear mxed effects model (1). Secton 3 presents a comparson of the dfferent approxmatons based on real and smulated data. Our conclusons and suestons for further nvestaton are ven n secton 4. Approxmatons to the Lolkelhood In ths secton we descrbe four dfferent approxmatons to the lolkelhood of y n the nonlnear mxed effects model (1). We show that there exsts a close relaton between the Laplacan approxmaton, mportance sampln and a Gaussan quadrature rule centered around the condtonal modes of the random effects b..1 Alternatn Approxmaton Lndstrom and Bates (1990) propose an alternatn alorthm for estmatn the parameters n model (1). Based on the current estmates of D (the scaled varance-covarance matrx of the random effects), the condtonal modes of the random effects b and the condtonal estmates of the fxed effects are obtaned by mnmzn a penalzed nonlnear least squares (PNLS) objectve functon where [f (; b )] j = f MX =1 ky? f (; b )k + b T D?1 b j ; x j ; = 1; : : :; M; j = 1; : : :; n. To update the estmate of D at the wth teraton, Lndstrom and Bates use a frst order Taylor expanson of the model functon around the current estmates of and the condtonal modes of the (w) (w) random effects b, whch we wll denote by b and b respectvely. Lettn T ; d X b;b T ; and bw (w) = y? f (b (w) ; b b (w) ) + c X (w) b (w) (w) + Z b bb (w) (3) 3
4 the approxmate lolkelhood used for the estmaton of D s? `A ; ; D j y =? 1 MX =1 lo +? bw (w) (w) (w) I + Z b D Z b (w)? X c T T I + b Z (w) (4) ) T (w) (w) D Z b? X c?1 bw (w) Ths lolkelhood s dentcal to that of a lnear mxed effects (LME) model (Lard and Ware, 198) n whch the response vector s ven by bw (w) and the fxed and random effects desn matrces are ven by c X (w) and b Z (w). Usn the results n Lndstrom and Bates (1988) one can express the optmal values of and as functons of D and work wth the profle lolkelhood of D, reatly smplfyn the optmzaton problem. Lndstrom and Bates (1990) have also proposed an approxmate restrcted lolkelhood for the estmaton of D? `RA ; ; D j y =? 1 MX =1 lo (w) c X T (w) (w) I + Z b D Z b T cx (w)? + `A ; ; D j y (5) Ther estmaton alorthm alternates between the PNLS and LME steps untl some converence crteron s met. Such alternatn alorthms tend to be more effcent when the estmates of the varancecovarance components (D and ) are not hhly correlated wth the estmates of the fxed effects (). Pnhero and Bates (1993) have demonstrated that, n the lnear mxed effects model, the (restrcted) maxmum lkelhood estmates of D and are asymptotcally ndependent of the (restrcted) maxmum lkelhood estmates of. These results have not yet been extended to the nonlnear mxed effects model (1). It can be shown that the maxmum lkelhood estmate of and the condtonal modes of the random effects b correspondn to the approxmate lolkelhood (4) are the values obtaned n the frst teraton of the Gauss-Newton alorthm used to mnmze the PNLS objectve functon (3). Therefore, at the convered value of c D, the estmates of and b obtaned from the LME and PNLS steps concde. We wll use `A when comparn the dfferent approxmatons at the optmal values n secton 3, but we do note that n Lndstrom and Bates (1990) approxmaton (4) s used only to update the estmates of D and not for estmatn. 4
5 . Laplacan Approxmaton Laplacan approxmatons are frequently used n Bayesan nference to estmate marnal posteror denstes and predctve dstrbutons (Terney and Kadane, 1986; Leonard, Hsu and Tsu, 1989). These technques can also be used for the nteraton consdered here. The nteral that we want to estmate for the marnal dstrbuton of y n model (1) can be wrtten as p(y j ; D; ) = Z??(n +q)= jdj?1= exp?(; D; y ; b )= db ; where (; D; y ; b ) = ky? f (; b )k + b T D?1 b Let bb = b b (; D; y ) = ar mn (; D; y ; b ) b 0 (; D; y ; b ) D; y ; b 00 (; D; y ; b ) (; D; y ; T and consder a second order Taylor expanson of around b b (; D; y ; b ) ' ; D; y ; b b + 1 h b? b b T 00 ; D; y ; b b hb? b b (6) where the lnear term of the approxmaton vanshes snce 0 (; D; y ; b b ) = 0. approxmaton s defned as p? y j ; D; '? "?N= jdj?m= exp? 1 Z? ( q= exp =??N= jdj?m= M Y =1? 1 MX =1 MX =1 ; D; y ; b b # h b? b b T ; D; y ; b b?1= exp The Laplacan ; D; y ; b b h b? b b ) db h? ; D; y ; b b = where N = P M =1 n. Now we consder an approxmaton to 00 smlar to the one used n Gauss-Newton optmzaton. 5
6 We have 00 ; D; y ; b b f (; T b=bb At b b, the contrbuton f (; b T that b )=@b T b=bb 00 ; D; y ; b b h y? f (; b b ) (; b T b ) + b=bb h y? f (; b b ) s usually nelble compared to (; b )=@b j b=bb (Bates and Watts, 1980) so we use the approxmaton ' G (; D; y ) (; b T (; b ) + b=bb Ths has the advantae of requrn only the frst order partal dervatves of the model functon wth respect to the random effects, whch are usually avalable from the estmaton of b b. Ths estmaton of bb s a penalzed least squares problem, for whch standard and relable code s avalable. The modfed Laplacan approxmaton to the lolkelhood of model (1) s then ven by? `LA ; D; j y = (? 1 N lo? + M lo jdj + MX =1 lo [G (; D; y )] +? M X =1 ; D; y ; b b ) (7) Snce b b does not depend upon, for ven and D the maxmum lkelhood estmate of (based upon `LA ) s b = b (; D; y) = MX =1 ; D; y ; b b =N We can profle `LA on to reduce the dmenson of the optmzaton problem, obtann `LAp =? 1 ( ) N 1 + lo () + lo? b + M lo jdj + M X =1 lo [G (; D; y )] (8) We note that f f s lnear n b then the modfed Laplacan approxmaton s exact because the second order Taylor expanson n (6) s exact when f (; b) = f () + Z () b.. There does not yet seem to be a strahtforward eneralzaton of the concept of restrcted maxmum lkelhood (Harvlle, 1974) to nonlnear mxed effects models. The dffculty s that restrcted maxmum lkelhood depends heavly upon the lnearty of the fxed effects n the model functon, whch does not occur n nonlnear models. Lndstrom and Bates (1990) crcumvented that problem by usn an 6
7 approxmaton to the model functon f n whch the fxed effects occur lnearly. Ths cannot be done for the Laplacan approxmaton, unless we consder yet another Taylor expanson of the model functon, what would lead us back to somethn very smlar to Lndstrom and Bates approach. We wll return to ths topc later n secton 4..3 Importance Sampln Importance sampln provdes a smple and effcent way of performn Monte Carlo nteraton. The crtcal step for the success of ths method s the choce of an mportance dstrbuton from whch the sample s drawn and the mportance wehts calculated. Ideally ths dstrbuton corresponds to the densty that we are tryn to nterate, but n practce one uses an easly sampled approxmaton. For the nonlnear mxed effects model the functon that we want to nterate s, up to a multplcatve constant, equal to exp [? (; D; y ; b ) = ]. As shown n subsecton., by takn a second order Taylor expanson of (; D; y ; b ) around b b the nterand s, up to a multplcatve constant, approxmately equal to a N bb ; [G(; D; y )]?1 densty. Ths ves us a natural choce for the mportance dstrbuton. Let N IS denote the number of mportance samples to be drawn. In practce one such sample can be enerated by selectn a vector z wth dstrbuton N (0; I) and calculatn the sample of random effects as b = b b + [G (; D; y )]?1= z, where [G (; D; y )]?1= denotes the nverse of the Cholesky factor of G (; D; y ). The mportance sampln approxmaton to the lolkelhood of y s then defned as? " `IS ; D; j y =? 1? N lo # MX + M lo jdj + lo jg (; D; y )j =1 8 9 MX < NX IS h = + lo exp? ; D; y : ; b j = + kz j k = =N IS ; =1 j=1 (9) Note that we cannot n eneral obtan a closed form expresson for the MLE of for fxed and D, so that profln on s no loner reasonable. As n the modfed Laplacan approxmaton, mportance sampln ves exact results when the 7
8 model functon s lnear n b because n ths case? p(y j b ; ; D; ) p(b ) = p y j ; D; N bb ; [G (; D; y?1 )] so that the mportance wehts are equal to p (y j ; D; )..4 Gaussan quadrature Gaussan quadrature s used to approxmate nterals of functons wth respect to a ven kernel by a wehted averae of the nterand evaluated at pre-determned abscssas. The wehts and abscssas used n Gaussan quadrature rules for the most common kernels can be obtaned from the tables of Abramowtz and Steun (1964) or by usn an alorthm proposed by Golub (1973) (see also Golub and Welsch (1969)). Gaussan quadrature rules for multple nterals are known to be numercally complex (Davs and Rabnowtz, 1984), but usn the structure of the nterand n the nonlnear mxed effects model we can transform the problem nto successve applcatons of smple one dmensonal Gaussan quadrature rules. Lettn z j ; w j j = 1; : : :; N GQ denote respectvely the abscssas and the wehts for the (one dmensonal) Gaussan quadrature rule wth N GQ ponts based on the N (0; 1) kernel, we et Z ( )?q= jdj?1= exp h? ky? f (; b )k = exp?b T D?1 b= Z = ' where z j 1;:::;j q = NX GQ ()?q= exp j 1=1 NX GQ j q=1 exp? y? f? y? f ; D T= z = ; D T= z j 1;:::;j q = exp? kz k = q Y k=1 w jk db (10) z j 1 ; : : :; z j q T. The correspondn approxmaton to the lolkelhood functon s dz? `GQ ; D; j y = (11) 8 MX < NX GQ 9? N lo( )= + lo exp? y :? f ; D T= z j Y q = = w jk ; =1 j k=1 where j = (j 1 ; : : :; j q ) T. The Gaussan quadrature rule n ths case can be vewed as a determnstc verson of Monte Carlo nteraton n whch random samples of b are enerated from the N (0; D) dstrbuton. The samples 8
9 (z j ) and the wehts (w j) are fxed beforehand, whle n Monte Carlo nteraton they are left to random choce. Snce mportance sampln tends to be much more effcent than smple Monte Carlo nteraton, we also consder the equvalent of mportance sampln n the Gaussan quadrature context, whch we wll denote by adaptve Gaussan quadrature. In ths approach the rd of abscssas n the b scale s centered around the condtonal modes b b rather than 0, as n (10). Another modfcaton s the use of G (; D; y ) nstead of D n the scaln of the z. The adaptve Gaussan quadrature s then ven by Z ( )?q= jdj?1= exp h? ky? f (; b )k = exp?b T D?1 b= Z = ' ()?q= jg (; D; y ) Dj?1= exp + kz k = exp? kz k = dz NX GQ j 1=1 NX GQ j q=1 exp?? db n; D; y ;b b + [G (; D; y )]?1= z o = n; D; y ; b b + [G (; D; y )]?1= z j 1;:::;j q o The correspondn approxmaton to the lolkelhood s then? `AGQ ; D; j y =? MX X + =1 lo 4 NGQ j " N lo? + M lo jdj + MX =1 = + lo jg (; D; y )j n exp? ; D; y ;b b + [G (; D; y )]?1= z j # o = + z Y q j 1;:::;j q = = (1) z j = qy The adaptve Gaussan quadrature approxmaton very closely resembles that obtaned for mportance sampln. The basc dfference s that the former uses fxed abscssas and wehts, whle the latter allows them to be determned by a pseudo-random mechansm. It s also nterestn to note that the one pont (.e. N GQ = 1) adaptve Gaussan quadrature approxmaton s smply the modfed Laplacan approxmaton (8), snce n ths case z 1 = 0 and w 1 = 1. The adaptve Gaussan quadrature also ves the exact lolkelhood when the model functon s lnear n b, but that s not true n eneral for the Gaussan quadrature approxmaton (10). Lke the mportance sampln approxmaton, the Gaussan quadrature approxmaton cannot be profled on to reduce the dmensonalty of the optmzaton problem. k=1 w jk 3 5 k=1 w jk 9
10 3 Comparn the Approxmatons In ths secton we present a comparson of the dfferent approxmatons to the lolkelhood of model (1) descrbed n secton. Two real data examples, the orane trees and Theophyllne data sets, and smulaton results are used to compare the statstcal and computatonal aspects of the varous approxmatons. 3.1 Orane Trees The data are presented on Fure 1 and consst of seven measurements of the trunk crcumference (n mllmeters) on each of fve orane trees, taken over a perod of 1600 days. These data were ornally presented n Draper and Smth (1981, p. 54) and were descrbed n Lndstrom and Bates (1990). Tree crcumference Days Fure 1: Trunk crcumference (n mllmeters) of fve orane trees: Data and ndvdual ftted curves from maxmum lkelhood estmaton usn the exact lolkelhood. The dashed lne represents the mean curve. The lostc model y = 1 = f1 + exp [? (t? ) = 3 ] seems to ft the data well. Lndstrom and Bates (1990) concluded n ther analyss that only the asymptotc crcumference 1 needs a random 10
11 effect to account for tree to tree varaton and suested the follown nonlnear mxed effects model y j = 1 + b exp [? (t j? ) = 3 ] + " j (13) where y j represents the jth crcumference measurement on the th tree, t j represents the day correspondn to the jth measurement on the th tree, b 1 ; = 1; : : :; 5 are..d. N (0; D), and " j ; = 1; : : :; 5; j = 1; : : :; 7 are..d. N (0; ) and ndependent of the b 1. Note that the snle random effect occurs lnearly n (13) and therefore the modfed Laplacan (8), the mportance sampln (9), and the adaptve Gaussan quadrature (1) approxmatons are all exact. Table 1 presents the results of estmaton usn the alternatn approxmaton, Gaussan quadrature wth 10 and 00 abscssas, and the exact lolkelhood. Snce only the alternatn approxmaton provdes a verson of restrcted maxmum lolkelhood, we wll just consder maxmum lkelhood estmaton n ths and the next subsecton. The subscrpt on Gaussan refers to the number of abscssas used n the approxmaton and the scalar L s p D, the square root of the scaled varance of the random effects. In eneral ths s a matrx but there s only one random effect here. Table 1: Estmaton Results Orane Trees Data Approxmaton lo(l) 1 3 lo( ) ` Alternatn Gaussan Gaussan Exact The estmaton results n Table 1 ndcate that the dfferent approxmatons produce smlar fts. The Gaussan approxmaton wth only 10 abscssas ves the worst approxmaton, n terms of the value of the lolkelhood, but even that s not far from the exact value. The Gaussan quadrature wth 00 abscssas s almost dentcal to the exact lolkelhood. The alternatn approxmaton s also very close to the exact value. Another mportant ssue reardn the dfferent approxmatons s how well they behave n a nehborhood of the optmal value, snce ths behavor s often used to assess the varablty of maxmum lkelhood estmates. Fure dsplays the profle traces and contours (Bates and Watts, 1988) for the exact lolkelhood and the alternatn approxmaton. Ths plot could not be obtaned for the Gaussan approxmaton because the objectve functon presented several local optma durn the profln alo- 11
12 rthm. We beleve that ths s related to the fact that the Gaussan approxmaton s centered at b = 0 and not at the condtonal modes of the random effects, where the nterand n () takes ts hhest values lo(l) lo(s) beta beta beta Fure : Profle traces and profle contour plots for the orane trees data based on the exact lolkelhood (sold lne) and the alternatn approxmaton (dashed lne). Plots below the daonal are n the ornal scale and plots above the daonal are n the zeta scale (Bates and Watts, 1988). Interpolated contours correspond approxmately to jont confdence levels of 68%, 87%, and 95%. It can be seen from Fure that the alternatn method ves a ood approxmaton to the lolkelhood n a nehborhood of the optmal values. It s nterestn to note that the profle traces for the varance-covarance components (D and ) and the fxed effects () meet almost perpendcularly. Ths ndcates a local uncorrelaton between the varance-covarance components and the fxed effects, whch explans why the alternatn method was so successful n approxmatn the lolkelhood. The same pattern was observed n several other data sets that we have analyzed, leadn us to conjecture that the asymptotc uncorrelaton between the estmators of the varance-covarance components and the fxed 1
13 effects verfed n the lnear mxed effects model also holds, at least approxmately, for the nonlnear mxed effects model. To compare the computatonal effcency of the dfferent approxmatons we consder the number of functon evaluatons needed untl converence. For the alternatn approxmaton there are two dfferent functons ben evaluated durn the teratons: the objectve functon (3) wthn the PNLS step and the approxmate lolkelhood `A (4) wthn the LME step. We wll use here the total number of evaluatons of ether (3) or `A, multpled by the number of clusters. For the other approxmatons we wll use the total number of calls to (; D; y ; b ). Even thouh the number of functon evaluatons used for the alternatn approxmaton s not drectly comparable to the number of functon evaluatons of the remann approxmatons, t ves a ood dea of the relatve computatonal effcency of ths alorthm. Table presents the number of functon evaluatons for the dfferent approxmatons n the orane trees example. The Gaussan quadrature approxmatons are consderably less effcent than ether the alternatn approxmaton or the exact lolkelhood. As expected the alternatn approxmaton s the most computatonally effcent. Table : Number of Functon Evaluatons to Converence Orane Trees Data Approxmaton Functon Evaluatons Alternatn 00 Exact 40 Gaussan 10 8,150 Gaussan , Theophyllne Knetcs The data consdered here are courtesy of Dr. Robert A. Upton of the Unversty of Calforna, San Francsco. Theophyllne was admnstered orally to 1 subjects whose serum concentratons were measured at 11 tmes over the next 5 hours. Ths s an example of a laboratory pharmacoknetc study characterzed by many observatons on a moderate number of ndvduals (clusters). Fure 3 dsplays the data and the ndvdual fts obtaned throuh maxmum lkelhood usn the adaptve Gaussan approxmaton wth 10 abscssas. A common model for such data s a frst order compartment model wth absorpton n a perpheral 13
14 Tme (hrs) Concentraton (m/l) a a a a a a a a a a a b b b b b b b b b b b c c c c c c c c c c c d d d d d d d d d d d e e e e e e e e e e e f f f f f f f f f f f h hh h h h h h h h h j j j j j j j j j j j k k k k k k k k k k k l l l l l l l l l l l Fure 3: Theophyllne concentratons (n m/l) of twelve patents: Data and ndvdual ftted curves from maxmum lkelhood estmaton usn the adaptve Gaussan approxmaton. compartment C t = DKk a Cl(k a? K) [exp (?Kt)? exp (?k at)] (14) where C t s the observed concentraton at tme t (m/l), t s the tme (hr), D s the dose (m/k), Cl s the clearance (L/k), K s the elmnaton rate constant (1/hr), and k a s the absorpton rate constant (1/hr). In order to ensure postvty of the rate constants and the clearance, the loarthms of these quanttes were used n the ft. Analyss of the Theophyllne data usn model (14) ndcated that only lo(cl) and lo(k a ) needed random effects to account for the patent-to-patent varablty. The nonlnear mxed effects model used for the Theophyllne data s C t = D exp [? ( 1 + b 1 ) + ( + b ) + 3 ] exp ( + b )? exp ( 3 ) fexp [? exp ( 3 ) t]? exp [? exp ( + b ) t] (15) Table 3 presents the estmaton results from the varous approxmatons to the lolkelhood. Only maxmum lkelhood estmaton s consdered. The subscrpts on Gaussan and on Adap. Gaussan 14
15 refer to the number of abscssas used n the Gaussan and adaptve Gaussan approxmatons, whle the subscrpt on Imp. Sampln refers to the number of mportance samples used n ths approxmaton. L denotes the vector wth elements ven by the upper tranular half of the Cholesky decomposton of D, stacked by columns. Table 3: Estmaton Results Theophyllne Data Approxmaton lo(l 1 ) L lo(l 3 ) 1 3 lo( ) ` Alternatn Laplacan Imp. Sampln Gaussan Gaussan Gaussan Adap. Gaussan Adap. Gaussan We can see from Table 3 that the alternatn approxmaton, the Laplacan approxmaton, the mportance sampln approxmaton, and the adaptve Gaussan approxmaton all ve smlar estmaton results. The Gaussan approxmaton only approaches the other approxmatons when the number of abscssas s ncreased consderably. Note that the actual number of ponts used n the rd that defnes the Gaussan approxmaton for ths example s the square of the number of abscssas. The adaptve Gaussan approxmatons for 1 (Laplacan), 5, and 10 abscssas ve smlar results, ndcatn that just a few ponts are needed for ths approxmaton to be accurate. The mportance sampln approxmaton caused some numercal dffcultes for the optmzaton alorthm (the ms() functon n S (Chambers and Haste, 199)) used to obtan the maxmum lkelhood estmates, snce the stochastc varablty assocated wth dfferent mportance samples overwhelmed the numercal varablty of the lolkelhood for small chanes n the parameter values (used to calculate numercal dervatves). We ended up havn to keep the random number enerator seed fxed durn the optmzaton process, thus usn the same mportance samples throuhout the calculatons. Snce the results obtaned usn mportance sampln were very smlar to those of the adaptve Gaussan approxmaton, we concluded that the latter s to be preferred for ts reater smplcty and computatonal effcency. Table 4 ves the number of functon evaluatons untl converence for the dfferent approxmatons. The alternatn approxmaton s the most effcent, followed by the Laplacan and adaptve Gaussan approxmatons. Gaussan quadrature wth 5 abscssas s effcent compared to the adaptve Gaussan, but s qute naccurate. The more relable Gaussan approxmaton wth 100 abscssas takes about
16 tmes more functon evaluatons than the adaptve Gaussan wth 10 abscssas. The mportance sampln approxmaton had the worst performance n terms of functon evaluatons. Table 4: Number of Functon Evaluatons to Converence Theophyllne Data Approxmaton Functon Evaluatons Alternatn 1,51 Laplacan 7,683 Adap. Gaussan 5 30,00 Adap. Gaussan 10 96,784 Gaussan 5 47,700 Gaussan ,000 Gaussan ,00,000 Imp. Sampln ,11,84 Next we consder the approxmatons n a nehborhood of the optmal value. We wll restrct ourselves here to the alternatn, the Laplacan, and the adaptve Gaussan approxmaton, as the Gaussan approxmaton for a moderate number of abscssas s not relable, and both the Gaussan approxmaton wth a larer number of abscssas and the mportance sampln approxmaton are very neffcent computatonally and ve results qute smlar to the adaptve Gaussan approxmaton. We used fve abscssas for the adaptve Gaussan quadrature, as ths ves rouhly the same precson as the ten-abscssa quadrature rule. The alternatn approxmaton ves results very smlar to the adaptve Gaussan quadrature. As n the orane trees example, the profle traces of the varance-covarance components and the fxed effects meet almost perpendcularly, ndcatn a local uncorrelaton between these estmates. The Laplacan and the adaptve Gaussan approxmatons ve vrtually dentcal plots (not ncluded here). Ths suests there s lttle to be aned by ncreasn the number of abscssas past one n the quadrature rule. The major an n precson s obtaned by centern the rd at the condtonal modes and scaln t usn the approxmate Hessan. 3.3 Smulaton Results In ths secton we nclude a comparson of the approxmatons to the lolkelhood n model (1) usn smulaton. We restrct ourselves to the alternatn, the Laplacan, and the (fve-abscssa) adaptve Gaussan approxmatons as these seem to be more accurate and/or more effcent than the Gaussan and the mportance sampln approxmatons. Two models were used n the smulaton analyss: a 16
17 lo(l1) L lo(l3) lo(s) lo(cl) lo(ka) lo(k) Fure 4: Profle traces and profle contour plots for the Theophyllne data based on the adaptve Gaussan approxmaton wth 5 abscssas (sold lne) and the alternatn approxmaton (dashed lne). Plots below the daonal are n the ornal scale and plots above the daonal are n the zeta scale (Bates and Watts, 1988). Interpolated contours correspond approxmately to jont confdence levels of 68%, 87%, and 95%. lostc model smlar to the one used for the orane trees data and a frst order open compartment model smlar to the one used for the Theophyllne example. For both models 1000 samples were enerated and maxmum lkelhood (ML) estmates based on the dfferent approxmatons obtaned. For the alternatn approxmaton, restrcted maxmum lkelhood (RML) estmates were also obtaned. 17
18 3.3.1 Lostc Model A lostc model smlar to (13), but wth two random effects nstead of one, was used to enerate the data. The model s ven by y j = 1 + b exp f? [t j? ( + b )] = 3 + " j; = 1; : : :; M; j = 1; : : :; n (16) where the b are..d. N (0; D), and the " j are..d. N (0; ) and ndependent of the b. We used M = 15, n = 10; ; : : :; 15, = 5, = (00; 700; 350) T, and D = 6 4 4?? 5 Table 5 summarzes the smulaton results for the varance-covarance components (MSE denotes the mean square error of the estmators). The dfferent approxmatons to the lolkelhood ve smlar smulaton results for all the parameters nvolved. The cluster specfc varance ( ) s estmated wth more relatve precson than the elements of the scaled varance-covarance matrx of the random effects (D). Ths s probably because the precson of the estmate of (as well as the estmates of ) s determned by the total number of observatons, whle the precson of the estmates of D s determned by the number of clusters. We can also see a tendency for the restrcted maxmum lkelhood to ve postvely based estmates of D 11 and D, whle the other approxmatons ve neatvely based estmates. The ratonale for restrcted maxmum lkelhood s to reduce bas n estmatn varance components. It does not seem to do so n ths case; t just chanes ts drecton. Table 5: Smulaton results for the varance-covarance components n the lostc model D 11 D 1 Approxmaton Mean Bas MSE Mean Bas MSE Alternatn RML Alternatn ML Laplacan Adap. Gaussan D Approxmaton Mean Bas MSE Mean Bas MSE Alternatn RML Alternatn ML Laplacan Adap. Gaussan Fure 5 presents the scatter plots of the varance-covarance component ( and D) estmates for 18
19 the alternatn RML, the alternatn ML, and the Laplacan approxmatons versus the adaptve Gaussan approxmaton. We see that, except for the alternatn RML approxmaton, all methods lead to very smlar estmates. In eneral the alternatn RML approxmaton ves larer values for the estmates of the varance components (especally D 11 and D ) than the other methods. The hher mean square error for D 1 from the alternatn ML and RML methods s vsble n the plot, as each of the panels comparn these estmates to those from the adaptve Gaussan method has a vertcal clump of ponts at the true value. s s s Alternatn - RML Alternatn - ML Laplacan D D D11 Alternatn - RML Alternatn - ML Laplacan D D D1 Alternatn - RML Alternatn - ML Laplacan Alternatn - RML D Alternatn - ML D Laplacan D Fure 5: Scatter plots of varance-covarance components estmates for the alternatn (RML and ML), Laplacan, and adaptve Gaussan approxmatons n the lostc model (16). The dashed lnes ndcate the true values of the parameters. Table 6 presents the smulaton results for the fxed effects estmates. The results are very smlar for all approxmatons consdered. We also note that the relatve varablty of the fxed effects estmates 19
20 s much smaller than those of the estmates of the elements of D. There s very lttle, f any, bas n the fxed effects estmates. Table 6: Smulaton results for the fxed effects n the lostc model 1 3 Approxmaton Mean Bas MSE Mean Bas MSE Mean Bas MSE Alternatn RML Alternatn ML Laplacan Adap. Gaussan Fure 6 presents the scatter plots of the fxed effects estmates for the alternatn RML, alternatn ML, and Laplacan approxmatons versus the adaptve Gaussan approxmaton. Aan we observe a stron areement n the estmates obtaned throuh the varous approxmatons. The alternatn approxmatons tend to ve estmates slhtly smaller than the Laplacan and adaptve Gaussan, but the dfferences are mnor Frst Order Compartment Model The model used n the smulaton s dentcal to (15). As n the Theophyllne example we set M = 1 and n = 11; = 1; : : :; 1. The parameter values used were = 0:5, = (?3:0; 0:5;?:5) T, and D = 6 4 0: Table 7 summarzes the smulaton results for the varance-covarance components estmates. As n the lostc model analyss, we observe that the elements of D are estmated wth less relatve precson than. The alternatn ML, Laplacan, and adaptve Gaussan approxmatons seem to lead to slhtly downward based estmates of D 11 and D, whle the alternatn RML approxmaton appears to ve unbased estmates (thus achevn ts man purpose). Note however that the unbasedeness of the RML estmates does not translate nto smaller mean square error all four estmaton methods lead to smlar MSE, for all parameters. Fure 7 presents the scatter plots of the varance-covarance estmates for the alternatn RML, alternatn ML, and Laplacan approxmatons versus the adaptve Gaussan approxmaton. The alternatn RML approxmaton tends to ve larer values for D 11 and D, and larer absolute values for D 1, whle the remann approxmatons lead to very smlar estmates. There was one sample for whch the alternatn approxmatons apparently convered to a dfferent soluton than the Laplacan and 0
21 Alternatn - RML Alternatn - RML beta beta Alternatn - ML Alternatn - ML beta beta Laplacan Laplacan beta beta beta beta beta3 Alternatn - RML Alternatn - ML Laplacan Fure 6: Scatter plots of fxed effects estmates for the alternatn (RML and ML), Laplacan, and adaptve Gaussan approxmatons n the lostc model (16). The dashed lnes ndcate the true values of the parameters. adaptve Gaussan. Overall there were no major dfferences between the approxmatons n estmatn the varance-covarance components. Table 6 ves the smulaton results for the fxed effects estmates. All four approxmatons ve vrtually dentcal results for the estmaton of the fxed effects. They all show very lttle bas and smaller relatve varablty when compared to the estmates of the varance-covarance components. The scatter plots of the fxed effects estmates, not ncluded here, show practcally dentcal results for the alternatn RML and ML, the Laplacan, and the adaptve Gaussan approxmatons. 1
22 Table 7: Smulaton results for the varance-covarance components n the frst order compartment model D 11 D 1 Approxmaton Mean Bas MSE Mean Bas MSE Alternatn RML Alternatn ML Laplacan Adap. Gaussan D Approxmaton Mean Bas MSE Mean Bas MSE Alternatn RML Alternatn ML Laplacan Adap. Gaussan Table 8: Smulaton results for the fxed effects n the frst order compartment model 1 3 Approxmaton Mean Bas MSE Mean Bas MSE Mean Bas MSE Alternatn RML Alternatn ML Laplacan Adap. Gaussan Conclusons The results of secton 3 ndcate that the alternatn approxmaton (4) to the lolkelhood functon n the nonlnear mxed effects model (1) proposed by Lndstrom and Bates (1990) ves accurate and relable estmaton results. The man advantaes of ths approxmaton are ts computatonal effcency (allown the use of lnear mxed effects technques to estmate the scaled varance-covarance matrx of the random effects D) and the avalablty of a restrcted lkelhood verson of t, whch s not yet defned for other approxmatons/estmaton methods. Wth reard to the restrcted maxmum lkelhood estmaton thouh, the results of secton 3 suest that the bas correcton ablty of ths method depends on the nonlnear model that s ben consdered: RML estmaton acheved ts purpose for the frst order compartment model (15), but t ncreased the bas n the lostc model (16). More research s needed n ths area. Snce t s smpler computatonally the alternatn approxmaton should be used to provde startn values for the more accurate approxmatons (e.. Laplacan and adaptve Gaussan) f they are preferred. The Gaussan quadrature approxmaton (11) only seems to ve accurate results for lare number of
23 s s s Alternatn - RML D11 Alternatn - ML D11 Laplacan D11 Alternatn - RML Alternatn - ML Laplacan D1 D1 D1 Alternatn - RML Alternatn - ML Laplacan D D D Alternatn - RML Alternatn - ML Laplacan Fure 7: Scatter plots of varance-covarance components estmates for the alternatn (RML and ML), Laplacan, and adaptve Gaussan approxmatons n the frst order compartment model (15). The dashed lnes ndcate the true values of the parameters. abscssa (> 100), what makes t very neffcent computatonally. The basc problem s that t centers the rd of abscssas at 0 (the expected value of the random effects) and scales t accordn to D, whle the hhest values of the nterand n () are concentrated around the posteror modes of the random effects (b b) and scaled accordn to 00 ; D; y; b b. The advantaes of ths approxmaton are that t does not requre the estmaton of the posteror modes of the random effects at each teraton and t admts closed form partal dervatves wth respect to the parameters of nterest (; D; and ), provded these are avalable for the model functon f (Davdan and Gallant, 199). We feel that these advantaes do not compensate for the naccuracy or computatonal neffcency of the Gaussan approxmaton. The mportance sampln approxmaton (9) ves relable estmaton results, comparable to those of 3
24 the adaptve Gaussan and Laplacan approxmatons, but s consderably less effcent computatonally than these approxmatons. Also, the stochastc varablty assocated wth the dfferent mportance samples may overwhelm the numercal varablty of the lolkelhood for small chanes n the parameter values, makn t dffcult to calculate numercal dervatves. The man advantae of the mportance sampln approxmaton s ts versatlty n handln dstrbutons other than the normal, for both the random effects and the cluster-specfc error term (). For example t would be rather strahtforward to adapt the mportance sampln nteraton to handle a multvarate t dstrbuton for the random effects, but that would not be a trval task for ether the alternatn, the Laplacan, or the adaptve Gaussan approxmatons. Wakefeld et al. (1994) use the smlar property of Gbbs sampler methods to check for outlers n nonlnear mxed effects models. If one s wlln to stck wth the normal dstrbuton for b and n the nonlnear mxed effects model (1) then the mportance sampln approxmaton s not the most effcent choce. Of all approxmatons consdered here, the Laplacan and adaptve Gaussan approxmatons probably ve the best mx of effcency and accuracy. The former can be rearded as a partcular case of the latter, where just one abscssa s used. Both approxmatons (and the mportance sampln approxmaton as well) ve the exact lolkelhood when the model functon f n (1) s a lnear functon of the random effects. In the examples that we analyzed not much was aned by on from a one-pont adaptve Gaussan quadrature (Laplacan) approxmaton to approxmatons wth a larer number of abscssas. It appears that the major an n adaptve Gaussan approxmatons s related to the centern and scaln of the abscssas. Increasn the number of ponts n the evaluaton rd only ves marnal mprovement. The Laplacan approxmaton has the addtonal advantae over the adaptve Gaussan approxmaton wth more than one abscssa of allown profln of the lolkelhood over, thus reducn the dmensonalty of the optmzaton problem. For statstcal analyss purpose we would recommend usn a hybrd scheme n whch the alternatn alorthm would be used to et ood ntal values for the more refned Laplacan approxmaton to the lolkelhood of model (1). Ths way the computatonal effcency of the alternatn alorthm would be combned wth the reater accuracy of the Laplacan approxmaton. Acknowledment Ths research was partally supported by Coordenação de Aperfeçoamento de Pessoal de Nível Superor, Brazl and NSF Grant DMS
25 References Abramowtz, M. and Steun, I. A. (1964). Handbook of Mathematcal Functons wth Formulas, Graphs, and mathematcal Tables, Dover, New York. Bates, D. M. and Watts, D. G. (1980). Relatve curvature measures of nonlnearty, Journal of the Royal Statstcal Socety, Ser. B 4: 1 5. Bates, D. M. and Watts, D. G. (1988). Nonlnear Reresson Analyss and Its Applcatons, Wley, New York. Chambers, J. M. and Haste, T. J. (eds) (199). Statstcal Models n S, Wadsworth, Belmont, CA. Davdan, M. and Gallant, A. R. (199). Smooth nonparametrc maxmum lkelhood estmaton for populaton pharmacoknetcs, wth applcaton to qundne, Journal of Pharmacoknetcs and Bopharmaceutcs 0: Davs, P. J. and Rabnowtz, P. (1984). Methods of Numercal Interaton, second edn, Academc Press, New York. Draper, N. R. and Smth, H. (1981). Appled Reresson Analyss, nd edn, Wley, New York. Geweke, J. (1989). Bayesan nference n econometrc models usn Monte Carlo nteraton, Econometrca 57: Golub, G. H. (1973). Some modfed matrx eenvalue problems, SIAM Revew 15: Golub, G. H. and Welsch, J. H. (1969). Calculaton of Gaussan quadrature rules, Math. Comp. 3: Harvlle, D. A. (1974). Bayesan nference for varance components usn only error contrasts, Bometrka 61: Lard, N. M. and Ware, J. H. (198). Random-effects models for lontudnal data, Bometrcs 38: Leonard, T., Hsu, J. S. J. and Tsu, K. W. (1989). Bayesan marnal nference, Journal of the Amercan Statstcal Assocaton 84:
26 Lndstrom, M. J. and Bates, D. M. (1988). Newton-Raphson and EM alorthms for lnear mxed-effects models for repeated-measures data, Journal of the Amercan Statstcal Assocaton 83: Lndstrom, M. J. and Bates, D. M. (1990). Nonlnear mxed effects models for repeated measures data, Bometrcs 46: Mallet, A., Mentre, F., Stemer, J.-L. and Lokek, F. (1988). Nonparametrc maxmum lkelhood estmaton for populaton pharmacoknetcs, wth applcatons to Cyclosporne, J. Pharmacokn. Bopharm. 16: Pnhero, J. C. and Bates, D. M. (1993). Asymptotc propertes of maxmum leklhood estmates n the eneral lnear mxed effects model. Submtted to Annals of Statstcs. Shener, L. B. and Beal, S. L. (1980). Evaluaton of methods for estmatn populaton pharmacoknetc parameters. I. mchaels-menten model: Routne clncal pharmacoknetc data, Journal of Pharmacoknetcs and Bopharmaceutcs 8(6): Terney, L. and Kadane, J. B. (1986). Accurate approxmatons for posteror moments and denstes, Journal of the Amercan Statstcal Assocaton 81(393): Vonesh, E. F. and Carter, R. L. (199). Mxed-effects nonlnear reresson for unbalanced repeated measures, Bometrcs 48: Wakefeld, J. C., Smth, A. F. M., Racne-Poon, A. and Gelfand, A. E. (1994). Bayesan analyss of lnear and nonlnear populaton models usn the Gbbs sampler, Appled Statstcs. Accepted for publcaton. 6
Population Design in Nonlinear Mixed Effects Multiple Response Models: extension of PFIM and evaluation by simulation with NONMEM and MONOLIX
Populaton Desgn n Nonlnear Mxed Effects Multple Response Models: extenson of PFIM and evaluaton by smulaton wth NONMEM and MONOLIX May 4th 007 Carolne Bazzol, Sylve Retout, France Mentré Inserm U738 Unversty
More informationFlux-Uncertainty from Aperture Photometry. F. Masci, version 1.0, 10/14/2008
Flux-Uncertanty from Aperture Photometry F. Masc, verson 1.0, 10/14/008 1. Summary We derve a eneral formula for the nose varance n the flux of a source estmated from aperture photometry. The 1-σ uncertanty
More informationOn an Extension of Stochastic Approximation EM Algorithm for Incomplete Data Problems. Vahid Tadayon 1
On an Extenson of Stochastc Approxmaton EM Algorthm for Incomplete Data Problems Vahd Tadayon Abstract: The Stochastc Approxmaton EM (SAEM algorthm, a varant stochastc approxmaton of EM, s a versatle tool
More informationParametric fractional imputation for missing data analysis. Jae Kwang Kim Survey Working Group Seminar March 29, 2010
Parametrc fractonal mputaton for mssng data analyss Jae Kwang Km Survey Workng Group Semnar March 29, 2010 1 Outlne Introducton Proposed method Fractonal mputaton Approxmaton Varance estmaton Multple mputaton
More informationGlobal Sensitivity. Tuesday 20 th February, 2018
Global Senstvty Tuesday 2 th February, 28 ) Local Senstvty Most senstvty analyses [] are based on local estmates of senstvty, typcally by expandng the response n a Taylor seres about some specfc values
More informationLecture Notes on Linear Regression
Lecture Notes on Lnear Regresson Feng L fl@sdueducn Shandong Unversty, Chna Lnear Regresson Problem In regresson problem, we am at predct a contnuous target value gven an nput feature vector We assume
More informationBayesian decision theory. Nuno Vasconcelos ECE Department, UCSD
Bayesan decson theory Nuno Vasconcelos ECE Department UCSD Notaton the notaton n DHS s qute sloppy e.. show that error error z z dz really not clear what ths means we wll use the follown notaton subscrpts
More informationx = , so that calculated
Stat 4, secton Sngle Factor ANOVA notes by Tm Plachowsk n chapter 8 we conducted hypothess tests n whch we compared a sngle sample s mean or proporton to some hypotheszed value Chapter 9 expanded ths to
More informationComposite Hypotheses testing
Composte ypotheses testng In many hypothess testng problems there are many possble dstrbutons that can occur under each of the hypotheses. The output of the source s a set of parameters (ponts n a parameter
More informationBayesian decision theory. Nuno Vasconcelos ECE Department, UCSD
Bayesan decson theory Nuno Vasconcelos ECE Department, UCSD Bayesan decson theory recall that we have state of the world observatons decson functon L[,y] loss of predctn y wth the epected value of the
More informationLINEAR REGRESSION ANALYSIS. MODULE IX Lecture Multicollinearity
LINEAR REGRESSION ANALYSIS MODULE IX Lecture - 30 Multcollnearty Dr. Shalabh Department of Mathematcs and Statstcs Indan Insttute of Technology Kanpur 2 Remedes for multcollnearty Varous technques have
More informationInner Product. Euclidean Space. Orthonormal Basis. Orthogonal
Inner Product Defnton 1 () A Eucldean space s a fnte-dmensonal vector space over the reals R, wth an nner product,. Defnton 2 (Inner Product) An nner product, on a real vector space X s a symmetrc, blnear,
More informationModule 3 LOSSY IMAGE COMPRESSION SYSTEMS. Version 2 ECE IIT, Kharagpur
Module 3 LOSSY IMAGE COMPRESSION SYSTEMS Verson ECE IIT, Kharagpur Lesson 6 Theory of Quantzaton Verson ECE IIT, Kharagpur Instructonal Objectves At the end of ths lesson, the students should be able to:
More informationEcon107 Applied Econometrics Topic 3: Classical Model (Studenmund, Chapter 4)
I. Classcal Assumptons Econ7 Appled Econometrcs Topc 3: Classcal Model (Studenmund, Chapter 4) We have defned OLS and studed some algebrac propertes of OLS. In ths topc we wll study statstcal propertes
More informationLinear Regression Analysis: Terminology and Notation
ECON 35* -- Secton : Basc Concepts of Regresson Analyss (Page ) Lnear Regresson Analyss: Termnology and Notaton Consder the generc verson of the smple (two-varable) lnear regresson model. It s represented
More informationNumerical Heat and Mass Transfer
Master degree n Mechancal Engneerng Numercal Heat and Mass Transfer 06-Fnte-Dfference Method (One-dmensonal, steady state heat conducton) Fausto Arpno f.arpno@uncas.t Introducton Why we use models and
More informationANSWERS. Problem 1. and the moment generating function (mgf) by. defined for any real t. Use this to show that E( U) var( U)
Econ 413 Exam 13 H ANSWERS Settet er nndelt 9 deloppgaver, A,B,C, som alle anbefales å telle lkt for å gøre det ltt lettere å stå. Svar er gtt . Unfortunately, there s a prntng error n the hnt of
More informationComputation of Higher Order Moments from Two Multinomial Overdispersion Likelihood Models
Computaton of Hgher Order Moments from Two Multnomal Overdsperson Lkelhood Models BY J. T. NEWCOMER, N. K. NEERCHAL Department of Mathematcs and Statstcs, Unversty of Maryland, Baltmore County, Baltmore,
More information2016 Wiley. Study Session 2: Ethical and Professional Standards Application
6 Wley Study Sesson : Ethcal and Professonal Standards Applcaton LESSON : CORRECTION ANALYSIS Readng 9: Correlaton and Regresson LOS 9a: Calculate and nterpret a sample covarance and a sample correlaton
More informationLimited Dependent Variables and Panel Data. Tibor Hanappi
Lmted Dependent Varables and Panel Data Tbor Hanapp 30.06.2010 Lmted Dependent Varables Dscrete: Varables that can take onl a countable number of values Censored/Truncated: Data ponts n some specfc range
More informationThe Multiple Classical Linear Regression Model (CLRM): Specification and Assumptions. 1. Introduction
ECONOMICS 5* -- NOTE (Summary) ECON 5* -- NOTE The Multple Classcal Lnear Regresson Model (CLRM): Specfcaton and Assumptons. Introducton CLRM stands for the Classcal Lnear Regresson Model. The CLRM s also
More informationA Comparative Study for Estimation Parameters in Panel Data Model
A Comparatve Study for Estmaton Parameters n Panel Data Model Ahmed H. Youssef and Mohamed R. Abonazel hs paper examnes the panel data models when the regresson coeffcents are fxed random and mxed and
More informationRESIDUALS AND INFLUENCE IN NONLINEAR REGRESSION FOR REPEATED MEASUREMENT DATA
Operatons Research and Applcatons : An Internatonal Journal (ORAJ), Vol.4, No.3/4, November 17 RESIDUALS AND INFLUENCE IN NONLINEAR REGRESSION FOR REPEAED MEASUREMEN DAA Munsr Al, Yu Feng, Al choo, Zamr
More informationLinear Approximation with Regularization and Moving Least Squares
Lnear Approxmaton wth Regularzaton and Movng Least Squares Igor Grešovn May 007 Revson 4.6 (Revson : March 004). 5 4 3 0.5 3 3.5 4 Contents: Lnear Fttng...4. Weghted Least Squares n Functon Approxmaton...
More informationChapter 5 Multilevel Models
Chapter 5 Multlevel Models 5.1 Cross-sectonal multlevel models 5.1.1 Two-level models 5.1.2 Multple level models 5.1.3 Multple level modelng n other felds 5.2 Longtudnal multlevel models 5.2.1 Two-level
More informationChapter 9: Statistical Inference and the Relationship between Two Variables
Chapter 9: Statstcal Inference and the Relatonshp between Two Varables Key Words The Regresson Model The Sample Regresson Equaton The Pearson Correlaton Coeffcent Learnng Outcomes After studyng ths chapter,
More informationFEATURE ANALYSIS ON QUEUE LENGTH OF ASYMMETRIC TWO-QUEUE POLLING SYSTEM WITH GATED SERVICES *
Journal of Theoretcal and Appled Informaton Technoloy 0 th January 03. Vol. 4 No. 005-03 JATIT & LLS. All rhts reserved. ISSN: 99-845 www.att.or E-ISSN: 8-395 FEATURE ANALYSIS ON QUEUE LENGTH OF ASYMMETRIC
More informationOn Outlier Robust Small Area Mean Estimate Based on Prediction of Empirical Distribution Function
On Outler Robust Small Area Mean Estmate Based on Predcton of Emprcal Dstrbuton Functon Payam Mokhtaran Natonal Insttute of Appled Statstcs Research Australa Unversty of Wollongong Small Area Estmaton
More informationMarkov Chain Monte Carlo (MCMC), Gibbs Sampling, Metropolis Algorithms, and Simulated Annealing Bioinformatics Course Supplement
Markov Chan Monte Carlo MCMC, Gbbs Samplng, Metropols Algorthms, and Smulated Annealng 2001 Bonformatcs Course Supplement SNU Bontellgence Lab http://bsnuackr/ Outlne! Markov Chan Monte Carlo MCMC! Metropols-Hastngs
More informationStatistics for Economics & Business
Statstcs for Economcs & Busness Smple Lnear Regresson Learnng Objectves In ths chapter, you learn: How to use regresson analyss to predct the value of a dependent varable based on an ndependent varable
More informationChapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems
Numercal Analyss by Dr. Anta Pal Assstant Professor Department of Mathematcs Natonal Insttute of Technology Durgapur Durgapur-713209 emal: anta.bue@gmal.com 1 . Chapter 5 Soluton of System of Lnear Equatons
More informationMACHINE APPLIED MACHINE LEARNING LEARNING. Gaussian Mixture Regression
11 MACHINE APPLIED MACHINE LEARNING LEARNING MACHINE LEARNING Gaussan Mture Regresson 22 MACHINE APPLIED MACHINE LEARNING LEARNING Bref summary of last week s lecture 33 MACHINE APPLIED MACHINE LEARNING
More information3.1 Expectation of Functions of Several Random Variables. )' be a k-dimensional discrete or continuous random vector, with joint PMF p (, E X E X1 E X
Statstcs 1: Probablty Theory II 37 3 EPECTATION OF SEVERAL RANDOM VARIABLES As n Probablty Theory I, the nterest n most stuatons les not on the actual dstrbuton of a random vector, but rather on a number
More informationLOW BIAS INTEGRATED PATH ESTIMATORS. James M. Calvin
Proceedngs of the 007 Wnter Smulaton Conference S G Henderson, B Bller, M-H Hseh, J Shortle, J D Tew, and R R Barton, eds LOW BIAS INTEGRATED PATH ESTIMATORS James M Calvn Department of Computer Scence
More informationEfficient Algorithms for Robust Estimation in Linear Mixed-Effects Models Using the Multivariate t-distribution
Effcent Algorthms for Robust Estmaton n Lnear Mxed-Effects Models Usng the Multvarate t-dstrbuton José C. Pnhero, Chuanha Lu Bell Laboratores Lucent Technologes Murray Hll, NJ 07974 and Yngnan Wu Department
More information28. SIMPLE LINEAR REGRESSION III
8. SIMPLE LINEAR REGRESSION III Ftted Values and Resduals US Domestc Beers: Calores vs. % Alcohol To each observed x, there corresponds a y-value on the ftted lne, y ˆ = βˆ + βˆ x. The are called ftted
More informationImproved conditional approximations of the population Fisher Information Matrix
Improved condtonal approxmatons of the populaton Fsher Informaton Matrx Joakm Nyberg, Sebastan Ueckert and Andrew C. Hooker Pharmacometrcs research group Department of Pharmaceutcal Boscences Uppsala Unversty
More informationChapter 11: Simple Linear Regression and Correlation
Chapter 11: Smple Lnear Regresson and Correlaton 11-1 Emprcal Models 11-2 Smple Lnear Regresson 11-3 Propertes of the Least Squares Estmators 11-4 Hypothess Test n Smple Lnear Regresson 11-4.1 Use of t-tests
More informationLINEAR REGRESSION ANALYSIS. MODULE IX Lecture Multicollinearity
LINEAR REGRESSION ANALYSIS MODULE IX Lecture - 31 Multcollnearty Dr. Shalabh Department of Mathematcs and Statstcs Indan Insttute of Technology Kanpur 6. Rdge regresson The OLSE s the best lnear unbased
More information1. Inference on Regression Parameters a. Finding Mean, s.d and covariance amongst estimates. 2. Confidence Intervals and Working Hotelling Bands
Content. Inference on Regresson Parameters a. Fndng Mean, s.d and covarance amongst estmates.. Confdence Intervals and Workng Hotellng Bands 3. Cochran s Theorem 4. General Lnear Testng 5. Measures of
More informationNETWORK IDENTIFICATION AND TIME SERIES ANALYIS USING S-SYSTEMS
NETWORK IDENTIFICATION AND TIME SERIES ANALYIS USING S-SYSTEMS SIREN RØST VEFLINGSTAD Norwean Unversty of Lfe Scences Research done n collaboraton wth JONAS ALMEIDA and EBERHARD VOIT durn a vst at the
More informationPinheiro, J. C. and Bates, D. M. (1995). Mixed-eects methods and classes for s and splus, Technical
Pnhero, J. C. and Bates, D. M. (1995). Mxed-eects methods and classes for s and splus, Techncal Report 89, Department of Bostatstcs, Unversty of Wsconsn { Madson. Pottho, R. F. and Roy, S. N. (1964). A
More informationSome Comments on Accelerating Convergence of Iterative Sequences Using Direct Inversion of the Iterative Subspace (DIIS)
Some Comments on Acceleratng Convergence of Iteratve Sequences Usng Drect Inverson of the Iteratve Subspace (DIIS) C. Davd Sherrll School of Chemstry and Bochemstry Georga Insttute of Technology May 1998
More information18. SIMPLE LINEAR REGRESSION III
8. SIMPLE LINEAR REGRESSION III US Domestc Beers: Calores vs. % Alcohol Ftted Values and Resduals To each observed x, there corresponds a y-value on the ftted lne, y ˆ ˆ = α + x. The are called ftted values.
More informationSingular Value Decomposition: Theory and Applications
Sngular Value Decomposton: Theory and Applcatons Danel Khashab Sprng 2015 Last Update: March 2, 2015 1 Introducton A = UDV where columns of U and V are orthonormal and matrx D s dagonal wth postve real
More informationDepartment of Quantitative Methods & Information Systems. Time Series and Their Components QMIS 320. Chapter 6
Department of Quanttatve Methods & Informaton Systems Tme Seres and Ther Components QMIS 30 Chapter 6 Fall 00 Dr. Mohammad Zanal These sldes were modfed from ther orgnal source for educatonal purpose only.
More informationLossy Compression. Compromise accuracy of reconstruction for increased compression.
Lossy Compresson Compromse accuracy of reconstructon for ncreased compresson. The reconstructon s usually vsbly ndstngushable from the orgnal mage. Typcally, one can get up to 0:1 compresson wth almost
More informationLecture 21: Numerical methods for pricing American type derivatives
Lecture 21: Numercal methods for prcng Amercan type dervatves Xaoguang Wang STAT 598W Aprl 10th, 2014 (STAT 598W) Lecture 21 1 / 26 Outlne 1 Fnte Dfference Method Explct Method Penalty Method (STAT 598W)
More informationNUMERICAL DIFFERENTIATION
NUMERICAL DIFFERENTIATION 1 Introducton Dfferentaton s a method to compute the rate at whch a dependent output y changes wth respect to the change n the ndependent nput x. Ths rate of change s called the
More informationwhere v means the change in velocity, and t is the
1 PHYS:100 LECTURE 4 MECHANICS (3) Ths lecture covers the eneral case of moton wth constant acceleraton and free fall (whch s one of the more mportant examples of moton wth constant acceleraton) n a more
More informationLecture 3 Stat102, Spring 2007
Lecture 3 Stat0, Sprng 007 Chapter 3. 3.: Introducton to regresson analyss Lnear regresson as a descrptve technque The least-squares equatons Chapter 3.3 Samplng dstrbuton of b 0, b. Contnued n net lecture
More informationDr. Shalabh Department of Mathematics and Statistics Indian Institute of Technology Kanpur
Analyss of Varance and Desgn of Experment-I MODULE VIII LECTURE - 34 ANALYSIS OF VARIANCE IN RANDOM-EFFECTS MODEL AND MIXED-EFFECTS EFFECTS MODEL Dr Shalabh Department of Mathematcs and Statstcs Indan
More informationA Hybrid Variational Iteration Method for Blasius Equation
Avalable at http://pvamu.edu/aam Appl. Appl. Math. ISSN: 1932-9466 Vol. 10, Issue 1 (June 2015), pp. 223-229 Applcatons and Appled Mathematcs: An Internatonal Journal (AAM) A Hybrd Varatonal Iteraton Method
More informationPhase I Monitoring of Nonlinear Profiles
Phase I Montorng of Nonlnear Profles James D. Wllams Wllam H. Woodall Jeffrey B. Brch May, 003 J.D. Wllams, Bll Woodall, Jeff Brch, Vrgna Tech 003 Qualty & Productvty Research Conference, Yorktown Heghts,
More informationCOMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS
Avalable onlne at http://sck.org J. Math. Comput. Sc. 3 (3), No., 6-3 ISSN: 97-537 COMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS
More informationStatistical Evaluation of WATFLOOD
tatstcal Evaluaton of WATFLD By: Angela MacLean, Dept. of Cvl & Envronmental Engneerng, Unversty of Waterloo, n. ctober, 005 The statstcs program assocated wth WATFLD uses spl.csv fle that s produced wth
More informationDiscussion of Extensions of the Gauss-Markov Theorem to the Case of Stochastic Regression Coefficients Ed Stanek
Dscusson of Extensons of the Gauss-arkov Theorem to the Case of Stochastc Regresson Coeffcents Ed Stanek Introducton Pfeffermann (984 dscusses extensons to the Gauss-arkov Theorem n settngs where regresson
More informationEstimation: Part 2. Chapter GREG estimation
Chapter 9 Estmaton: Part 2 9. GREG estmaton In Chapter 8, we have seen that the regresson estmator s an effcent estmator when there s a lnear relatonshp between y and x. In ths chapter, we generalzed the
More informationPsychology 282 Lecture #24 Outline Regression Diagnostics: Outliers
Psychology 282 Lecture #24 Outlne Regresson Dagnostcs: Outlers In an earler lecture we studed the statstcal assumptons underlyng the regresson model, ncludng the followng ponts: Formal statement of assumptons.
More informationChapter 13: Multiple Regression
Chapter 13: Multple Regresson 13.1 Developng the multple-regresson Model The general model can be descrbed as: It smplfes for two ndependent varables: The sample ft parameter b 0, b 1, and b are used to
More informationComparison of Regression Lines
STATGRAPHICS Rev. 9/13/2013 Comparson of Regresson Lnes Summary... 1 Data Input... 3 Analyss Summary... 4 Plot of Ftted Model... 6 Condtonal Sums of Squares... 6 Analyss Optons... 7 Forecasts... 8 Confdence
More informationImplicit Integration Henyey Method
Implct Integraton Henyey Method In realstc stellar evoluton codes nstead of a drect ntegraton usng for example the Runge-Kutta method one employs an teratve mplct technque. Ths s because the structure
More informationStatistics Chapter 4
Statstcs Chapter 4 "There are three knds of les: les, damned les, and statstcs." Benjamn Dsrael, 1895 (Brtsh statesman) Gaussan Dstrbuton, 4-1 If a measurement s repeated many tmes a statstcal treatment
More informationSTAT 511 FINAL EXAM NAME Spring 2001
STAT 5 FINAL EXAM NAME Sprng Instructons: Ths s a closed book exam. No notes or books are allowed. ou may use a calculator but you are not allowed to store notes or formulas n the calculator. Please wrte
More informationIII. Econometric Methodology Regression Analysis
Page Econ07 Appled Econometrcs Topc : An Overvew of Regresson Analyss (Studenmund, Chapter ) I. The Nature and Scope of Econometrcs. Lot s of defntons of econometrcs. Nobel Prze Commttee Paul Samuelson,
More informationANOMALIES OF THE MAGNITUDE OF THE BIAS OF THE MAXIMUM LIKELIHOOD ESTIMATOR OF THE REGRESSION SLOPE
P a g e ANOMALIES OF THE MAGNITUDE OF THE BIAS OF THE MAXIMUM LIKELIHOOD ESTIMATOR OF THE REGRESSION SLOPE Darmud O Drscoll ¹, Donald E. Ramrez ² ¹ Head of Department of Mathematcs and Computer Studes
More informationMarkov Chain Monte Carlo Lecture 6
where (x 1,..., x N ) X N, N s called the populaton sze, f(x) f (x) for at least one {1, 2,..., N}, and those dfferent from f(x) are called the tral dstrbutons n terms of mportance samplng. Dfferent ways
More informationSee Book Chapter 11 2 nd Edition (Chapter 10 1 st Edition)
Count Data Models See Book Chapter 11 2 nd Edton (Chapter 10 1 st Edton) Count data consst of non-negatve nteger values Examples: number of drver route changes per week, the number of trp departure changes
More informationLab 4: Two-level Random Intercept Model
BIO 656 Lab4 009 Lab 4: Two-level Random Intercept Model Data: Peak expratory flow rate (pefr) measured twce, usng two dfferent nstruments, for 17 subjects. (from Chapter 1 of Multlevel and Longtudnal
More informationStatistics for Managers Using Microsoft Excel/SPSS Chapter 13 The Simple Linear Regression Model and Correlation
Statstcs for Managers Usng Mcrosoft Excel/SPSS Chapter 13 The Smple Lnear Regresson Model and Correlaton 1999 Prentce-Hall, Inc. Chap. 13-1 Chapter Topcs Types of Regresson Models Determnng the Smple Lnear
More informationS Advanced Digital Communication (4 cr) Targets today
S.72-3320 Advanced Dtal Communcaton (4 cr) Convolutonal Codes Tarets today Why to apply convolutonal codn? Defnn convolutonal codes Practcal encodn crcuts Defnn qualty of convolutonal codes Decodn prncples
More informationThe Ordinary Least Squares (OLS) Estimator
The Ordnary Least Squares (OLS) Estmator 1 Regresson Analyss Regresson Analyss: a statstcal technque for nvestgatng and modelng the relatonshp between varables. Applcatons: Engneerng, the physcal and chemcal
More informationTime-Varying Coefficient Model with Linear Smoothing Function for Longitudinal Data in Clinical Trial
Tme-Varyng Coeffcent Model wth Lnear Smoothng Functon for Longtudnal Data n Clncal Tral Masanor Ito, Toshhro Msum and Hdek Hrooka Bostatstcs Group, Data Scence Dept., Astellas Pharma Inc. Introducton In
More informationTesting for seasonal unit roots in heterogeneous panels
Testng for seasonal unt roots n heterogeneous panels Jesus Otero * Facultad de Economía Unversdad del Rosaro, Colomba Jeremy Smth Department of Economcs Unversty of arwck Monca Gulett Aston Busness School
More informationECONOMICS 351*-A Mid-Term Exam -- Fall Term 2000 Page 1 of 13 pages. QUEEN'S UNIVERSITY AT KINGSTON Department of Economics
ECOOMICS 35*-A Md-Term Exam -- Fall Term 000 Page of 3 pages QUEE'S UIVERSITY AT KIGSTO Department of Economcs ECOOMICS 35* - Secton A Introductory Econometrcs Fall Term 000 MID-TERM EAM ASWERS MG Abbott
More informationx i1 =1 for all i (the constant ).
Chapter 5 The Multple Regresson Model Consder an economc model where the dependent varable s a functon of K explanatory varables. The economc model has the form: y = f ( x,x,..., ) xk Approxmate ths by
More informationSIMPLE LINEAR REGRESSION
Smple Lnear Regresson and Correlaton Introducton Prevousl, our attenton has been focused on one varable whch we desgnated b x. Frequentl, t s desrable to learn somethng about the relatonshp between two
More information4 Analysis of Variance (ANOVA) 5 ANOVA. 5.1 Introduction. 5.2 Fixed Effects ANOVA
4 Analyss of Varance (ANOVA) 5 ANOVA 51 Introducton ANOVA ANOVA s a way to estmate and test the means of multple populatons We wll start wth one-way ANOVA If the populatons ncluded n the study are selected
More informationAn (almost) unbiased estimator for the S-Gini index
An (almost unbased estmator for the S-Gn ndex Thomas Demuynck February 25, 2009 Abstract Ths note provdes an unbased estmator for the absolute S-Gn and an almost unbased estmator for the relatve S-Gn for
More informationBOOTSTRAP METHOD FOR TESTING OF EQUALITY OF SEVERAL MEANS. M. Krishna Reddy, B. Naveen Kumar and Y. Ramu
BOOTSTRAP METHOD FOR TESTING OF EQUALITY OF SEVERAL MEANS M. Krshna Reddy, B. Naveen Kumar and Y. Ramu Department of Statstcs, Osmana Unversty, Hyderabad -500 007, Inda. nanbyrozu@gmal.com, ramu0@gmal.com
More informationThe Order Relation and Trace Inequalities for. Hermitian Operators
Internatonal Mathematcal Forum, Vol 3, 08, no, 507-57 HIKARI Ltd, wwwm-hkarcom https://doorg/0988/mf088055 The Order Relaton and Trace Inequaltes for Hermtan Operators Y Huang School of Informaton Scence
More informationComputing MLE Bias Empirically
Computng MLE Bas Emprcally Kar Wa Lm Australan atonal Unversty January 3, 27 Abstract Ths note studes the bas arses from the MLE estmate of the rate parameter and the mean parameter of an exponental dstrbuton.
More informationLecture 6: Introduction to Linear Regression
Lecture 6: Introducton to Lnear Regresson An Manchakul amancha@jhsph.edu 24 Aprl 27 Lnear regresson: man dea Lnear regresson can be used to study an outcome as a lnear functon of a predctor Example: 6
More informationLecture 3. Ax x i a i. i i
18.409 The Behavor of Algorthms n Practce 2/14/2 Lecturer: Dan Spelman Lecture 3 Scrbe: Arvnd Sankar 1 Largest sngular value In order to bound the condton number, we need an upper bound on the largest
More informationEconomics 130. Lecture 4 Simple Linear Regression Continued
Economcs 130 Lecture 4 Contnued Readngs for Week 4 Text, Chapter and 3. We contnue wth addressng our second ssue + add n how we evaluate these relatonshps: Where do we get data to do ths analyss? How do
More informationDr. Shalabh Department of Mathematics and Statistics Indian Institute of Technology Kanpur
Analyss of Varance and Desgn of Experment-I MODULE VII LECTURE - 3 ANALYSIS OF COVARIANCE Dr Shalabh Department of Mathematcs and Statstcs Indan Insttute of Technology Kanpur Any scentfc experment s performed
More informationParametric fractional imputation for missing data analysis
Secton on Survey Research Methods JSM 2008 Parametrc fractonal mputaton for mssng data analyss Jae Kwang Km Wayne Fuller Abstract Under a parametrc model for mssng data, the EM algorthm s a popular tool
More information2E Pattern Recognition Solutions to Introduction to Pattern Recognition, Chapter 2: Bayesian pattern classification
E395 - Pattern Recognton Solutons to Introducton to Pattern Recognton, Chapter : Bayesan pattern classfcaton Preface Ths document s a soluton manual for selected exercses from Introducton to Pattern Recognton
More informationLecture 16 Statistical Analysis in Biomaterials Research (Part II)
3.051J/0.340J 1 Lecture 16 Statstcal Analyss n Bomaterals Research (Part II) C. F Dstrbuton Allows comparson of varablty of behavor between populatons usng test of hypothess: σ x = σ x amed for Brtsh statstcan
More information( ) ( ) ( ) ( ) STOCHASTIC SIMULATION FOR BLOCKED DATA. Monte Carlo simulation Rejection sampling Importance sampling Markov chain Monte Carlo
SOCHASIC SIMULAIO FOR BLOCKED DAA Stochastc System Analyss and Bayesan Model Updatng Monte Carlo smulaton Rejecton samplng Importance samplng Markov chan Monte Carlo Monte Carlo smulaton Introducton: If
More informationA Note on Test of Homogeneity Against Umbrella Scale Alternative Based on U-Statistics
J Stat Appl Pro No 3 93- () 93 NSP Journal of Statstcs Applcatons & Probablty --- An Internatonal Journal @ NSP Natural Scences Publshng Cor A Note on Test of Homogenety Aganst Umbrella Scale Alternatve
More informationOn the Influential Points in the Functional Circular Relationship Models
On the Influental Ponts n the Functonal Crcular Relatonshp Models Department of Mathematcs, Faculty of Scence Al-Azhar Unversty-Gaza, Gaza, Palestne alzad33@yahoo.com Abstract If the nterest s to calbrate
More informationPredictive Analytics : QM901.1x Prof U Dinesh Kumar, IIMB. All Rights Reserved, Indian Institute of Management Bangalore
Sesson Outlne Introducton to classfcaton problems and dscrete choce models. Introducton to Logstcs Regresson. Logstc functon and Logt functon. Maxmum Lkelhood Estmator (MLE) for estmaton of LR parameters.
More informationSupplementary Notes for Chapter 9 Mixture Thermodynamics
Supplementary Notes for Chapter 9 Mxture Thermodynamcs Key ponts Nne major topcs of Chapter 9 are revewed below: 1. Notaton and operatonal equatons for mxtures 2. PVTN EOSs for mxtures 3. General effects
More informationProperties of Least Squares
Week 3 3.1 Smple Lnear Regresson Model 3. Propertes of Least Squares Estmators Y Y β 1 + β X + u weekly famly expendtures X weekly famly ncome For a gven level of x, the expected level of food expendtures
More informationHow its computed. y outcome data λ parameters hyperparameters. where P denotes the Laplace approximation. k i k k. Andrew B Lawson 2013
Andrew Lawson MUSC INLA INLA s a relatvely new tool that can be used to approxmate posteror dstrbutons n Bayesan models INLA stands for ntegrated Nested Laplace Approxmaton The approxmaton has been known
More informationJoint Statistical Meetings - Biopharmaceutical Section
Iteratve Ch-Square Test for Equvalence of Multple Treatment Groups Te-Hua Ng*, U.S. Food and Drug Admnstraton 1401 Rockvlle Pke, #200S, HFM-217, Rockvlle, MD 20852-1448 Key Words: Equvalence Testng; Actve
More informationA New Refinement of Jacobi Method for Solution of Linear System Equations AX=b
Int J Contemp Math Scences, Vol 3, 28, no 17, 819-827 A New Refnement of Jacob Method for Soluton of Lnear System Equatons AX=b F Naem Dafchah Department of Mathematcs, Faculty of Scences Unversty of Gulan,
More informationEcon Statistical Properties of the OLS estimator. Sanjaya DeSilva
Econ 39 - Statstcal Propertes of the OLS estmator Sanjaya DeSlva September, 008 1 Overvew Recall that the true regresson model s Y = β 0 + β 1 X + u (1) Applyng the OLS method to a sample of data, we estmate
More informationAppendix B: Resampling Algorithms
407 Appendx B: Resamplng Algorthms A common problem of all partcle flters s the degeneracy of weghts, whch conssts of the unbounded ncrease of the varance of the mportance weghts ω [ ] of the partcles
More information