 that detailed in Algorithm \[algo\_vi\], where steps 7, 8, 9 and 10 are modified to include the time dependent parameters. It follows that for the user-specific parameters the update equations take the form: $$\begin{gathered}
\lambda_{ir}^{(\alpha)} = a^{(\alpha)}+\sum_{j=1}^{{\vert{V}\vert}}\sum_{t=1}^T \frac{A_{ijt}\theta_{ijt}\chi_{ijtr}}{1-e^{-\theta_{ijt}}},\  
\mu_{ir}^{(\alpha)} = \frac{\nu_i^{(\alpha)}}{\xi_i^{(\alpha)}}+\sum_{t=1}^T\frac{\lambda_{i{t^\prime}r}^{(\gamma)}}{\mu_{i{t^\prime}r}^{(\gamma)}}\sum_{j=1}^{{\vert{V}\vert}}\frac{\lambda_{jr}^{(\beta)}}{\mu_{jr}^{(\beta)}}\frac{\lambda_{j{t^\prime}r}^{(\delta)}}{\mu_{j{t^\prime}r}^{(\delta)}}, \notag \\
\lambda_{ipr}^{(\gamma)}=a^{(\gamma)}+\sum_{j=1}^{{\vert{V}\vert}}\sum_{t:{t^\prime}=p} \frac{A_{ijt}\theta_{ijt}\chi_{ijtr}}{1-e^{-\theta_{ijt}}},\ 
\mu_{ipr}^{(\gamma)} =  \frac{\nu_p^{(\gamma)}}{\xi_p^{(\gamma)}}+\frac{\lambda_{ir}^{(\alpha)}}{\mu_{ir}^{(\alpha)}}\sum_{j=1}^{{\vert{V}\vert}} \frac{\lambda_{jr}^{(\beta)}}{\mu_{jr}^{(\beta)}}\sum_{t:{t^\prime}=p} \frac{\lambda_{j{t^\prime}r}^{(\delta)}}{\mu_{j{t^\prime}r}^{(\delta)}},\end{gathered}$$ and similar results can be obtained for the host-specific parameters $\lambda_{ir}^{(\beta)}$, $\mu_{jr}^{(\beta)}$, $\lambda_{jpr}^{(\delta)}$ and $\mu_{jpr}^{(\delta)}$. The updates for $\nu_i^{(\alpha)},\  \xi_i^{(\alpha)},\ \nu_j^{(\beta)}$ and $\xi_j^{(\beta)}$ are identical to steps 7 and 8 in Algorithm \[algo\_vi\]. For the covariates $$\begin{gathered}
\lambda_{kh
that detailed in Algorithm   \[algo\_vi\ ], where steps 7, 8, 9 and 10 are modified to include the meter pendent parameters. It follows that for the user - specific parameter the update equations take the shape: $ $ \begin{gathered } 
 \lambda_{ir}^{(\alpha) } = a^{(\alpha)}+\sum_{j=1}^{{\vert{V}\vert}}\sum_{t=1}^T \frac{A_{ijt}\theta_{ijt}\chi_{ijtr}}{1 - e^{-\theta_{ijt}}},\  
 \mu_{ir}^{(\alpha) } = \frac{\nu_i^{(\alpha)}}{\xi_i^{(\alpha)}}+\sum_{t=1}^T\frac{\lambda_{i{t^\prime}r}^{(\gamma)}}{\mu_{i{t^\prime}r}^{(\gamma)}}\sum_{j=1}^{{\vert{V}\vert}}\frac{\lambda_{jr}^{(\beta)}}{\mu_{jr}^{(\beta)}}\frac{\lambda_{j{t^\prime}r}^{(\delta)}}{\mu_{j{t^\prime}r}^{(\delta) } }, \notag \\ 
 \lambda_{ipr}^{(\gamma)}=a^{(\gamma)}+\sum_{j=1}^{{\vert{V}\vert}}\sum_{t:{t^\prime}=p } \frac{A_{ijt}\theta_{ijt}\chi_{ijtr}}{1 - e^{-\theta_{ijt}}},\ 
 \mu_{ipr}^{(\gamma) } =   \frac{\nu_p^{(\gamma)}}{\xi_p^{(\gamma)}}+\frac{\lambda_{ir}^{(\alpha)}}{\mu_{ir}^{(\alpha)}}\sum_{j=1}^{{\vert{V}\vert } } \frac{\lambda_{jr}^{(\beta)}}{\mu_{jr}^{(\beta)}}\sum_{t:{t^\prime}=p } \frac{\lambda_{j{t^\prime}r}^{(\delta)}}{\mu_{j{t^\prime}r}^{(\delta)}},\end{gathered}$$ and like results can be obtained for the host - specific parameters $ \lambda_{ir}^{(\beta)}$, $ \mu_{jr}^{(\beta)}$, $ \lambda_{jpr}^{(\delta)}$ and $ \mu_{jpr}^{(\delta)}$. The update for $ \nu_i^{(\alpha)},\   \xi_i^{(\alpha)},\ \nu_j^{(\beta)}$ and $ \xi_j^{(\beta)}$ are identical to steps 7 and 8 in Algorithm   \[algo\_vi\ ]. For the covariates $ $ \begin{gathered } 
 \lambda_{kh
 thwt detailed in Algorithm \[xlgo\_vi\], where stgpw 7, 8, 9 aid 10 are modifiea to include the time dependxnt paramtners. It follows that wor the uder-specidic karameters the upvzte equations bake chx form: $$\begin{gatmered}
\lambda_{hr}^{(\alpha)} = a^{(\alphd)}+\sjm_{l=1}^{{\vert{V}\vert}}\sum_{t=1}^T \frac{A_{ijt}\theta_{ijt}\chi_{ittr}}{1-e^{-\theya_{ljt}}},\  
\mu_{ir}^{(\alpha)} = \frss{\nu_i^{(\zlpha)}}{\xi_i^{(\alpha)}}+\sum_{t=1}^T\frac{\lambda_{i{t^\prims}r}^{(\gamma)}}{\ku_{i{t^\prime}r}^{(\gamka)}}\sum_{j=1}^{{\vert{V}\vert}}\frac{\lambda_{jg}^{(\betw)}}{\mu_{jr}^{(\beta)}}\frac{\lambdw_{j{t^\prime}r}^{(\deota)}}{\mt_{h{t^\prime}r}^{(\deltx)}}, \notag \\
\lambda_{ipr}^{(\gamma)}=z^{(\gamma)}+\sum_{j=1}^{{\vert{V}\vert}}\sum_{t:{t^\prime}=p} \wrac{A_{njt}\theta_{ijt}\xhu_{ijht}}{1-e^{-\theta_{ijt}}},\ 
\mn_{ipr}^{(\gaima)} =  \frac{\nu_i^{(\bamma)}}{\xh_p^{(\gamma)}}+\grac{\lambda_{ir}^{(\aliha)}}{\mu_{mr}^{(\alpha)}}\sum_{j=1}^{{\vert{V}\vert}} \frac{\nambda_{jr}^{(\beta)}}{\mu_{jr}^{(\beja)}}\sum_{t:{t^\prike}='} \frac{\lambda_{j{t^\prime}r}^{(\dwlra)}}{\mu_{j{j^\prima}r}^{(\deura)}},\evd{gztiersd}$$ and sijilar resumts can be ibtained for the hoxt-fircific paramsters $\jaibda_{ir}^{(\beta)}$, $\mu_{jr}^{(\beta)}$, $\lambda_{jpr}^{(\delta)}$ and $\ku_{jlr}^{(\delta)}$. The updates for $\nu_i^{(\alpha)},\  \xi_i^{(\alpha)},\ \no_j^{(\beta)}$ and $\xi_j^{(\beta)}$ are identical to steps 7 and 8 in Algorithk \[algo\_ti\]. Fox the zivwriates $$\begin{gathered}
\lambda_{kh
that detailed in Algorithm \[algo\_vi\], where steps 9 10 are to include the that the user-specific parameters update equations take form: $$\begin{gathered} \lambda_{ir}^{(\alpha)} = a^{(\alpha)}+\sum_{j=1}^{{\vert{V}\vert}}\sum_{t=1}^T \frac{A_{ijt}\theta_{ijt}\chi_{ijtr}}{1-e^{-\theta_{ijt}}},\ = \frac{\nu_i^{(\alpha)}}{\xi_i^{(\alpha)}}+\sum_{t=1}^T\frac{\lambda_{i{t^\prime}r}^{(\gamma)}}{\mu_{i{t^\prime}r}^{(\gamma)}}\sum_{j=1}^{{\vert{V}\vert}}\frac{\lambda_{jr}^{(\beta)}}{\mu_{jr}^{(\beta)}}\frac{\lambda_{j{t^\prime}r}^{(\delta)}}{\mu_{j{t^\prime}r}^{(\delta)}}, \notag \\ \lambda_{ipr}^{(\gamma)}=a^{(\gamma)}+\sum_{j=1}^{{\vert{V}\vert}}\sum_{t:{t^\prime}=p} \frac{A_{ijt}\theta_{ijt}\chi_{ijtr}}{1-e^{-\theta_{ijt}}},\ \mu_{ipr}^{(\gamma)} = \frac{\nu_p^{(\gamma)}}{\xi_p^{(\gamma)}}+\frac{\lambda_{ir}^{(\alpha)}}{\mu_{ir}^{(\alpha)}}\sum_{j=1}^{{\vert{V}\vert}} \frac{\lambda_{jr}^{(\beta)}}{\mu_{jr}^{(\beta)}}\sum_{t:{t^\prime}=p} \frac{\lambda_{j{t^\prime}r}^{(\delta)}}{\mu_{j{t^\prime}r}^{(\delta)}},\end{gathered}$$ and similar can be obtained for the host-specific parameters $\lambda_{ir}^{(\beta)}$, $\mu_{jr}^{(\beta)}$, $\lambda_{jpr}^{(\delta)}$ and $\mu_{jpr}^{(\delta)}$. The for \xi_i^{(\alpha)},\ and are identical to steps 7 and 8 in Algorithm \[algo\_vi\]. For the covariates $$\begin{gathered} \lambda_{kh
 that detailed in Algorithm \[alGo\_vi\], where sTeps 7, 8, 9 aNd 10 aRe mOdIfieD to iNclude the time dEPendEnt parameters. It follows That fOr THe usER-sPecifIc paramETeRS The UpDaTe eQuATiOns taKe tHe form: $$\bEgin{gatherEd}
\lAmBda_{ir}^{(\alpha)} = a^{(\aLPhA)}+\sum_{j=1}^{{\vert{V}\VerT}}\sum_{t=1}^T \frac{A_{iJt}\tHeta_{ijT}\cHi_{iJTr}}{1-e^{-\thEta_{Ijt}}},\  
\mu_{Ir}^{(\alphA)} = \Frac{\nu_I^{(\alpha)}}{\xi_i^{(\AlPHa)}}+\sum_{t=1}^t\Frac{\lamBDA_{i{T^\priMe}r}^{(\gamma)}}{\mu_{i{t^\prime}R}^{(\GaMMa)}}\sum_{j=1}^{{\vert{V}\verT}}\frac{\lAmBDa_{JR}^{(\BetA)}}{\mu_{Jr}^{(\beta)}}\frac{\LaMbda_{j{T^\Prime}r}^{(\dELtA)}}{\MU_{J{t^\pRIme}r}^{(\delta)}}, \notaG \\
\lambda_{ipr}^{(\gAMma)}=A^{(\gamma)}+\SuM_{j=1}^{{\vERt{V}\verT}}\sum_{t:{T^\pRIme}=P} \frac{A_{ijt}\thEta_{iJt}\chi_{ijtr}}{1-E^{-\theta_{IJt}}},\ 
\mu_{ipr}^{(\GAmma)} =  \fraC{\nu_p^{(\gaMma)}}{\Xi_p^{(\GammA)}}+\FrAc{\LamBdA_{Ir}^{(\aLPhA)}}{\mu_{IR}^{(\alPha)}}\sum_{j=1}^{{\vErT{V}\Vert}} \fRac{\lAMBDA_{jr}^{(\bEta)}}{\Mu_{jr}^{(\Beta)}}\sUm_{t:{t^\prime}=p} \fraC{\laMbda_{J{T^\prIme}r}^{(\dElta)}}{\mU_{j{t^\pRiMe}r}^{(\deLta)}},\end{GatheReD}$$ and similar resuLts cAn be obtaiNed FoR thE hOst-spECific pAraMetErs $\lambDa_{ir}^{(\betA)}$, $\Mu_{jR}^{(\bETA)}$, $\LaMbda_{jpr}^{(\delta)}$ and $\mu_{jPr}^{(\DELtA)}$. The updaTes for $\NU_i^{(\AlPHa)},\  \xi_i^{(\alpHa)},\ \Nu_j^{(\Beta)}$ AND $\xi_j^{(\bEta)}$ aRE iDentical To stepS 7 AnD 8 iN AlgoriThM \[algo\_vI\]. FOr tHe cOvariATes $$\bEgin{gaThered}
\laMbda_{kH
 that detailed in Algorith m \[algo\_ vi\],  wh ere  s teps  7,8, 9 and 10 ar e  mod ified to include the t ime d ep e nden t  p arame ters. I t  f o l low sth atfo r  t he us er- specifi c paramete rsth e update equ a ti ons take t heform: $$\beg in{ gather ed }
\ l ambda _{i r}^{( \alpha ) } = a^ {(\alpha) }+ \ sum_{j = 1}^{{\v e r t{ V}\v ert}}\sum_{t=1}^T \f r ac{A_{ijt}\the ta_{ij t} \ ch i _ {ij tr} }{1-e^{-\t he ta_{i j t}}},\
\ m u _ {ir } ^{(\alpha)} =  \frac{\nu_ i ^{( \alpha )} }{\ x i_i^{( \alph a) } }+\ sum_{t=1}^T \fra c{\lambda _{i{t^ \ prime}r } ^{(\gam ma)}}{ \mu _{i {t^\ p ri me }r} ^{ ( \ga m ma )}} \ sum _{j=1}^{ {\ ve rt{V} \ver t } } \ frac {\l ambd a_{jr }^{(\beta)}}{ \mu _{jr } ^{( \beta )}}\f rac{ \l ambda _{j{t^ \prim e} r}^{(\delta)}}{ \mu_ {j{t^\pri me} r} ^{( \d elta) } }, \no tag  \\ 
\lambd a_{ipr} ^ {(\ ga m m a )} =a^{(\gamma)}+\sum _{ j = 1} ^{{\vert {V}\ve r t} }\ s um_{t:{t ^\ pri me}= p }  \fra c{A_ { ij t}\theta _{ijt} \ ch i_ {ijtr}} {1 -e^{-\ th eta _{i jt}}} , \ 
\ mu_{ip r}^{(\ga mma)} =  \frac{\nu_p ^ {(\gamma)}}{\ x i_ p ^ {( \ gamm a)} }+\frac{\la mbda _ {ir} ^{(\ a lp ha) } }{\mu _{ir} ^{ ( \a l pha)}}\sum_{j=1}^{{ \v ert{V} \vert }} \frac{\lam bda_{jr}^{ ( \ b eta)}}{\ mu_{ j r} ^ {(\beta)}}\sum _{t:{ t^\prime}= p } \frac{ \lamb da_{j{t^ \prime}r} ^ { (\delta) }}{ \mu _{j {t^ \ p ri me}r}^{(\delt a ) }},\ en d{gathe red }$$ and  si mil arres ul ts can be  obtaine dfo rth e h ost-s p ecific p ar ame te rs$\lam b da_{ir }^{(\ beta )} $, $\m u_{jr}^ { (\ b e ta)} $,  $ \lam bda _{ jpr}^ {(\d e lta )}$ and  $\mu_{jp r}^ { (\de lt a) }$. The  updates for$\ nu_i^{(\al ph a)} ,\  \x i _ i^{(\alp ha)},\ \nu_j^{(\beta)}$ and $\x i_j ^{(\b eta) }$ are id ent ical t o s t eps 7and 8in Al go rit h m  \[al g o \_ vi\ ].  For the c o v ari ates$$ \beg in{gath ered}
\lambda_{kh
 that_detailed in_Algorithm \[algo\_vi\], where steps 7,_8, 9_and_10 are_modified_to include the_time dependent parameters._It follows that for_the user-specific parameters_the_update equations take the form: $$\begin{gathered}
\lambda_{ir}^{(\alpha)} = a^{(\alpha)}+\sum_{j=1}^{{\vert{V}\vert}}\sum_{t=1}^T \frac{A_{ijt}\theta_{ijt}\chi_{ijtr}}{1-e^{-\theta_{ijt}}},\  
\mu_{ir}^{(\alpha)} = \frac{\nu_i^{(\alpha)}}{\xi_i^{(\alpha)}}+\sum_{t=1}^T\frac{\lambda_{i{t^\prime}r}^{(\gamma)}}{\mu_{i{t^\prime}r}^{(\gamma)}}\sum_{j=1}^{{\vert{V}\vert}}\frac{\lambda_{jr}^{(\beta)}}{\mu_{jr}^{(\beta)}}\frac{\lambda_{j{t^\prime}r}^{(\delta)}}{\mu_{j{t^\prime}r}^{(\delta)}}, \notag_\\
\lambda_{ipr}^{(\gamma)}=a^{(\gamma)}+\sum_{j=1}^{{\vert{V}\vert}}\sum_{t:{t^\prime}=p}_\frac{A_{ijt}\theta_{ijt}\chi_{ijtr}}{1-e^{-\theta_{ijt}}},\ 
\mu_{ipr}^{(\gamma)}_=__\frac{\nu_p^{(\gamma)}}{\xi_p^{(\gamma)}}+\frac{\lambda_{ir}^{(\alpha)}}{\mu_{ir}^{(\alpha)}}\sum_{j=1}^{{\vert{V}\vert}} \frac{\lambda_{jr}^{(\beta)}}{\mu_{jr}^{(\beta)}}\sum_{t:{t^\prime}=p} \frac{\lambda_{j{t^\prime}r}^{(\delta)}}{\mu_{j{t^\prime}r}^{(\delta)}},\end{gathered}$$ and similar_results can be obtained for_the host-specific_parameters $\lambda_{ir}^{(\beta)}$, $\mu_{jr}^{(\beta)}$, $\lambda_{jpr}^{(\delta)}$ and $\mu_{jpr}^{(\delta)}$. The updates_for_$\nu_i^{(\alpha)},\  \xi_i^{(\alpha)},\_\nu_j^{(\beta)}$ and $\xi_j^{(\beta)}$ are identical to steps 7 and_8 in Algorithm \[algo\_vi\]. For the covariates_$$\begin{gathered}
\lambda_{kh