Statistics and Data Analysis
|
|
- Lewis Mathews
- 5 years ago
- Views:
Transcription
1 Simulation of Propenity Scoring Method Dee H. Wu, Ph.D, David M. Thompon, Ph.D., David Bard, Ph.D. Univerity of Oklahoma Health Science Center, Oklahoma City, OK ABSTRACT In certain clinical trial or obervational tudie, proper random aignment of treatment and control group i not alway poible, o that election bia may become an iue. Recent effort to addre iue of nonrandom aignment, including a cla of method known a Propenity Scoring, are alternative to reduce bia in the etimation of treatment effect when aignment i not random. The original technique wa decribed by Roenbaum and Rubin in However, we have found that much of current literature decribe the technique uing a large number of mathematical contruct. The multiplicity of mathematical approache reduce the technique acceibility to uer from different interdiciplinary department. Failure to undertand the tatitical method deter proper uage. We will firt decribe the technique through diagram and imulation to improve teaching and undertanding the material intuitively in SAS. We hope that the introduction permit reader to enhance their own undertanding of the relevant mathematical formulim preented in the literature. INTRODUCTION On returning from a national SAS uer group meeting, we preented the Propenity Scoring Method to our local tudent/faculty SAS group. It wa apparent from the group repone that thoe who did not attend the national meeting were only able to follow part of the complex preentation. We decided to create a imulation and viual preentation of the method in SAS for teaching and to put the problem in the hand of tudent. Some Are highly correlated With treatment aignment Some Are highly correlated With treatment aignment For Simplicity in Drawing (we actually Keep multiple in the problem) A vector that i highly correlated with treatment aignment A vector that i highly correlated with treatment aignment e i diregarded Jut howing two compoite Vector a et which highly impact The treatment election and one that doen t Increae in Probability Of treatment Now generate a propenity Score (a calar function) Figure 1: Overview of the propenity coring method in picture in a nuthell. I II III IV V e I I I II III IV V Group (match ) Now We can Do our analyi On group 1
2 Our firt goal wa to decribe the mathematic purely through image (ee figure 1). The method multiple tep make propenity coring method (PSM) difficult to grap. PSM involve projection of vector, the aumption that a vector of i highly correlated with treatment aignment, the ue of logitic regreion, ome algorithm for matching or tratification (or claification uch a deign tree), and the contruction of general linear model. Only when we created thee picture did intuition about the method became more apparent to u. We generated a tet et baed on 4 example model that we deigned and elected to how the technique trength and weaknee. The imulation procedure i decribed below Method: 1. Generate four imulation et (with X and output Y) 1 A Figure 2: Hitogram of four exemplary ditribution. Although we produced thee ditribution with N=5000 for diplay, we imulated tet cae with N=500 to create larger enitivity in the etimate. The parameter for the four cae were imulated in Matlab but can be eaily run in SAS: cae 1: x1 = GenerateNHumpRand(x,[-4 4],[0 1 0],[0 0 0 ], [1 1 1]); cae 2: x1 = GenerateNHumpRand(x,[-4 4],[ ],[-1 0 1], [ ]); cae 3: x1 = GenerateNHumpRand(x,[-4 4],[ ],[-1 0 1], [ ]); cae 4: x1 = GenerateNHumpRand(x,[-4 4],[ ],[ ], [ ]); 2
3 The input to the function GenerateNHumpRand are: The firt vector [-4 4] i the permiible range from which the imulation draw the x value from a tandard normal ditribution. The econd vector ize indicate the ditribution number of hump, and the fractional multiplier pecify the ize of the obervation The third vector pecifie location for the mean of the ditribution, The fourth vector pecifie tandard deviation. 2. Create 3 (x2, x3, x4) that are uncorrelated with x1, which wa generated above. Although the vector x1 may in reality be a combination of other vector, we reduced x1 to a ingle one for demontration implicity.!"## $$%%"$ &!&"$"&&$'!'%$%%&% ($($($($($ $%%"$$!&)$! $!$'! ($$" "&($$%($ $ &!&"$!&)$$%$ '$ ') ($$" "&$'&'"$%'($!$"&&$! $%$*$ %)$!)%%% ($($$'&'"$'!($ "$'!'%$!$ '*$ %)$$'' ($$%$%'$'!$) ##$%%&% $'!$ ')$!)%%%$''$ ($($($($$) 3. Calculate a propenity core baed a linear combination (weight um to 1) of the X1 to X4, along with a mall random normal error term The variable 'propenity core wa generated by a linear combination of the but could have been generated by any function. propenity = 0.79*x1+0.15*x2+0.05*x *x4+0.3*normal(0,1); 4. Aign to Treatment group baed on the propenity core. Note we ue run a tandard random normal probability with mean 0 and td=1 to aign treatment group (i.e. x1 ha a large influence on etting the treatment group). Note, we ued the invere Probit tranformation to control the treatment/control ratio (proportion) and for generating the proper threhold for claifying the treatment. Notice z1 ha a large influence on the treatment claification. The core we et up typically had range of -2 to Create obervation vector Y for each of the tet cae 3
4 Y= a1*x1 + a2*x2 + a3*x3 + a4*x4 + b*t + c1*t*x1 + c2*t*x2 + c3*t*x3 + c4*t*x *Norm(0,1) where, a1=0.28,a2=0.13,a3=0.23,a4=0.352; b=0.4; c1=c2=c3=0.2,c4=0.02 y y x1 Treat 0 1 x2 Treat 0 1 Figure 3: The treatment et and control group are plotted againt x1. In many application, everal may be aociated with treatment aignment, not jut a ingle one like x1. x1 can repreent a modeled linear combination of. Note that, Figure 3 illutrate heterogeneity of regreion that would challenge conventional ANCOVA method but may not affect PSM. Reult and Dicuion: Ditribut ion Decription Niave Quintile Method Full Model Error Quint Error Full Delta Quint Delta Full - 1 1lump % -6.3% lump % 17.8% lumpbroad % -5.5% aymetric % 30.0% Table 1: We compare the Propenity Quintile method v the Full Model Note that in cae where the propenity core behave well and cover a wide range of value, the propenity quintile method outperformed the Full regreion model: +!",,,!," However, when the propenity core are more aymmetrically ditributed or mixed, the propenity method wa le ueful. A we taught ourelve, the efficacy of the propenity coring technique depended on it ability to create a well-behaved function, and on the matching technique. Matching can be improved with better claification cheme; we performed a imple quintile tratification. Concluion: Our approach to demontrating the Propenity Scoring technique by pictorial decription over mathematical formulation wa better received by faculty and tudent, helped our undertanding, 4
5 and built collaboration. Thi demontration provided a tool that tudent and faculty can ue to better undertand the propenity technique. Thi work can be extended to the multiple varietie of dicriminant/claification cheme available for matching. Regardle, uer of the technique hould beware of it pitfall and challenge once they have undertood the bae method. The Propenity Method in Parimoniou Mathematical Formulim. Now we return to the mathematic, a extracted from our undertanding and from an online earch for the topic. Y i = α + τw + X β + ε i ' i i α, the grand mean, i the average treatment effect (in condition that aignment to the treatment or control group i random) τ i the treatment effect W i = 1 for individual in the treated group. W i = 0 for thoe in the control group. Y i i the oberved outcome of individual i X i include a et of oberved characteritic (), ome of which affection election into treatment. ε i i the error term denoting unoberved characteritic. Z i i an oberved variable that affect election into the treatment. Matching i neceary becaue of dependence between ε i and W i that i due to a et of oberved, X i, that are aociated with election into the treatment. Matching pecifically et: E(ε i W i, X i, Z i ) = E(ε it X i, Z i ). We want to etimate W i Z i (the probability of Wi given Zi) by contructing an independent vector that i orthogonal to it and i therefore free from election or other bia. The matching condition can be decribed a: (Y i1, Y i0 ) W i Z i and provide a vector from which we can calculate unbiaed etimate of treatment effect. 5
6 Appendix: -.#/0! 123 & #9 :+ :8;:/7:/;9; 7!'%$ '&&'$!"'!$)!($ <"'"!""$%&)%!$%'%%'% 4"''&"$%))) % *:/ 97:<7 $!!"! $)!$)!"%)!$!) :8;4#=::7:/;9; #/"!)'$'"'" &"'&$")&!%)!'$'!($ 4%$&!%$&!'$"!$ :8;4#===::7:/;9; #/" $") )!$!% &"!$)($ 4%$&!%$&!'$"!$ : +<+<9 =#$)!"!!!">$)!"%)! $')($ #/* $ ")" ))>$&%'")%* $'($ #/ *$&!'!'%'%'>$%')%*$ '($ #/!*$!"" ')>$!!!*$!"($ #/"*$"!&%!'>$)'%) *)$!($ #/$>$$$ 4*$!''%%'' >$!!!%!*!$&$ 4$>$$$ -.#/0! 44++# #+-03?3@ 6
7 ;2A92+!",,,!,"BC 123 & #9 :+ :8;:/7:/;9; 7'%%%$&''!'&"$))))")%$%($ <"'&"$%)))"&$ %! 4"''&"$%))) % *:/ 97:<7 $)'' &)$)"! $!' '$!) :8;4#=::7:/;9; 4 &$&)&"& &$&)&"&)%$($ $'%"'!$'%"'!%%'$%"($!"$!!''"&!"$!!''"& $!($!%$'" %$'" "$!($ "!&$'")!&$'")% $"%($,4$''%$''%!$&$!,4&$"!& &$"!& "'$"($!,4$!)%$!)%$&$&)"% ",4$!!''$!!'' $! $ )" :8;4#===::7:/;9; 4 $ ") $ ")"$!$ $'% % $'% % %$!$!!&$&")!&$&") "!$"%($!$ $ $$&! "!)$"%!) "!)$"%!) "%%$"($,4$"&)!$"&)!$!$&),4&$'!) &$'!) "'$)&($!,4$'!'$'!'$&$&))' ",4$!!''$!!'' $! $ )" : +<+<9 =#$! &&'"& >$&)!"$!%($ 4*$!&"&) >$')!&%"*!$)$ Table 2: The reult from the 1 lump model (ditribution 1) with the full model (including convariate and treatment effect 7
8 % Dee Wu, David Thompon % Univerity of Oklahoma Health Science Center %macro r; %do i=2 %to 2; /*create working SAS data tep from alump tab-delimited text file*/ %let indata="f:\wu\latet\propenity&i.lump.txt"; data in; infile &indata firtob=2 dlm='09'x; input num propenity x1 x2 x3 x4 a1 a2 a3 a4 b c1 c2 c3 c4 d1 ort1 ort2 ort3 Treat y T2 y2 Y3; goption reet=all; title "Ditribution of Y veru confounder x1 for dataet &i"; proc gplot data=in; ymbol v=plu i=rl; plot y*x1=treat; goption reet=all; title "Ditribution of Y veru confounder x2 for dataet &i"; proc gplot data=in; ymbol v=plu i=rl; plot y*x2=treat; goption reet=all; proc corr data=in; var y x1 x2 x3 x4 propenity; proc univariate data=in; title "Ditribution of Y veru confounder x1 for dataet &i"; var y; hitogram y; /*Thi tep calculate the mean and SD of the outcome Y for the two oberved treatment group. While our initial attempt at dicerning true group memberhip in a mixture model ued thi information, it' le important here.*/ proc mean data=in mean td; cla treat; var y; od output ummary=treatmean; /*Logitic regreion obtain PROPENSITY SCORES in the form of predicted probabilitie that an obervation i aigned to treatment group 1, conditional on it oberved. The regreion function conider jut the (including potential confounder), not the outcome Y.*/ proc logitic data=in; model treat (event='1')= x1 x2 x3 x4; output out=predprob p=phat; 8
9 /*Obervation can be grouped on the bai of the predicted probabilitie that are output from the previou tep. One approach imply dichotomize the grouping. Becaue phat=p(treat=1), predgroup i 1 if phat i high (.5 or higher), and 2 if phat i le than.5*/ data predict1; et predprob; predgroup=2-(phat ge.5); /*A more flexible approach group obervation on the bai of quantile.*/ proc rank data=predprob group=5 out=rank; rank pquintile; var phat; proc freq data=rank; table pquintile*treat; proc ort data=rank; by decending treat; /*Checking aociation among and quantile aignment.*/ title('model x1 x2 x3 x4 = pquintile treat pquintile*treat'); cla pquintile treat; model x1 x2 x3 x4 = pquintile treat pquintile*treat; /*THIS INTERACTION MODEL PRODUCES DATA S1 WITH BETA ESTIMATE*/ title('model y = pquintile treat pquintile*treat'); cla pquintile treat; model y = pquintile treat pquintile*treat / olution; od output parameteretimate=1 (where=(ubtr(parameter,1,5)='treat')); /*THIS NO INTERACTION MODEL THAT ADJUSTS FOR THE PROPENSITY SCORE'S QUINTILE PRODUCES DATA S2 WITH BETA ESTIMATE*/ title('model y = pquintile treat'); cla pquintile treat; model y = pquintile treat / olution; od output parameteretimate=2 (where=(ubtr(parameter,1,5)='treat')); /*THIS NAIVE MODEL MAKES NO ADJUSTMENT AND RECORDS BETA ESTIMATE IN DATA S3*/ title('model y=treat'); cla treat; model y = treat / olution; od output parameteretimate=3 (where=(ubtr(parameter,1,5)='treat')); 9
10 /*WEIGHTING THE NAIVE MODEL BY THE PROPENSITY SCORE*/ weight phat; cla treat; model y=treat / olution; od output parameteretimate=4 (where=(ubtr(parameter,1,5)='treat')); /*TREATING THE PROPENSITY SCORE AS A SURROGATE FOR INFORMATION ON THE COVARIATES*/ cla treat ; model y=treat phat / olution; od output parameteretimate=5 (where=(ubtr(parameter,1,5)='treat')); /*A final method ue IPTW etimator by contructing invere weight from original predicted probabilitie, then uing thee in another logitic regreion.*/ data weight; et predprob; if treat=1 then w=1/phat; if treat=0 then w=1/(1-phat); /*and run the regreion model with IPTW taking place of that might confound expoure of interet.*/ proc glm data=weight order=data; cla treat; model y = treat w / olution; od output parameteretimate=6 (where=(ubtr(parameter,1,5)='treat')); /*INSTEAD OF USING THE IPTW ESTIMATOR AS A SURROGATE FOR THE COVARIATES, USE IT AS A WEIGHT.*/ proc glm data=weight order=data; cla treat; WEIGHT W; model y = treat / olution; od output parameteretimate=7 (where=(ubtr(parameter,1,5)='treat')); /*THIS STEP COLLECTS THE BETA ESTIMATES IN A SINGLE DATASET and applie a format to identify the model that produced the etimate. */ proc format; value methf 1='model y = pquintile treat pquintile*treat' 2='model y = pquintile treat' 3='unadjuted model y=treat' 4='unadjuted model weighted by propenity core' 5='model y=treat phat' 6='IPTW etimator a urrogate for model y=treat w' 7='IPTW etimator a weight in model y=treat' ; data reult; et S7; 10
11 if etimate ne 0; method+1; format method methf.; drop dependent biaed; %end; %mend r; %r; quit; 11
Suggested Answers To Exercises. estimates variability in a sampling distribution of random means. About 68% of means fall
Beyond Significance Teting ( nd Edition), Rex B. Kline Suggeted Anwer To Exercie Chapter. The tatitic meaure variability among core at the cae level. In a normal ditribution, about 68% of the core fall
More informationIf Y is normally Distributed, then and 2 Y Y 10. σ σ
ull Hypothei Significance Teting V. APS 50 Lecture ote. B. Dudek. ot for General Ditribution. Cla Member Uage Only. Chi-Square and F-Ditribution, and Diperion Tet Recall from Chapter 4 material on: ( )
More informationOptimal Coordination of Samples in Business Surveys
Paper preented at the ICES-III, June 8-, 007, Montreal, Quebec, Canada Optimal Coordination of Sample in Buine Survey enka Mach, Ioana Şchiopu-Kratina, Philip T Rei, Jean-Marc Fillion Statitic Canada New
More informationSocial Studies 201 Notes for November 14, 2003
1 Social Studie 201 Note for November 14, 2003 Etimation of a mean, mall ample ize Section 8.4, p. 501. When a reearcher ha only a mall ample ize available, the central limit theorem doe not apply to the
More informationSocial Studies 201 Notes for March 18, 2005
1 Social Studie 201 Note for March 18, 2005 Etimation of a mean, mall ample ize Section 8.4, p. 501. When a reearcher ha only a mall ample ize available, the central limit theorem doe not apply to the
More informationLecture 4 Topic 3: General linear models (GLMs), the fundamentals of the analysis of variance (ANOVA), and completely randomized designs (CRDs)
Lecture 4 Topic 3: General linear model (GLM), the fundamental of the analyi of variance (ANOVA), and completely randomized deign (CRD) The general linear model One population: An obervation i explained
More informationStratified Analysis of Probabilities of Causation
Stratified Analyi of Probabilitie of Cauation Manabu Kuroki Sytem Innovation Dept. Oaka Univerity Toyonaka, Oaka, Japan mkuroki@igmath.e.oaka-u.ac.jp Zhihong Cai Biotatitic Dept. Kyoto Univerity Sakyo-ku,
More informationMINITAB Stat Lab 3
MINITAB Stat 20080 Lab 3. Statitical Inference In the previou lab we explained how to make prediction from a imple linear regreion model and alo examined the relationhip between the repone and predictor
More informationCHAPTER 4 DESIGN OF STATE FEEDBACK CONTROLLERS AND STATE OBSERVERS USING REDUCED ORDER MODEL
98 CHAPTER DESIGN OF STATE FEEDBACK CONTROLLERS AND STATE OBSERVERS USING REDUCED ORDER MODEL INTRODUCTION The deign of ytem uing tate pace model for the deign i called a modern control deign and it i
More informationStandard Guide for Conducting Ruggedness Tests 1
Deignation: E 69 89 (Reapproved 996) Standard Guide for Conducting Ruggedne Tet AMERICA SOCIETY FOR TESTIG AD MATERIALS 00 Barr Harbor Dr., Wet Conhohocken, PA 948 Reprinted from the Annual Book of ASTM
More informationComparing Means: t-tests for Two Independent Samples
Comparing ean: t-tet for Two Independent Sample Independent-eaure Deign t-tet for Two Independent Sample Allow reearcher to evaluate the mean difference between two population uing data from two eparate
More informationGain and Phase Margins Based Delay Dependent Stability Analysis of Two- Area LFC System with Communication Delays
Gain and Phae Margin Baed Delay Dependent Stability Analyi of Two- Area LFC Sytem with Communication Delay Şahin Sönmez and Saffet Ayaun Department of Electrical Engineering, Niğde Ömer Halidemir Univerity,
More informationPARAMETERS OF DISPERSION FOR ON-TIME PERFORMANCE OF POSTAL ITEMS WITHIN TRANSIT TIMES MEASUREMENT SYSTEM FOR POSTAL SERVICES
PARAMETERS OF DISPERSION FOR ON-TIME PERFORMANCE OF POSTAL ITEMS WITHIN TRANSIT TIMES MEASUREMENT SYSTEM FOR POSTAL SERVICES Daniel Salava Kateřina Pojkarová Libor Švadlenka Abtract The paper i focued
More informationAlternate Dispersion Measures in Replicated Factorial Experiments
Alternate Diperion Meaure in Replicated Factorial Experiment Neal A. Mackertich The Raytheon Company, Sudbury MA 02421 Jame C. Benneyan Northeatern Univerity, Boton MA 02115 Peter D. Krau The Raytheon
More informationClustering Methods without Given Number of Clusters
Clutering Method without Given Number of Cluter Peng Xu, Fei Liu Introduction A we now, mean method i a very effective algorithm of clutering. It mot powerful feature i the calability and implicity. However,
More informationSIMPLE LINEAR REGRESSION
SIMPLE LINEAR REGRESSION In linear regreion, we conider the frequency ditribution of one variable (Y) at each of everal level of a econd variable (). Y i known a the dependent variable. The variable for
More informationChapter 12 Simple Linear Regression
Chapter 1 Simple Linear Regreion Introduction Exam Score v. Hour Studied Scenario Regreion Analyi ued to quantify the relation between (or more) variable o you can predict the value of one variable baed
More information1. The F-test for Equality of Two Variances
. The F-tet for Equality of Two Variance Previouly we've learned how to tet whether two population mean are equal, uing data from two independent ample. We can alo tet whether two population variance are
More informationSource slideplayer.com/fundamentals of Analytical Chemistry, F.J. Holler, S.R.Crouch. Chapter 6: Random Errors in Chemical Analysis
Source lideplayer.com/fundamental of Analytical Chemitry, F.J. Holler, S.R.Crouch Chapter 6: Random Error in Chemical Analyi Random error are preent in every meaurement no matter how careful the experimenter.
More informationSuggestions - Problem Set (a) Show the discriminant condition (1) takes the form. ln ln, # # R R
Suggetion - Problem Set 3 4.2 (a) Show the dicriminant condition (1) take the form x D Ð.. Ñ. D.. D. ln ln, a deired. We then replace the quantitie. 3ß D3 by their etimate to get the proper form for thi
More informationStatistical Analysis Using Combined Data Sources Ray Chambers Centre for Statistical and Survey Methodology University of Wollongong
Statitical Analyi Uing Combined Data Source Ray Chamber Centre for Statitical and Survey Methodology Univerity of Wollongong JPSM Preentation, Univerity of Maryland, April 7, 2011 1 What Are The Data?
More informationinto a discrete time function. Recall that the table of Laplace/z-transforms is constructed by (i) selecting to get
Lecture 25 Introduction to Some Matlab c2d Code in Relation to Sampled Sytem here are many way to convert a continuou time function, { h( t) ; t [0, )} into a dicrete time function { h ( k) ; k {0,,, }}
More informationMultipurpose Small Area Estimation
Multipurpoe Small Area Etimation Hukum Chandra Univerity of Southampton, U.K. Ray Chamber Univerity of Wollongong, Autralia Weighting and Small Area Etimation Sample urvey are generally multivariate, in
More informationBy Xiaoquan Wen and Matthew Stephens University of Michigan and University of Chicago
Submitted to the Annal of Applied Statitic SUPPLEMENTARY APPENDIX TO BAYESIAN METHODS FOR GENETIC ASSOCIATION ANALYSIS WITH HETEROGENEOUS SUBGROUPS: FROM META-ANALYSES TO GENE-ENVIRONMENT INTERACTIONS
More informationMULTI-LAYERED LOSSY FINITE LENGTH DIELECTRIC CYLINDIRICAL MODEL OF MAN AT OBLIQUE INCIDENCE
Proceeding 3rd Annual Conference IEEE/EMBS Oct.5-8, 1, Itanbul, TURKEY MULTI-LAYERED LOSSY FINITE LENGTH DIELECTRIC CYLINDIRICAL MODEL OF MAN AT OBLIQUE INCIDENCE S.S. Şeker, B. Yeldiren Boğaziçi Univerity,
More information7.2 INVERSE TRANSFORMS AND TRANSFORMS OF DERIVATIVES 281
72 INVERSE TRANSFORMS AND TRANSFORMS OF DERIVATIVES 28 and i 2 Show how Euler formula (page 33) can then be ued to deduce the reult a ( a) 2 b 2 {e at co bt} {e at in bt} b ( a) 2 b 2 5 Under what condition
More informationEE Control Systems LECTURE 14
Updated: Tueday, March 3, 999 EE 434 - Control Sytem LECTURE 4 Copyright FL Lewi 999 All right reerved ROOT LOCUS DESIGN TECHNIQUE Suppoe the cloed-loop tranfer function depend on a deign parameter k We
More informationEvolutionary Algorithms Based Fixed Order Robust Controller Design and Robustness Performance Analysis
Proceeding of 01 4th International Conference on Machine Learning and Computing IPCSIT vol. 5 (01) (01) IACSIT Pre, Singapore Evolutionary Algorithm Baed Fixed Order Robut Controller Deign and Robutne
More informationAfter the invention of the steam engine in the late 1700s by the Scottish engineer
Introduction to Statitic 22 After the invention of the team engine in the late 1700 by the Scottih engineer Jame Watt, the production of machine-made good became widepread during the 1800. However, it
More informationREPRESENTATION OF ALGEBRAIC STRUCTURES BY BOOLEAN FUNCTIONS. Logic and Applications 2015 (LAP 2015) September 21-25, 2015, Dubrovnik, Croatia
REPRESENTATION OF ALGEBRAIC STRUCTURES BY BOOLEAN FUNCTIONS SMILE MARKOVSKI Faculty of Computer Science and Engineering, S Ciryl and Methodiu Univerity in Skopje, MACEDONIA mile.markovki@gmail.com Logic
More informationPreemptive scheduling on a small number of hierarchical machines
Available online at www.ciencedirect.com Information and Computation 06 (008) 60 619 www.elevier.com/locate/ic Preemptive cheduling on a mall number of hierarchical machine György Dóa a, Leah Eptein b,
More informationEE 4443/5329. LAB 3: Control of Industrial Systems. Simulation and Hardware Control (PID Design) The Inverted Pendulum. (ECP Systems-Model: 505)
EE 4443/5329 LAB 3: Control of Indutrial Sytem Simulation and Hardware Control (PID Deign) The Inverted Pendulum (ECP Sytem-Model: 505) Compiled by: Nitin Swamy Email: nwamy@lakehore.uta.edu Email: okuljaca@lakehore.uta.edu
More informationConfusion matrices. True / False positives / negatives. INF 4300 Classification III Anne Solberg The agenda today: E.g., testing for cancer
INF 4300 Claification III Anne Solberg 29.10.14 The agenda today: More on etimating claifier accuracy Cure of dimenionality knn-claification K-mean clutering x i feature vector for pixel i i- The cla label
More informationNonlinear Single-Particle Dynamics in High Energy Accelerators
Nonlinear Single-Particle Dynamic in High Energy Accelerator Part 6: Canonical Perturbation Theory Nonlinear Single-Particle Dynamic in High Energy Accelerator Thi coure conit of eight lecture: 1. Introduction
More informationAssessing Hidden Bias in the Estimation of Causal Effect in Longitudinal Data by. Using a Matching Estimator with Rosenbaum Bounds.
Aeing Hidden Bia in the Etimation of Caual Effect in Longitudinal Data by Uing a Matching Etimator with Roenbaum Bound Tim Futing Liao Univerity of Illinoi Paper Prepared for the BHPS Conference June 30-July
More informationAPPLICATION OF THE SINGLE IMPACT MICROINDENTATION FOR NON- DESTRUCTIVE TESTING OF THE FRACTURE TOUGHNESS OF NONMETALLIC AND POLYMERIC MATERIALS
APPLICATION OF THE SINGLE IMPACT MICROINDENTATION FOR NON- DESTRUCTIVE TESTING OF THE FRACTURE TOUGHNESS OF NONMETALLIC AND POLYMERIC MATERIALS REN A. P. INSTITUTE OF APPLIED PHYSICS OF THE NATIONAL ACADEMY
More informationμ + = σ = D 4 σ = D 3 σ = σ = All units in parts (a) and (b) are in V. (1) x chart: Center = μ = 0.75 UCL =
Our online Tutor are available 4*7 to provide Help with Proce control ytem Homework/Aignment or a long term Graduate/Undergraduate Proce control ytem Project. Our Tutor being experienced and proficient
More informationDesign spacecraft external surfaces to ensure 95 percent probability of no mission-critical failures from particle impact.
PREFERRED RELIABILITY PAGE 1 OF 6 PRACTICES METEOROIDS & SPACE DEBRIS Practice: Deign pacecraft external urface to enure 95 percent probability of no miion-critical failure from particle impact. Benefit:
More informationAssessing the Discriminatory Power of Credit Scores under Censoring
Aeing the Dicriminatory Power of Credit Score under Cenoring Holger Kraft, Gerald Kroiandt, Marlene Müller Fraunhofer Intitut für Techno- und Wirtchaftmathematik (ITWM) Thi verion: Augut 27, 2003 Abtract:
More informationZ a>2 s 1n = X L - m. X L = m + Z a>2 s 1n X L = The decision rule for this one-tail test is
M09_BERE8380_12_OM_C09.QD 2/21/11 3:44 PM Page 1 9.6 The Power of a Tet 9.6 The Power of a Tet 1 Section 9.1 defined Type I and Type II error and their aociated rik. Recall that a repreent the probability
More informationIntegration Testing - Looking for a Solution to Testing Concurrent Components and Systems
Integration Teting - Looking for a Solution to Teting Concurrent Component and Sytem Integration Teting 23. TAV 1 Andrej Pietchker, PhD Siemen AG, CT SE 1 andrej.pietchker@iemen.com ph: +49 (89) 636 55130
More informationWhat lies between Δx E, which represents the steam valve, and ΔP M, which is the mechanical power into the synchronous machine?
A 2.0 Introduction In the lat et of note, we developed a model of the peed governing mechanim, which i given below: xˆ K ( Pˆ ˆ) E () In thee note, we want to extend thi model o that it relate the actual
More informationAdvanced D-Partitioning Analysis and its Comparison with the Kharitonov s Theorem Assessment
Journal of Multidiciplinary Engineering Science and Technology (JMEST) ISSN: 59- Vol. Iue, January - 5 Advanced D-Partitioning Analyi and it Comparion with the haritonov Theorem Aement amen M. Yanev Profeor,
More informationARTICLE Overcoming the Winner s Curse: Estimating Penetrance Parameters from Case-Control Data
ARTICLE Overcoming the Winner Cure: Etimating Penetrance Parameter from Cae-Control Data Sebatian Zöllner and Jonathan K. Pritchard Genomewide aociation tudie are now a widely ued approach in the earch
More informationCodes Correcting Two Deletions
1 Code Correcting Two Deletion Ryan Gabry and Frederic Sala Spawar Sytem Center Univerity of California, Lo Angele ryan.gabry@navy.mil fredala@ucla.edu Abtract In thi work, we invetigate the problem of
More informationA BATCH-ARRIVAL QUEUE WITH MULTIPLE SERVERS AND FUZZY PARAMETERS: PARAMETRIC PROGRAMMING APPROACH
Mathematical and Computational Application Vol. 11 No. pp. 181-191 006. Aociation for Scientific Reearch A BATCH-ARRIVA QEE WITH MTIPE SERVERS AND FZZY PARAMETERS: PARAMETRIC PROGRAMMING APPROACH Jau-Chuan
More informationLecture 7: Testing Distributions
CSE 5: Sublinear (and Streaming) Algorithm Spring 014 Lecture 7: Teting Ditribution April 1, 014 Lecturer: Paul Beame Scribe: Paul Beame 1 Teting Uniformity of Ditribution We return today to property teting
More informationMolecular Dynamics Simulations of Nonequilibrium Effects Associated with Thermally Activated Exothermic Reactions
Original Paper orma, 5, 9 7, Molecular Dynamic Simulation of Nonequilibrium Effect ociated with Thermally ctivated Exothermic Reaction Jerzy GORECKI and Joanna Natalia GORECK Intitute of Phyical Chemitry,
More informationProblem Set 8 Solutions
Deign and Analyi of Algorithm April 29, 2015 Maachuett Intitute of Technology 6.046J/18.410J Prof. Erik Demaine, Srini Devada, and Nancy Lynch Problem Set 8 Solution Problem Set 8 Solution Thi problem
More informationCHAPTER 6. Estimation
CHAPTER 6 Etimation Definition. Statitical inference i the procedure by which we reach a concluion about a population on the bai of information contained in a ample drawn from that population. Definition.
More informationQuantifying And Specifying The Dynamic Response Of Flowmeters
White Paper Quantifying And Specifying The Dynamic Repone Of Flowmeter DP Flow ABSTRACT The dynamic repone characteritic of flowmeter are often incompletely or incorrectly pecified. Thi i often the reult
More informationCHAPTER 8 OBSERVER BASED REDUCED ORDER CONTROLLER DESIGN FOR LARGE SCALE LINEAR DISCRETE-TIME CONTROL SYSTEMS
CHAPTER 8 OBSERVER BASED REDUCED ORDER CONTROLLER DESIGN FOR LARGE SCALE LINEAR DISCRETE-TIME CONTROL SYSTEMS 8.1 INTRODUCTION 8.2 REDUCED ORDER MODEL DESIGN FOR LINEAR DISCRETE-TIME CONTROL SYSTEMS 8.3
More informationTheoretical Population Biology
Theoretical Population Biology 8 () 3 Content lit available at ScienceDirect Theoretical Population Biology journal homepage: www.elevier.com/locate/tpb On the ize ditribution of private microatellite
More informationASSESSING EXPECTED ACCURACY OF PROBE VEHICLE TRAVEL TIME REPORTS
ASSESSING EXPECTED ACCURACY OF PROBE VEHICLE TRAVEL TIME REPORTS By Bruce Hellinga, 1 P.E., and Liping Fu 2 (Reviewed by the Urban Tranportation Diviion) ABSTRACT: The ue of probe vehicle to provide etimate
More informationRecent progress in fire-structure analysis
EJSE Special Iue: Selected Key Note paper from MDCMS 1 1t International Conference on Modern Deign, Contruction and Maintenance of Structure - Hanoi, Vietnam, December 2007 Recent progre in fire-tructure
More informationObserving Condensations in Atomic Fermi Gases
Oberving Condenation in Atomic Fermi Gae (Term Eay for 498ESM, Spring 2004) Ruqing Xu Department of Phyic, UIUC (May 6, 2004) Abtract Oberving condenation in a ga of fermion ha been another intereting
More informationA Simplified Methodology for the Synthesis of Adaptive Flight Control Systems
A Simplified Methodology for the Synthei of Adaptive Flight Control Sytem J.ROUSHANIAN, F.NADJAFI Department of Mechanical Engineering KNT Univerity of Technology 3Mirdamad St. Tehran IRAN Abtract- A implified
More information[Saxena, 2(9): September, 2013] ISSN: Impact Factor: INTERNATIONAL JOURNAL OF ENGINEERING SCIENCES & RESEARCH TECHNOLOGY
[Saena, (9): September, 0] ISSN: 77-9655 Impact Factor:.85 IJESRT INTERNATIONAL JOURNAL OF ENGINEERING SCIENCES & RESEARCH TECHNOLOGY Contant Stre Accelerated Life Teting Uing Rayleigh Geometric Proce
More informationPredicting Random Effects in Group Randomized Trials. Edward J. Stanek III. PredictingRandomEffects1-stanek.doc 10/29/2003 8:55 AM
Predicting Random Effect in Group Randomized Trial Edward J. Stanek III Abtract The main objective in a group randomized trial i a comparion between treatment. Often, there i interet in the expected repone
More informationMicroblog Hot Spot Mining Based on PAM Probabilistic Topic Model
MATEC Web of Conference 22, 01062 ( 2015) DOI: 10.1051/ matecconf/ 2015220106 2 C Owned by the author, publihed by EDP Science, 2015 Microblog Hot Spot Mining Baed on PAM Probabilitic Topic Model Yaxin
More informationRELIABILITY OF REPAIRABLE k out of n: F SYSTEM HAVING DISCRETE REPAIR AND FAILURE TIMES DISTRIBUTIONS
www.arpapre.com/volume/vol29iue1/ijrras_29_1_01.pdf RELIABILITY OF REPAIRABLE k out of n: F SYSTEM HAVING DISCRETE REPAIR AND FAILURE TIMES DISTRIBUTIONS Sevcan Demir Atalay 1,* & Özge Elmataş Gültekin
More informationMathematical modeling of control systems. Laith Batarseh. Mathematical modeling of control systems
Chapter two Laith Batareh Mathematical modeling The dynamic of many ytem, whether they are mechanical, electrical, thermal, economic, biological, and o on, may be decribed in term of differential equation
More informationRepresentation of a Group of Three-phase Induction Motors Using Per Unit Aggregation Model A.Kunakorn and T.Banyatnopparat
epreentation of a Group of Three-phae Induction Motor Uing Per Unit Aggregation Model A.Kunakorn and T.Banyatnopparat Abtract--Thi paper preent a per unit gregation model for repreenting a group of three-phae
More informationResearch Article Reliability of Foundation Pile Based on Settlement and a Parameter Sensitivity Analysis
Mathematical Problem in Engineering Volume 2016, Article ID 1659549, 7 page http://dxdoiorg/101155/2016/1659549 Reearch Article Reliability of Foundation Pile Baed on Settlement and a Parameter Senitivity
More informationThe variance theory of the mirror effect in recognition memory
Pychonomic Bulletin & Review 001, 8 (3), 408-438 The variance theory of the mirror effect in recognition memory SVERKER SIKSTRÖM Stockholm Univerity, Stockholm, Sweden The mirror effect refer to a rather
More informationSTOCHASTIC GENERALIZED TRANSPORTATION PROBLEM WITH DISCRETE DISTRIBUTION OF DEMAND
OPERATIONS RESEARCH AND DECISIONS No. 4 203 DOI: 0.5277/ord30402 Marcin ANHOLCER STOCHASTIC GENERALIZED TRANSPORTATION PROBLEM WITH DISCRETE DISTRIBUTION OF DEMAND The generalized tranportation problem
More informationStochastic Optimization with Inequality Constraints Using Simultaneous Perturbations and Penalty Functions
Stochatic Optimization with Inequality Contraint Uing Simultaneou Perturbation and Penalty Function I-Jeng Wang* and Jame C. Spall** The John Hopkin Univerity Applied Phyic Laboratory 11100 John Hopkin
More informationAnnex-A: RTTOV9 Cloud validation
RTTOV-91 Science and Validation Plan Annex-A: RTTOV9 Cloud validation Author O Embury C J Merchant The Univerity of Edinburgh Intitute for Atmo. & Environ. Science Crew Building King Building Edinburgh
More informationGiven the following circuit with unknown initial capacitor voltage v(0): X(s) Immediately, we know that the transfer function H(s) is
EE 4G Note: Chapter 6 Intructor: Cheung More about ZSR and ZIR. Finding unknown initial condition: Given the following circuit with unknown initial capacitor voltage v0: F v0/ / Input xt 0Ω Output yt -
More informationLecture 10 Filtering: Applied Concepts
Lecture Filtering: Applied Concept In the previou two lecture, you have learned about finite-impule-repone (FIR) and infinite-impule-repone (IIR) filter. In thee lecture, we introduced the concept of filtering
More informationThe machines in the exercise work as follows:
Tik-79.148 Spring 2001 Introduction to Theoretical Computer Science Tutorial 9 Solution to Demontration Exercie 4. Contructing a complex Turing machine can be very laboriou. With the help of machine chema
More informationALLOCATING BANDWIDTH FOR BURSTY CONNECTIONS
SIAM J. COMPUT. Vol. 30, No. 1, pp. 191 217 c 2000 Society for Indutrial and Applied Mathematic ALLOCATING BANDWIDTH FOR BURSTY CONNECTIONS JON KLEINBERG, YUVAL RABANI, AND ÉVA TARDOS Abtract. In thi paper,
More informationA FUNCTIONAL BAYESIAN METHOD FOR THE SOLUTION OF INVERSE PROBLEMS WITH SPATIO-TEMPORAL PARAMETERS AUTHORS: CORRESPONDENCE: ABSTRACT
A FUNCTIONAL BAYESIAN METHOD FOR THE SOLUTION OF INVERSE PROBLEMS WITH SPATIO-TEMPORAL PARAMETERS AUTHORS: Zenon Medina-Cetina International Centre for Geohazard / Norwegian Geotechnical Intitute Roger
More informationAcceptance sampling uses sampling procedure to determine whether to
DOI: 0.545/mji.203.20 Bayeian Repetitive Deferred Sampling Plan Indexed Through Relative Slope K.K. Sureh, S. Umamahewari and K. Pradeepa Veerakumari Department of Statitic, Bharathiar Univerity, Coimbatore,
More informationSMALL-SIGNAL STABILITY ASSESSMENT OF THE EUROPEAN POWER SYSTEM BASED ON ADVANCED NEURAL NETWORK METHOD
SMALL-SIGNAL STABILITY ASSESSMENT OF THE EUROPEAN POWER SYSTEM BASED ON ADVANCED NEURAL NETWORK METHOD S.P. Teeuwen, I. Erlich U. Bachmann Univerity of Duiburg, Germany Department of Electrical Power Sytem
More informationAutomatic Control Systems. Part III: Root Locus Technique
www.pdhcenter.com PDH Coure E40 www.pdhonline.org Automatic Control Sytem Part III: Root Locu Technique By Shih-Min Hu, Ph.D., P.E. Page of 30 www.pdhcenter.com PDH Coure E40 www.pdhonline.org VI. Root
More informationPerformance Measures for BSkSP-3 with BMChSP-1 as a reference plan
International Journal of Advanced Scientific and Technical Reearch Iue 7 volume 4 July-Aug 2017 Available online on http://www.rpublication.com/ijt/index.html ISSN 2249-9954 Performance Meaure for BSkSP-3
More informationExtending MFM Function Ontology for Representing Separation and Conversion in Process Plant
Downloaded from orbit.dtu.dk on: Oct 05, 2018 Extending MFM Function Ontology for Repreenting Separation and Converion in Proce Plant Zhang, Xinxin; Lind, Morten; Jørgenen, Sten Bay; Wu, Jing; Karnati,
More informationEXTENDED STABILITY MARGINS ON CONTROLLER DESIGN FOR NONLINEAR INPUT DELAY SYSTEMS. Otto J. Roesch, Hubert Roth, Asif Iqbal
EXTENDED STABILITY MARGINS ON CONTROLLER DESIGN FOR NONLINEAR INPUT DELAY SYSTEMS Otto J. Roech, Hubert Roth, Aif Iqbal Intitute of Automatic Control Engineering Univerity Siegen, Germany {otto.roech,
More informationOne Class of Splitting Iterative Schemes
One Cla of Splitting Iterative Scheme v Ciegi and V. Pakalnytė Vilniu Gedimina Technical Univerity Saulėtekio al. 11, 2054, Vilniu, Lithuania rc@fm.vtu.lt Abtract. Thi paper deal with the tability analyi
More informationTRANSITION PROBABILITY MATRIX OF BRIDGE MEMBERS DAMAGE RATING
TRANSITION PROBABILITY MATRIX OF BRIDGE MEMBERS DAMAGE RATING Hirohi Sato and Ryoji Hagiwara 2 Abtract Bridge member damage characteritic were tudied uing the inpection record. Damage can be claified into
More informationMassachusetts Institute of Technology Dynamics and Control II
I E Maachuett Intitute of Technology Department of Mechanical Engineering 2.004 Dynamic and Control II Laboratory Seion 5: Elimination of Steady-State Error Uing Integral Control Action 1 Laboratory Objective:
More informationControl Systems Analysis and Design by the Root-Locus Method
6 Control Sytem Analyi and Deign by the Root-Locu Method 6 1 INTRODUCTION The baic characteritic of the tranient repone of a cloed-loop ytem i cloely related to the location of the cloed-loop pole. If
More informationA Constraint Propagation Algorithm for Determining the Stability Margin. The paper addresses the stability margin assessment for linear systems
A Contraint Propagation Algorithm for Determining the Stability Margin of Linear Parameter Circuit and Sytem Lubomir Kolev and Simona Filipova-Petrakieva Abtract The paper addree the tability margin aement
More informationA NEW LOAD MODEL OF THE PEDESTRIANS LATERAL ACTION
A NEW LOAD MODEL OF THE PEDESTRIANS LATERAL ACTION Fiammetta VENUTI PhD Politecnico di Torino Torino, IT Luca Bruno Aociate Profeor Politecnico di Torino Torino, IT Summary Thi paper propoe a new load
More informationAn estimation approach for autotuning of event-based PI control systems
Acta de la XXXIX Jornada de Automática, Badajoz, 5-7 de Septiembre de 08 An etimation approach for autotuning of event-baed PI control ytem Joé Sánchez Moreno, María Guinaldo Loada, Sebatián Dormido Departamento
More informationA simple construction procedure for resolvable incomplete block designs. for any number of treatments
A imple contruction procedure for reolvable incomplete block deign for any number of treatment By MEENA KHARE and W. T. FEDERER Biometric Unit, Cornell Univerity, Ithaca, New York BU-666-M June, 1979 SUMMARY.
More informationInteraction of Pile-Soil-Pile in Battered Pile Groups under Statically Lateral Load
Interaction of Pile-Soil-Pile in Battered Pile Group under Statically Lateral Load H. Ghaemadeh 1*, M. Alibeikloo 2 1- Aitant Profeor, K. N. Tooi Univerity of Technology 2- M.Sc. Student, K. N. Tooi Univerity
More informationGNSS Solutions: What is the carrier phase measurement? How is it generated in GNSS receivers? Simply put, the carrier phase
GNSS Solution: Carrier phae and it meaurement for GNSS GNSS Solution i a regular column featuring quetion and anwer about technical apect of GNSS. Reader are invited to end their quetion to the columnit,
More informationS_LOOP: SINGLE-LOOP FEEDBACK CONTROL SYSTEM ANALYSIS
S_LOOP: SINGLE-LOOP FEEDBACK CONTROL SYSTEM ANALYSIS by Michelle Gretzinger, Daniel Zyngier and Thoma Marlin INTRODUCTION One of the challenge to the engineer learning proce control i relating theoretical
More informationCHEAP CONTROL PERFORMANCE LIMITATIONS OF INPUT CONSTRAINED LINEAR SYSTEMS
Copyright 22 IFAC 5th Triennial World Congre, Barcelona, Spain CHEAP CONTROL PERFORMANCE LIMITATIONS OF INPUT CONSTRAINED LINEAR SYSTEMS Tritan Pérez Graham C. Goodwin Maria M. Serón Department of Electrical
More informationA Bluffer s Guide to... Sphericity
A Bluffer Guide to Sphericity Andy Field Univerity of Suex The ue of repeated meaure, where the ame ubject are teted under a number of condition, ha numerou practical and tatitical benefit. For one thing
More informationTests of Statistical Hypotheses with Respect to a Fuzzy Set
Modern pplied cience; Vol 8, No 1; 014 IN 1913-1844 E-IN 1913-185 Publihed by Canadian Center of cience and Education Tet of tatitical Hypothee with Repect to a uzzy et P Pandian 1 & D Kalpanapriya 1 1
More informationOn the Isomorphism of Fractional Factorial Designs 1
journal of complexity 17, 8697 (2001) doi:10.1006jcom.2000.0569, available online at http:www.idealibrary.com on On the Iomorphim of Fractional Factorial Deign 1 Chang-Xing Ma Department of Statitic, Nankai
More informationUnavoidable Cycles in Polynomial-Based Time-Invariant LDPC Convolutional Codes
European Wirele, April 7-9,, Vienna, Autria ISBN 978--87-4-9 VE VERLAG GMBH Unavoidable Cycle in Polynomial-Baed Time-Invariant LPC Convolutional Code Hua Zhou and Norbert Goertz Intitute of Telecommunication
More informationPOWER SYSTEM SMALL SIGNAL STABILITY ANALYSIS BASED ON TEST SIGNAL
POWE YEM MALL INAL ABILIY ANALYI BAE ON E INAL Zheng Xu, Wei hao, Changchun Zhou Zheang Univerity, Hangzhou, 37 PChina Email: hvdc@ceezueducn Abtract - In thi paper, a method baed on ome tet ignal (et
More informationEC381/MN308 Probability and Some Statistics. Lecture 7 - Outline. Chapter Cumulative Distribution Function (CDF) Continuous Random Variables
EC38/MN38 Probability and Some Statitic Yanni Pachalidi yannip@bu.edu, http://ionia.bu.edu/ Lecture 7 - Outline. Continuou Random Variable Dept. of Manufacturing Engineering Dept. of Electrical and Computer
More informationUnified Correlation between SPT-N and Shear Wave Velocity for all Soil Types
6 th International Conference on Earthquake Geotechnical Engineering 1-4 ovember 15 Chritchurch, ew Zealand Unified Correlation between SPT- and Shear Wave Velocity for all Soil Type C.-C. Tai 1 and T.
More informationStatistical Moments of Polynomial Dimensional Decomposition
Statitical Moment of Polynomial Dimenional Decompoition Sharif Rahman, M.ASCE 1 Abtract: Thi technical note preent explicit formula for calculating the repone moment of tochatic ytem by polynomial dimenional
More informationWhite Rose Research Online URL for this paper: Version: Accepted Version
Thi i a repoitory copy of Identification of nonlinear ytem with non-peritent excitation uing an iterative forward orthogonal leat quare regreion algorithm. White Roe Reearch Online URL for thi paper: http://eprint.whiteroe.ac.uk/107314/
More information