comparecc function

Compare colour classes of an RCOX model

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

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.

Returns

A list with entries: - tab: A data frame with the test results - cc1, cc2: Lists of colour classes

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.

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")

See Also

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

Author(s)

Søren Højsgaard, sorenh@math.aau.dk