Let M_{i}(t) be some summary of risk factors at time t for participant i from the longitudinal risk factor model. The hazard sub-model is h_{i}(t \vert M_{i}(t)) = h_{0}\exp\{\alpha^{'T}w_{i} + \alpha \int_{s=0}^{t} m_{i}(s) I(m_{i}(s) > \gamma)ds where \gamma is a risk factor threshold and linea…

Right. Here’s some stan code: functions{ real indicator(real x, real c){ real ans = x > c ? 1 : 0; return(ans); } real phi( real x, real x_c, array[] real theta, array[] real x_r, array[] int x_i ) { real ans = x^2 * exp(-x^2/2) * indicator(x, 0); return(ans); } …

Area above threshold feature

General

maxbiostat May 3, 2023, 8:07pm 3

Hi,

This thread might be of value. Just beware that the behaviour of integrate_1d might have changed in recent releases.

1 Like

Topic		Replies	Views
Joint Model Specification Modeling	3	745	November 8, 2019
Little question about the error of int ( In a large longitudinal model) General bugs , rstan , integrator	2	287	May 23, 2023
A matrix of parameters in Stan Integral Modeling fitting-issues	41	1434	November 25, 2019
Calculating hazard ratios and interpreting differences in survival and between a multi-level covariate from rstanarm::stan_surv rstanarm techniques , interpret-results , survival	3	1272	May 25, 2021
Estimate an indicator variable within a linear regression Modeling	8	1480	June 15, 2018

Area above threshold feature

Related Topics