A MATLAB Implementation of a Halton Sequence-Based GHK Simulator for Multinomial Probit Models *
|
|
- Annis Wade
- 5 years ago
- Views:
Transcription
1 A MATLAB Implementation of a Halton Sequene-Based GHK Simulato fo Multinomial Pobit Models * Jose Miguel M. Abito Depatment of Eonomis National Uvesity of Singapoe Apil 005 * I would like to thank Kut Shmidheiny fo aknowledging my equest fo his MATLAB Eonometis Toolbox and Kenneth Tain fo his online ouse in UC Bekeley. Nonetheless, all eos ae mine.
2 Intodution The populaity of logit-based models is due to the availability of losed-fom expessions fo thei hoie pobabilities. Unfotunately, Logit-based models have impotant limitations patiulaly the inability to handle andom taste vaiation, geneal foms of substitution pattens, and epeated hoie situations with seially oelated eos (Tain, 003. To emedy these defiiees, the Multinomial Pobit (MNP model is often used. Howeve losed-fom hoie pobabilities fo MNP ae not eadily available. Numeial integation methods an be utilized to solve these integals. Howeve, most numeial integation methods ae only feasible and patial to use when the dimension of integation is small enough (< 5 beause of limitations in omputing powe. Theefoe simulation methods an be employed to addess this omputational limitation. This is the pime motivation fo studying simulation-based methods applied to MNP and othe Disete Choie models with analytially intatable hoie pobabilities. In using simulation to estimate the hoie pobabilities of a high-dimensional MNP model, the fist step is to find an appopiate algoithm to use that would povide an estimate with mimal vaiane elative to the tue value of the integal. In geneal, when estimating an integal, vaiane edution tehques suh as Antitheti Monte Calo and Impotane Sampling ae used to intodue some negative oelation o othe e popeties among the pseudo-andom daws fo the simulation. Cude o Non-weighted Monte Calo tehques that ely solely on pseudo-andom geneatos often pefom pooly due to thei low equidistibution popety (Sando and Andas, 004 and high disepany (Bandimate, 00. With this in mind, diffeent simulatos have been used in the MNP ontext. Among these diffeent methods, the GHK (Geweke, 989, 99; Hajivassiliou and MFadden, 998; Keane, 990, 994 simulato is the most widely utilized (Tain, 003 and has been studied and obseved to pefom bette than othe simulatos (Geweke et al., 994; Boh-Supan and Hajivassiliou, 993; Hajivassilious et al., 996. Aside fom Monte Calo based vaiane edution tehques, Quasi-Monte Calo methods an be employed to impove estimation. The main idea behind this is that sine the geneation of pseudo-andom numbes is atually detemisti, the eseahe an mapulate and devise altenative detemisti sequenes of numbes that display bette equidistibution and low disepany popeties. Sando and Tain (00 obseve that (t, m, s-nets and Halton daws pefom bette than independent daws in an appliation to maximum simulated likelihood estimation of a mixed logit model. Moe eently, Sando and Andas (004 show that quasi-monte Calo samples and samples based on othogonal aays have bette pefomane in geneal, ompaed to Cude Monte Calo and Antitheti Monte Calo samples when used with GHK simulatos Howeve, the authos do not examine the quality of these GHK simulatos by applying it in a MNP model and examing whethe estimates using GHK simulatos with diffeent samples diffe sigfiantly. Thus this pape fouses on the patiality of using simulatos based on a patiula Quasi-Monte Calo method.
3 This pape povides a simple and patial inopoation of a quasi-monte Calo sample into Kut Shmidheiny s MNP outine in MATLAB (based on the toolbox used in Shmidheiny [003]. The quasi-monte Calo sample used is the Halton sequene whih is the simplest. Speifially, a andomized Halton sequene is inopoated in the GHK simulato of the toolbox. Though moe sophistiated and effetive andomization tehques fo Halton sequenes ae available (Tuffin, 996; Wang and Hikenell, 000, a simple andomization poedue suggested by Tain (003 is used. Lastly, ompaisons between a pseudo-andom geneato based GHK simulato and a andomized Halton sequene-based GHK simulato is aied out using the tavel hoie data studied in Geene (003 and Geene and Henshe (997. Both GHK simulatos povide simila estimates even with only few simulation uns. This shows that fistly, the GHK simulato itself pefoms well enough even with few simulation uns (given that the numbe of obsevations is suffiiently lage, and seondly, only maginal impovements, if any, an be ahieved with a (ude andomized Halton sequene-based GHK simulato that tend to expend moe estimation time. The next setion povides a bief bakgound on how the GHK simulato is onstuted in geneal settings. In Setion 3, a non-tehal disussion of Halton sequenes is given. Setion 4 biefly disusses the data used followed by the main esults of ou expeiment. The last setion onludes. GHK Simulato The GHK Simulato is based on Impotane Sampling with a eusively tunated multivaiate nomal as its Impotane pdf. The basi idea is to dietly appoximate a etangula pobability (Gouieoux and Monfot, 996. As mentioned, the GHK simulato is the most popula simulato fo MNP (Tain, 003 and has been obseved to pefom bette than othe simulatos (Geweke et al., 994; Boh-Supan and Hajivassiliou, 993; Hajivassilious et al., 996. Conside a geneal model with J altenatives. The utility of peson n fom hoosing altenative j is given by U = + ε ( nj V nj nj fo j =,, J. Afte applying the nomalization (fo identifiation method suggested by Tain (003 and expessing ( as diffeenes against altenative i, we have the following: The main issue is the estimation of hoie pobabilities using simulation. Issues egading the full estimation of the MNP model using simulation (whethe and how to implement Maximum Simulated Likelihood Estimation o Method of Simulated Moments ae not dealt with in this shot pape. 3
4 U = + ε ( nji V nji nji whee ε (,..., = ε ε nji ' N(0, Ω i and is (J- x. The pobability that peson n hooses altenative i will then be given as P = P( V + ε < 0 j i (. nji nji Note that via Choleski deomposition, Ω i = L i L i whee 3 L i = M J 0 M J L L O L 0 0 M JJ and onsequently ( an be eintepeted as U = V + Lη ( i n whee U = ( U n i,..., U ', V = ( Vn i,..., V ' and η n = ( η n,..., η nj '. Note also that η nj nji nji N(0,. The above tansfomation allows us to poeed with simulations using independent daws. onveently as Using Bayes Rule eusively, the hoie pobability fo altenative i an then be ewitten P Vn i = P( η < ( V + η Vn i P( η < η < ( Vn3i + 3η + 3η ( V + P( η3 < η < L 33 η η < V (. The GHK simulato is then alulated as follows (Tain, 003: Vn i Vn i. Calulate P( η < = Φ(.. Daw a value η, denoted as η fom a tunated standad nomal tunated at V n i. 4 3 L i is (J- x (J-. 4 Refe to Tain (003 on how to obtain a daw fom a tunated standad nomal distibution. It is impotant to note that in dawing fom a tunated standad nomal, U(0, daws ae used. 4
5 ( V + η ( V + η 3. Calulate P( η < η = η = Φ(. 4. Daw a value η, denoted as η fom a tunated standad nomal tunated at ( Vn i + η. 5. Calulate ( Vn3i + 3η + 3η ( Vn3i + 3η + 3η P( η3 < η = η η = η = Φ( Repeat steps -5 fo all altenatives exept i. 7. The simulated pobability fo the th daw of η, η J is given as ( Vn i ( V + η ( Vn3i + 3η + 3η P = Φ( Φ( Φ( L Repeat steps -7 fo =,, R whee R is the numbe of epliations fo the GHK simulato. ( R ( 9. Finally, the simulated hoie pobability is P = P. R = This algoithm is implemented in Shmidheiny s MATLAB toolbox. 3 Halton Sequene This setion povides a non-tehal intodution to Halton Sequenes. Fo a moe tehal teatment, onsult the efeenes ited in Sando and Andas (004. Conside a one-dimensional Halton sequene in base. Divide the ut line into equal patitions oesponding to the base hosen (i.e. ut in half. The fist element of the sequene is -. Add - to 0 and 0.5 to get the seond and thid elements of the sequene espetively. Ou sequene now is {0.5, 0.5, 0.75}. Add -3 to zeo and eah element of the sequene in pope ode. Ou sequene is now {0.5, 0.5, 0.75, 0.5, 0.65, 0.375, 0.875}. Repeat the poess until one has the desied numbe of elements. The main ationale fo hoosing a Halton sequene ove a simple pseudo-andom geneato is that we want ou daws to mimi a ufom distibution as muh as possible. A possible way to do this is to intodue some negative oelation between suessive daws so that eah side of the distibution is epesented almost equipopotionally. This is the same woking idea fo using antitheti Monte Calo tehques. With the Halton sequene, eah suessive daw jumps fom one side of the distibution to the othe. One measue of pefomane fo pseudo-andom o detemisti samples that appoximate the ontinuous ufom distibution is Disepany. Mathematially, Disepany is defined as n D( x,..., x = sup Sn ( G nxl x X I m X m 5
6 whee I m is an m-dimensional hypeube, G X is a subspae of I m, and S n (G X is a funtion that ounts the numbe of sample points ontained in eah subspae. Lowe disepany means that ou sample points ae moe dispesed in I m whih is ou desied haateisti. Conside a ompaison of the pseudo-andom geneato of MATLAB (and( and a (andomized Halton sequene in a two-dimensional spae: Figue a: Pseudo-Random Geneato Figue b: (Randomized Halton Sequene By visual inspetion, Figue a has highe disepany ompaed to Figue b. In othe wods, the Halton sequene is moe dispesed aoss the egion ompaed to the simple pseudo-andom geneato and theefoe we expet that the fome will pefom bette than the latte in estimation via simulation. 4 Data and Results Shmidheiny s MATLAB toolbox ontains the data set and applies the MNP outine using the data set as a demo. The data set oiginally omes fom Geene and Henshe (997 and is losely studied in 5 The two-dimensional andomized Halton sequene is of base and 3. See the appendix fo the MATLAB ode. 6
7 Geene (003. The data set ontains 840 obsevations on 4 tavel modes fo 0 individuals. The 4 modes onsist of Ai, Tain, Bus o Ca. The hoie attibutes in the model ae g (genealized ost measue, ttime (teminal waiting time, hin (household inome, and hinai (household inome inteated with ai dummy. Estimation is pefomed using maximum simulated likelihood. The MSLE esults ae given below: Table : MSLE Results A B C D E Constants Ai ( ( ( ( ( Tain *** *** ***.6474***.57948*** ( (0.363 ( ( ( Bus 0.906*** *** *** *** *** (0.988 (0.775 ( (0.538 ( Slope Coeffiients g *** *** *** *** *** ( ( ( ( ( ttime ** ** ** *** *** ( ( ( ( ( hinai ** * ** * ( ( ( ( ( # of Obsevations (Individuals Type of daws and( and( Halton and( Halton # of GHK Repliations Substitution Patten Unestited, Nomalized Unestited, Nomalized Unestited, Nomalized Unestited, Nomalized Unestited, Nomalized Computational Time (in minutes Sigfiane Levels: % (***, 5% (**, 0% (* Note: Nomalization algotihm fo unestited substitution patten is based on Tain (003 and is implemented in Shmidheiny s toolbox 7
8 Models A, B and D ae based on the pseudo-andom geneato and( and Model A is used as the benhmak model. Models C and E utilize daws fom andomized Halton sequenes of k- x dimension whee k is the numbe of altenatives (oesponds to the numbe of bases used and is the numbe of epliations (oesponds to the numbe of elements in eah base. With a few mino exeptions, all of the models have simila estimation esults whih an mean two things. Fist, the (odinay GHK simulato is obust enough to handle extemely small epliations and exhibits the same pefomane as a Halton sequene-based GHK simulato. Seond, ou expeiment is not entiely onlusive fom a theoetial point of view sine thee seems to be some divegene in tems of ou objetive and the measue of pefomane used in the analysis and ompaison of the two samples. Speifially, ou objetive is to see how well the diffeent samples atually estimate the hoie pobabilities but the yadstik of pefomane is based on estimated oeffiients fo the hoie attibutes. We ae atually looking at how well these samples yield an estimate via MSLE whih is affeted by the numbe of obsevations used in the study and not simply on the numbe of epliations. With a lage enough sampled population, the estimate will be athe peise even if the numbe of epliations is small (Tain, 003. Nonetheless fom a patial point of view, if the eseahe has a suffiiently lage numbe of obsevations and the only issue is the numbe of epliations to be used, then an odinay GHK simulato pefoms as well as a Halton-sequene based GHK simulato 6, as mentioned peviously. 5 Conlusion As we have seen in the pevious setion, a GHK simulato pefoms as well as a Halton sequenebased GHK simulato, assuming the eseahe has enough obsevations. Though a andomized Halton sequene is easy to inopoate in any pogamming outine fo MNP models (see the Appendix, a onsideable amount of omputing time is expended duing estimation and thus may not be woth it (again, given that we have a suffiiently lage numbe of obsevations. Nevetheless, othe quasi-monte Calo samples and samples based on othogonal aays an be inopoated. How these samples atually pefom in patial appliations still emains to be studied. 6 Given that the numbe of obsevations is suffiiently lage, using a Halton sequene-based GHK simulato has only maginal impovements with espet to the estimates but takes onsideably moe time to pefom the estimation. Theefoe, on these gounds and as explained above, an odinay GHK simulato would be most likely pefeed in patial appliations. 8
9 Appendix The andomized Halton sequene geneato an be easily inopoated into the GHK simulato by eplaing the funtion that alls and( with a funtion that gives us a andomized Halton sequene of size k- x. 7 The MATLAB ode fo geneating andomized Halton sequenes is as follows: funtion = RandHaltVe(m,n % Retuns a matix of andom Halton daws of size m x n (n daws fom m bases = []; vp = GeneatePime(m; fo i = :m (i,: = RandomHalton(n,vP(i; end funtion = RandomHalton(m,b % Randomizes Halton daws with base b and m elements Halton_nonR = GetHalton(m,b; u = and(,; % The following ensues that ou andomized Halton element is within % (0,: fo i = : m if Halton_nonR(i + u < ; t(i = Halton_nonR(i + u; else t(i = Halton_nonR(i + u - ; end end =t'; whee GetHalton(m,b is a funtion that povides a Halton sequene in base b and with size m x (Bandimate, 00. The funtion RandomHalton(m,bis based on a andomizing poedue suggested by Kenneth Tain in his online ouse in UC Bekeley 8. 7 This ous when we daw a sample point fom a tunated standad nomal
10 Refeenes Bosh-Supan, A. and V. Hajivassiliou (993, Smooth Unbiased Multivaiate Pobability Simulation fo Maximum Likelihood Estimation of Limited Dependent Vaiable Models, Jounal of Eonometis 58, Bandimate, P. (00, Numeial Methods in Finane: A MATLAB -Based Intodution, Wiley & Sons, NY. Geweke, J. (989, Bayesian infeene in eonometi models using Monte Calo integation, Eonometia 57, Geweke, J. (99, Effiient simulation fom the multivaiate nomal and Student-t distibutions subjet to linea onstaints, in E. M. Keamidas, ed., Compute Siene and Statistis: Poeedings of the Twenty-Thid Symposium on the Intefae, pp Intefae Foundation of Noth Ameia In., Faifax. Geweke, J., M. Keane, and D. Runkle (994, Altenative Computational Appoahes to Infeene in the Multinomial Pobit Model, Review of Eonomis and Statistis 76, Gouieoux, C. and A. Monfot (996, Simulation-Based Eonometi Methods, Oxfod Uvesity Pess, NY. Geene, W (003, Eonometi Analysis, 5 th Ed., Pentie Hall, NJ. Geene, W. and D. Henshe (997, Multinomial Logit and Disete Choie Models, in Geene, W., LIMDEP Vesion 7.0 Use s Manual, Revised, Eonometi Softwae In., NY. Hajivassiliou, V. and D. MFadden (998, The method of simulated soes fo the estimation of LDV models, Eonometia 66, Hajivassiliou, V., D. MFadden and P. Ruud (996, Simulation of Multivaiate Nomal Retangle Pobabilities and Thei Deivatives: Theoetial and Computational Results, Jounal of Eonometis 7,
11 Sando, Z. and P. Andas (004, Altenative Sampling Methods fo Estimating Multivaiate Nomal Pobabilities, Jounal of Eonometis 0, Sando, Z. and K. Tain (00, Quasi-Random Simulation of Disete Choie Models, mimeo. Shmidheiny, K. (003, Inome Segegation and Loal Pogessive Taxation: Empiial Evidene fom Switzeland, Uvesity of Ben, mimeo. Tain, K. (003, Disete Choie Methods with Simulation, Cambidge Uvesity Pess, NY. Tuffin, B. (996, On the Use of Low-Disepany Sequenes in Monte Calo Methods, Monte Calo Methods and Appliations, Wang, X. and F. J. Hikenell (000, Randomized Halton Sequenes, Mathematial and Compute Modelling 3,
COMPARING MORE THAN TWO POPULATION MEANS: AN ANALYSIS OF VARIANCE
COMPARING MORE THAN TWO POPULATION MEANS: AN ANALYSIS OF VARIANCE To see how the piniple behind the analysis of vaiane method woks, let us onside the following simple expeiment. The means ( 1 and ) of
More informationSKP-2 ALGORITHM: ON FORMING PART AND MACHINE CLUSTERS SEPARATELY
Poeedings of the 1998 Paifi Confeene on Manufatuing, August 18-20, 1998, Bisbane, Queensland, Austalia SKP-2 ALGORITHM: ON FORMING PART AND MACHINE CLUSTERS SEPARATELY Susanto,S., Kennedy,R.D. and Pie,
More informationConfidence Intervals for the Squared Multiple Semipartial Correlation Coefficient. James Algina. University of Florida. H. J.
Eet Size Conidene Inteval 1 Conidene Intevals o the Squaed Multiple Semipatial Coelation Coeiient by James Algina Univesity o Floida H. J. Keselman Univesity o Manitoba all D. Penield Univesity o Miami
More informationChapter 4. Sampling of Continuous-Time Signals
Chapte 4 Sampling of Continuous-Time Signals 1 Intodution Disete-time signals most ommonly ou as epesentations of sampled ontinuous-time signals. Unde easonable onstaints, a ontinuous-time signal an be
More informationGeneralized Vapor Pressure Prediction Consistent with Cubic Equations of State
Genealized Vapo Pessue Pedition Consistent with Cubi Equations of State Laua L. Petasky and Mihael J. Misovih, Hope College, Holland, MI Intodution Equations of state may be used to alulate pue omponent
More informationSAMPLE LABORATORY SESSION FOR JAVA MODULE B. Calculations for Sample Cross-Section 2
SAMPLE LABORATORY SESSION FOR JAVA MODULE B Calulations fo Sample Coss-Setion. Use Input. Setion Popeties The popeties of Sample Coss-Setion ae shown in Figue and ae summaized below. Figue : Popeties of
More informationAnisotropic 2-D Wavelet Packets and Rectangular Tiling: Theory and Algorithms
Anisotopi -D Wavelet Pakets and Retangula Tiling: Theoy and Algoithms Dan Xu and Minh N. Do Depatment of Eletial and Compute Engineeing and Bekman Institute Univesity of Illinois at Ubana-Champaign Email:
More informationCenter for Advanced Studies in Measurement and Assessment. CASMA Research Report. Using G Theory to Examine Confounded Effects: The Problem of One
Cente fo Advaned Studies in Measuement and Assessment CASMA Reseah Repot Numbe 51 Using G Theoy to Examine Confounded Effets: The Poblem of One Robet L. Bennan 1 Januay 017 1 Robet L. Bennan is E. F. Lindquist
More informationOn the indirect e ect in the Stokes±Helmert method of geoid determination
Jounal of Geodesy (1999) 7: 87±9 On the indiet e et in the Stokes±Helmet method of geoid detemination L. E. SjoÈ beg, H. Nahavandhi oyal Institute of Tehnology, Depatment of Geodesy and Photogammety, S-100
More informationPhotographing a time interval
Potogaping a time inteval Benad Rotenstein and Ioan Damian Politennia Univesity of imisoaa Depatment of Pysis imisoaa Romania benad_otenstein@yaoo.om ijdamian@yaoo.om Abstat A metod of measuing time intevals
More informationMultiple Experts with Binary Features
Multiple Expets with Binay Featues Ye Jin & Lingen Zhang Decembe 9, 2010 1 Intoduction Ou intuition fo the poect comes fom the pape Supevised Leaning fom Multiple Expets: Whom to tust when eveyone lies
More informationNumerical Modeling in Biomedical Systems
Numeial Modeling in Biomedial Systems BME 15:35 Letue 7 9/6/17 Nonlinea Systems Dunn Chapte 5 Nonlinea equations Root finding Baketing methods Open methods Gaphial Bisetion False Position Newton s method
More informationReview for the previous lecture
Review fo the pevious letue Definition: sample spae, event, opeations (union, intesetion, omplementay), disjoint, paiwise disjoint Theoem: ommutatitivity, assoiativity, distibution law, DeMogan s law Pobability
More informationThe Substring Search Problem
The Substing Seach Poblem One algoithm which is used in a vaiety of applications is the family of substing seach algoithms. These algoithms allow a use to detemine if, given two chaacte stings, one is
More informationPhysics 211: Newton s Second Law
Physics 211: Newton s Second Law Reading Assignment: Chapte 5, Sections 5-9 Chapte 6, Section 2-3 Si Isaac Newton Bon: Januay 4, 1643 Died: Mach 31, 1727 Intoduction: Kinematics is the study of how objects
More informationLab #4: Newton s Second Law
Lab #4: Newton s Second Law Si Isaac Newton Reading Assignment: bon: Januay 4, 1643 Chapte 5 died: Mach 31, 1727 Chapte 9, Section 9-7 Intoduction: Potait of Isaac Newton by Si Godfey Knelle http://www.newton.cam.ac.uk/at/potait.html
More informationIMPLEMENTATION OF MUR S ABSORBING BOUNDARIES WITH PERIODIC STRUCTURES TO SPEED UP THE DESIGN PROCESS USING FINITE-DIFFERENCE TIME-DOMAIN METHOD
Pogess In Eletomagnetis Reseah, PIER 58, 101 114, 006 IMPLEMENTATION OF MUR S ABSORBING BOUNDARIES WITH PERIODIC STRUCTURES TO SPEED UP THE DESIGN PROCESS USING FINITE-DIFFERENCE TIME-DOMAIN METHOD G.
More informationCentral Coverage Bayes Prediction Intervals for the Generalized Pareto Distribution
Statistics Reseach Lettes Vol. Iss., Novembe Cental Coveage Bayes Pediction Intevals fo the Genealized Paeto Distibution Gyan Pakash Depatment of Community Medicine S. N. Medical College, Aga, U. P., India
More informationLOCAL VARIABILITY MODELING FOR TEXT-INDEPENDENT SPEAKER VERIFICATION
Odyssey 2014: he Speake and Language eognition Wokshop 16-19 June 2014, Joensuu, Finland LOAL VAIABILIY MODELING FO EX-INDEPENDEN SPEAKE VEIFIAION Liping hen 1, Kong Aik Lee 2, Bin Ma 2, Wu Guo 1, Haizhou
More informationNon-Ideal Gas Behavior P.V.T Relationships for Liquid and Solid:
hemodynamis Non-Ideal Gas Behavio.. Relationships fo Liquid and Solid: An equation of state may be solved fo any one of the thee quantities, o as a funtion of the othe two. If is onsideed a funtion of
More informationShrinkage Testimator in Gamma Type-II Censored Data under LINEX Loss Function
Open Jounal of Statistis, 03, 3, 45-57 http://dx.doi.og/0.436/ojs.03.3408 Published Online August 03 (http://www.sip.og/jounal/ojs) Shinkage Testimato in Gamma Type-II Censoed Data unde LINEX Loss Funtion
More informationAutodesk Robot Structural Analysis Professional - Verification Manual for Italian Codes
Autodesk Robot Stutual Analysis Pofessional VERIFICATIO MAUAL FOR ITALIA CODES Mah 2014 ITRODUCTIO... 1 COCRETE... 2 1. DM 9/1/96 RC COLUMS... 3 VERIFICATIO EXAMPLE 1 - COLUM SUBJECTED TO AXIAL LOAD AD
More informationRed Shift and Blue Shift: A realistic approach
Red Shift and Blue Shift: A ealisti appoah Benhad Rothenstein Politehnia Uniesity of Timisoaa, Physis Dept., Timisoaa, Romania E-mail: benhad_othenstein@yahoo.om Coina Nafonita Politehnia Uniesity of Timisoaa,
More informationDARK MATTER AND THE DYNAMICS OF GALAXIES: A NEWTONIAN APPROACH 1. INTRODUCTION
DARK MATTER AND THE DYNAMICS OF GALAXIES: A NEWTONIAN APPROACH Mugu B. RĂUŢ Coesponding autho: Mugu RĂUŢ, E-mail: m_b_aut@yahoo.om Abstat In this pape I popose a oetion to the well-known Newtonian gavitational
More information3.1 Random variables
3 Chapte III Random Vaiables 3 Random vaiables A sample space S may be difficult to descibe if the elements of S ae not numbes discuss how we can use a ule by which an element s of S may be associated
More informationFrom E.G. Haug Escape Velocity To the Golden Ratio at the Black Hole. Branko Zivlak, Novi Sad, May 2018
Fom E.G. Haug Esape eloity To the Golden Ratio at the Blak Hole Banko Zivlak, bzivlak@gmail.om Novi Sad, May 018 Abstat Esape veloity fom the E.G. Haug has been heked. It is ompaed with obital veloity
More informationOBSTACLE DETECTION USING RING BEAM SYSTEM
OBSTACLE DETECTION USING RING BEAM SYSTEM M. Hiaki, K. Takamasu and S. Ozono Depatment of Peision Engineeing, The Univesity of Tokyo 7-3-1 Hongo, Bunkyo-ku, Tokyo, Japan Abstat: In this pape, we popose
More informationSpecial Relativity in Acoustic and Electromagnetic Waves Without Phase Invariance and Lorentz Transformations 1. Introduction n k.
Speial Relativit in Aousti and Eletomagneti Waves Without Phase Invaiane and Loentz Tansfomations Benhad Rothenstein bothenstein@gmail.om Abstat. Tansfomation equations fo the phsial quantities intodued
More informationarxiv: v2 [physics.data-an] 15 Jul 2015
Limitation of the Least Squae Method in the Evaluation of Dimension of Factal Bownian Motions BINGQIANG QIAO,, SIMING LIU, OUDUN ZENG, XIANG LI, and BENZONG DAI Depatment of Physics, Yunnan Univesity,
More information1. Review of Probability.
1. Review of Pobability. What is pobability? Pefom an expeiment. The esult is not pedictable. One of finitely many possibilities R 1, R 2,, R k can occu. Some ae pehaps moe likely than othes. We assign
More informationOn the Poisson Approximation to the Negative Hypergeometric Distribution
BULLETIN of the Malaysian Mathematical Sciences Society http://mathusmmy/bulletin Bull Malays Math Sci Soc (2) 34(2) (2011), 331 336 On the Poisson Appoximation to the Negative Hypegeometic Distibution
More informationEncapsulation theory: the transformation equations of absolute information hiding.
1 Encapsulation theoy: the tansfomation equations of absolute infomation hiding. Edmund Kiwan * www.edmundkiwan.com Abstact This pape descibes how the potential coupling of a set vaies as the set is tansfomed,
More informationA STUDY OF HAMMING CODES AS ERROR CORRECTING CODES
AGU Intenational Jounal of Science and Technology A STUDY OF HAMMING CODES AS ERROR CORRECTING CODES Ritu Ahuja Depatment of Mathematics Khalsa College fo Women, Civil Lines, Ludhiana-141001, Punjab, (India)
More informationVision Sensor. Vision. (Phase 1) pre-shaping. Actuator. Tactile Sensor. Vision. (Phase 2) shaping. Actuator. Tactile Sensor.
Optimal Gasping using Visual and Tatile Feedbak Akio NAMIKI Masatoshi ISHIKAWA Depatment of Mathematial Engineeing and Infomation Physis Univesity of Tokyo Tokyo 3, Japan namik@k.t.u-tokyo.a.jp Abstat
More informationCorrespondence Analysis & Related Methods
Coespondene Analysis & Related Methods Oveview of CA and basi geometi onepts espondents, all eades of a etain newspape, osstabulated aoding to thei eduation goup and level of eading of the newspape Mihael
More informationValue Distribution of L-Functions with Rational Moving Targets
Advanes in Pue Mathematis 3 3 79-73 Pulished Online Deeme 3 (http://wwwsipog/jounal/apm) http://dxdoiog/436/apm33998 Value Distiution of -Funtions wh ational Moving agets Matthew Cadwell Zhuan Ye * ntelligent
More informationPROCESSING QUALITY OF NEW POTATO CULTIVARS FOLLOWING PROLONGED STORAGE
PRESSIG QUALITY F EW PTAT ULTIVARS FLLWIG PRLGED STRAGE Glenn E. Vogt J.R. Simplot o., aldwell, ID and Ei P. Eldedge and linton. Shok Malheu Expeiment Station egon State Univesity ntaio, R Intodution The
More informationThermodynamic derivations of conditions for chemical equilibrium and of Onsager reciprocal relations for chemical reactors
JOURNAL OF CHEMICAL PHYSICS VOLUME 121, NUMBER 6 8 AUGUST 2004 Themodynami deivations of onditions fo hemial equilibium and of Onsage eipoal elations fo hemial eatos Gian Paolo Beetta a) Univesità di Besia,
More informationTHEORETICAL AND EXPERIMENTAL STUDY ON DROPWISE CONDENSATION IN PLATE HEAT EXCHANGERS
Abstat THEORETICAL AND EXPERIMENTAL STUDY ON DROPWISE CONDENSATION IN PLATE HEAT EXCHANGERS V. Bendt, S. Zunft and H. Mülle-Steinhagen Geman Aeospae Cente (DLR), Stuttgat, Gemany This pape desibes the
More informationAlternative Tests for the Poisson Distribution
Chiang Mai J Sci 015; 4() : 774-78 http://epgsciencecmuacth/ejounal/ Contibuted Pape Altenative Tests fo the Poisson Distibution Manad Khamkong*[a] and Pachitjianut Siipanich [b] [a] Depatment of Statistics,
More informationINTEGRATION OF THE SELF-ORGANIZING MAP AND NEURAL GAS WITH MULTIDIMENSIONAL SCALING
ISSN 39 4X INFORMATION TECHNOLOGY AND CONTROL, 0, Vol.40, No. INTEGRATION OF THE SELF-ORGANIZING MAP AND NEURAL GAS WITH MULTIDIMENSIONAL SCALING Olga Kuasova, Alma Molytė Vilnius Univesity, Institute
More informationMULTILAYER PERCEPTRONS
Last updated: Nov 26, 2012 MULTILAYER PERCEPTRONS Outline 2 Combining Linea Classifies Leaning Paametes Outline 3 Combining Linea Classifies Leaning Paametes Implementing Logical Relations 4 AND and OR
More informationOn the integration of the equations of hydrodynamics
Uebe die Integation de hydodynamischen Gleichungen J f eine u angew Math 56 (859) -0 On the integation of the equations of hydodynamics (By A Clebsch at Calsuhe) Tanslated by D H Delphenich In a pevious
More informationGraphs of Sine and Cosine Functions
Gaphs of Sine and Cosine Functions In pevious sections, we defined the tigonometic o cicula functions in tems of the movement of a point aound the cicumfeence of a unit cicle, o the angle fomed by the
More informationStanford University CS259Q: Quantum Computing Handout 8 Luca Trevisan October 18, 2012
Stanfod Univesity CS59Q: Quantum Computing Handout 8 Luca Tevisan Octobe 8, 0 Lectue 8 In which we use the quantum Fouie tansfom to solve the peiod-finding poblem. The Peiod Finding Poblem Let f : {0,...,
More informationQIP Course 10: Quantum Factorization Algorithm (Part 3)
QIP Couse 10: Quantum Factoization Algoithm (Pat 3 Ryutaoh Matsumoto Nagoya Univesity, Japan Send you comments to yutaoh.matsumoto@nagoya-u.jp Septembe 2018 @ Tokyo Tech. Matsumoto (Nagoya U. QIP Couse
More informationCOMPARATIVE ANALYSIS OF LDPC AND BCH CODES ERROR-CORRECTING CAPABILITIES
5 UDC 621.391 COMPARAIVE ANALYSIS OF LDPC AND BCH CODES ERROR-CORRECING CAPABILIIES Leonid O. Uyvsky, Sehii O. Osyphuk eleommuniation Netwoks Depatment Institute of eleommuniation Systems National ehnial
More informationEstimation of the Correlation Coefficient for a Bivariate Normal Distribution with Missing Data
Kasetsat J. (Nat. Sci. 45 : 736-74 ( Estimation of the Coelation Coefficient fo a Bivaiate Nomal Distibution with Missing Data Juthaphon Sinsomboonthong* ABSTRACT This study poposes an estimato of the
More informationChem 453/544 Fall /08/03. Exam #1 Solutions
Chem 453/544 Fall 3 /8/3 Exam # Solutions. ( points) Use the genealized compessibility diagam povided on the last page to estimate ove what ange of pessues A at oom tempeatue confoms to the ideal gas law
More information1 Fundamental Solutions to the Wave Equation
1 Fundamental Solutions to the Wave Equation Physial insight in the sound geneation mehanism an be gained by onsideing simple analytial solutions to the wave equation. One example is to onside aousti adiation
More informationDissolution of Solid Particles in Liquids: A Shrinking Core Model
Wold Aademy of Siene, Engineeing and Tehnology 5 9 Dissolution of Solid Patiles in Liquids: A Shining oe Model Wei-Lun Hsu, Mon-Jyh Lin, and Jyh-Ping Hsu Astat The dissolution of spheial patiles in liquids
More information6 PROBABILITY GENERATING FUNCTIONS
6 PROBABILITY GENERATING FUNCTIONS Cetain deivations pesented in this couse have been somewhat heavy on algeba. Fo example, detemining the expectation of the Binomial distibution (page 5.1 tuned out to
More informationPlanck Quantization of Newton and Einstein Gravitation
Plank Quantization of Newton and Einstein Gavitation Espen Gaade Haug Nowegian Univesity of Life Sienes Mah 0, 06 Abstat In this pape we ewite the gavitational onstant based on its elationship with the
More informationFractional Zero Forcing via Three-color Forcing Games
Factional Zeo Focing via Thee-colo Focing Games Leslie Hogben Kevin F. Palmowski David E. Robeson Michael Young May 13, 2015 Abstact An -fold analogue of the positive semidefinite zeo focing pocess that
More informationAn analytic calculation method on air gap flux in permanent magnet. brushless DC motor with ironless rotor
Intenational Confeene on Enegy and Envionmental Potetion ICEEP 6 An analyti alulation method on ai gap flux in pemanent magnet bushless DC moto with ionless oto Xinghua Wang,Yaolong Sheng andshugang Zhao,,
More informationDiscrete-Time Immersion and Invariance Adaptive Control of a Slider-crank Mechanism
Pepints of the 9th Wold Congess he Intenational Fedeation of Automati Contol Cape own South Afia August -9 0 Disete-ime Immesion and Invaiane Adaptive Contol of a Slide-an Mehanism Yapa Yalçın Dept of
More informationA Relativistic BCS Theory of S-Wave Superconductivity
A Relativisti BCS Theoy of S-Wave Supeondutivity Damien BERTRAND Jan GOVAERTS John MENDY UC-FYNU / Supeondutivity Goup ouvain-la-neuve - June 4 Outline Ginzbug-andau fundamentals Oiginal theoy Covaiant
More informationSurveillance Points in High Dimensional Spaces
Société de Calcul Mathématique SA Tools fo decision help since 995 Suveillance Points in High Dimensional Spaces by Benad Beauzamy Januay 06 Abstact Let us conside any compute softwae, elying upon a lage
More informationarxiv: v2 [astro-ph] 16 May 2008
New Anomalies in Cosmic Micowave Backgound Anisotopy: Violation of the Isotopic Gaussian Hypothesis in Low-l Modes Shi Chun, Su and M.-C., Chu Depatment of Physics and Institute of Theoetical Physics,
More informationTemporal-Difference Learning
.997 Decision-Making in Lage-Scale Systems Mach 17 MIT, Sping 004 Handout #17 Lectue Note 13 1 Tempoal-Diffeence Leaning We now conside the poblem of computing an appopiate paamete, so that, given an appoximation
More informationClassical Worm algorithms (WA)
Classical Wom algoithms (WA) WA was oiginally intoduced fo quantum statistical models by Pokof ev, Svistunov and Tupitsyn (997), and late genealized to classical models by Pokof ev and Svistunov (200).
More informationReflectance spectra for Si
Refletane speta fo Si Notie R and ε i and ε show onsideable stutues in the fom of peas and shouldes. These stutues aise fom the optial tansitions between alene bands to the ondution bands. 16 Miosopi Theoy:
More informationMultiple Criteria Secretary Problem: A New Approach
J. Stat. Appl. Po. 3, o., 9-38 (04 9 Jounal of Statistics Applications & Pobability An Intenational Jounal http://dx.doi.og/0.785/jsap/0303 Multiple Citeia Secetay Poblem: A ew Appoach Alaka Padhye, and
More informationExtra Examples for Chapter 1
Exta Examples fo Chapte 1 Example 1: Conenti ylinde visomete is a devie used to measue the visosity of liquids. A liquid of unknown visosity is filling the small gap between two onenti ylindes, one is
More informationThe Research of AQI Index Changing Regularity Mainly in Tianjin Ziyu Guo
nd Intenational Confeene on Eduation Tehnology, Management and Humanities Siene (ETMHS 06 The Reseah of AQI Index Changing Regulaity Mainly in Tianjin Ziyu Guo Shool of Institute of Eletial and Eletoni
More informationCSCE 478/878 Lecture 4: Experimental Design and Analysis. Stephen Scott. 3 Building a tree on the training set Introduction. Outline.
In Homewok, you ae (supposedly) Choosing a data set 2 Extacting a test set of size > 3 3 Building a tee on the taining set 4 Testing on the test set 5 Repoting the accuacy (Adapted fom Ethem Alpaydin and
More informationFZX: Personal Lecture Notes from Daniel W. Koon St. Lawrence University Physics Department CHAPTER 7
FZX: Pesonal Lectue Notes fom Daniel W. Koon St. Lawence Univesity Physics Depatment CHAPTER 7 Please epot any glitches, bugs o eos to the autho: dkoon at stlawu.edu. 7. Momentum and Impulse Impulse page
More informationPearson s Chi-Square Test Modifications for Comparison of Unweighted and Weighted Histograms and Two Weighted Histograms
Peason s Chi-Squae Test Modifications fo Compaison of Unweighted and Weighted Histogams and Two Weighted Histogams Univesity of Akueyi, Bogi, v/noduslód, IS-6 Akueyi, Iceland E-mail: nikolai@unak.is Two
More informationworking pages for Paul Richards class notes; do not copy or circulate without permission from PGR 2004/11/3 10:50
woking pages fo Paul Richads class notes; do not copy o ciculate without pemission fom PGR 2004/11/3 10:50 CHAPTER7 Solid angle, 3D integals, Gauss s Theoem, and a Delta Function We define the solid angle,
More informationMethod for Approximating Irrational Numbers
Method fo Appoximating Iational Numbes Eic Reichwein Depatment of Physics Univesity of Califonia, Santa Cuz June 6, 0 Abstact I will put foth an algoithm fo poducing inceasingly accuate ational appoximations
More informationIntroduction to Mathematical Statistics Robert V. Hogg Joeseph McKean Allen T. Craig Seventh Edition
Intoduction to Mathematical Statistics Robet V. Hogg Joeseph McKean Allen T. Caig Seventh Edition Peason Education Limited Edinbugh Gate Halow Essex CM2 2JE England and Associated Companies thoughout the
More informationTOTAL VARIANCE AS AN EXACT ANALYSIS OF THE SAMPLE VARIANCE*
29th Annual Peise Time and Time Inteval (PTTI) Meeting V TOTAL VARIANCE AS AN EXACT ANALYSIS OF THE SAMPLE VARIANCE* b I F Donald B. Peival Applied Physis Laboatoy Univesity of Washington Box 355640 Seattle,
More informationAPPENDIX D COMPRESSIBILITY FACTOR EQUATIONS D.1 THE REDLICH KWONG EQUATION
AENDIX D COMRESSIBILIY FACOR EQUAIONS D.1 HE REDLICH KWONG EQUAION he Redlih-Kwong equation is atually an equation of state. It was fomulated by Otto Redlih and Joseph N. S. Kwong in 1949 [Chemial Review
More informationThe Persistence of Most Probable Explanations in Bayesian Networks
EAI 2014 T. Schaub et al. (Eds.) 2014 The Authos and IOS Pess. This aticle is published online with Open Access by IOS Pess and distibuted unde the tems of the eative ommons Attibution Non-ommecial License.
More informationInformation Retrieval Advanced IR models. Luca Bondi
Advanced IR models Luca Bondi Advanced IR models 2 (LSI) Pobabilistic Latent Semantic Analysis (plsa) Vecto Space Model 3 Stating point: Vecto Space Model Documents and queies epesented as vectos in the
More informationPROBLEM SET #1 SOLUTIONS by Robert A. DiStasio Jr.
POBLM S # SOLUIONS by obet A. DiStasio J. Q. he Bon-Oppenheime appoximation is the standad way of appoximating the gound state of a molecula system. Wite down the conditions that detemine the tonic and
More informationdp p v= = ON SHOCK WAVES AT LARGE DISTANCES FROM THE PLACE OF THEIR ORIGIN By Lev D. Landau J. Phys. U.S.S.R. 9, 496 (1945).
ON SHOCK WAVES AT LARGE DISTANCES FROM THE PLACE OF THEIR ORIGIN By Lev D. Landau J. Phys. U.S.S.R. 9, 496 (1945). It is shown that at lage distanes fom the body, moving with a. veloity exeeding that of
More informationF-IF Logistic Growth Model, Abstract Version
F-IF Logistic Gowth Model, Abstact Vesion Alignments to Content Standads: F-IFB4 Task An impotant example of a model often used in biology o ecology to model population gowth is called the logistic gowth
More informationNumerical Integration
MCEN 473/573 Chapte 0 Numeical Integation Fall, 2006 Textbook, 0.4 and 0.5 Isopaametic Fomula Numeical Integation [] e [ ] T k = h B [ D][ B] e B Jdsdt In pactice, the element stiffness is calculated numeically.
More information1) (A B) = A B ( ) 2) A B = A. i) A A = φ i j. ii) Additional Important Properties of Sets. De Morgan s Theorems :
Additional Impotant Popeties of Sets De Mogan s Theoems : A A S S Φ, Φ S _ ( A ) A ) (A B) A B ( ) 2) A B A B Cadinality of A, A, is defined as the numbe of elements in the set A. {a,b,c} 3, { }, while
More informationTime Dilation in Gravity Wells
Time Dilation in Gavity Wells By Rihad R. Shiffman Digital Gaphis Asso. 038 Dunkik Ave. L.A., Ca. 9005 s@isi.edu This doument disusses the geneal elativisti effet of time dilation aused by a spheially
More information10/04/18. P [P(x)] 1 negl(n).
Mastemath, Sping 208 Into to Lattice lgs & Cypto Lectue 0 0/04/8 Lectues: D. Dadush, L. Ducas Scibe: K. de Boe Intoduction In this lectue, we will teat two main pats. Duing the fist pat we continue the
More informationA near-perfect invisibility cloak constructed with homogeneous materials
A nea-pefet invisibility loak onstuted with homogeneous mateials Wei Li, Jianguo Guan, * Zhigang Sun, Wei Wang, and Qingjie Zhang 1 State Key Lab of Advaned Tehnology fo Mateials Synthesis and Poessing,
More information1 Fundamental Solutions to the Wave Equation
1 Fundamental Solutions to the Wave Equation Physial insight in the sound geneation mehanism an be gained by onsideing simple analytial solutions to the wave equation One example is to onside aousti adiation
More informationA Bijective Approach to the Permutational Power of a Priority Queue
A Bijective Appoach to the Pemutational Powe of a Pioity Queue Ia M. Gessel Kuang-Yeh Wang Depatment of Mathematics Bandeis Univesity Waltham, MA 02254-9110 Abstact A pioity queue tansfoms an input pemutation
More informationPrediction and Correlations of Residual Entropy of Superheated Vapor for Pure Compounds
Iaqi Jounal of Chemial and etoleum Engineeing Iaqi Jounal of Chemial and etoleum Engineeing Vol. No. (June 0) - 7 ISSN: 997-4884 Univesity of Baghdad College of Engineeing edition and Coelations of esidual
More informationLecture 8 - Gauss s Law
Lectue 8 - Gauss s Law A Puzzle... Example Calculate the potential enegy, pe ion, fo an infinite 1D ionic cystal with sepaation a; that is, a ow of equally spaced chages of magnitude e and altenating sign.
More informationGENLOG Multinomial Loglinear and Logit Models
GENLOG Multinomial Loglinea and Logit Models Notation This chapte descibes the algoithms used to calculate maximum-likelihood estimates fo the multinomial loglinea model and the multinomial logit model.
More informationA Multivariate Normal Law for Turing s Formulae
A Multivaiate Nomal Law fo Tuing s Fomulae Zhiyi Zhang Depatment of Mathematics and Statistics Univesity of Noth Caolina at Chalotte Chalotte, NC 28223 Abstact This pape establishes a sufficient condition
More informationLET a random variable x follows the two - parameter
INTERNATIONAL JOURNAL OF MATHEMATICS AND SCIENTIFIC COMPUTING ISSN: 2231-5330, VOL. 5, NO. 1, 2015 19 Shinkage Bayesian Appoach in Item - Failue Gamma Data In Pesence of Pio Point Guess Value Gyan Pakash
More informationSuppose you have a bank account that earns interest at rate r, and you have made an initial deposit of X 0
IOECONOMIC MODEL OF A FISHERY (ontinued) Dynami Maximum Eonomi Yield In ou deivation of maximum eonomi yield (MEY) we examined a system at equilibium and ou analysis made no distintion between pofits in
More informationANA BERRIZBEITIA, LUIS A. MEDINA, ALEXANDER C. MOLL, VICTOR H. MOLL, AND LAINE NOBLE
THE p-adic VALUATION OF STIRLING NUMBERS ANA BERRIZBEITIA, LUIS A. MEDINA, ALEXANDER C. MOLL, VICTOR H. MOLL, AND LAINE NOBLE Abstact. Let p > 2 be a pime. The p-adic valuation of Stiling numbes of the
More informationASTR415: Problem Set #6
ASTR45: Poblem Set #6 Cuan D. Muhlbege Univesity of Mayland (Dated: May 7, 27) Using existing implementations of the leapfog and Runge-Kutta methods fo solving coupled odinay diffeential equations, seveal
More informationTo Feel a Force Chapter 7 Static equilibrium - torque and friction
To eel a oce Chapte 7 Chapte 7: Static fiction, toque and static equilibium A. Review of foce vectos Between the eath and a small mass, gavitational foces of equal magnitude and opposite diection act on
More informationAVS fiziks. Institute for NET/JRF, GATE, IIT-JAM, JEST, TIFR and GRE in PHYSICAL SCIENCES
ELECTROMAGNETIC THEORY SOLUTIONS GATE- Q. An insulating sphee of adius a aies a hage density a os ; a. The leading ode tem fo the eleti field at a distane d, fa away fom the hage distibution, is popotional
More informationE(r,t) = e 3. r 3. (b) Show that the transverse current, J t,is 3n(n e 3 ) e 3
Polem Set 3 (Jakson 6.20).. An example of the pesevation of ausality and finite speed of popagation in spite of the use of the Coulomg gauge is affoded y a unit stength dipole soue that is flashed on and
More informationLight Time Delay and Apparent Position
Light Time Delay and ppaent Position nalytical Gaphics, Inc. www.agi.com info@agi.com 610.981.8000 800.220.4785 Contents Intoduction... 3 Computing Light Time Delay... 3 Tansmission fom to... 4 Reception
More informationPower transformations in correspondence analysis
Powe tansfomations in oespondene analysis Mihael Geenae Depatment of Eonomis and Business Univesitat Pompeu Faba Ramon Tias Fagas, 5-7 08005 Baelona SPAIN E-mail: mihael@upfes Abstat: Powe tansfomations
More information1 Explicit Explore or Exploit (E 3 ) Algorithm
2.997 Decision-Making in Lage-Scale Systems Mach 3 MIT, Sping 2004 Handout #2 Lectue Note 9 Explicit Exploe o Exploit (E 3 ) Algoithm Last lectue, we studied the Q-leaning algoithm: [ ] Q t+ (x t, a t
More informationHOW TO TEACH THE FUNDAMENTALS OF INFORMATION SCIENCE, CODING, DECODING AND NUMBER SYSTEMS?
6th INTERNATIONAL MULTIDISCIPLINARY CONFERENCE HOW TO TEACH THE FUNDAMENTALS OF INFORMATION SCIENCE, CODING, DECODING AND NUMBER SYSTEMS? Cecília Sitkuné Göömbei College of Nyíegyháza Hungay Abstact: The
More information