Compare colour classes of an RCOX model

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

1
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, sorenh@agrsci.dk

See Also

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

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
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")