Spatial Extreme Value Using Bayesian Hierarchical Model For Precipitation Return Levels Prediction
|
|
- Ambrose Booker
- 5 years ago
- Views:
Transcription
1 PROCEEDING OF 3 RD INTERNATIONAL CONFERENCE ON RESEARCH, IMPLEMENTATION AND EDUCATION OF MATHEMATICS AND SCIENCE YOGYAKARTA, 6 7 MAY 06 Spatal Extreme Value Usng Bayesan Herarchcal Model For Precptaton Return Levels Predcton Indra Tsan Hazhah, Sutkno, Dedy Dw Prastyo 3 Dept. of Statstcs, Insttut Teknolog Sepuluh Nopember, Surabaya, Indonesa Dept. of Statstcs, Insttut Teknolog Sepuluh Nopember, Surabaya, Indonesa 3 Dept. of Statstcs, Insttut Teknolog Sepuluh Nopember, Surabaya, Indonesa e-mal: ndratsan@gmal.com Abstract Extreme precptaton s a rare natural phenomenon categorzed as one of extreme clmates ndcator. It leads to natural dsasters such as floods and landsldes. Predcton of precptaton return level,.e. quantle whch has probablty /m of beng exceeded n m perod, become qute mportant to reduce the negatve mpact caused by ths extreme event. The lnk between observaton recorded on a partcular tme frame and quanttes of longer tme scales such as return level s provded by Extreme Value Theory (EVT) commonly used to learn the behavors of extreme events. Gven the observatons are recorded from several locatons, the extreme events at dfferent locatons are drven by geographcal and clmatologcally factors. Unfortunately, the data of these factors are not always avalable. In ths study, the spatal Bayesan herarchcal model (BHM) was employed to update the nformaton n the lkelhood that s not fully descrbed by those unobservable covarates. The proposed method was appled to predct the extreme precptaton return level n Lamongan dstrct, East Java, Indonesa. The Peak Over Threshold (POT) scheme was used to obtan extreme observatons. The pror dstrbuton was employed to update the lkelhood of Generalzed Pareto Dstrbuton as an asymptotc dstrbuton of exceedng resulted from POT procedure. The emprcal results showed that the return level got hgher for longer perods. Keywords: Extreme value theory, Peaks Over Threshold, Bayesan herarchcal model, spatal, return level. M 09 I. INTRODUCTION Extreme value theory (EVT) s a method that s commonly used to analyze extreme natural events. Such extreme events assocated wth the locaton are called spatal extreme value. The ntensve ran that usually happens durng the rany season can cause floods and landsldes. These dsasters are examples of negatves effect caused by extreme events. The studes to learn these extreme phenomena are qute mportant to predct and antcpate ther occurrence n the future. The EVT s used to analyze the heavy-taled dstrbuton of the extreme events. There are two approaches used EVT: () Block Maxma (BM) approach that has the Generalzed Extreme Value (GEV) dstrbuton and () Peaks over Threshold (POT) approach that follows the Generalzed Pareto Dstrbuton (GPD) dstrbuton. These research employed POT scheme. The characterstcs of the GPD are expressed by the three parameters of EVT,.e. threshold ( u ), scale ( ), and shape ( ). One of the most mportant thngs n the EVT s to calculate the value of return level,.e. the probablty of occurrence of a partcular level of extreme events n the comng perod. Knowledge of the return level s mportant for dsaster mtgaton and for the preparaton of long-term programs []. Ths paper dscusses the procedures to obtan the parameters estmates GPD and to calculate the return level usng Bayesan Herarchcal Model (BHM). The BHM s chosen because t has some advantages, for nstance, t able to accommodate geography and clmatology at each level of the herarchy structure. Ths paper s organzed as follows. Ths secton s followed by a secton for a method that brefly dscuss POT and BHM. The next part s used to descrbe the data and research methodology. The last part gves concluson and suggestons. M - 5
2 ISBN II. LITERATURE REVIEW In ths part s dscussed how the procedure was done to model the spatal extreme value by usng the POT approach that follows GPD combned wth BHM. A. Peaks Over Threshold (POT) POT s one of the methods to dentfy the extreme value by usng the so-called threshold (u) as a reference value. The data above the threshold wll be dentfed as extreme values. Fgure shows how to collect extreme data usng POT. The values of, and are larger than the threshold (u) therefore these sx data are consdered as extreme values that wll be used n further analyss u F FIGURE. ILLUSTRATION OF EXTREME DATA COLLECTION USING POT The hgher threshold s the hgher probablty that the extremes wll approach the GPD []. The probablty densty functon (pdf) of GPD s formulated as follow: f x u x u, 0 exp x u, 0 () where 0 x u f 0 and u x u u f 0 parameter and shape parameter (the tal ndex), respectvely.. The and n () represent scale B. Parameter Estmaton of Generalzed Pareto Dstrbuton Parameter estmaton of GPD can be done usng several ways. One of them commonly used s a Maxmum Lkelhood Estmaton (MLE) that maxmzes the lkelhood functon of y, y,, y n, where y x u. Therefore, the equaton () s rewrtten n () as follows: f,, y The log-lkelhood functon of () s expressed as: y, 0 exp y, 0 n y ln L, y, y,, y n ln ln n () (3) M - 5
3 PROCEEDING OF 3 RD INTERNATIONAL CONFERENCE ON RESEARCH, IMPLEMENTATION AND EDUCATION OF MATHEMATICS AND SCIENCE YOGYAKARTA, 6 7 MAY 06 The frst dervatve of log-lkelhood wth resct to each parameters are: ln L n y n y ln y (4) ln L n y n y The soluton for shape parameter estmate usng MLE s not close form. The Newton-Raphson was employed to solve ths ssue. C. Determnaton of the Threshold The extreme values were obtaned once the threshold (u) was determned as a reference value. The determnaton of the threshold ams to fnd the optmum balance n order to obtan the mnmum error. In ths paper, the Mean Resdual Lfe Plot (MRLP) was used to determne the threshold wth the followng steps.. Frst, makng the MRP wth the coordnates of the ponts based on the followng equaton, n u u, x u : u x n maks u (6) where n u s the number of observatons above the threshold. The lnear functon of u s: u uu u 0 0 E X u X u where s extreme threshold value and s the scale parameter changes. (5), (7). Chosen value threshold pont when the plot began to form a lnear pattern around the value of u. D. Autocorrelaton Functon (ACF) In ths paper, the ACF was employed to detect f the extreme values larger than the threshold has a stochastc nature. Stochastc nature s random nature that can only be explaned by a dstrbuton [3]. Therefore, when the ACF s not sgnfcant or statonary, t s necessary to declusterng data that take the extreme values of the hghest n the range (r) determned n advance usng the extremal ndex at (8) [4]. N T mn, f max T : N N N T E. Clusterng the Locaton N T N N T T mn, f max T : N Groupng of locatons used to obtan shape parameter that estmated from statonary process. In ths paper, f the slope values that are not statonary, then the groupngs of locatons wth non-herarchcal clusterng method,.e. k-mean cluster, was employed. The number of groups was determned usng the average value of slhouette coeffcent calculated as n (9) [5]: S j B j k O j O j B j k mn, (9) max,, (8) M - 53
4 ISBN where denotes the average dstance the members j wth other members from dfferent cluster and descrbes the average dstance between the member j to other members n the same cluster. F. Bayesan Herarchcal Model (BHM) The purpose of usng BHM s to accommodate geographcal and clmatologcal nformaton nto the dstrbuton of parameters and to update the parameters through the data that were known as a posteror dstrbuton. The posteror dstrbuton was determned by the followng Bayes theorem. z k, j, c k, j y k, j, (0) p θ, θ Z s p Z s θ p θ θ p 3 θ, () where c k, js the constant transformaton of resultng the cluster, z k, j, s the result of the transformaton n the cluster k, staton j, and observaton. The s jont probablty (or lkelhood), s condtonal where ϕ, and are scale and shape parameters n lkelhood. The and are parameters n layer process. pror dstrbuton, and s pror dstrbuton. The parameters n () are θ, and θ, III. DATA AND RESEARCH METHODOLOGY Ths secton dscusses how to obtan estmates of the parameters and return level for the modelng of extreme precptaton events usng BHM. The softwares used n ths work are openbugs, Mcrosoft Excel, and R. The proposed method was appled to model ranfall data observed n sx locatons n Lamongan dstrct, East Java, Indonesa. A. Data and Pre-Processng Pre-processng the data used to dentfy the mssng value, outlers, and the observatons that do not ft the requrement. The Mcrosoft Excel was used for pre-processng and open-source software R was used to calculate descrptve statstcs such as hstograms and normalty plot, see [6]. TABLE. DATA STRUCTURE Year Month Staton S S S Observaton u v u v u v 98 x, x, x, 98 x, x, x, x,n x,n x,n B. Samplng Extreme Value of the Data. Installng EVA packages n R used to analyze extreme value.. Actvatng extreme toolkt and create MRLP for a specfed nterval. 3. Identfyng the mean resdual lfe plot, fnd upper and lower bounds that are stable on the plot. Use the lower and upper lmts of the modfed parameters to make the plot of scale and shape parameters. Next, select the approprate threshold values based on the results of the modfed parameter plot. 4. Collectng extreme data above the threshold value. 5. The data obtaned were tested wth ACF plot. When there s no sgnfcant lag (statonary), then test the GPD usng Anderson-Darlng test [3]. If t s not statonary yet, then the declusterng must be done. 6. A groupng staton usng the k-mean method and fll t wth the value Slhouette as n (9) to obtan a proper number of the cluster. C. Parameter Estmaton wth BHM. Complng the Drected Acyclc Graph (DAG) to see the relatonshp between the data, parameters of the model, and the pror dstrbuton usng openbugs software [7].. Determnng the ntal value for each parameter to be estmated. M - 54
5 Autocorrelaton Shape Mean Excess Modfed Scale PROCEEDING OF 3 RD INTERNATIONAL CONFERENCE ON RESEARCH, IMPLEMENTATION AND EDUCATION OF MATHEMATICS AND SCIENCE YOGYAKARTA, 6 7 MAY Determnng pror for,,,0,, and as the parameters n layer process and lkelhood. D. Applcaton to Real Data The steps descrbed before were appled to ranfall data n Lamongan dstrct observed n seven locatons. The data used s daly ranfall data span from 98 to 03. The frst step done s determnng the threshold value. Input the data nto Mcrosoft Excel CSV (MC-DOS) and comple t nto extreme toolkt already actvated by typng lbrary (extremes) n R, as shown Fgure. The MRLP, modfed shape, scale and all about extreme value usng extremes toolkt were produced as dsplayed n Fgure 3. FIGURE. EXTREME TOOLKIT Mean Resdual Lfe Plot: runnng V Threshold u FIGURE 3. MRLP (LEFT) AND MODIFIED SCALE AND SHAPE PARAMETER PLOT (RIGHT) As shown n Fgure 3, the threshold lmt (u) 60 began to show nstablty. Creatng the ACF plot and conductng Anderson Darlng s test (usng EVA package n R) determne f the extreme data comples the dstrbuton,.e. GPD. If the test ACF nsgnfcant and Anderson darlng s test show that the GPD s satsfed, then the next step s dong declusterng to take samples agan. But, f the test does not meet the GPD, t s necessary to re-determne the threshold. Threshold Autocorrelaton Functon for V6 (wth 5% sgnfcance lmts for the autocorrelatons) Lag LOCATIONS STATISTIC P-VALUE REMARKS BABAT NOT REJECT H0 BLAWI NOT REJECT H0 KARANGBINANGUN NOT REJECT H0 KEDUNGPRING NOT REJECT H0 LAMONGAN NOT REJECT H0 PANGKATREJO NOT REJECT H0 SUKODADI NOT REJECT H0 FIGURE 4. ACF PLOT (LEFT) AND THE RESULT OF ANDERSON DARLING TEST (RIGHT) M - 55
6 ISBN The ACF plot for Blaw staton s dsplayed n Fgure 3 (left). The plot shows that the extreme observatons are ndependent. The Anderson Darlng test wth null hypothess assumes that the extreme observaton follows GPD. Ths test was employed usng R software wth command gpdad(y) from EVA package. The observaton n all locatons follows GPD as P-values are larger than Type-I error 5%. The next step s estmatng the scale and shape parameter by means of MLE. The results are summarzed n Table. TABLE. ESTIMATES OF PARAMETERS LOCATIONS Scale Shape BABAT BLAWI KARANGBINANGUN KEDUNGPRING LAMONGAN PANGKATREJO SUKODADI Once the parameters estmators were obtaned, the next step s applyng BHM (n ths work usng openbugs software) to the data. The condtonal pror s updated by nvolvng the locaton nformaton as follows: s MVN, θ θ μ Σ p ~ () 0 lats longs (3) exp( s s ' ) (4) The mean accommodate the locaton coordnate represented by lattude and longtude whereas the varance-covarance matrx s accommodate the dstances among locatons. The equaton (3) and (4) are the mean and varance of the process layer () to fnd the posteror dstrbuton n(). Fgure 5 shows the posteror densty of scale parameters. FIGURE 5. POSTERIOR DENSITY OF SCALE PARAMETERS TABLE 3. PARAMETERS ESTIMATOR OBTAINED FROM BHM APPROACH Param. mean Sd MC_error val.5pc medan val97.5pc Start sample Beta Beta alpha alpha alpha scale[] scale[] scale[3] scale[4] M - 56
7 PROCEEDING OF 3 RD INTERNATIONAL CONFERENCE ON RESEARCH, IMPLEMENTATION AND EDUCATION OF MATHEMATICS AND SCIENCE YOGYAKARTA, 6 7 MAY 06 scale[5] scale[6] scale[7] Table 3 shows that the parameters 0, and are not sgnfcant because there s zero value wthn credble nterval whle the others are sgnfcant. The sgnfcant parameters can be used n the model. Thus, t can be concluded that the lattude and longtude varable are sutable for spatal models but not sutable as an effect of predctor models. The model of layer process wth equal to zero s expressed n (5): s MVN 0, θ θ Σ p ~ (5) Havng obtaned a sgnfcant model, the next step s to estmate the scale and shape parameter as reported n Table 4. TABLE 4. ESTIMATOR OF SCALE PARAMETERS AND RETURN LEVEL OBTAINED FROM BHM APPROACH LOCATIONS Return Level (Years) n (x>u) u scale Shape Pr(X>u) BABAT BLAWI KARANGBINANGUN KEDUNGPRING LAMONGAN PANGKATREJO SUKODADI The predctons of return level for 5, 0, 30, 40 and 00 years were summarzed n Table 4. In 08 (fve years return level), the predcted heavest ranfall daly n Babat Staton s 30.9 mm. The return level ncreased for twenty years become mm contnuously ncreased for longer perods. Ths also happened n other sx locatons. IV. CONCLUSION AND SUGGESTIONS Ths paper showed the procedure for obtanng estmates of GPD parameters and calculatng the return level usng a Bayesan Herarchcal models appled to ranfall data n Lamongan. The locaton nformaton represented by longtude and lattude varable were consdered as an nput varable n the model. The emprcal result shows that these two varables are not sgnfcant n the mean model. The return level got hgher for longer perods n all statons. These emprcal fndngs suggest that the future research should explore other clmatology varables as nput n the model. ACKNOWLEDGMENT The authors would lke to thank BMKG Surabaya for gvng permsson to use the data employed n ths paper. REFERENCES [] Cooley, D., Nychka, D., and Naveau, P. Bayesan Spatal Modelng of Extreme Precptaton Return Levels. New York : Journal of the Amerca Assocaton, 007. Vol. 0. [] Gll, M., and Kellez, E. An Applcaton of Extreme Value Theory for Measurng Rsk. Amsterdam: Elseber Scence, 003. [3] We, W. Tme Seres Analyss: Unvarate and Multvarate. New York: Addson-Wesley Publshng Co. 006 [4] Ferro, C., and Seger, J. Inference for Clusters of Extreme Values. Journal R. Stat. Socety, 65, [5] Handoyo, R., Ruman, and Nasuton, S. M. Orygnaly n Indonesa Perbandngan Metode Clusterng Menggunakan Metode Sngle Lnkage dan K - Means Pada Pengelompokan Dokumen. JSM STMIK Mkrosk, 5 (), [6] Evertt, Bran S., and Hothorn, Torsten. A handbook os Statstcal Analyses Usng R. London: Chapman & Hall/CRC, 006. M - 57
8 ISBN [7] Ntzoufras, I. Bayesan Modelng Usng WnBUGS. Greece: WILEY, 009. M - 58
CS 534: Computer Vision Model Fitting
CS 534: Computer Vson Model Fttng Sprng 004 Ahmed Elgammal Dept of Computer Scence CS 534 Model Fttng - 1 Outlnes Model fttng s mportant Least-squares fttng Maxmum lkelhood estmaton MAP estmaton Robust
More informationS1 Note. Basis functions.
S1 Note. Bass functons. Contents Types of bass functons...1 The Fourer bass...2 B-splne bass...3 Power and type I error rates wth dfferent numbers of bass functons...4 Table S1. Smulaton results of type
More informationEmpirical Distributions of Parameter Estimates. in Binary Logistic Regression Using Bootstrap
Int. Journal of Math. Analyss, Vol. 8, 4, no. 5, 7-7 HIKARI Ltd, www.m-hkar.com http://dx.do.org/.988/jma.4.494 Emprcal Dstrbutons of Parameter Estmates n Bnary Logstc Regresson Usng Bootstrap Anwar Ftranto*
More informationParallelism for Nested Loops with Non-uniform and Flow Dependences
Parallelsm for Nested Loops wth Non-unform and Flow Dependences Sam-Jn Jeong Dept. of Informaton & Communcaton Engneerng, Cheonan Unversty, 5, Anseo-dong, Cheonan, Chungnam, 330-80, Korea. seong@cheonan.ac.kr
More informationSLAM Summer School 2006 Practical 2: SLAM using Monocular Vision
SLAM Summer School 2006 Practcal 2: SLAM usng Monocular Vson Javer Cvera, Unversty of Zaragoza Andrew J. Davson, Imperal College London J.M.M Montel, Unversty of Zaragoza. josemar@unzar.es, jcvera@unzar.es,
More informationNAG Fortran Library Chapter Introduction. G10 Smoothing in Statistics
Introducton G10 NAG Fortran Lbrary Chapter Introducton G10 Smoothng n Statstcs Contents 1 Scope of the Chapter... 2 2 Background to the Problems... 2 2.1 Smoothng Methods... 2 2.2 Smoothng Splnes and Regresson
More informationy and the total sum of
Lnear regresson Testng for non-lnearty In analytcal chemstry, lnear regresson s commonly used n the constructon of calbraton functons requred for analytcal technques such as gas chromatography, atomc absorpton
More informationProper Choice of Data Used for the Estimation of Datum Transformation Parameters
Proper Choce of Data Used for the Estmaton of Datum Transformaton Parameters Hakan S. KUTOGLU, Turkey Key words: Coordnate systems; transformaton; estmaton, relablty. SUMMARY Advances n technologes and
More informationUnsupervised Learning and Clustering
Unsupervsed Learnng and Clusterng Why consder unlabeled samples?. Collectng and labelng large set of samples s costly Gettng recorded speech s free, labelng s tme consumng 2. Classfer could be desgned
More informationA Semi-parametric Regression Model to Estimate Variability of NO 2
Envronment and Polluton; Vol. 2, No. 1; 2013 ISSN 1927-0909 E-ISSN 1927-0917 Publshed by Canadan Center of Scence and Educaton A Sem-parametrc Regresson Model to Estmate Varablty of NO 2 Meczysław Szyszkowcz
More informationStatistical analysis on mean rainfall and mean temperature via functional data analysis technique
Statstcal analyss on mean ranfall and mean temperature va functonal data analyss technque Jamaludn Suhala Ctaton: AIP Conference Proceedngs 1613, 368 (01); do: 10.1063/1.89361 Vew onlne: http://dx.do.org/10.1063/1.89361
More informationUnsupervised Learning
Pattern Recognton Lecture 8 Outlne Introducton Unsupervsed Learnng Parametrc VS Non-Parametrc Approach Mxture of Denstes Maxmum-Lkelhood Estmates Clusterng Prof. Danel Yeung School of Computer Scence and
More informationOutline. Type of Machine Learning. Examples of Application. Unsupervised Learning
Outlne Artfcal Intellgence and ts applcatons Lecture 8 Unsupervsed Learnng Professor Danel Yeung danyeung@eee.org Dr. Patrck Chan patrckchan@eee.org South Chna Unversty of Technology, Chna Introducton
More informationX- Chart Using ANOM Approach
ISSN 1684-8403 Journal of Statstcs Volume 17, 010, pp. 3-3 Abstract X- Chart Usng ANOM Approach Gullapall Chakravarth 1 and Chaluvad Venkateswara Rao Control lmts for ndvdual measurements (X) chart are
More informationSupport Vector Machines
/9/207 MIST.6060 Busness Intellgence and Data Mnng What are Support Vector Machnes? Support Vector Machnes Support Vector Machnes (SVMs) are supervsed learnng technques that analyze data and recognze patterns.
More informationModeling Local Uncertainty accounting for Uncertainty in the Data
Modelng Local Uncertanty accontng for Uncertanty n the Data Olena Babak and Clayton V Detsch Consder the problem of estmaton at an nsampled locaton sng srrondng samples The standard approach to ths problem
More informationFEATURE EXTRACTION. Dr. K.Vijayarekha. Associate Dean School of Electrical and Electronics Engineering SASTRA University, Thanjavur
FEATURE EXTRACTION Dr. K.Vjayarekha Assocate Dean School of Electrcal and Electroncs Engneerng SASTRA Unversty, Thanjavur613 41 Jont Intatve of IITs and IISc Funded by MHRD Page 1 of 8 Table of Contents
More informationAnalysis of Malaysian Wind Direction Data Using ORIANA
Modern Appled Scence March, 29 Analyss of Malaysan Wnd Drecton Data Usng ORIANA St Fatmah Hassan (Correspondng author) Centre for Foundaton Studes n Scence Unversty of Malaya, 63 Kuala Lumpur, Malaysa
More informationSolitary and Traveling Wave Solutions to a Model. of Long Range Diffusion Involving Flux with. Stability Analysis
Internatonal Mathematcal Forum, Vol. 6,, no. 7, 8 Soltary and Travelng Wave Solutons to a Model of Long Range ffuson Involvng Flux wth Stablty Analyss Manar A. Al-Qudah Math epartment, Rabgh Faculty of
More informationA Binarization Algorithm specialized on Document Images and Photos
A Bnarzaton Algorthm specalzed on Document mages and Photos Ergna Kavalleratou Dept. of nformaton and Communcaton Systems Engneerng Unversty of the Aegean kavalleratou@aegean.gr Abstract n ths paper, a
More informationFeature Reduction and Selection
Feature Reducton and Selecton Dr. Shuang LIANG School of Software Engneerng TongJ Unversty Fall, 2012 Today s Topcs Introducton Problems of Dmensonalty Feature Reducton Statstc methods Prncpal Components
More informationAn Entropy-Based Approach to Integrated Information Needs Assessment
Dstrbuton Statement A: Approved for publc release; dstrbuton s unlmted. An Entropy-Based Approach to ntegrated nformaton Needs Assessment June 8, 2004 Wllam J. Farrell Lockheed Martn Advanced Technology
More informationA New Approach For the Ranking of Fuzzy Sets With Different Heights
New pproach For the ankng of Fuzzy Sets Wth Dfferent Heghts Pushpnder Sngh School of Mathematcs Computer pplcatons Thapar Unversty, Patala-7 00 Inda pushpndersnl@gmalcom STCT ankng of fuzzy sets plays
More informationSolutions to Programming Assignment Five Interpolation and Numerical Differentiation
College of Engneerng and Coputer Scence Mechancal Engneerng Departent Mechancal Engneerng 309 Nuercal Analyss of Engneerng Systes Sprng 04 Nuber: 537 Instructor: Larry Caretto Solutons to Prograng Assgnent
More informationTECHNIQUE OF FORMATION HOMOGENEOUS SAMPLE SAME OBJECTS. Muradaliyev A.Z.
TECHNIQUE OF FORMATION HOMOGENEOUS SAMPLE SAME OBJECTS Muradalyev AZ Azerbajan Scentfc-Research and Desgn-Prospectng Insttute of Energetc AZ1012, Ave HZardab-94 E-mal:aydn_murad@yahoocom Importance of
More informationSmoothing Spline ANOVA for variable screening
Smoothng Splne ANOVA for varable screenng a useful tool for metamodels tranng and mult-objectve optmzaton L. Rcco, E. Rgon, A. Turco Outlne RSM Introducton Possble couplng Test case MOO MOO wth Game Theory
More information2x x l. Module 3: Element Properties Lecture 4: Lagrange and Serendipity Elements
Module 3: Element Propertes Lecture : Lagrange and Serendpty Elements 5 In last lecture note, the nterpolaton functons are derved on the bass of assumed polynomal from Pascal s trangle for the fled varable.
More informationSubspace clustering. Clustering. Fundamental to all clustering techniques is the choice of distance measure between data points;
Subspace clusterng Clusterng Fundamental to all clusterng technques s the choce of dstance measure between data ponts; D q ( ) ( ) 2 x x = x x, j k = 1 k jk Squared Eucldean dstance Assumpton: All features
More informationWishing you all a Total Quality New Year!
Total Qualty Management and Sx Sgma Post Graduate Program 214-15 Sesson 4 Vnay Kumar Kalakband Assstant Professor Operatons & Systems Area 1 Wshng you all a Total Qualty New Year! Hope you acheve Sx sgma
More informationHermite Splines in Lie Groups as Products of Geodesics
Hermte Splnes n Le Groups as Products of Geodescs Ethan Eade Updated May 28, 2017 1 Introducton 1.1 Goal Ths document defnes a curve n the Le group G parametrzed by tme and by structural parameters n the
More informationImplementation Naïve Bayes Algorithm for Student Classification Based on Graduation Status
Internatonal Journal of Appled Busness and Informaton Systems ISSN: 2597-8993 Vol 1, No 2, September 2017, pp. 6-12 6 Implementaton Naïve Bayes Algorthm for Student Classfcaton Based on Graduaton Status
More informationSome Advanced SPC Tools 1. Cumulative Sum Control (Cusum) Chart For the data shown in Table 9-1, the x chart can be generated.
Some Advanced SP Tools 1. umulatve Sum ontrol (usum) hart For the data shown n Table 9-1, the x chart can be generated. However, the shft taken place at sample #21 s not apparent. 92 For ths set samples,
More informationA Statistical Model Selection Strategy Applied to Neural Networks
A Statstcal Model Selecton Strategy Appled to Neural Networks Joaquín Pzarro Elsa Guerrero Pedro L. Galndo joaqun.pzarro@uca.es elsa.guerrero@uca.es pedro.galndo@uca.es Dpto Lenguajes y Sstemas Informátcos
More informationExercises (Part 4) Introduction to R UCLA/CCPR. John Fox, February 2005
Exercses (Part 4) Introducton to R UCLA/CCPR John Fox, February 2005 1. A challengng problem: Iterated weghted least squares (IWLS) s a standard method of fttng generalzed lnear models to data. As descrbed
More informationSynthesizer 1.0. User s Guide. A Varying Coefficient Meta. nalytic Tool. Z. Krizan Employing Microsoft Excel 2007
Syntheszer 1.0 A Varyng Coeffcent Meta Meta-Analytc nalytc Tool Employng Mcrosoft Excel 007.38.17.5 User s Gude Z. Krzan 009 Table of Contents 1. Introducton and Acknowledgments 3. Operatonal Functons
More information6.854 Advanced Algorithms Petar Maymounkov Problem Set 11 (November 23, 2005) With: Benjamin Rossman, Oren Weimann, and Pouya Kheradpour
6.854 Advanced Algorthms Petar Maymounkov Problem Set 11 (November 23, 2005) Wth: Benjamn Rossman, Oren Wemann, and Pouya Kheradpour Problem 1. We reduce vertex cover to MAX-SAT wth weghts, such that the
More informationA Robust Method for Estimating the Fundamental Matrix
Proc. VIIth Dgtal Image Computng: Technques and Applcatons, Sun C., Talbot H., Ourseln S. and Adraansen T. (Eds.), 0- Dec. 003, Sydney A Robust Method for Estmatng the Fundamental Matrx C.L. Feng and Y.S.
More informationClassifier Selection Based on Data Complexity Measures *
Classfer Selecton Based on Data Complexty Measures * Edth Hernández-Reyes, J.A. Carrasco-Ochoa, and J.Fco. Martínez-Trndad Natonal Insttute for Astrophyscs, Optcs and Electroncs, Lus Enrque Erro No.1 Sta.
More informationFINDING IMPORTANT NODES IN SOCIAL NETWORKS BASED ON MODIFIED PAGERANK
FINDING IMPORTANT NODES IN SOCIAL NETWORKS BASED ON MODIFIED PAGERANK L-qng Qu, Yong-quan Lang 2, Jng-Chen 3, 2 College of Informaton Scence and Technology, Shandong Unversty of Scence and Technology,
More informationUSING GRAPHING SKILLS
Name: BOLOGY: Date: _ Class: USNG GRAPHNG SKLLS NTRODUCTON: Recorded data can be plotted on a graph. A graph s a pctoral representaton of nformaton recorded n a data table. t s used to show a relatonshp
More informationCluster Analysis of Electrical Behavior
Journal of Computer and Communcatons, 205, 3, 88-93 Publshed Onlne May 205 n ScRes. http://www.scrp.org/ournal/cc http://dx.do.org/0.4236/cc.205.350 Cluster Analyss of Electrcal Behavor Ln Lu Ln Lu, School
More informationTN348: Openlab Module - Colocalization
TN348: Openlab Module - Colocalzaton Topc The Colocalzaton module provdes the faclty to vsualze and quantfy colocalzaton between pars of mages. The Colocalzaton wndow contans a prevew of the two mages
More informationFusion Performance Model for Distributed Tracking and Classification
Fuson Performance Model for Dstrbuted rackng and Classfcaton K.C. Chang and Yng Song Dept. of SEOR, School of I&E George Mason Unversty FAIRFAX, VA kchang@gmu.edu Martn Lggns Verdan Systems Dvson, Inc.
More informationMachine Learning 9. week
Machne Learnng 9. week Mappng Concept Radal Bass Functons (RBF) RBF Networks 1 Mappng It s probably the best scenaro for the classfcaton of two dataset s to separate them lnearly. As you see n the below
More informationA Bootstrap Approach to Robust Regression
Internatonal Journal of Appled Scence and Technology Vol. No. 9; November A Bootstrap Approach to Robust Regresson Dr. Hamadu Dallah Department of Actuaral Scence and Insurance Unversty of Lagos Akoka,
More informationDetermining the Optimal Bandwidth Based on Multi-criterion Fusion
Proceedngs of 01 4th Internatonal Conference on Machne Learnng and Computng IPCSIT vol. 5 (01) (01) IACSIT Press, Sngapore Determnng the Optmal Bandwdth Based on Mult-crteron Fuson Ha-L Lang 1+, Xan-Mn
More informationLearning the Kernel Parameters in Kernel Minimum Distance Classifier
Learnng the Kernel Parameters n Kernel Mnmum Dstance Classfer Daoqang Zhang 1,, Songcan Chen and Zh-Hua Zhou 1* 1 Natonal Laboratory for Novel Software Technology Nanjng Unversty, Nanjng 193, Chna Department
More informationLife Tables (Times) Summary. Sample StatFolio: lifetable times.sgp
Lfe Tables (Tmes) Summary... 1 Data Input... 2 Analyss Summary... 3 Survval Functon... 5 Log Survval Functon... 6 Cumulatve Hazard Functon... 7 Percentles... 7 Group Comparsons... 8 Summary The Lfe Tables
More informationEXST7034 Regression Techniques Geaghan Logistic regression Diagnostics Page 1
Logstc regresson Dagnostcs Page 1 1 **********************************************************; 2 *** Logstc Regresson - Dsease outbreak example ***; 3 *** NKNW table 14.3 (Appendx C3) ***; 4 *** Study
More informationHierarchical clustering for gene expression data analysis
Herarchcal clusterng for gene expresson data analyss Gorgo Valentn e-mal: valentn@ds.unm.t Clusterng of Mcroarray Data. Clusterng of gene expresson profles (rows) => dscovery of co-regulated and functonally
More informationThe BGLR (Bayesian Generalized Linear Regression) R- Package. Gustavo de los Campos, Amit Pataki & Paulino Pérez. (August- 2013)
Bostatstcs Department Bayesan Generalzed Lnear Regresson (BGLR) The BGLR (Bayesan Generalzed Lnear Regresson) R- Package By Gustavo de los Campos, Amt Patak & Paulno Pérez (August- 03) (contact: gdeloscampos@gmal.com
More informationLecture 4: Principal components
/3/6 Lecture 4: Prncpal components 3..6 Multvarate lnear regresson MLR s optmal for the estmaton data...but poor for handlng collnear data Covarance matrx s not nvertble (large condton number) Robustness
More informationBioTechnology. An Indian Journal FULL PAPER. Trade Science Inc.
[Type text] [Type text] [Type text] ISSN : 0974-74 Volume 0 Issue BoTechnology 04 An Indan Journal FULL PAPER BTAIJ 0() 04 [684-689] Revew on Chna s sports ndustry fnancng market based on market -orented
More informationBayesian Approach for Fatigue Life Prediction from Field Inspection
Bayesan Approach for Fatgue Lfe Predcton from Feld Inspecton Dawn An, and Jooho Cho School of Aerospace & Mechancal Engneerng, Korea Aerospace Unversty skal@nate.com, jhcho@kau.ac.kr Nam H. Km, and Srram
More informationA CLASS OF TRANSFORMED EFFICIENT RATIO ESTIMATORS OF FINITE POPULATION MEAN. Department of Statistics, Islamia College, Peshawar, Pakistan 2
Pa. J. Statst. 5 Vol. 3(4), 353-36 A CLASS OF TRANSFORMED EFFICIENT RATIO ESTIMATORS OF FINITE POPULATION MEAN Sajjad Ahmad Khan, Hameed Al, Sadaf Manzoor and Alamgr Department of Statstcs, Islama College,
More informationParameter estimation for incomplete bivariate longitudinal data in clinical trials
Parameter estmaton for ncomplete bvarate longtudnal data n clncal trals Naum M. Khutoryansky Novo Nordsk Pharmaceutcals, Inc., Prnceton, NJ ABSTRACT Bvarate models are useful when analyzng longtudnal data
More informationBiostatistics 615/815
The E-M Algorthm Bostatstcs 615/815 Lecture 17 Last Lecture: The Smplex Method General method for optmzaton Makes few assumptons about functon Crawls towards mnmum Some recommendatons Multple startng ponts
More informationUB at GeoCLEF Department of Geography Abstract
UB at GeoCLEF 2006 Mguel E. Ruz (1), Stuart Shapro (2), June Abbas (1), Slva B. Southwck (1) and Davd Mark (3) State Unversty of New York at Buffalo (1) Department of Lbrary and Informaton Studes (2) Department
More informationWhy visualisation? IRDS: Visualization. Univariate data. Visualisations that we won t be interested in. Graphics provide little additional information
Why vsualsaton? IRDS: Vsualzaton Charles Sutton Unversty of Ednburgh Goal : Have a data set that I want to understand. Ths s called exploratory data analyss. Today s lecture. Goal II: Want to dsplay data
More informationPage 0 of 0 SPATIAL INTERPOLATION METHODS
Page 0 of 0 SPATIAL INTERPOLATION METHODS 2018 1. Introducton Spatal nterpolaton s the procedure to predct the value of attrbutes at unobserved ponts wthn a study regon usng exstng observatons (Waters,
More informationMathematics 256 a course in differential equations for engineering students
Mathematcs 56 a course n dfferental equatons for engneerng students Chapter 5. More effcent methods of numercal soluton Euler s method s qute neffcent. Because the error s essentally proportonal to the
More informationParallel Numerics. 1 Preconditioning & Iterative Solvers (From 2016)
Technsche Unverstät München WSe 6/7 Insttut für Informatk Prof. Dr. Thomas Huckle Dpl.-Math. Benjamn Uekermann Parallel Numercs Exercse : Prevous Exam Questons Precondtonng & Iteratve Solvers (From 6)
More informationA Similarity-Based Prognostics Approach for Remaining Useful Life Estimation of Engineered Systems
2008 INTERNATIONAL CONFERENCE ON PROGNOSTICS AND HEALTH MANAGEMENT A Smlarty-Based Prognostcs Approach for Remanng Useful Lfe Estmaton of Engneered Systems Tany Wang, Janbo Yu, Davd Segel, and Jay Lee
More informationTHE THEORY OF REGIONALIZED VARIABLES
CHAPTER 4 THE THEORY OF REGIONALIZED VARIABLES 4.1 Introducton It s ponted out by Armstrong (1998 : 16) that Matheron (1963b), realzng the sgnfcance of the spatal aspect of geostatstcal data, coned the
More informationNUMERICAL SOLVING OPTIMAL CONTROL PROBLEMS BY THE METHOD OF VARIATIONS
ARPN Journal of Engneerng and Appled Scences 006-017 Asan Research Publshng Network (ARPN). All rghts reserved. NUMERICAL SOLVING OPTIMAL CONTROL PROBLEMS BY THE METHOD OF VARIATIONS Igor Grgoryev, Svetlana
More informationResearch on Categorization of Animation Effect Based on Data Mining
MATEC Web of Conferences 22, 0102 0 ( 2015) DOI: 10.1051/ matecconf/ 2015220102 0 C Owned by the authors, publshed by EDP Scences, 2015 Research on Categorzaton of Anmaton Effect Based on Data Mnng Na
More informationEstimating Regression Coefficients using Weighted Bootstrap with Probability
Norazan M R, Habshah Md, A H M R Imon Estmatng Regresson Coeffcents usng Weghted Bootstrap wth Probablty NORAZAN M R, HABSHAH MIDI AND A H M R IMON Faculty of Computer and Mathematcal Scences, Unversty
More informationDetection of hand grasping an object from complex background based on machine learning co-occurrence of local image feature
Detecton of hand graspng an object from complex background based on machne learnng co-occurrence of local mage feature Shnya Moroka, Yasuhro Hramoto, Nobutaka Shmada, Tadash Matsuo, Yoshak Shra Rtsumekan
More informationC2 Training: June 8 9, Combining effect sizes across studies. Create a set of independent effect sizes. Introduction to meta-analysis
C2 Tranng: June 8 9, 2010 Introducton to meta-analyss The Campbell Collaboraton www.campbellcollaboraton.org Combnng effect szes across studes Compute effect szes wthn each study Create a set of ndependent
More informationHelp for Time-Resolved Analysis TRI2 version 2.4 P Barber,
Help for Tme-Resolved Analyss TRI2 verson 2.4 P Barber, 22.01.10 Introducton Tme-resolved Analyss (TRA) becomes avalable under the processng menu once you have loaded and selected an mage that contans
More informationEXTENDED BIC CRITERION FOR MODEL SELECTION
IDIAP RESEARCH REPORT EXTEDED BIC CRITERIO FOR ODEL SELECTIO Itshak Lapdot Andrew orrs IDIAP-RR-0-4 Dalle olle Insttute for Perceptual Artfcal Intellgence P.O.Box 59 artgny Valas Swtzerland phone +4 7
More informationAdaptive Transfer Learning
Adaptve Transfer Learnng Bn Cao, Snno Jaln Pan, Yu Zhang, Dt-Yan Yeung, Qang Yang Hong Kong Unversty of Scence and Technology Clear Water Bay, Kowloon, Hong Kong {caobn,snnopan,zhangyu,dyyeung,qyang}@cse.ust.hk
More informationAnalysis of Continuous Beams in General
Analyss of Contnuous Beams n General Contnuous beams consdered here are prsmatc, rgdly connected to each beam segment and supported at varous ponts along the beam. onts are selected at ponts of support,
More informationSkew Angle Estimation and Correction of Hand Written, Textual and Large areas of Non-Textual Document Images: A Novel Approach
Angle Estmaton and Correcton of Hand Wrtten, Textual and Large areas of Non-Textual Document Images: A Novel Approach D.R.Ramesh Babu Pyush M Kumat Mahesh D Dhannawat PES Insttute of Technology Research
More informationExtension of the Application of Conway-Maxwell-Poisson Models: Analyzing Traffic Crash Data Exhibiting Under-Dispersion
Extenson of the Applcaton of Conway-Maxwell-Posson Models: Analyzng Traffc Crash Data Exhbtng Under-Dsperson Domnque Lord 1 Assstant Professor Zachry Department of Cvl Engneerng Texas A&M Unversty 3136
More informationContent Based Image Retrieval Using 2-D Discrete Wavelet with Texture Feature with Different Classifiers
IOSR Journal of Electroncs and Communcaton Engneerng (IOSR-JECE) e-issn: 78-834,p- ISSN: 78-8735.Volume 9, Issue, Ver. IV (Mar - Apr. 04), PP 0-07 Content Based Image Retreval Usng -D Dscrete Wavelet wth
More information7/12/2016. GROUP ANALYSIS Martin M. Monti UCLA Psychology AGGREGATING MULTIPLE SUBJECTS VARIANCE AT THE GROUP LEVEL
GROUP ANALYSIS Martn M. Mont UCLA Psychology NITP AGGREGATING MULTIPLE SUBJECTS When we conduct mult-subject analyss we are tryng to understand whether an effect s sgnfcant across a group of people. Whether
More informationProgramming in Fortran 90 : 2017/2018
Programmng n Fortran 90 : 2017/2018 Programmng n Fortran 90 : 2017/2018 Exercse 1 : Evaluaton of functon dependng on nput Wrte a program who evaluate the functon f (x,y) for any two user specfed values
More informationA Multivariate Analysis of Static Code Attributes for Defect Prediction
Research Paper) A Multvarate Analyss of Statc Code Attrbutes for Defect Predcton Burak Turhan, Ayşe Bener Department of Computer Engneerng, Bogazc Unversty 3434, Bebek, Istanbul, Turkey {turhanb, bener}@boun.edu.tr
More informationA Comparative Study for Outlier Detection Techniques in Data Mining
A Comparatve Study for Outler Detecton Technques n Data Mnng Zurana Abu Bakar, Rosmayat Mohemad, Akbar Ahmad Department of Computer Scence Faculty of Scence and Technology Unversty College of Scence and
More informationModeling Waveform Shapes with Random Effects Segmental Hidden Markov Models
Modelng Waveform Shapes wth Random Effects Segmental Hdden Markov Models Seyoung Km, Padhrac Smyth Department of Computer Scence Unversty of Calforna, Irvne CA 9697-345 {sykm,smyth}@cs.uc.edu Abstract
More informationPositive Semi-definite Programming Localization in Wireless Sensor Networks
Postve Sem-defnte Programmng Localzaton n Wreless Sensor etworks Shengdong Xe 1,, Jn Wang, Aqun Hu 1, Yunl Gu, Jang Xu, 1 School of Informaton Scence and Engneerng, Southeast Unversty, 10096, anjng Computer
More informationREFRACTIVE INDEX SELECTION FOR POWDER MIXTURES
REFRACTIVE INDEX SELECTION FOR POWDER MIXTURES Laser dffracton s one of the most wdely used methods for partcle sze analyss of mcron and submcron sze powders and dspersons. It s quck and easy and provdes
More informationAPPLICATION OF A COMPUTATIONALLY EFFICIENT GEOSTATISTICAL APPROACH TO CHARACTERIZING VARIABLY SPACED WATER-TABLE DATA
RFr"W/FZD JAN 2 4 1995 OST control # 1385 John J Q U ~ M Argonne Natonal Laboratory Argonne, L 60439 Tel: 708-252-5357, Fax: 708-252-3 611 APPLCATON OF A COMPUTATONALLY EFFCENT GEOSTATSTCAL APPROACH TO
More information3. CR parameters and Multi-Objective Fitness Function
3 CR parameters and Mult-objectve Ftness Functon 41 3. CR parameters and Mult-Objectve Ftness Functon 3.1. Introducton Cogntve rados dynamcally confgure the wreless communcaton system, whch takes beneft
More informationFuzzy Filtering Algorithms for Image Processing: Performance Evaluation of Various Approaches
Proceedngs of the Internatonal Conference on Cognton and Recognton Fuzzy Flterng Algorthms for Image Processng: Performance Evaluaton of Varous Approaches Rajoo Pandey and Umesh Ghanekar Department of
More informationUSING LINEAR REGRESSION FOR THE AUTOMATION OF SUPERVISED CLASSIFICATION IN MULTITEMPORAL IMAGES
USING LINEAR REGRESSION FOR THE AUTOMATION OF SUPERVISED CLASSIFICATION IN MULTITEMPORAL IMAGES 1 Fetosa, R.Q., 2 Merelles, M.S.P., 3 Blos, P. A. 1,3 Dept. of Electrcal Engneerng ; Catholc Unversty of
More informationUser Authentication Based On Behavioral Mouse Dynamics Biometrics
User Authentcaton Based On Behavoral Mouse Dynamcs Bometrcs Chee-Hyung Yoon Danel Donghyun Km Department of Computer Scence Department of Computer Scence Stanford Unversty Stanford Unversty Stanford, CA
More informationSteps for Computing the Dissimilarity, Entropy, Herfindahl-Hirschman and. Accessibility (Gravity with Competition) Indices
Steps for Computng the Dssmlarty, Entropy, Herfndahl-Hrschman and Accessblty (Gravty wth Competton) Indces I. Dssmlarty Index Measurement: The followng formula can be used to measure the evenness between
More informationImproved Methods for Lithography Model Calibration
Improved Methods for Lthography Model Calbraton Chrs Mack www.lthoguru.com, Austn, Texas Abstract Lthography models, ncludng rgorous frst prncple models and fast approxmate models used for OPC, requre
More informationJournal of Chemical and Pharmaceutical Research, 2014, 6(6): Research Article
Avalable onlne www.jocpr.com Journal of Chemcal and Pharmaceutcal Research, 2014, 6(6):2512-2520 Research Artcle ISSN : 0975-7384 CODEN(USA) : JCPRC5 Communty detecton model based on ncremental EM clusterng
More informationA mathematical programming approach to the analysis, design and scheduling of offshore oilfields
17 th European Symposum on Computer Aded Process Engneerng ESCAPE17 V. Plesu and P.S. Agach (Edtors) 2007 Elsever B.V. All rghts reserved. 1 A mathematcal programmng approach to the analyss, desgn and
More informationA Fast Visual Tracking Algorithm Based on Circle Pixels Matching
A Fast Vsual Trackng Algorthm Based on Crcle Pxels Matchng Zhqang Hou hou_zhq@sohu.com Chongzhao Han czhan@mal.xjtu.edu.cn Ln Zheng Abstract: A fast vsual trackng algorthm based on crcle pxels matchng
More informationUnderstanding K-Means Non-hierarchical Clustering
SUNY Albany - Techncal Report 0- Understandng K-Means Non-herarchcal Clusterng Ian Davdson State Unversty of New York, 1400 Washngton Ave., Albany, 105. DAVIDSON@CS.ALBANY.EDU Abstract The K-means algorthm
More informationMachine Learning. Topic 6: Clustering
Machne Learnng Topc 6: lusterng lusterng Groupng data nto (hopefully useful) sets. Thngs on the left Thngs on the rght Applcatons of lusterng Hypothess Generaton lusters mght suggest natural groups. Hypothess
More informationA Clustering Algorithm for Chinese Adjectives and Nouns 1
Clusterng lgorthm for Chnese dectves and ouns Yang Wen, Chunfa Yuan, Changnng Huang 2 State Key aboratory of Intellgent Technology and System Deptartment of Computer Scence & Technology, Tsnghua Unversty,
More informationRStudio for Data Management,
K23166 2015/1/9 17:35 page 4 #4 Usng R and RStudo for Data Management, Statstcal Analyss and Graphcs Second Edton Ncholas J. Horton Department of Mathematcs and Statstcs Amherst College Massachusetts,
More information5.0 Quality Assurance
5.0 Dr. Fred Omega Garces Analytcal Chemstry 25 Natural Scence, Mramar College Bascs of s what we do to get the rght answer for our purpose QA s planned and refers to planned and systematc producton processes
More informationRobust data analysis in innovation project portfolio management
MATEC Web of Conferences 70, 007 (08) SPbWOSCE-07 https://do.org/0.05/matecconf/0870007 Robust data analyss n nnovaton project portfolo management Bors Ttarenko,*, Amr Hasnaou, Roman Ttarenko 3 and Llya
More informationCluster-Based Profile Monitoring in Phase I Analysis. Yajuan Chen. Doctor of Philosophy In Statistics
Cluster-Based Profle Montorng n Phase I Analyss Yajuan Chen Dssertaton submtted to the faculty of the Vrgna Polytechnc Insttute and State Unversty n partal fulfllment of the requrements for the degree
More information