score_copula {CopulaCenR}

Generalized score test for covariate effect(s)

Description

Generalized score test on covariate effect(s) under a fitted copula model.

Usage

score_copula(object, var_score)

Arguments

object

The output object from the main functions (rc_par_copula, rc_spCox_copula, ic_spTran_copula, ic_par_copula) under the null hypothesis

var_score

the list of covariates to be tested by the score test

Value

the score statistics, p value

Examples

# Score test for "rs2284665" in AREDS data
# fit a Copula2-semiparametric model under NULL
data(AREDS)
copula2_sp_null <- ic_spTran_copula(data = AREDS, copula = "Copula2",
                   l = 0, u = 15, m = 3, r = 3,
                   var_list = c("ENROLLAGE","SevScaleBL"))
score_copula(object = copula2_sp_null, var_score = "rs2284665")