Description Usage Arguments Details Value Author(s) References See Also Examples
The function implements the MCMC algorithm with data augmentation to estimate the parameters in the zero-inflated Poisson model while correcting for the measurement error arising from the responses. The function returns the trace of the sampled parameters in each interaction. To obtain the summary estimation, use summary
().
1 2 3 4 5 6 7 8 | ZIPInt (Ystar, Covarmainphi, Covarmainmu, Covarplus, Covarminus,
Ystarval, Yval, Covarvalmainphi, Covarvalmainmu, Covarvalplus, Covarvalminus,
betaphi, betamu, alphaplus, alphaminus,
Uibound = c(7,11),
priorgamma, priormu, priorSigma,
propsigmaphi, propsigmamu = propsigmaphi,
propsigmaplus = propsigmaphi, propsigmaminus = propsigmaphi,
seed = 1, nmcmc = 500)
|
Ystar |
a count vector of length n specifying the error-prone response in the zero-inflated Poisson model. |
Covarmainphi |
a n \times p_1 dimensional data.frame or matrix of the covariate data with respect to the probability component of the zero-inflated Poisson model |
Covarmainmu |
a n \times p_2 dimensional data.frame or matrix of the covariate data with respect to the mean component of the zero-inflated Poisson model |
Covarplus |
a n \times q_1 dimensional data.frame or matrix of the covariate data for the measurement error model of the add-in error process |
Covarminus |
a n \times q_2 dimensional data.frame or matrix of the covariate data for the measurement error model of the leave-out error process |
Ystarval |
a count vector of length m specifying the error-prone response in the validation data. |
Yval |
a count vector of length m specifying the precisely measured response in the validation data. |
Covarvalmainphi |
a m \times p_1 dimensional data.frame or matrix of the covariate for validation data with respect to the probability component of ZIP model |
Covarvalmainmu |
a m \times p_2 dimensional data.frame or matrix of the covariate for validation data with respect to the mean component of ZIP model |
Covarvalplus |
a m \times q_1 dimensional data.frame or matrix of the covariate for validation data of the add-in error process |
Covarvalminus |
a m \times q_2 dimensional data.frame or matrix of the covariate for validation data of the leave-out error process |
betaphi |
a vector of length p_1 specifying the initial values of the parameters in the probability component of the zero-inflated Poisson model |
betamu |
a vector of length p_2 specifying the initial values of the parameters in the probability component of the zero-inflated Poisson model |
alphaplus |
a vector of length q_1 specifying the initial values of the parameters for the measurement error model of the add-in error process |
alphaminus |
a vector of length q_2 specifying the initial values of the parameters in the probability component of the leave-out error process |
Uibound |
a vector of length 2 specifying the maximum number of the count in the inverse sampling method |
priorgamma |
a vector of length 2 specifying the two parameters of gamma prior |
priormu |
a vector of length q_2 specifying the mean vector of the normal prior for the measurement error model of the leave-out error process |
priorSigma |
a vector of length q_2 specifying the standard erros of the normal prior for the measurement error model of the leave-out error process |
propsigmaphi |
a vector of length p_1 specifying the standard error of the Gaussian proposal distribution for the parameters corresponds to the probability component. |
propsigmamu |
a vector of length p_2 specifying the standard error of the Gaussian proposal distribution for the parameters corresponds to the mean component. |
propsigmaplus |
a vector of length q_1 specifying the standard error of the Gaussian proposal distribution for the parameters corresponds to the add-in error process. |
propsigmaminus |
a vector of length q_1 specifying the standard error of the Gaussian proposal distribution for the parameters corresponds to the leave-out error process. |
seed |
a numeric value specifying the seed for random generator |
nmcmc |
a integer specify the number of the generation of MCMC algorithm |
Comparing to the ZIPExt function for the external validation study, this function has an addition component – covariates of response model in the validation data. Here, the argument “Covarvalmainphi” and “Covarvalmainmu” are new for the sceanrio with external validation.
BayesResults |
the list of trace of generated parameters for each component of the models. Data frame "betaphi_trace" corresponds to the probability component of ZIP response model; "betamu_trace" refers to the mean component of the ZIP response model. Data frames "alphaplus_trace" and "alphaminus_trace", respectively, correspond to the add-in error and leave-out error process in the measruement error model. |
Qihuang Zhang and Grace Y. Yi
Zhang, Qihuang. "Inference Methods for Noisy Correlated Responses with Measurement Error." (2020).
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 | data(datasimInt)
set.seed(0)
result <- ZIPInt(Ystar = datasimInt$main$Ystar,
Covarmainphi = datasimInt$main[,c("intercept","X1")],
Covarmainmu = datasimInt$main[,c("intercept","X2")],
Covarplus = datasimInt$main[,c("intercept", "Zplus")],
Covarminus = datasimInt$main[,c("intercept", "Zminus")],
Ystarval = datasimInt$val$Ystar,
Yval = datasimInt$val$Y,
Covarvalmainphi = datasimInt$val[,c("intercept","X1")],
Covarvalmainmu = datasimInt$val[,c("intercept","X2")],
Covarvalplus = datasimInt$val[,c("intercept", "Zplus")],
Covarvalminus = datasimInt$val[,c("intercept", "Zminus")],
betamu = c(-0.5,0.5), betaphi = c(0.5,0),
alphaplus = c(0,0), alphaminus= c(0,0),
priorgamma = c(1,1), priormu = c(0,0), priorSigma = c(1,1),
propsigmaphi = c(0.05,0.05), nmcmc = 10)
summary(result)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.