aco2arm | R Documentation |
This function estimates parameters of proportional hazards model with gene-treatment interaction. It employs case-cohort estimation incorporating the case-only estimators. The method was published in Dai et al. (2015) Biometrics.
aco2arm(data, svtime, event, treatment, BaselineMarker, subcohort=NULL, esttype = NULL, weight=NULL, extra=NULL)
data |
A data frame used to access the following data. |
svtime |
A character string of column name, corresponds to one column of the data frame, which is used to store the failure time variable (numeric). |
event |
A character string of column name, corresponds to one column of the data frame, which is used to store the indicator of failure event (1: failure, 0: not failure). |
treatment |
A character string of column name, corresponds to one column of the data frame, which is used to store the binary vector of treatment variable (1: treatment, 0: placebo). |
BaselineMarker |
A character string of column name, corresponds to one column of the data frame, which is used to store a vector of biomarker. |
subcohort |
A character string of column name, corresponds to one column of the data frame, which is used to store the indicator of sub-cohort in the case-cohort sampling (1: sample belong to the sub-cohort, 0: not belong to the sub-cohort). In case-control sampling, this variable is set to be NULL. |
esttype |
The option of estimation methods (1: Self-Prentice estimator, 0: Lin-Ying estimator). |
weight |
If the genotype data are obtained through case-control sampling, weight is a vector of sampling weights (inverse of sampling probability) corresponding to rows of data. If the genotype data are obtained through case-cohort sampling, weight is NULL. If a vector of weights have been supplied by user, then esttype is automatically set to 0: Lin-Ying estimator. |
extra |
A string vector of column name(s), corresponds to more or more column(s) of the data frame, which is/are used to store the extra baseline covariate(s) to be adjusted for in addition to treatment and biomarker. |
The function returns estimates of the proportional hazards model, and variance of the estimates. The method was published in Dai et al. (2016) Biometrics.
A list of estimates and variance of the estimates.
Estimate |
A data frame of beta(Estimated parameter), stder(Standard error),and pVal(p value) |
Covariance |
covariance data frame of genotype,treatment,and interaction |
James Y. Dai
J. Y. Dai, X. C. Zhang,C. Y. Wang, and C. Kooperberg. Augmented case-only designs for randomized clinical trials with failure time endpoints. Biometrics, DOI: 10.1111/biom.12392, 2016.
aco1arm
## Load the example data data(acodata) ## Case-cohort + case-only estimators rfit1 <- aco2arm(data=acodata, svtime="vacc1_evinf", event="f_evinf", treatment="f_treat", BaselineMarker="fcgr2a.3", subcohort="subcoh", esttype=1, weight=NULL, extra=c("f_agele30","f_hsv_2","f_ad5gt18","f_crcm","any_drug", "num_male_part_cat","uias","uras")) rfit1
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.