ModelSelection {gRcox}R Documentation

Stepwise selection of RCOX models

Description

~~ A concise (1-5 lines) description of what the function does. ~~

Usage

stepJoinVCC(m,  alpha = 0.05, details=TRUE,
                 criterion = c("test", "aic", "bic"), statistic =
                 c("wald", "dev"), display = FALSE)
stepJoinECC(m,  alpha = 0.05, details=TRUE,
                 criterion = c("test", "aic", "bic"), statistic =
                 c("wald", "dev"), display = FALSE)
stepJoinCC(m, type = "ecc", alpha = 0.05, details=TRUE,
                 criterion = c("test", "aic", "bic"), statistic =
                 c("wald", "dev"), display = FALSE)
stepSplitVCC(m,  alpha = 0.05,details=TRUE,
                 criterion = c("test", "aic", "bic"), display = FALSE)
stepSplitECC(m,  alpha = 0.05,details=TRUE,
                 criterion = c("test", "aic", "bic"), display = FALSE)
stepSplitCC(m, type = "ecc", alpha = 0.05,details=TRUE,
                 criterion = c("test", "aic", "bic"), display = FALSE)
stepDropECC(m, alpha = 0.05, statistic = c("wald", "dev"),
                 criterion = c("test", "aic", "bic"), details = TRUE,  display = FALSE)
stepAddECC(m, alpha = 0.05, criterion = c("test", "aic", "bic"), details = TRUE, display = FALSE)

Arguments

m An RCOX model object
alpha Significance level
type Either 'ecc' (edge colour class) or 'vcc' (vertex colour class)
details Show details of the tests
statistic Use Wald or likelihood ratio test
criterion Either 'test', 'aic' or 'bic'
display If TRUE and if the Rgraphviz package is loaded then the models are displayed graphically.

Details

stepSplitVCC and stepSplitECC are simply wrappers for stepSplitCC and likewise for stepJoinVCC, stepJoinECC and stepJoinCC.

Value

An RCOX model object.

Author(s)

Søren Højsgaard, sorenh@agrsci.dk

Examples

data(math)
bufly0  <- rcox(vcc=list(~me, ~ve, ~al, ~an, ~st), 
             ecc=list(~me:ve, ~me:al, ~ve:al, ~al:an, ~al:st, ~an:st), 
             data=math)

bufly01<-stepJoinVCC(bufly0,statistic='dev')

bufly02<-stepJoinECC(bufly01, details=FALSE, statistic='dev')

bufly2  <- rcox(vcc=list(~al, ~me+st, ~ve+an), 
              ecc=list(~me:ve+me:al+ve:al, ~al:an+al:st+an:st), 
              data=math)
stepSplitECC(bufly2)

[Package gRcox version 0.1 Index]