comparecc: Compare colour classes of an RCOX model

In gRc: Inference in Graphical Gaussian Models with Edge and Vertex Symmetries

Description

A general function for pairwise comparisons of colour classes in an RCOX model, i.e. for testing whether the corresponding parameters are significantly different

Usage

comparecc(object, cc1 = NULL, cc2 = NULL, type = "ecc", stat = "wald", details = 1)

Arguments

object	An RCOX model, an object of class 'rcox'
cc1, cc2	Lists of colour classes of type 'type', see 'details' for an explanation of the defaults.
type	Either "ecc" for edge colour classes or "vcc" for vertex colour classes
stat	Base the comparison on either "wald" for a Wald statistic or "dev" for a deviance statistic
details	Control the amount of output created.

Details

All colour classes specified in cc1 are compared with all those given in cc2 (duplicate entries are not compared). If cc2=NULL (the default) then all colour classes specified in cc1 are compared with all colour classes in the model except those specified in cc1. If cc1=NULL (the default) and cc2=NULL then all pairwise comparisons are made.

Value

A list with entries:

tab	A data frame with the test results
cc1, cc2	Lists of colour classes

Author(s)

S<f8>ren H<f8>jsgaard, [e-mail]

See Also

add1.rcox, drop1.rcox, stepadd1, stepdrop1, join1, split1, stepjoin1, stepsplit1

Examples

data (math)

gm  = ~al:an:st
vcc = list(~me+st, ~ve+an, ~al)
ecc = list(~me:ve+me:al, ~ve:al+al:st)

m1 <- rcox(gm=gm, vcc=vcc, ecc=ecc, data=math)
m1

comparecc(m1, type="vcc")
comparecc(m1, type="ecc")

gRc documentation built on May 2, 2019, 5:22 p.m.