A modified hyperbolic secant distribution

Similar documents
Revisiting what you have learned in Advanced Mathematical Analysis

Section 2: The Z-Transform

Week 06 Discussion Suppose a discrete random variable X has the following probability distribution: f ( 0 ) = 8

On the Derivatives of Bessel and Modified Bessel Functions with Respect to the Order and the Argument

The Laplace Transform

Laplace Transform. National Chiao Tung University Chun-Jen Tsai 10/19/2011

Bicomplex Version of Laplace Transform

Songklanakarin Journal of Science and Technology SJST R1 Thongchan. A Modified Hyperbolic Secant Distribution

Fourier Series and Parseval s Relation Çağatay Candan Dec. 22, 2013

3.4 Repeated Roots; Reduction of Order

The model proposed by Vasicek in 1977 is a yield-based one-factor equilibrium model given by the dynamic

A Study on the Nature of an Additive Outlier in ARMA(1,1) Models

Inverse Fourier Transform. Properties of Continuous time Fourier Transform. Review. Linearity. Reading Assignment Oppenheim Sec pp.289.

K x,y f x dx is called the integral transform of f(x). The function

Single Correct Type. cos z + k, then the value of k equals. dx = 2 dz. (a) 1 (b) 0 (c)1 (d) 2 (code-v2t3paq10) l (c) ( l ) x.

Fourier. Continuous time. Review. with period T, x t. Inverse Fourier F Transform. x t. Transform. j t

Math 266, Practice Midterm Exam 2

Systems of First Order Linear Differential Equations

A Study of the Solutions of the Lotka Volterra. Prey Predator System Using Perturbation. Technique

Systems of First Order Linear Differential Equations

Relation between Fourier Series and Transform

Available at Vol. 1, No. 1 (2006) pp

Applied Statistics and Probability for Engineers, 6 th edition October 17, 2016

Chapter 5 The Laplace Transform. x(t) input y(t) output Dynamic System

1 Finite Automata and Regular Expressions

LAPLACE TRANSFORMS AND THEIR APPLICATIONS

Engine Thrust. From momentum conservation

THE ECONOMETRIC MODELING OF A SYSTEM OF THREE RANDOM VARIABLES WITH THE β DEPENDENCE

Analytical Solution of A Differential Equation that Predicts the Weather Condition by Lorenz Equations Using Homotopy Perturbation Method

16.512, Rocket Propulsion Prof. Manuel Martinez-Sanchez Lecture 3: Ideal Nozzle Fluid Mechanics

Lecture 21 : Graphene Bandstructure

SOLUTIONS. 1. Consider two continuous random variables X and Y with joint p.d.f. f ( x, y ) = = = 15. Stepanov Dalpiaz

A Condition for Stability in an SIR Age Structured Disease Model with Decreasing Survival Rate

Elementary Differential Equations and Boundary Value Problems

CSE 245: Computer Aided Circuit Simulation and Verification

Chapter 3: Fourier Representation of Signals and LTI Systems. Chih-Wei Liu

Stability of time-varying linear system

General Article Application of differential equation in L-R and C-R circuit analysis by classical method. Abstract

4.1 The Uniform Distribution Def n: A c.r.v. X has a continuous uniform distribution on [a, b] when its pdf is = 1 a x b

ELECTRIC VELOCITY SERVO REGULATION

TOPIC 5: INTEGRATION

Boyce/DiPrima 9 th ed, Ch 7.8: Repeated Eigenvalues

An Indian Journal FULL PAPER. Trade Science Inc. A stage-structured model of a single-species with density-dependent and birth pulses ABSTRACT

Midterm exam 2, April 7, 2009 (solutions)

Lecture 11 Waves in Periodic Potentials Today: Questions you should be able to address after today s lecture:

The Procedure Abstraction Part II: Symbol Tables and Activation Records

Option markets and the stochastic behavior of commodity prices 1

More on FT. Lecture 10 4CT.5 3CT.3-5,7,8. BME 333 Biomedical Signals and Systems - J.Schesser

On Ψ-Conditional Asymptotic Stability of First Order Non-Linear Matrix Lyapunov Systems

Advanced Engineering Mathematics, K.A. Stroud, Dexter J. Booth Engineering Mathematics, H.K. Dass Higher Engineering Mathematics, Dr. B.S.

