spatialIC: PH model for spatial general interval-censored data
In PICBayes: Bayesian Models for Partly Interval-Censored Data

Description Usage Arguments Details Value Author(s) References

Fit a Bayesian semiparametric PH model with spatial frailty for spatially dependent general interval-censored data.

1
2
3

spatialIC(L, R, y, xcov, IC, scale.designX, scaled, area, binary, I, C, nn, 
order, knots, grids, a_eta, b_eta, a_ga, b_ga, a_lamb, b_lamb, beta_iter, 
phi_iter, beta_cand, beta_sig0, x_user, total, burnin, thin, conf.int, seed)

`L`	The vector of left endpoints of the observed time intervals.
`R`	The vector of right endponts of the observed time intervals.
`y`	The vector of censoring indicator: 0=left-censored, 1=interval-censored, 2=right-censored, 3=exact.
`xcov`	The covariate matrix for the p predictors.
`IC`	The vector of general interval-censored indicator: 1=general interval-censored, 0=exact.
`scale.designX`	The TRUE or FALSE indicator of whether or not to scale the design matrix X.
`scaled`	The vector indicating whether each covariate is to be scaled: 1=to be scaled, 0=not.
`area`	The vector of area ID.
`I`	The number of areas.
`C`	The adjacency matrix.
`nn`	The vector of number of neighbors for each area.
`binary`	The vector indicating whether each covariate is binary.
`order`	The degree of basis I-splines: 1=linear, 2=quadratic, 3=cubic, etc.
`knots`	A sequence of knots to define the basis I-splines.
`grids`	A sequence of points at which baseline survival function is to be estimated.
`a_eta`	The shape parameter of Gamma prior for `gamma_l`.
`b_eta`	The rate parameter of Gamma prior for `gamma_l`.
`a_ga`	The shape parameter of Gamma prior for `e^{beta_r}`.
`b_ga`	The rate parameter of Gamma prior for `e^{beta_r}`.
`a_lamb`	The shape parameter of Gamma prior for spatial precision `lambda`.
`b_lamb`	The rate parameter of Gamma prior for spatial precision `lambda`.
`beta_iter`	The number of initial iterations in the Metropolis-Hastings sampling for `beta_r`.
`phi_iter`	The number of initial iterations in the Metropolis-Hastings sampling for `phi_i`.
`beta_cand`	The sd of the proposal normal distribution in the MH sampling for `beta_r`.
`beta_sig0`	The sd of the prior normal distribution for `beta_r`.
`x_user`	The user-specified covariate vector at which to estimate survival function(s).
`total`	The number of total iterations.
`burnin`	The number of burnin.
`thin`	The frequency of thinning.
`conf.int`	The confidence level of the CI for `beta_r`.
`seed`	A user-specified random seed.

The baseline cumulative hazard is approximated by a linear combination of I-splines:

sum_{l=1}^{K}(gamma_l*b_l(t)).

For a binary prdictor, we sample e^{beta_r}, with Gamma prior.

The regression coefficient beta_r for a continuous predictor is sampled using MH algorithm. During the initial beta_iter iterations, sd of the proposal distribution is beta_cand. Afterwards, proposal sd is set to be the sd of available MCMC draws.

a list containing the following elements:

`N`	The sample size.
`parbeta`	A `total` by `p` matrix of MCMC draws of `beta_r`, r=1, ..., p.
`parsurv0`	A `total` by `length(grids)` matrix, each row contains the baseline survival at `grids` from one iteration.
`parsurv`	A `total` by `length(grids)*G` matrix, each row contains the survival at `grids` from one iteration. G is the number of sets of user-specified covariate values.
`parphi`	A `total` by `I` matrix of MCMC draws of `phi_i`, i=1,...,I.
`parlamb`	A `total` by 1 matrix of MCMC draws of `lambda`.
`coef`	A vector of regression coefficient estimates.
`coef_ssd`	A vector of sample standard deviations of regression coefficient estimates.
`coef_ci`	The credible intervals for the regression coefficients.
`S0_m`	The estimated baseline survival at `grids`.
`S_m`	The estimated survival at `grids` with user-specified covariate values `x_user`.
`grids`	The sequance of points where baseline survival functions is estimated.
`DIC`	Deviance information criterion.
`NLLK`	Negative log pseudo-marginal likelihood.

Chun Pan

Pan, C. and Cai, B. (2020). A Bayesian model for spatial partly interval-censored data. Communications in Statistics - Simulation and Computation, DOI: 10.1080/03610918.2020.1839497.

PICBayes documentation built on Aug. 5, 2021, 9:06 a.m.