Bayesian Survival Analysis without Covariates Using Optimization and Simulation Tools
|
|
- Anthony Eaton
- 5 years ago
- Views:
Transcription
1 ISSN(Onlne) : ISSN (Prnt) : Internatonal Journal of Innovatve Research n Scence, Engneerng and Technology (An ISO 3297: 2007 Certfed Organzaton) Bayesan Survval Analyss wthout Covarates Usng Optmzaton and Smulaton Tools Yasmn Khan 1, Athar Al Khan 2 Research Scholar, Department of Statstcs & Operatons Research, Algarh Muslm Unversty, Algarh, Uttar Pradesh, Inda 1 Professor, Department of Statstcs & Operatons Research, Algarh Muslm Unversty, Algarh, Uttar Pradesh, Inda 2 ABSTRACT: Webull and exponental dstrbuton s one of the most mportant and flexble dstrbutons n survval analyses. In ths paper, Bayesan regresson analyss wthout covarate wth censorng mechansm s carred out for ntrauterne devce (IUD) data problem. Throughout the Bayesan approach s mplemented usng R and approprate llustratons are made. KEYWORDS: Bayesan Inference, Rght censorng, LaplaceApproxmaton, Survval functon. I. INTRODUCTION The analyss of survval data s a major focus of the statstcs busness. The man theme of ths chapter s the analyss of survval data usng two parametrc models, namely Webull and exponental. These models are very flexble and have been found to provde a good descrpton of many types of tme-to-event data. These are the dstrbutons whch occupy a central role because of ther demonstrated usefulness n a wde range of stuatons. There are many potental lfe tme models but these models are used qute effectvely to analyze skewed data sets and gve best data ft. Webull dstrbuton has two parameter, shape and scale. Its densty, survval and hazard functons, respectvely are: a y f ( y) b b y exp b a1 a a y S( y) exp b Where a > 0 and b > 0 are the shape and scale parameters respectvely. Correspondng probablty densty, survval and hazard functon of exponental model are: 1 f ( t) exp b y S( t) exp b One mportant feature of survval data s the presence of censorng, whch create specal problems n the analyss of the survval data. Lfetme data are censored when the exact falure tme for a specfc tral s unknown. When analyzng censored data, Bayesan methods have an mportant advantage over classcal methods. From a classcal perspectve, confdence nterval and other nferental statements must be made wth respect to repeated samplng of the data. From Bayesan perspectve, only the observed censorng pattern s relevant. There are several categores of censorng but n ths chapter we wll dscuss only rght censorng mechansm. y b Copyrght to IJIRSET DOI: /IJIRSET
2 ISSN(Onlne) : ISSN (Prnt) : Internatonal Journal of Innovatve Research n Scence, Engneerng and Technology Lkelhood functon for rght censored data s, L (An ISO 3297: 2007 Certfed Organzaton) n n Pr( y, ) 1 1 [ f ( y )] [ S( y )] Where s an ndcator varable whch takes value 1 f observaton s censored and 0 otherwse. Secton III begns wth a bref dscusson of the Bayesan analyss of ntercept model based on the Webull and then exponental dstrbuton. There has been growth n the development and applcaton of the Bayesan nference. The Bayesan nference enables us to ft very complex model that cannot be ft by alternatve frequentst methods. To ft the Bayesan models, one needs a statstcal computng envronment. An envronment that meets these requrements s the R (R Development Core Team [1]) software. In ths paper, an attempt has been made to llustrate the Bayesan modellng by usng the R language. The package used n the paper to get the posteror summary s LaplacesDemon Statstcat [2] whch reveals the conceptual smplcty of the Bayesan approach for survval data analyss. Ths package has several optmzaton and smulaton algorthms. The default optmzaton algorthm s LBFGS (Broyden-Fletcher-Goldfarb- Shanno). Also a very popular algorthm,.e. the Nelder-Mead [3], s a dervatve-free, drect search method whch s effcent n small-dmensonal problems. The LaplacesDemon offers numerous MCMC algorthms for smulaton n Bayesan nference, and are, Random-Walk Metropols (RWM), Metropols-wthn-Gbbs (MWG), Delayed Rejecton Metropols (DRM), etc. Ths package provdes a complete Bayesan envronment to the user. Secton VII provdes a model compatblty study based on DIC and devance crteron so that the choce of these two models can be justfed for the data under consderaton. Fnally, n the last secton a bref dscusson and concluson s gven. 1 II. RELATED WORK A sgnfcant amount of work has been done for estmatng the lnear regresson model for censored data khan and khan[4], Khan and khan[5], Khan et al[6], Collet[7] and Puja et al. [8] dscuss the Bayesan survval analyss of head and neck cancer data and conclude that whch therapy has gve better performance for patents. In ths paper all analyses and computaton were undertaken usng LaplacesDemon package avalable n R software. III. BAYESIAN FITTING OF INTERCEPT MODEL WITH CENSORING Here we want to make nferences about the response and ntercept. In ths secton, we wll consder these two models, Webull and exponental models for the data whch s dscuss here are the tme to dscontnuaton of the use of an IUD. The data n Table 1 refer to the number of weeks from the commencement of use of a partcular type of ntrauterne devce (IUD), known as Mult load 250, untl dscontnuaton because of menstrual bleedng problems. Data are gven for 18 women, all of whom were aged between 18 and 35 years and who had experenced two prevous pregnances. Dscontnuaton tmes that are censored are labelled wth an astersk. Table I: Tme n weeks to dscontnuaton of the use of an IUD * 18* 19 23* * 54* 56* * * 107* The ntercept model s n the form of y e y ~ Webull( a, b) Copyrght to IJIRSET DOI: /IJIRSET
3 ISSN(Onlne) : ISSN (Prnt) : Internatonal Journal of Innovatve Research n Scence, Engneerng and Technology Where, a and b are shape and scale parameter and (An ISO 3297: 2007 Certfed Organzaton) X, log( b) X ~ N(0,1000) ~ halfcauchy(25) IV. FITTING OF WEIBULL MODEL Webull dstrbuton has two parameters, shape and scale. It can be denoted as thus, the lkelhood s p( y a, b) whch mples log-lkelhood as n 1 y ~ W ( a, b) a y b b a1 y exp b n a1 a a y y log p( y a, b) log exp 1 b b b Bayesan fttng of Webull model for ths data can be done n R by usng functon LaplaceApproxmaton and then wth LaplacesDemon. Its fttng ncludes codes for creaton of data and defnton of model. R codes to ft Webull model s beng descrbed below. A.CREATION OF DATA LaplaceApproxmaton functon requres data that s specfed n a lst. Though most R functons use data n the form of a data frame, Laplace's Demon uses one or more numerc matrces n a lst. It s much faster to process a numerc matrx than a data frame n teratve estmaton. For the above data of 18 patents wth dscontnuaton tme has gven the survval tmes n weeks as a response varable. Snce ntercept s the only term n the model, a vector of 1s s nserted nto desgn matrx X. Thus, X=1 ndcates only column of 1s n the desgn matrx. lbrary(laplacesdemon) y<-c(10,13,18,19,23,30,36,38,54,56,59,75,93,97,104,107,107,107) censor<-c(1,0,0,1,0,1,1,0,0,0,1,1,1,1,0,1,0,0) N<-18 J<-1 X<-matrx(1,nrow=length(y)) mon.names<-c("lp","shape") parm.names<-as.parm.names(lst(beta=rep(0,j),log.shape=0)) MyData<-lst(J=J,X=X,mon.names=mon.names,parm.names=parm.names,y=y) There are J=1 ndependent varables (response varable). The R codes defned above must have a name specfed for each parameter n the vector parm.names, and parameter names must be ncluded wth the data n a lst called as.parm.names. The user must specfy the number of observatons n the data as ether a scalar n or N. If these are not found by the LaplaceApproxmaton or LaplacesDemon functons, then t wll attempt to determne sample sze as the number of rows n y or Y. B. INITIAL VALUES Laplace's Demon requres a vector of ntal values for the parameters. Each ntal value s a startng pont for the estmaton of a parameter. When all ntal values are set to zero, Laplace's Demon wll optmze ntal values usng a ht-and-run algorthm wth adaptve length n the LaplaceApproxmaton functon. Intal.Values <- c(rep(0,j), log(1)) a Copyrght to IJIRSET DOI: /IJIRSET
4 ISSN(Onlne) : ISSN (Prnt) : Internatonal Journal of Innovatve Research n Scence, Engneerng and Technology (An ISO 3297: 2007 Certfed Organzaton) C. MODEL SPECIFICATION Model<-functon(parm,Data) { beta<-parm[1:data$j] shape<-exp(parm[data$j+1]) beta.pror<-sum(dnorm(beta,0,1000,log=t)) shape.pror<-dhalfcauchy(shape,25,log=t) mu<-tcrossprod(beta,data$x) scale<-exp(mu) w1<-log(scale)+log(shape)+(shape-1)*log(y)-scale*y^shape s1<--scale*y^shape LL<-censor*w1+(1-censor)*s1 LL<-sum(LL) LP<-LL+beta.pror+shape.pror Modelout<-lst(LP=LP,Dev=-2*LL,Montor=c(LP,shape),yhat=mu,parm=parm) return(modelout) } D. FITTING OF MODEL WITH LAPLACE APPROXIMATION M1<-LaplaceApproxmaton(Model,Intal.Values,Data=MyData,Sample=1000,Iteratons=10000) The output summary s tabulated n Table 2 and Table 3. Table 2: Approxmated posteror summary of IUD data usng LaplaceApproxmaton functon wth posteror mode, posteror sd and ther quantles. Mode SD LB UB beta Log.shape Table 3: Smuated posteror summary of IUD data usng samplng mportance resamplng n LaplaceApproxmaton functon wth posteror mean, posteror sd and ther quantles. Mean SD LB Medan UB beta Log.shape V. MEDIAN SURVIVAL TIME Medan survval tme s the tme beyond whch 50% of the ndvduals n the populaton under study are expected to survve, and s gven by that value t(50) whch s such that S(t(50)) = 0.5. As we know Survval tme dstrbuton s always postvely skew, the medan s the preferred summary measure of the locaton of the dstrbuton. Once the survvor functon has been estmated, then t s easy to obtan an estmate of the medan survval tme. The estmated medan survval tme, {t(50)}, s defned to be the smallest observed survval tme for whch the value of the estmated survvor functon s less than 0.5. t ( 50) mn t S( t ) 0.5, Where t s the observed survval tme for the th ndvdual, 1,2,... n. snce the estmated survvor functon only changes at a death tme, ths s equvalent to the defnton where j t t ( 50) mn S( t ) 0.5, t s the jth ordered death tme, j 1,2,... r. j Copyrght to IJIRSET DOI: /IJIRSET
5 ISSN(Onlne) : ISSN (Prnt) : Internatonal Journal of Innovatve Research n Scence, Engneerng and Technology (An ISO 3297: 2007 Certfed Organzaton) Fg 1. Medan survval tme and other percentle of IUD data. The correspondng nterval for the true medan dscontnuaton tme s (53.81,126.69), so that there s 95% chance that the nterval from 54 days to 127 days ncludes the true value of the medan dscontnuaton tme. VI. SIMULATION STUDY: INDEPENDENT METROPOLIS ALGORITHM For the purpose of smulaton from jont posteror dstrbuton, ndependent Metropols algorthm wll be performed. Now, n ths secton we have to explore IUD data usng functon LaplacesDemon. It maxmzes the logarthm of the unnormalzed jont posteror densty wth MCMC and provdes samples of the margnal posteror dstrbutons, devance, and other montored varables. In LaplacesDemon functon there s an argument called Algorthm, here the algorthm used for smulaton from jont posteror dstrbuton s ndependent-metropols algorthm. Multvarate normal has been treated as a proposal dstrbuton q ( ). Here, the proposal dstrbuton does not depend on the prevous state of the chan. The IM algorthm s effcent when the proposal s a good approxmaton of the target posteror dstrbuton. Good ndependent proposal denstes can be based on LaplaceApproxmaton (Terney and Kadane[9] Terney et al.[10] and Erkanl[11]). Thus, a generally successful proposal can be obtaned by a multvarate normal dstrbuton wth mean equal to the posteror mode and precson matrx 2 log p( y) H ( ) j that s, mnus the second dervatve matrx of the log-posteror densty Copyrght to IJIRSET DOI: /IJIRSET
6 ISSN(Onlne) : ISSN (Prnt) : Internatonal Journal of Innovatve Research n Scence, Engneerng and Technology (An ISO 3297: 2007 Certfed Organzaton) log p( y) log p( y ) log p( ) evaluated at the posteror mode. Consequently, an effcent proposal s gven by, 1 q ( ) N, H ( ) The acceptance probablty, when proposng a transton from to, s gven by p( y) q( ) mn 1, p( y) q( ) Whch can be re-expressed as, w( ) mn 1,, w( ) Where, w( ) p( y)/ q( ) s the rato between the target and the proposal dstrbuton and s equvalent to the mportance weght used n mportance samplng Ntzroufras[12]. Ths theory s mplement n LaplacesDemon wth object name M2, Intal.Values<-as.ntal.values(M1) M2<-LaplacesDemon(Model,Data=MyData,Intal.Values,Covar=M1$Covar,Algorthm="IM", Iteratons=10000,Status=F,Specs=lst(mu=M1$Summary1[1:length(Intal.Values),1])) The output obtaned from M2 object s summarzed n Table 4. Table 4 contans posteror mean, posteror sd and respectve quantles. The margnal posteror densty plots of both parameter.e. scale and shape parameters from LaplaceApproxmaton and LaplacesDemon functons are reported s Fgure 2. Table 4: Smulated posteror summary usng LaplacesDemon functonwth posteror mean and quantles. Mean SD LB Medan UB beta Posteror densty of scale Posteror densty of shape Densty LaplaceApproxmaton LaplacesDemon Densty LaplaceApproxmaton LaplacesDemon N = 1000 Bandwdth = N = 1000 Bandwdth = 0.25 Log.shape Fg 2: Margnal posteror densty plots of scale (left) and shape (rght) parameter. Sold lne depcts the margnal posteror densty obtan by LaplaceApproxmaton method wth posteror mode and dotted lne s posteror densty plot obtan by smulaton Copyrght to IJIRSET DOI: /IJIRSET
7 ISSN(Onlne) : ISSN (Prnt) : Internatonal Journal of Innovatve Research n Scence, Engneerng and Technology (An ISO 3297: 2007 Certfed Organzaton) usng LaplacesDemon functon wth posteror mean Smlarly the sold and dotted lne of shape parameter n rght, depcts the margnal posteror densty wth posteror mode 1.68 amd posteror mean 1.66, respectvely. VI. FITTING OF EXPONENTIAL MODEL R-codes for the fttng exponental dstrbuton s beng descrbe below, lbrary(laplacesdemon) y<-c(10,13,18,19,23,30,36,38,54,56,59,75,93,97,104,107,107,107) censor<-c(1,0,0,1,0,1,1,0,0,0,1,1,1,1,0,1,0,0) N<-18 J<-2 X<-matrx(1,nrow=length(y)) mon.names<-c("lp") parm.names<-as.parm.names(lst(beta=rep(0,j))) MyData<-lst(J=J,X=X,mon.names=mon.names,parm.names=parm.names,y=y) Intal.Values <- c(rep(0,j)) Model<-functon(parm,Data) { beta<-parm[1:data$j] beta.pror<-sum(dnorm(beta,0,1000,log=t)) mu<-tcrossprod(beta,data$x) scale<-exp(mu) w1<-log(scale)+log(1)+(1-1)*log(y)-scale*y^1 s1<--scale*y^1 LL<-censor*w1+(1-censor)*s1 LL<-sum(LL) LP<-LL+beta.pror Modelout<-lst(LP=LP,Dev=-2*LL,Montor=c(LP),yhat=mu,parm=parm) return(modelout) } M2<-LaplaceApproxmaton(Model,Intal.Values,Data=MyData,Sample=10000,Iteratons=10000) Table 5: Approxmated posteror summary of IUD data usng LaplaceApproxmaton functon wth posteror mode, posteror sd and ther quantles. Mode SD LB UB beta[1] beta[2] Table 6: Smuated posteror summary of IUD data usng LaplacesDemon functon wth posteror mean, posteror sd and ther quantles. Mean SD LB Medan UB beta[1] beta[2] Copyrght to IJIRSET DOI: /IJIRSET
8 ISSN(Onlne) : ISSN (Prnt) : Internatonal Journal of Innovatve Research n Scence, Engneerng and Technology (An ISO 3297: 2007 Certfed Organzaton) Posteror densty of scale Densty LaplaceApproxmaton LaplacesDemon N = Bandwdth = Fg 3: Margnal posteror densty of scale parameter. Sold lne ( ) s the densty obtan as a result of LaplaceApproxmaton functon and dotted lne (-----) s the densty obtan by LaplacesDemon functon. The posteror mode s and posteror mean s also evdent from Table 5 and Table 6 after dong exp(beta[1]). A. MEDIAN SURVIVAL TIME OF EXPONENTIAL DISTRIBUTION. Medan of exponental dstrbuton s gven as, t med log b Densty plots of medan dscontnuaton survval tmes and also percentle s reported n Fgure 4. Fg 4. Medan survval tme and other percentle of IUD data under the assumpton of exponental dstrbuton. The 90th percentle of the dstrbuton of dscontnuaton tmes s 287 days. Ths means that on the assumpton that the rsk of dscontnuaton the use of an IUD s ndependent of tme and 90% of women wll have a dscontnuaton tme less than 287 days. Copyrght to IJIRSET DOI: /IJIRSET
9 ISSN(Onlne) : ISSN (Prnt) : Internatonal Journal of Innovatve Research n Scence, Engneerng and Technology (An ISO 3297: 2007 Certfed Organzaton) VII. MODEL COMPARISON In ths secton, a goodness-of-ft crteron tests would be appled n order to verfy whch dstrbuton fts better for ths data. To compare the two models; namely exponental and Webull dstrbuton, the model selecton crteron preferred by the Bayesans and lkelhoodsts are devance and devance nformaton crteron (DIC, Spegelhalter et al.[13], s a model assessment tool). A smaller DIC and devance ndcates a better ft to the data set. Table 7: Model comparson of Webull and exponental model for IUD data. Both devance and DIC crteron support Webull dstrbuton s a better choce as compared to exponental dstrbuton. Model Devance DIC Webull Exponental VII. CONCLUSION In trals nvolvng contraceptves, preventon of pregnancy s an obvous crteron for acceptablty. The IUD data whch was dscussed n the whole chapter s the tme orgn corresponds to the frst day n whch woman uses IUD, and the end pont s the dscontnuaton because of bleedng problems. Some women n the study ceased usng IUD because of the desre of pregnancy, or because they had no further need for a contraceptve, whle others were smply lost to follow-up. The data was fully analyzed n Bayesan framework. Fg 5: Estmated survval curve for Webull and exponental dstrbuton showng ther medan survval tme. The IUD data had been analyzed both by both Webull and exponental models. The dscontnuaton tmes when modeled under an exponental dstrbuton, the medan dscontnuaton tme s days whch could also be seen n above Fgure, and an estmate of the 90 percentle of the dstrbuton of dscontnuaton tmes s 287 days as reported n Fgure 4. Ths means that on the assumpton that the rsk of dscontnuaton the use of an IUD s ndependent of tme and 90% of women wll have a dscontnuaton tme less than 287 days. On the other hand, under the assumpton of Copyrght to IJIRSET DOI: /IJIRSET
10 ISSN(Onlne) : ISSN (Prnt) : Internatonal Journal of Innovatve Research n Scence, Engneerng and Technology (An ISO 3297: 2007 Certfed Organzaton) Webull dstrbuton, the estmated medan dscontnuaton tme for IUD data s days and the 90 percentle s 162 days as reported n Fgure 1. So 90%of women wll have a dscontnuaton tme less than 162 days. Hence t s very clear that the medan s therefore estmated more precsely when the dscontnuaton tme are assumed to have a Webull dstrbuton. It would also be evdent from Table 7. Table 7 shows the devance and DIC of both models. Here the devance of Webull model s 100 whereas exponental model have 207, whch prove that Webull model s sutable for ths data and exponental model does not provde an acceptable ft to ths data. REFERENCES 1. R Development Core Team, R: A Language and Envronment for Statstcal Computng, R Foundaton for Statstcal Computng, Venna, Austra, ISBN , URL Statstcat, LLC, LaplacesDemon: Complete envronment for Bayesan Inference, R package verson , URL com, J. A. Nelder, and R. Mead, A smplex method for functon mnmzaton, The Computer Journal Vol. 7, , Y. Khan, and A. A. Khan, Bayesan regresson analyss of Webull model wth R and Bugs, Statstcal Methodologes and Applcatons, Vol. 2, , Y. Khan, and A. A. Khan, Bayesan analyss of Webull and log normal survval models wth censorng mechansm, Internatonal Journal of Appled Mathematcs, Vol. 26, , Y. Khan, M. T. Alhtar, R. Shehla, and A. A. Khan, Bayesan modellng of forestry data wth R usng optmzaton and smulaton tools, Journal of Appled Analyss and computaton, Vol. 5, 38-51, D. Collet, Modellng Survval Data n Medcal Research, London: Chapman & Hall, M. Puja, K. S. Puneet, R. S. Sngh, and S. K. Upadhyay, Bayesan Survval Analyss of Head and Neck Cancer Data Usng Lognormal model, Communcatons n Statstcs-Smulaton and computaton, L. Terney, and J. B. Kadane, Accurate approxmatons for posteror moments and margnal denstes", Journal of the Amercan Statstcal Assocaton, Vol. 81, 82-86, L. Terney, R. E. Kass, and J. B. Kadane, Fully exponental Laplace approxmatons to expectatons and varances of non-postve functons, Journal of the Amercan Statstcal Assocaton, Vol. 84, No. 407, pp , Erkanl, Laplace approxmatons for posteror expectaton when the model occurs at the boundary of the parameter space, Journal of the Amercan Statstcal Assocaton, Vol. 89, , Ntzoufras, Bayesan Modellng usng WnBugs, John Wley & Sons, D. Spegelhalter, N. Best, B. Carln, and A. van der Lnde, Bayesan measures of model complexty and ft, Journal of the Royal Statstcal Socety, B Vol. 64, , Copyrght to IJIRSET DOI: /IJIRSET
A Bayesian Approach to Survival Analysis of Inverse Gaussian Model with Laplace Approximation
Internatonal Journal of Statstcs and Applcatons 016, 6(6): 391-398 DOI: 10.593/j.statstcs.0160606.08 A Bayesan Approach to Survval Analyss of Inverse Gaussan Model wth Laplace Approxmaton Najrullah Khan
More informationNon-Mixture Cure Model for Interval Censored Data: Simulation Study ABSTRACT
Malaysan Journal of Mathematcal Scences 8(S): 37-44 (2014) Specal Issue: Internatonal Conference on Mathematcal Scences and Statstcs 2013 (ICMSS2013) MALAYSIAN JOURNAL OF MATHEMATICAL SCIENCES Journal
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 informationComparison of the Population Variance Estimators. of 2-Parameter Exponential Distribution Based on. Multiple Criteria Decision Making Method
Appled Mathematcal Scences, Vol. 7, 0, no. 47, 07-0 HIARI Ltd, www.m-hkar.com Comparson of the Populaton Varance Estmators of -Parameter Exponental Dstrbuton Based on Multple Crtera Decson Makng Method
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 informationStatistical inference for generalized Pareto distribution based on progressive Type-II censored data with random removals
Internatonal Journal of Scentfc World, 2 1) 2014) 1-9 c Scence Publshng Corporaton www.scencepubco.com/ndex.php/ijsw do: 10.14419/jsw.v21.1780 Research Paper Statstcal nference for generalzed Pareto dstrbuton
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 informationRELIABILITY ASSESSMENT
CHAPTER Rsk Analyss n Engneerng and Economcs RELIABILITY ASSESSMENT A. J. Clark School of Engneerng Department of Cvl and Envronmental Engneerng 4a CHAPMAN HALL/CRC Rsk Analyss for Engneerng Department
More informationChange Point Estimation for Pareto Type-II Model
Journal of Modern Appled Statstcal Methods Volume 3 Issue Artcle 5--04 Change Pont Estmaton for Pareto Type-II Model Gyan Prakash S. N. Medcal College, Agra, U. P., Inda, ggyanj@yahoo.com Follow ths and
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 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 informationUncertainty as the Overlap of Alternate Conditional Distributions
Uncertanty as the Overlap of Alternate Condtonal Dstrbutons Olena Babak and Clayton V. Deutsch Centre for Computatonal Geostatstcs Department of Cvl & Envronmental Engneerng Unversty of Alberta An mportant
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 informationHere is the rationale: If X and y have a strong positive relationship to one another, then ( x x) will tend to be positive when ( y y)
Secton 1.5 Correlaton In the prevous sectons, we looked at regresson and the value r was a measurement of how much of the varaton n y can be attrbuted to the lnear relatonshp between y and x. In ths secton,
More informationInternational Journal of Industrial Engineering Computations
Internatonal Journal of Industral Engneerng Computatons 4 (03) 47 46 Contents lsts avalable at GrowngScence Internatonal Journal of Industral Engneerng Computatons homepage: www.growngscence.com/ec Modelng
More informationBayesian predictive Configural Frequency Analysis
Psychologcal Test and Assessment Modelng, Volume 54, 2012 (3), 285-292 Bayesan predctve Confgural Frequency Analyss Eduardo Gutérrez-Peña 1 Abstract Confgural Frequency Analyss s a method for cell-wse
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 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 informationInterval Estimation of Stress-Strength Reliability for a General Exponential Form Distribution with Different Unknown Parameters
Internatonal Journal of Statstcs and Probablty; Vol. 6, No. 6; November 17 ISSN 197-73 E-ISSN 197-74 Publshed by Canadan Center of Scence and Educaton Interval Estmaton of Stress-Strength Relablty for
More informationDouble Acceptance Sampling Plan for Time Truncated Life Tests Based on Transmuted Generalized Inverse Weibull Distribution
J. Stat. Appl. Pro. 6, No. 1, 1-6 2017 1 Journal of Statstcs Applcatons & Probablty An Internatonal Journal http://dx.do.org/10.18576/jsap/060101 Double Acceptance Samplng Plan for Tme Truncated Lfe Tests
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 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 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 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 informationParameters Estimation of the Modified Weibull Distribution Based on Type I Censored Samples
Appled Mathematcal Scences, Vol. 5, 011, no. 59, 899-917 Parameters Estmaton of the Modfed Webull Dstrbuton Based on Type I Censored Samples Soufane Gasm École Supereure des Scences et Technques de Tuns
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 informationInfluence Diagnostics on Competing Risks Using Cox s Model with Censored Data. Jalan Gombak, 53100, Kuala Lumpur, Malaysia.
Proceedngs of the 8th WSEAS Internatonal Conference on APPLIED MAHEMAICS, enerfe, Span, December 16-18, 5 (pp14-138) Influence Dagnostcs on Competng Rsks Usng Cox s Model wth Censored Data F. A. M. Elfak
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 informationChapter 20 Duration Analysis
Chapter 20 Duraton Analyss Duraton: tme elapsed untl a certan event occurs (weeks unemployed, months spent on welfare). Survval analyss: duraton of nterest s survval tme of a subject, begn n an ntal state
More informationThis column is a continuation of our previous column
Comparson of Goodness of Ft Statstcs for Lnear Regresson, Part II The authors contnue ther dscusson of the correlaton coeffcent n developng a calbraton for quanttatve analyss. Jerome Workman Jr. and Howard
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 informationA Robust Method for Calculating the Correlation Coefficient
A Robust Method for Calculatng the Correlaton Coeffcent E.B. Nven and C. V. Deutsch Relatonshps between prmary and secondary data are frequently quantfed usng the correlaton coeffcent; however, the tradtonal
More informationBayesian Planning of Hit-Miss Inspection Tests
Bayesan Plannng of Ht-Mss Inspecton Tests Yew-Meng Koh a and Wllam Q Meeker a a Center for Nondestructve Evaluaton, Department of Statstcs, Iowa State Unversty, Ames, Iowa 5000 Abstract Although some useful
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 informationNegative Binomial Regression
STATGRAPHICS Rev. 9/16/2013 Negatve Bnomal Regresson Summary... 1 Data Input... 3 Statstcal Model... 3 Analyss Summary... 4 Analyss Optons... 7 Plot of Ftted Model... 8 Observed Versus Predcted... 10 Predctons...
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 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 informationDurban Watson for Testing the Lack-of-Fit of Polynomial Regression Models without Replications
Durban Watson for Testng the Lack-of-Ft of Polynomal Regresson Models wthout Replcatons Ruba A. Alyaf, Maha A. Omar, Abdullah A. Al-Shha ralyaf@ksu.edu.sa, maomar@ksu.edu.sa, aalshha@ksu.edu.sa Department
More informationSTK4080/9080 Survival and event history analysis
SK48/98 Survval and event hstory analyss Lecture 7: Regresson modellng Relatve rsk regresson Regresson models Assume that we have a sample of n ndvduals, and let N (t) count the observed occurrences of
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 informationCredit Card Pricing and Impact of Adverse Selection
Credt Card Prcng and Impact of Adverse Selecton Bo Huang and Lyn C. Thomas Unversty of Southampton Contents Background Aucton model of credt card solctaton - Errors n probablty of beng Good - Errors n
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 informationKernel Methods and SVMs Extension
Kernel Methods and SVMs Extenson The purpose of ths document s to revew materal covered n Machne Learnng 1 Supervsed Learnng regardng support vector machnes (SVMs). Ths document also provdes a general
More informationUsing T.O.M to Estimate Parameter of distributions that have not Single Exponential Family
IOSR Journal of Mathematcs IOSR-JM) ISSN: 2278-5728. Volume 3, Issue 3 Sep-Oct. 202), PP 44-48 www.osrjournals.org Usng T.O.M to Estmate Parameter of dstrbutons that have not Sngle Exponental Famly Jubran
More informationInterval Estimation in the Classical Normal Linear Regression Model. 1. Introduction
ECONOMICS 35* -- NOTE 7 ECON 35* -- NOTE 7 Interval Estmaton n the Classcal Normal Lnear Regresson Model Ths note outlnes the basc elements of nterval estmaton n the Classcal Normal Lnear Regresson Model
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 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 informationSimulated Power of the Discrete Cramér-von Mises Goodness-of-Fit Tests
Smulated of the Cramér-von Mses Goodness-of-Ft Tests Steele, M., Chaselng, J. and 3 Hurst, C. School of Mathematcal and Physcal Scences, James Cook Unversty, Australan School of Envronmental Studes, Grffth
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 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 informationMaximum Likelihood Estimation of Binary Dependent Variables Models: Probit and Logit. 1. General Formulation of Binary Dependent Variables Models
ECO 452 -- OE 4: Probt and Logt Models ECO 452 -- OE 4 Maxmum Lkelhood Estmaton of Bnary Dependent Varables Models: Probt and Logt hs note demonstrates how to formulate bnary dependent varables models
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 informationQUASI-LIKELIHOOD APPROACH TO RATER AGREEMENT PLUS LINEAR BY LINEAR ASSOCIATION MODEL FOR ORDINAL CONTINGENCY TABLES
Journal of Statstcs: Advances n Theory and Applcatons Volume 6, Number, 26, Pages -5 Avalable at http://scentfcadvances.co.n DOI: http://dx.do.org/.8642/jsata_72683 QUASI-LIKELIHOOD APPROACH TO RATER AGREEMENT
More information4.3 Poisson Regression
of teratvely reweghted least squares regressons (the IRLS algorthm). We do wthout gvng further detals, but nstead focus on the practcal applcaton. > glm(survval~log(weght)+age, famly="bnomal", data=baby)
More informationConjugacy and the Exponential Family
CS281B/Stat241B: Advanced Topcs n Learnng & Decson Makng Conjugacy and the Exponental Famly Lecturer: Mchael I. Jordan Scrbes: Bran Mlch 1 Conjugacy In the prevous lecture, we saw conjugate prors for the
More informationCalculating CLs Limits. Abstract
DØNote 4492 Calculatng CLs Lmts Harrson B. Prosper Florda State Unversty, Tallahassee, Florda 32306 (Dated: June 8, 2004) Abstract Ths note suggests how the calculaton of lmts based on the CLs method mght
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 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 informationStat 543 Exam 2 Spring 2016
Stat 543 Exam 2 Sprng 2016 I have nether gven nor receved unauthorzed assstance on ths exam. Name Sgned Date Name Prnted Ths Exam conssts of 11 questons. Do at least 10 of the 11 parts of the man exam.
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 informationUncertainty and auto-correlation in. Measurement
Uncertanty and auto-correlaton n arxv:1707.03276v2 [physcs.data-an] 30 Dec 2017 Measurement Markus Schebl Federal Offce of Metrology and Surveyng (BEV), 1160 Venna, Austra E-mal: markus.schebl@bev.gv.at
More informationProportional Hazards Regression Model for Time to Event Breast Cancer Data: A Bayesian Approach
IJIRST Internatonal Journal for Innovatve Research n Scence & Technology Volume 3 Issue 08 January 2017 ISSN (onlne): 2349-6010 Proportonal Hazards Regresson Model for Tme to Event Breast Cancer Data:
More informationDegradation Data Analysis Using Wiener Process and MCMC Approach
Engneerng Letters 5:3 EL_5_3_0 Degradaton Data Analyss Usng Wener Process and MCMC Approach Chunpng L Hubng Hao Abstract Tradtonal relablty assessment methods are based on lfetme data. However the lfetme
More informationChapter 8 Indicator Variables
Chapter 8 Indcator Varables In general, e explanatory varables n any regresson analyss are assumed to be quanttatve n nature. For example, e varables lke temperature, dstance, age etc. are quanttatve n
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 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 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 informationStat 543 Exam 2 Spring 2016
Stat 543 Exam 2 Sprng 206 I have nether gven nor receved unauthorzed assstance on ths exam. Name Sgned Date Name Prnted Ths Exam conssts of questons. Do at least 0 of the parts of the man exam. I wll score
More informationStatistics II Final Exam 26/6/18
Statstcs II Fnal Exam 26/6/18 Academc Year 2017/18 Solutons Exam duraton: 2 h 30 mn 1. (3 ponts) A town hall s conductng a study to determne the amount of leftover food produced by the restaurants n the
More informationSmall Area Interval Estimation
.. Small Area Interval Estmaton Partha Lahr Jont Program n Survey Methodology Unversty of Maryland, College Park (Based on jont work wth Masayo Yoshmor, Former JPSM Vstng PhD Student and Research Fellow
More informationA: Details of MCMC for Fitting DDP-GP Model
Supplement: Bayesan Nonparametrc Estmaton for Dynamc Treatment Regmes wth Sequental Transton Tmes A: Detals of MCMC for Fttng DDP-GP Model Summary of Model p(y k x k, F k ) = F k (y k x k ) F k DDP-GP
More informationExplaining the Stein Paradox
Explanng the Sten Paradox Kwong Hu Yung 1999/06/10 Abstract Ths report offers several ratonale for the Sten paradox. Sectons 1 and defnes the multvarate normal mean estmaton problem and ntroduces Sten
More informationBasically, if you have a dummy dependent variable you will be estimating a probability.
ECON 497: Lecture Notes 13 Page 1 of 1 Metropoltan State Unversty ECON 497: Research and Forecastng Lecture Notes 13 Dummy Dependent Varable Technques Studenmund Chapter 13 Bascally, f you have a dummy
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 informationAn adaptive SMC scheme for ABC. Bayesian Computation (ABC)
An adaptve SMC scheme for Approxmate Bayesan Computaton (ABC) (ont work wth Prof. Mke West) Department of Statstcal Scence - Duke Unversty Aprl/2011 Approxmate Bayesan Computaton (ABC) Problems n whch
More informationMATH 829: Introduction to Data Mining and Analysis The EM algorithm (part 2)
1/16 MATH 829: Introducton to Data Mnng and Analyss The EM algorthm (part 2) Domnque Gullot Departments of Mathematcal Scences Unversty of Delaware Aprl 20, 2016 Recall 2/16 We are gven ndependent observatons
More informationTHE ROYAL STATISTICAL SOCIETY 2006 EXAMINATIONS SOLUTIONS HIGHER CERTIFICATE
THE ROYAL STATISTICAL SOCIETY 6 EXAMINATIONS SOLUTIONS HIGHER CERTIFICATE PAPER I STATISTICAL THEORY The Socety provdes these solutons to assst canddates preparng for the eamnatons n future years and for
More informationHidden Markov Models & The Multivariate Gaussian (10/26/04)
CS281A/Stat241A: Statstcal Learnng Theory Hdden Markov Models & The Multvarate Gaussan (10/26/04) Lecturer: Mchael I. Jordan Scrbes: Jonathan W. Hu 1 Hdden Markov Models As a bref revew, hdden Markov models
More informationNotes on Frequency Estimation in Data Streams
Notes on Frequency Estmaton n Data Streams In (one of) the data streamng model(s), the data s a sequence of arrvals a 1, a 2,..., a m of the form a j = (, v) where s the dentty of the tem and belongs to
More informationStructure and Drive Paul A. Jensen Copyright July 20, 2003
Structure and Drve Paul A. Jensen Copyrght July 20, 2003 A system s made up of several operatons wth flow passng between them. The structure of the system descrbes the flow paths from nputs to outputs.
More informationComments on Detecting Outliers in Gamma Distribution by M. Jabbari Nooghabi et al. (2010)
Comments on Detectng Outlers n Gamma Dstrbuton by M. Jabbar Nooghab et al. (21) M. Magdalena Lucn Alejandro C. Frery September 17, 215 arxv:159.55v1 [stat.co] 16 Sep 215 Ths note shows that the results
More informationand V is a p p positive definite matrix. A normal-inverse-gamma distribution.
OSR Journal of athematcs (OSR-J) e-ssn: 78-578, p-ssn: 39-765X. Volume 3, ssue 3 Ver. V (ay - June 07), PP 68-7 www.osrjournals.org Comparng The Performance of Bayesan And Frequentst Analyss ethods of
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 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 informationBézier curves. Michael S. Floater. September 10, These notes provide an introduction to Bézier curves. i=0
Bézer curves Mchael S. Floater September 1, 215 These notes provde an ntroducton to Bézer curves. 1 Bernsten polynomals Recall that a real polynomal of a real varable x R, wth degree n, s a functon of
More informationA New Method for Estimating Overdispersion. David Fletcher and Peter Green Department of Mathematics and Statistics
A New Method for Estmatng Overdsperson Davd Fletcher and Peter Green Department of Mathematcs and Statstcs Byron Morgan Insttute of Mathematcs, Statstcs and Actuaral Scence Unversty of Kent, England Overvew
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 informationA Particle Filter Algorithm based on Mixing of Prior probability density and UKF as Generate Importance Function
Advanced Scence and Technology Letters, pp.83-87 http://dx.do.org/10.14257/astl.2014.53.20 A Partcle Flter Algorthm based on Mxng of Pror probablty densty and UKF as Generate Importance Functon Lu Lu 1,1,
More informationGaussian process classification: a message-passing viewpoint
Gaussan process classfcaton: a message-passng vewpont Flpe Rodrgues fmpr@de.uc.pt November 014 Abstract The goal of ths short paper s to provde a message-passng vewpont of the Expectaton Propagaton EP
More informationTime-Varying Systems and Computations Lecture 6
Tme-Varyng Systems and Computatons Lecture 6 Klaus Depold 14. Januar 2014 The Kalman Flter The Kalman estmaton flter attempts to estmate the actual state of an unknown dscrete dynamcal system, gven nosy
More informationChapter 14 Simple Linear Regression
Chapter 4 Smple Lnear Regresson Chapter 4 - Smple Lnear Regresson Manageral decsons often are based on the relatonshp between two or more varables. Regresson analss can be used to develop an equaton showng
More informationUSE OF DOUBLE SAMPLING SCHEME IN ESTIMATING THE MEAN OF STRATIFIED POPULATION UNDER NON-RESPONSE
STATISTICA, anno LXXV, n. 4, 015 USE OF DOUBLE SAMPLING SCHEME IN ESTIMATING THE MEAN OF STRATIFIED POPULATION UNDER NON-RESPONSE Manoj K. Chaudhary 1 Department of Statstcs, Banaras Hndu Unversty, Varanas,
More informationA New Evolutionary Computation Based Approach for Learning Bayesian Network
Avalable onlne at www.scencedrect.com Proceda Engneerng 15 (2011) 4026 4030 Advanced n Control Engneerng and Informaton Scence A New Evolutonary Computaton Based Approach for Learnng Bayesan Network Yungang
More informationFirst Year Examination Department of Statistics, University of Florida
Frst Year Examnaton Department of Statstcs, Unversty of Florda May 7, 010, 8:00 am - 1:00 noon Instructons: 1. You have four hours to answer questons n ths examnaton.. You must show your work to receve
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 informationCHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE
CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE Analytcal soluton s usually not possble when exctaton vares arbtrarly wth tme or f the system s nonlnear. Such problems can be solved by numercal tmesteppng
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 informationP R. Lecture 4. Theory and Applications of Pattern Recognition. Dept. of Electrical and Computer Engineering /
Theory and Applcatons of Pattern Recognton 003, Rob Polkar, Rowan Unversty, Glassboro, NJ Lecture 4 Bayes Classfcaton Rule Dept. of Electrcal and Computer Engneerng 0909.40.0 / 0909.504.04 Theory & Applcatons
More informationStatistical tables are provided Two Hours UNIVERSITY OF MANCHESTER. Date: Wednesday 4 th June 2008 Time: 1400 to 1600
Statstcal tables are provded Two Hours UNIVERSITY OF MNCHESTER Medcal Statstcs Date: Wednesday 4 th June 008 Tme: 1400 to 1600 MT3807 Electronc calculators may be used provded that they conform to Unversty
More informationMaximum Likelihood Estimation of Binary Dependent Variables Models: Probit and Logit. 1. General Formulation of Binary Dependent Variables Models
ECO 452 -- OE 4: Probt and Logt Models ECO 452 -- OE 4 Mamum Lkelhood Estmaton of Bnary Dependent Varables Models: Probt and Logt hs note demonstrates how to formulate bnary dependent varables models for
More information