3(8 ) (8 x x ) 3x x (8 )

A Tutorial of The Context Tree Weighting Method: Basic Properties

2.1. Differential Equations and Solutions #3, 4, 17, 20, 24, 35

CIVL 8/ D Boundary Value Problems - Rectangular Elements 1/7

VIBRATION ANALYSIS OF CURVED SINGLE-WALLED CARBON NANOTUBES EMBEDDED IN AN ELASTIC MEDIUM BASED ON NONLOCAL ELASTICITY

DEPARTMENT OF ELECTRICAL &ELECTRONICS ENGINEERING SIGNALS AND SYSTEMS. Assoc. Prof. Dr. Burak Kelleci. Spring 2018

Lecture 1: Numerical Integration The Trapezoidal and Simpson s Rule

UNSTEADY HEAT TRANSFER

2 T. or T. DSP First, 2/e. This Lecture: Lecture 7C Fourier Series Examples: Appendix C, Section C-2 Various Fourier Series

FREE VIBRATION AND BENDING ANALYSES OF CANTILEVER MICROTUBULES BASED ON NONLOCAL CONTINUUM MODEL

An Optimal Ordering Policy for Inventory Model with. Non-Instantaneous Deteriorating Items and. Stock-Dependent Demand

CS 688 Pattern Recognition. Linear Models for Classification

Inventory Management Model with Quadratic Demand, Variable Holding Cost with Salvage value

CHAPTER CHAPTER14. Expectations: The Basic Tools. Prepared by: Fernando Quijano and Yvonn Quijano

LGOVNATDEFUSAAD

Lecture 4: Laplace Transforms

Wave Equation (2 Week)

I) Title: Rational Expectations and Adaptive Learning. II) Contents: Introduction to Adaptive Learning

Digital Signal Processing. Digital Signal Processing READING ASSIGNMENTS. License Info for SPFirst Slides. Fourier Transform LECTURE OBJECTIVES

Contraction Mapping Principle Approach to Differential Equations

Lagrangian for RLC circuits using analogy with the classical mechanics concepts

3. Renewal Limit Theorems

EEE 303: Signals and Linear Systems

CPSC 211 Data Structures & Implementations (c) Texas A&M University [ 259] B-Trees

On the Speed of Heat Wave. Mihály Makai

Control Systems. Modelling Physical Systems. Assoc.Prof. Haluk Görgün. Gears DC Motors. Lecture #5. Control Systems. 10 March 2013

Reliability Analysis of a Bridge and Parallel Series Networks with Critical and Non- Critical Human Errors: A Block Diagram Approach.

A MATHEMATICAL MODEL FOR NATURAL COOLING OF A CUP OF TEA

Boyce/DiPrima 9 th ed, Ch 2.1: Linear Equations; Method of Integrating Factors

A Production Inventory Model for Different Classes of Demands with Constant Production Rate Considering the Product s Shelf-Life Finite

7.4 QUANTUM MECHANICAL TREATMENT OF FLUCTUATIONS *

Double Slits in Space and Time

Discussion 06 Solutions

ON NEW INEQUALITIES OF SIMPSON S TYPE FOR FUNCTIONS WHOSE SECOND DERIVATIVES ABSOLUTE VALUES ARE CONVEX

LaPlace Transform in Circuit Analysis

5.1-The Initial-Value Problems For Ordinary Differential Equations

A Simple Method for Determining the Manoeuvring Indices K and T from Zigzag Trial Data

Canadian Journal of Physics. Kantowski-Sachs modified holographic Ricci dark energy model in Saez-Ballester theory of gravitation

Logistic equation of Human population growth (generalization to the case of reactive environment).

UNIT #5 EXPONENTIAL AND LOGARITHMIC FUNCTIONS

Ratio-Product Type Exponential Estimator For Estimating Finite Population Mean Using Information On Auxiliary Attribute

Instructions for Section 1

e t dt e t dt = lim e t dt T (1 e T ) = 1

PHA Final Exam Fall On my honor, I have neither given nor received unauthorized aid in doing this assignment.

Spring 2006 Process Dynamics, Operations, and Control Lesson 2: Mathematics Review

Errata for Second Edition, First Printing

