join-split | R Documentation |
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.
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 |
A list with entries:
tab |
A data frame with the test results |
cc |
A list of colour classes |
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.
Søren Højsgaard, sorenh@math.aau.dk
rcox
, update
,
comparecc
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)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.