join1split1: Joining and splitting of colour classes in RCOX models

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

Description

Test for joining of two colour classes (of a specific type) by testing if their corresponding parameters are not significantly different. Split a colour class and test how much this changes the fit of the model.

Usage

join1(object, scope=NULL, type = "ecc", details = 1,stat = "wald")
split1(object, scope=NULL, type = "ecc", details = 1) 

Arguments

object	An RCOX model, an object of class RCOX
scope	A specification of colour classes which should be considered for joining/splitting. If NULL, then all colour classes are considered.
type	Either "ecc" for edge colour classes or "vcc" for vertex colour classes.
stat	Either "wald" for a Wald statistic or "dev" for deviance statistic.
details	Control the amount of output

Value

A list with entries:

tab	A data frame with the test results
cc	A list of colour classes

Note

Note that the keyword 'stat' is not available for split1 because this function expands the current and hence the Wald statistic is not available. Note also that join1 is simply a wrapper for comparecc applied to edge colour classes.

Author(s)

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

References

~put references to the literature/web site here ~

See Also

rcox, update, comparecc

Examples

data(math)
g1     <- ~me:ve:al+al:st:an
m1     <- rcox(gm=g1, data=math)
join1(m1)

gm  = ~al:an:st
vcc = list(~me+st, ~ve+an)
ecc = list(~me:ve+me:al, ~ve:al+al:st)
m2 <- rcox(gm=gm, vcc=vcc, ecc=ecc, data=math, type="rcon")
split1(m2)

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