Integration Continued. Integration by Parts Solving Definite Integrals: Area Under a Curve Improper Integrals

PHA First Exam Fall On my honor, I have neither given nor received unauthorized aid in doing this assignment.

Charging of capacitor through inductor and resistor

Chahrazed L Journal of Scientific and Engineering Research, 2018, 5(4): and

Transcription:

Songklnkrin J Sci Tchnol 39 (1 11-18 Jn - Fb 2017 hp://wwwsjspsuch Originl Aricl A modifid hyprbolic scn disribuion Pnu Thongchn nd Wini Bodhisuwn * Dprmn of Sisics Fculy of Scinc Kssr Univrsiy Chuchk Bngkok 10900 Thilnd Rcivd: 17 Dcmbr 2015; Rvisd: 4 Mrch 2016; Accpd: 24 Mrch 2016 Absrc Th ims of his ppr r o inroduc nd vlid nw disribuion which is rld o h hyprbolic funcion Th diffrnil quion is pplid o obin h survivl nd probbilisic funcions Th momn gnring funcion is providd by using h ingrl rnsform Th proposd disribuion is pplid o rl d ss I is dmonsrd h i cn b usd s n lrniv modl in vrious disciplins such s lcronics finncil whr nd rrivl ims Kywords: hyprbolic funcion hyprbolic scn disribuion scn funcion 1 Inroducion Th hyprbolic funcion is n imporn mhmicl funcion in rlion o rigonomric funcions xponnil funcions nd complx numbrs Th hyprbolic scn is pr of s of hyprbolic funcions which is dfind s 2 sch( x In 1934 Bn inroducd probbiliy x x ( disribuion rld o h hyprbolic scn funcion which is clld h hyprbolic scn disribuion (HSD (Bn 1934 This work ws hn xpndd by Tlcko (1956 who proposd h disribuion for finncil rurn modls Consqunly h HSD provids n opiml fi nd xhibis mor lpokurosis hn boh h norml nd logisic disribuions Howvr i is limid in is uiliy; i cnno k on vrious shps du o lck of flxibiliy in is prmrs hus nurl bhviors cnno b sufficinly xplind wih his modl Mny brnchs of nurl scincs mphsiz h sudy of phnomn such s h sprding of diss or growh of populion by looking rs of chng Diffrnil quion is n imporn chniqu usd o solv hs yps of problm Diffrnil quion is usful chniqu * Corrsponding uhor Emil ddrss: fsciwnb@kuch in xplining sisicl propris nd vil ool in proving or disproving sisic-bsd issus This modl cn b pplid o survivl nlysis lifim d nlysis nd rlibiliy nlysis This ppr inroducs chniqus usd o gnr nw disribuion h is consrucd by modifying h hyprbolic funcions ino mor flxibl modl using diffrnil quions W implmn h chniqu wih s ( S( b s im dfind ( S( S( ( S( S( 0 n 2 0 1 2 n n whr n blongs o {123 T} T / n nd S( S( S( k k k Th 1 procss cn b xprssd S( S( S( S( S 0 1 0 1 S( S( S S( S( S 2 1 2 k k 1 k Undr h ssumpion h is vry smll w obin S( S( S ds( k k 1 k Som propris of h survivl funcion S( r righ coninuous in Nvrhlss for coninuous survivl im T S( is coninuous wih non-incrsing funcion ( S( 1 0 nd lim S( 0 n n Th rs of h ppr is orgnizd s follows In Scion 2 nw fmily of disribuions is proposd from h hyprbolic funcions using diffrnil quions o dvlop his nw probbiliy funcion Also w will pply rl d ss o h proposd modl in Scion 3 bfor sring our conclusions in Scion 4

12 P Thongchn & W Bodhisuwn / Songklnkrin J Sci Tchnol 39 (1 11-18 2017 2 A Nw Disribuion This scion prsns n lrniv wy in which h nw probbiliy funcion could b dvlopd W pply h diffrnil quion o driv probbiliy dnsiy funcion by sing up h firs ordr diffrnil quion nd solv i o g survivl funcion W hn k h driviv of h survivl funcion Consqunly probbiliy dnsiy funcion ssocid wih survivl funcion will b obind Dfiniion 1: L T b rndom vribl on probbiliy spc ( F P wih probbiliy dnsiy funcion ( f ( ; disribuion funcion ( F ( ; nd survivl funcion ( S( ; 1 F ( ; nd ( f ( ; dfind by ( F ( ; or ( S ( ; Proposiion 1: L T b rndom vribl on probbiliy spc ( F P wih probbiliy dnsiy funcion f ( ; k nd T [0 which producs h quion blow: f ( k k k 2 2 (2 log[2 ] log[ ] 2 ( 2 whr is n iniil vlu nd 2 0 k 2 (2 1 log(2 log(1 Proof: W propos h nw disribuion by finding h soluion of h following diffrnil quion whn ( ( ki S S S ( 1 ( ( ( 1 [0 I 0 ohrwis Tking h firs driviv w will obin h probbiliy funcions s sn in Eq (1 As for h soluion funcion of h diffrnil quion w s U ( s n rbirry funcion hn muliplying Eq (2 by U ( w obin U ( ki U S ( ( ( U ( S( ( ( by driviv produc rul i is hn s o h following form U ( ki ( U ( S( ( ( ( whr U ( d Thrfor h rsul of diffrnil quion Eq (2 is (1 (2 (3 1 U ( ki S ( [ d C] ( ( U ( ( which is h survivl funcion In his cs U ( nd h soluion of h diffrnil quion wih iniil vlu problm S( 1 is ( 2 2 2 (2 k log(2 k log(2 2 S( 2 Now king h firs driviv h pdf of T is ( 2 2 2 ( (2 k log(2 k log(2 2 f ( 2 Corollry 1: k 2 2 ( L T b rndom vribl wih h pdf f ( ; b which could b xprssd s b 2 f ( ; b 2b (1 wih h locion prmr b nd iniil vlu c Proof: Similr o Eq (1 f ( 2 b 2b (1 solving h diffrnil quion xprssd s 2 ( S ( c 1 ( b ( b ( ' S b (4 obind by by ingring boh sids w obin h survivl funcion which is h rsul from solving Eq(5 cb b 2rcn( 2rcn( S( hn by king firs driviv h pdf is obind in h form b 2 f ( ; b 2b (1 whr T ( I is lso possibl o prsn h gnrl diffrnil quion s AkI 2(1 A ( AS ( ( ( ( b ( b ( ( ' S by imposing A 0 or 1 Th Eq(6 cn b rducd o Eq (2 nd Eq(5 rspcivly Following his logic S ( 1 nd 1 S ( c 1 r hrfor h iniil vlus of Eq (2 nd Eq (5 2 rspcivly (5 (6

P Thongchn & W Bodhisuwn / Songklnkrin J Sci Tchnol 39 (1 11-18 2017 13 Morovr h pdf drivd from h survivl funcion is probbiliy funcion which is ssifid h following propris; 1 Sing S( 1 nd S( c 1 r iniil vlu funcions 2 Th survivl funcion is monoonic dcrsing 3 lim S( 0 4 If h pdf corrspond o prmr spc hn i will b grr hn zro 21 Som propris of h nw disribuion Thr r mny mhods o solv h diffrnil quions An imporn mhod uss h Lplc rnsform which is rld o h momn gnring funcions Th Lplc rnsform of Eq (1 is givn by L( s s f ( d ; s 0 2 2 2 s k (2 k log[2 ] k log[ ] ( d( 2 ( 2 ( s 1 ( s 1 ( s1 2 k 2 k log[2 ] [ ( ( ( 2 (1 2 2 ( s 1 ( s 1 2 ( s 1 k (log[2 ] k ( d ( ( 2 (1 ( s 1 2 ( s 1 ( s1 2 k log[ ] ( ] d 2 ( s1 2 k log[ ] ( d 2 ( s 1 s 2 s ( s 1 2 k (log[2 ] k k log[ ] ( d ( ( ( d (7 2 2 (1 ( s 1 2 ( s 1 Focusing on h ls rm nd using h by pr ingrion chniqu w hv ( s 1 2 2 2 ( s1 ( s 1 k log[ ] k log[ ] k (8 ( d [ ( d] 2 2 2 2 ( s 1 ( s 1 (1 subsiu h ls rm of Eq (7 wih Eq (8 bcoms ( s 1 s 2 s k (log[2 ] k L( s ( d ( ( 2 2 (1 ( s 1 2 ( s 1 2 2 ( s 1 ( s 1 k log[ ] k [ ] ( d 2 2 ( s 1 ( s 1 (1 ( s1 s 2 s s 2 1 k (log[2 ] k k log[2 ] ( 1 ( d ( ( ( 2 ( s 1 (1 ( s 1 2 ( s 1 (( s 12 Th convrgnc of h ingrl pr of Eq (9 will b shown in Appndix Undr h condiion of 0 wih powr sris w rrrng Eq (9 o obin h following Lplc rnsformion k 1 L S d s ( s1 n 2( n ( ( 1 (( 1 ( ( s 1 n 0 ( s 1 s In ddiion h momn gnring funcion (mgf of is givn by M ( E( L( s for ll momns if h Lplc rnsform xiss (s shown in Appndix nd whr h rh momn is E( r ( 1 r L r (0 In priculr Th momn gnring funcion of Eq (4 is givn by (9

14 P Thongchn & W Bodhisuwn / Songklnkrin J Sci Tchnol 39 (1 11-18 2017 Figur 1 Som plos of h survivl funcion of h nw disribuion whr k = 0 1 2 6 from h boom lin rspcivly s M ( s f ( d s 2 d ( b ( b ( ( s 1 b ( s 1 b b 2 2 d d ; s 0 2 ( b b 2( b (1 ( 1 (1 ( 1 n sb 4 ( 1 (10 ( s 1 2n n 0 In h cs of Eq (4 whr symmric propris dfin h s of domin s ( hn M ( s cn b xprssd s Eq(10 Som plos of h survivl nd rlibiliy funcions rld o h proposd disribuion r prsnd in Figur 1 Th survivl funcion of Eq (4 r xhibid in Figur 2 In survivl nlysis h hzrd funcion is dfind s f ( h( whr f ( nd S( r h pdf nd survivl S( Figur 2 Som plos of h survivl funcion for h nw disribuion wih som prmr vlus of Eq (4 funcion rspcivly Th hzrd funcion of h proposd disribuion (Eq (1 cn b wrin s 2 k 2 2 2 h( 1 ( (2 k log[2 ] k log[ ] Som shps of h hzrd funcion of h proposd disribuion (Eq (1 wih som prmrs r shown in Figur 3 Figur 3 Som plos of hzrd funcions for h proposd disribuion wih prmr = 0 of Eq(1

P Thongchn & W Bodhisuwn / Songklnkrin J Sci Tchnol 39 (1 11-18 2017 15 Figur 4 Som plos of h hzrd funcions for h proposd disribuion wih som prmr vlus of Eq(4 Th hzrd funcion of Eq (4 cn b wrin s b 2 h( 2( b ( cb ( b (1 ( 2 rcn( 2rcn( whr svrl funcion wih som prmrs r illusrd in Figur 4 Figur 5 shows som shps of h pdf bsd on som slcd prmr vlus I illusrs h h nw proposd disribuion consiss of vrious shps Th symmric bhviors of Eq (4 r illusrd in Figur 6 3 Applicions Four d ss r fid wih h proposd disribuion Th firs xmpl dls wih h r of chng of filur ims of lcronic dvics rpord by Domm (2014 Th scond xmpl is h r of chng of soybn prics Chicgo Bord of Trd (CBOT from Jun 09 - Jun 14 fid wih h proposd disribuion Th hird nd fourh d ss r US July prcipiion (Top 8 soybns producion s from 1988-2015 nd h inr-rrivl ims (10 minus d s for crs (Lw 2015 rspcivly Figur 5 Som pdf plos of nw disribuion wih som prmr vlus of Eq (1

16 P Thongchn & W Bodhisuwn / Songklnkrin J Sci Tchnol 39 (1 11-18 2017 Figur 6 Som pdf plos of h nw disribuion wih som prmr vlus of Eq (4 Th simd prmrs could b crrid ou by mximum liklihood simion (MLE In his sudy w us bbml (Bolkr & R Tm 2014 pckg of R progrmming lngug (R Cor Tm 2014 o obin h prmr sims Th fiing disribuions for h four xmpls r vrifid s shown in Figur 7 In ddiion h simd prmrs using MLE r shown in Tbl 1 4 Conclusions Th im of his ppr is o propos n lrniv mhod o gnr nw survivl funcion W hv don so using h rlionship bwn diffrnil quion nd hyprbolic funcion nd obind h soluion of h mhod s nonincrsing funcion; w obin h survivl nd h probbiliy dnsiy funcions from hs chniqus Morovr Figur 7 Som fid disribuions wih min(

P Thongchn & W Bodhisuwn / Songklnkrin J Sci Tchnol 39 (1 11-18 2017 17 Tbl 1 Th simd prmrs using MLE D s Esimd prmrs 1 Th r of chng of filur ims of lcronic dvics 0 178 2 Th r of chng of soybn prics Chicgo Bord of Trd (CBOT from Jun 09 - Jun 14 0 3274479 10 14 3 July prcipiion from 1988-2015 0 316 4 Th inr-rrivl ims (10 minus d s for crs 0 2 k w orgniz h gnrl form of h modl Eq(6 which producs wo survivl nd h probbilisic funcions Th vrious grphicl syls of h rsul cn b usd o dmonsr is flxibiliy in nlyzing bhvior in rl d Acknowldgmns Th uhors would lik o hnk o Dprmn of Sisics Fculy of Scinc nd h Grdu School of Kssr Univrsiy Also w hnk o Chron Pokphnd Foods (CPF for supporing h firs uhor Rfrncs Bn W D (1934 Th probbiliy lw for h sum of n indpndn vribls Bullin of Amricn Mhmicl Sociy 40 284-290 Bolkr B & Tm R D C (2014 Tools for Gnrl Mximum Liklihood Esimion (R pckg vrsion 1017 Domm F & Condino F (2014 A nw clss of disribuion funcions for lifim d Rlibiliy Enginring nd Sysm Sfy 129 36-45 Fishr M J (2014 Gnrlizd hyprbolic scn disribuion Brlin Grmny: Springr Jffry A & Zwillingr D (2007 Tbl of ingrl sris nd produc Nw York NY: Acdmic Prss Jing R (2013 A nw bhub curv modl wih fini suppor Rlibiliy Enginring nd Sysm Sfy 119 44-51 Klbflisch J D & Prnic R L (1980 Th sisicl nlysis of filur d Nw York NY: John Wily & Son Lw A (2015 Simulion modling nd nlysis (McGrw- Hill sris in indusril nginring nd mngmn Nw York NY: McGrw-Hill Lwlss J F (1982 Sisicl modls nd mhods for lifim d Nw York NY: John Wily & Son Mnoukin E B & Ndu P (1988 A no on h hyprbolic scn disribuion Th Amricn Sisicin 42 77-79 R Cor Tm (2014 R: A lngug nd nvironmn for sisicl compuing Vinn Ausri: R Foundion for Sisicl Compuing Tlcko J (1956 Prk disribuions nd hir rol in h hory of winr s sochsic vribls Trbjos d Esisic 17 159-174 Wng F K (2000 A nw modl wih bhub-shpd filur r using n ddiiv burr xii disribuion Rlibiliy Enginring nd Sysm Sfy 70 305-312

18 P Thongchn & W Bodhisuwn / Songklnkrin J Sci Tchnol 39 (1 11-18 2017 APPENDIX A convrgnc will b provd ( s 1 s 2 s 1 k (log[2 ] k d 2 2 ( s 1 (1 ( s 1 2 ( s 1 ( 1 ( ( ( Th ingrl is impropr ingrl w us h comprison s o prov convrgnc by h following horm A comprison s horm Suppos h f nd g r coninuous funcions wih 0 f ( g( for hn 1 If g( d is convrgn hn f ( d is convrgn 2 If g( d is divrgn hn f ( d is divrgn Espcilly his impropr ingrl s shown nd ( s1 k ( d 2 (1 W s ( s 1 g( d k ( s1 g( d k = k ( s 1 ( s 1 k f ( d ( d 2 2 (1 whrs ( ( s 1 [ ] ; s 0 k ( s1 I is vrifid h g( d convrg o [ ] hrrfor ( s 1 g is grr hn f ( for ll s 0 ( s 1 k f ( d ( d 2 (1 is convrgn oo