 respectively.

### Transit parameters

Before analyzing the TTV signals, we revise the transit parameters of KOI-94c, KOI-94d, and KOI-94e obtained by the Kepler team (Table \[ppar\_prior\]) so that they are consistent with the light curves obtained in Section \[sec:dataprocess\]. Here we first use the parameters publicized by the Kepler team to phase-fold the observed transit light curves, and then refit those phase curves to obtain the revised transit parameters.

In the first step, we fit each of the detrended light curve centered at the transit (for $\sim 1.7$ times its duration) to obtain the times of transit centers $t_c$, using a Markov chain Monte Carlo (MCMC) algorithm. Here we use a light curve model by @2009ApJ...690....1O, and fix $a/R_{\star}$, $R_p/R_{\star}$, and $b$ to the values obtained by the Kepler team, assuming $e = 0$. We model the limb-darkening using a quadratic law in Eq.(\[qlimb\]) and adopt the limb-darkening coefficients $u_1$ and $u_2$ obtained by @2012ApJ...759L..36H (all these parameters are summarized in Table \[ppar\_prior\]). Since the detrend procedure in Section \[sec:dataprocess\] can remove only the out-of-transit outliers, we also exclude in-transit $5\sigma$ outliers of this fit, if any, and fit the light curve again. Using the series of $t_c$ obtained in this way, we construct the phase-folded transit light curve for each planet.

As the second step, we fit the resulting phase-folded transit light curve for $a/R_{\star}$, $R_p/R_{\star}$, $b$, $u_1$, and $u_2$ using the same light curve model as above. In this way, we obtain the revised values of the set of parameters shown in Table \[ppar\_fit\] and the corresponding best-fit light curves (Figures \[94c\_phase\] to \[94e\_phase\]).[^1] In this fit, all the parameters converge well in the case of KOI-94d. In contrast, $a/R_{\star}$ and $b$ of KOI-94c and KO
respectively. 

 # # # Transit parameters 

 Before analyzing the TTV signals, we retool the passage parameters of KOI-94c, KOI-94d, and KOI-94e obtained by the Kepler team (mesa \[ppar\_prior\ ]) so that they are coherent with the light curves obtain in Section \[sec: dataprocess\ ]. Here we first use the parameter publicized by the Kepler team to phase - fold the note transit light curves, and then refit those phase curves to obtain the revised theodolite parameters. 

 In the first step, we equip each of the detrended light curve centered at the transit (for $ \sim 1.7 $ multiplication its duration) to obtain the times of transit centers $ t_c$, using a Markov range Monte Carlo (MCMC) algorithm. Here we use a light curve model by @2009ApJ... 690.... 1O, and fix $ a / R_{\star}$, $ R_p / R_{\star}$, and $ b$ to the values prevail by the Kepler team, assuming $ e = 0$. We model the limb - darkening using a quadratic law in Eq.(\[qlimb\ ]) and adopt the limb - darkening coefficients $ u_1 $ and $ u_2 $ obtain by @2012ApJ... 759L.. 36H (all these parameters are summarized in Table \[ppar\_prior\ ]). Since the detrend routine in Section \[sec: dataprocess\ ] can absent only the out - of - transit outlier, we also exclude in - transit $ 5\sigma$ outliers of this paroxysm, if any, and fit the light curve again. Using the series of $ t_c$ obtained in this means, we construct the phase - folded transit light curvature for each planet. 

 As the second tone, we fit the resulting phase - folded transit light curvature for $ a / R_{\star}$, $ R_p / R_{\star}$, $ b$, $ u_1 $, and $ u_2 $ using the same light curvature model as above. In this way, we obtain the revised value of the set of parameter shown in Table \[ppar\_fit\ ] and the corresponding well - fit light curve (Figures \[94c\_phase\ ] to \[94e\_phase\]).[^1 ] In this paroxysm, all the parameters converge well in the casing of KOI-94d. In contrast, $ a / R_{\star}$ and $ b$ of KOI-94c and KO
 redpectively.

### Transit paramtters

