ANOTEONPARAMETERDIFFERENTIATIONOFMATRIXEXPONENTIALS,WITHAPPLICATIONSTOCONTINUOUS-TIMEMODELING(SHORTRUNNINGTITLE:PARAMETERDIFFERENTIATIONOFMATRIXEXPONENTIALS)HENGHSIUTSAIInstituteofStatisticalScience,AcademiaSinica,Taipei,Taiwan115,R.O.C.htsai@stat.sinica.edu.twK.S.CHANDepartmentofStatistics&ActuarialScience,UniversityofIowa,IowaCity,IA52242,U.S.A.kchan@stat.uiowa.eduDecember27,2001SummaryWehavederivedanewanalyticformulaforevaluatingthederivativesofamatrixexponential.Incontrasttosomeexistingmethods,eigenvaluesandeigenvectorsdonotappearexplicitlyintheformulae,althoughweshowthatanecessaryandsucientconditionforthevalidityoftheformulaeisthatthematrixhasdistincteigenvalues.Thenewformulaexpressesthederivativesofamatrixexponentialintermsofminors,polynomials,exponentialofthematrixaswellasmatrixinver-1sion,andhenceisalgebraicallymoremanageable.Forsparsematrices,theformulacanbefurthersimplied.Twoexamplesarediscussedinsomedetails.Forthecompanionmatrixofacontinuous-timeautoregressivemovingaverageprocess,thederivativesoftheexponentialofthecompanionmatrixcanbecomputedre-cursively.WeillustrateanuseoftheseformulaeintheconditionalleastsquareestimationofaCAR(p)modelthatleadstoanumericallystableestimationpro-cedure.Thesecondexampleconcernstheexponentialofthetridiagonaltransitionintensitymatrixofanite-state-spacecontinuous-timeMarkovchainwhosein-stantaneoustransitionsmustbebetweenadjacentstates.Keywords:Cayley-HamiltonTheorem;CARMAmodels;companionmatrix;nite-state-spacecontinuous-timeMarkovprocesses;maximumlikelihoodestimation;minimalpolynomial;tridiagonalmatrix.1IntroductionVariousmethodsofparameterdierentiationofamatrixexponentialhavebeenstudiedinstatisticalmechanicsandquantumtheory(see,e.g.,Wilcox,1967),aswellasinthemathematics,economicsandstatisticsliterature,seee.g.,JennrichandBright(1976),VanLoan(1978),KalbeischandLawless(1985),Graham(1986),HornandJohnson(1991),ChanandMunoz-Hernandez(1997),ChenandZadrozny(2001).Forcontinuous/discretestatespacemodelling(see,e.g.,Jazwin-ski,1970andSinger,1995),parameterdierentiationofamatrixexponentialisneededforcomputingtheanalyticalscorefunction.Forcontinuous-timeMarkovmodeling,ecientalgorithmforthecomputationofthetransitionprobabilityma-trixanditsderivativeswithrespecttothetransitionintensityparametersisneededformaximumlikelihoodestimation.Forexample,seeKalbeischandLawless(1985)foranapproachofanalyzingapanelofcategoricaldatabyassumingthatthedataareobtainedfromsamplingalatentcontinuous-timenite-state-spaceMarkovprocess.Weproposeinthisnoteanalternativemethodforcomputingthederivativesofamatrixexponential.Incontrasttosomeexistingmethods,eigenvaluesandeigenvectorsdonotappearexplicitlyintheformulae,althoughweshowthatanec-2essaryandsucientconditionforthevalidityoftheformulaeisthatthematrixhasdistincteigenvalues.Thenewformulaexpressesthederivativesofamatrixexponentialintermsofminors,polynomials,exponentialofthematrixaswellasmatrixinversion,andhenceisalgebraicallymoremanageable.Whenthema-trixhasrepeatedeigenvalues,itseemshardtoextendtheresults.Seetheendofsection3fordiscussion.Fortunately,inmoststatisticalapplicationsthatinvolvematrixexponentials,thedistincteigenvalueassumptionoftenholds.Forexam-ple,incontinuous-timeMarkovchainmodelling,formostmodelsofinterest,thetransitionintensitymatrixhasdistincteigenvaluesforalmostallparametervalues(see,e.g.,KalbeischandLawless,1985).Thisnoteisorganizedasfollows.Inx2,wederivethenewformulaforcomputingthederivativesofamatrixexponentialandanecessaryandsucientconditionforthevalidityoftheformula.Forsparsematrices,theformulamaybefurthersimplied.Twointerestingexamplesaretheexponentialofthecompanionmatrixarisingfromacontinuous-timeautoregressivemovingaverageprocessandthatofthetridiagonaltransitionintensitymatrixarisingfromacontinuous-timeMarkovchainwhoseinstantaneoustransitionsmustbejumpsbetweenadjacentcategories.Thesimpliedformulaeforthesetwoexamplesaregiveninx3.2MainresultsLetA=[aij]beappmatrixwhoseelementsarefunctionsof#=(#1;:::;#r)0.Byequation(2.1)ofWilcox(1967),wehavethat,fori=1;:::;r,@etA@#i=Zt0e(t u)A@A@#i!euAdu:(1)Alternatively,ifweassumeAhasdistincteigenvaluesd1;;dpandXistheppmatrixwhosejthcolumnisarighteigenvectorcorrespondingtodj,thenA=XDX 1,whereD=diag(d1;;dp).ThenetA=Xdiag(ed1t;;edpt)X 1,and@etA@#u=XVuX 1;u=1;;r;(2)whereVuisappmatrixwith(i;j)entryg(u)ij(edit edjt)=(di dj);i6=j;3g(u)iitedit;i=j;andg(u)ijisthe(i;j)entryinG(u)=X 1(@A=@u)X.SeeKalbeischandLawless(1985)fortheaboveformulaandrelateddiscussions.SeealsoJennrichandBright(1976)andChanandMunos-Hernandez(1997).WhenAhasrepeatedeigenvalues,ananalogousdecompositionofAtoJordancanonicalformispossible(seechapter4ofCoxandMiller,1965).ButaspointedoutbyKalbeischandLawless(1985),thisisrarelynecessary,sinceformostmodelsofinterestincontinuous-timeMarkovmodelling,Ahasdistincteigenvaluesforalmostallparameters.Oneofthemainresultsofthispaperistoderiveanotherclosedformsolutionfor@etA=@#i.Forr=1;:::;p,denertobeap1vectorwith1inpositionrand0elsewhere.For1i;jp,letBij=i0j,anddeneij=Zt0e(