caseonly {TwoPhaseInd}

A function to deal with case-only designs

Description

This function estimates parameters of case-only designs.

Usage

caseonly(data, treatment, BaselineMarker, extra = NULL, fraction = 0.5)

Arguments

data

A data frame used to access the following data.

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.

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 included in case-only regression. Note that extra covariates are not needed unless the interactions of treatment and extra coviarates are of interest.

fraction

The randomization fraction of active treatment assignment.

Details

This function estimates parameters of case-only designs. It estimates two parameters for "treatment effect when baselineMarker=0"" and treatment+baselineMarker interaction".

Value

For each paramter, it returns:

beta

Estimated parameter

stder

Standard error

pVal

p value

Author(s)

James Y. Dai

References

J. Y. Dai, S. S. Li, and P. B. Gilbert. Case-only methods for competing risks models with application to assessing differential vaccine efficacy by viral and host genetics. Biometrics, 15(1):196-203, 2014.

Examples

#form the data
data(acodata)
cdata=acodata[acodata[,2]==1,]
cfit=caseonly(data=cdata,
              treatment="f_treat",
              BaselineMarker="fcgr2a.3",
              extra=c("f_agele30","f_hsv_2","f_ad5gt18","f_crcm",
              "any_drug","num_male_part_cat","uias","uras"))
cfit