Before analyeibg the TTV sjgnals, wd revise the transit parametxrs if KOU-94c, KOI-94d, and KOI-94e obtaived by thv Kepler ream (Rable \[ppar\_'dior\]) so that fmey axe consistent wijh the light curves obtainad iu Section \[sec:dataprocess\]. Here we firft use yhf parameters poblicpzqd bg the Kepler team to phase-fold the observtd transit light cirves, and then refit those phade curves to obtaij the revisgs twqnsit paramegers.

In the first step, se fit each of the detrended lieht cbrve centergb at jhe transit (hor $\sii 1.7$ times its duratiot) to obyain the times of trqnsit centers $t_c$, usinj a Markov chain Monje Carlo (MWME) algorithm. Here we uwe a licht wurvd moael bb @2009AlJ...690....1O, anf fmx $a/R_{\star}$, $R_l/R_{\star}$, and $v$ to the values obtsigvc by the Kepmer tewm, assuming $e = 0$. We model the limb-darkeninc uaing a quadratic law in Eq.(\[qlimb\]) and adopt thg limb-darkqning coefficients $u_1$ and $u_2$ obtained by @2012ApJ...759L..36H (all dhese oarcneterr age summarized in Table \[ppar\_prior\]). Since the deedemd procedure in Fection \[sec:cahaltocess\] can remuve only fhe out-of-transit oktliers, we aoso excluqe im-transit $5\sigma$ outliers of rhis fit, if cny, and fit the light curve agaiu. Usinb the series of $t_c$ obtained nn thia way, we cojstruct tgd phase-folded trxnspt lhght curve for each planet.

Ws the serond xtep, we fit the rqsulting pjase-folded transit light curvg for $d/R_{\star}$, $R_p/R_{\dtar}$, $b$, $u_1$, and $u_2$ using the same lmjht curve modgl ds dbove. In this way, we obtaig the revised ralues oy the ret of parzmeters shown in Twble \[ppar\_fit\] djd the correvponding besr-fit light zjrves (Figures \[94v\_phase\] to \[94e\_phase\]).[^1] In this fit, all the iaramgtsrs converge weol un the case of LOI-94a. Ig bonvrast, $d/R_{\star}$ and $b$ of YOI-94z and YO
respectively. ### Transit parameters Before analyzing the we the transit of KOI-94c, KOI-94d, Kepler (Table \[ppar\_prior\]) so they are consistent the light curves obtained in Section Here we first use the parameters publicized by the Kepler team to phase-fold observed transit light curves, and then refit those phase curves to obtain the transit In first we fit each of the detrended light curve centered at the transit (for $\sim 1.7$ times duration) to obtain the times of transit centers using a Markov chain Carlo (MCMC) algorithm. Here we a curve model @2009ApJ...690....1O, fix $R_p/R_{\star}$, and $b$ the values obtained by the Kepler team, assuming $e = 0$. We model the limb-darkening using a law in adopt the coefficients and obtained by @2012ApJ...759L..36H parameters are summarized in Table \[ppar\_prior\]). procedure in Section \[sec:dataprocess\] can remove only the outliers, we exclude in-transit $5\sigma$ outliers of this if any, and fit the light curve again. the series of $t_c$ obtained in this way, we construct the phase-folded transit light curve planet. As the second we fit the phase-folded light for $R_p/R_{\star}$, $b$, and $u_2$ using the same light curve model as above. In way, we obtain the revised values of the set of in \[ppar\_fit\] and the best-fit light curves (Figures to In this fit, all converge in KOI-94d. contrast, and $b$ of KOI-94c KO
 respectively.

### Transit parameTers

Before AnalyZinG thE TtV siGnalS, we revise the trANsit Parameters of KOI-94c, KOI-94d, anD KOI-94e ObTAineD By The KePler teaM (taBLE \[ppAr\_PrIor\]) So THaT they Are ConsistEnt with the LigHt Curves obtainED iN Section \[seC:daTaprocess\]. HerE we First uSe The PArameTerS publIcized BY the KePler team tO pHAse-folD The obseRVEd TranSit light curves, and THeN Refit those phasE curveS tO ObTAIn tHe rEvised tranSiT paraMEters.

In THe FIRSt sTEp, we fit each of The detrendeD LigHt curvE cEntERed at tHe traNsIT (foR $\sim 1.7$ times itS durAtion) to obTain thE Times of TRansit cEnters $T_c$, uSinG a MaRKoV cHaiN MONte cArLo (McmC) aLgorithm. heRe We use A ligHT CURve mOdeL by @2009APJ...690....1O, anD fix $a/R_{\star}$, $R_p/R_{\StaR}$, and $B$ To tHe valUes obTainEd By the kepler Team, aSsUming $e = 0$. We model thE limB-darkeninG usInG a qUaDratiC Law in EQ.(\[qlImb\]) And adopT the limB-DarKeNING cOefficients $u_1$ and $u_2$ obTaINEd By @2012ApJ...759L..36H (aLl thesE PaRaMEters are SuMmaRizeD IN TablE \[ppaR\_PrIor\]). Since The detREnD pRocedurE iN SectiOn \[Sec:DatAprocESs\] caN removE only the Out-of-TRansit outliers, WE also exclude iN-TrANSiT $5\SigmA$ ouTliers of thiS fit, IF any, And fIT tHe lIGht cuRve agAiN. usINg the series of $t_c$ obtaInEd in thIs way, We construct thE phase-foldED TRansit liGht cURvE For each planet.

AS the sEcond step, wE Fit the reSultiNg phase-fOlded tranSIT light cuRve For $A/R_{\sTar}$, $r_P/r_{\sTar}$, $b$, $u_1$, and $u_2$ usinG THe saMe Light cuRve Model as AboVe. IN thIs wAy, We obtain tHe reviseD vAlUeS oF thE set oF ParameteRs ShoWn In TAble \[pPAr\_fit\] aNd the CorrEsPoNDinG best-fiT LiGHT curVeS (FIgurEs \[94c\_PhAse\] to \[94E\_phaSE\]).[^1] In This fit, All the parAmeTErs cOnVeRge well In the case of KOi-94d. in contrast, $A/R_{\StaR}$ and $b$ oF koI-94c and KO
 respectively.

### Transi t paramete rs

B efo rean alyz ingthe TTV signal s , we  revise the transit pa ramet er s  ofK OI -94c,  KOI-94 d ,a n d K OI -9 4eob t ai ned b y t he Kepl er team (T abl e\[ppar\_prio r \] ) so thatthe y are consis ten t with  t hel ightcur ves o btaine d  in Se ction \[s ec : datapr o cess\]. H er e we  first use the pa r am e ters publicize d by t he Ke p l ertea m to phase -f old t h e obser v ed t r ans i t light curve s, and then ref it tho se  ph a se cur ves t oo bta in the revi sedtransit p aramet e rs.

In the fir st ste p,wefite ac hofth e  de t re nde d  li ght curv ece ntere d at t h e  tra nsi t (f or $\ sim 1.7$ time s i ts d u rat ion)to ob tain  t he ti mes of  tran si t centers $t_c$ , us ing a Mar kov  c hai nMonte Carlo(MC MC)  algori thm. He r e w eu s e  a  light curve model  b y @2 009ApJ.. .690.. . .1 O, and fix$a /R_ {\st a r }$, $ R_p/ R _{ \star}$,  and $ b $to  the va lu es obt ai ned  by  theK eple r team , assumi ng $e = 0$. We model the limb-dark e ni n g  u s inga q uadratic la w in Eq.( \[ql i mb \]) and a doptth e  l i mb-darkening coeffi ci ents $ u_1$and $u_2$ obt ained by @ 2 0 1 2ApJ...7 59L. . 36 H  (all these pa ramet ers are su m marizedin Ta ble \[pp ar\_prior \ ] ). Since  th e d etr end p ro cedure in Sec t i on \ [s ec:data pro cess\]can  re mov e o nl y the out -of-tran si tou tl ier s, we also exc lu dein -tr ansit $5\sig ma$ o utli er so f t his fit ,  i f any,  a nd  fit  th elight  cur v e a gain. U sing theser i es o f$t _c$ obt ained in this  w ay, we con st ruc t thep h ase-fold ed transit light curvef or each  pl anet. 

As  the seco ndstep,wef it the  resul tingph ase - f olded t ra nsi tlight curv e for  $a/R _{ \sta r}$, $R _p/R_{\star}$, $b$ ,  $u _1$, and $u_2 $ u sing t he  sa m el igh tc urv e model as above.  In this w ay ,  w e obtain t h e r ev ised va lues of  thes et of p arameters  shown in  T able \ [pp ar\_fit\]and thecorrespon d ing b e st -fitlig ht cur ve s ( Figur es \[9 4 c\_ phase \] to\[ 94e\_p hase\ ]) .[^1] In  this fit, all the para meters  conv erg e well in  th e  ca se of KOI -94d . In contr ast , $ a/R_{ \st a r}$ a nd $ b $ofK OI-94 c an d  KO
 respectively.

###_Transit parameters

Before_analyzing the TTV signals,_we revise_the_transit parameters_of_KOI-94c, KOI-94d, and_KOI-94e obtained by_the Kepler team (Table_\[ppar\_prior\]) so that_they_are consistent with the light curves obtained in Section \[sec:dataprocess\]. Here we first use_the_parameters publicized_by_the_Kepler team to phase-fold the_observed transit light curves, and_then refit_those phase curves to obtain the revised transit_parameters.

In_the first step,_we fit each of the detrended light curve centered_at the transit (for $\sim 1.7$_times its duration)_to_obtain_the times of transit_centers $t_c$, using a Markov chain_Monte Carlo (MCMC) algorithm. Here we_use a light curve model by @2009ApJ...690....1O,_and fix $a/R_{\star}$, $R_p/R_{\star}$, and $b$_to the values obtained by_the Kepler_team, assuming $e = 0$._We model the_limb-darkening using_a quadratic law_in Eq.(\[qlimb\]) and adopt the limb-darkening_coefficients $u_1$ and_$u_2$ obtained by @2012ApJ...759L..36H (all these_parameters_are summarized in_Table_\[ppar\_prior\])._Since the_detrend procedure in_Section_\[sec:dataprocess\] can_remove_only the out-of-transit outliers, we also_exclude_in-transit $5\sigma$ outliers of this fit, if_any, and fit the_light_curve again. Using the_series of $t_c$ obtained in_this way, we construct the phase-folded_transit light_curve for_each planet.

As the second step, we fit the resulting phase-folded transit_light curve for $a/R_{\star}$, $R_p/R_{\star}$, $b$,_$u_1$, and $u_2$ using_the same_light_curve model as_above._In this_way, we obtain the revised values of_the set_of parameters shown in Table \[ppar\_fit\]_and the corresponding best-fit_light_curves (Figures \[94c\_phase\] to \[94e\_phase\]).[^1] In_this fit, all the parameters converge_well in the case of_KOI-94d._In_contrast, $a/R_{\star}$ and $b$ of_KOI-94c and KO