showDeviance {CoCoRaw} | R Documentation |
Print on standard output the test statistics (deviance and F-test statistics) with p-values based on likelihoods and dimensions of the two argument models. For discrete models the adjustment of the degrees of freedom is also reported.
showDeviance(model.1 = "current", model.2 = "base", data = NULL, object = .object.of.models(model.1, model.2, data = data, ...), ...)
model.1 |
See isSubmodel . |
model.2 |
See isSubmodel . |
data |
See exportCoCo . |
object |
See exportCoCo . |
... |
Additional arguments to generate the CoCo object
from the data argument. See propertyModel . |
The value TRUE is returned. If an error occurs then a warning is printed.
Jens Henrik Badsberg
Badsberg, J.H.: A guide to CoCo, JSS, 2001 ( http://www.jstatsoft.org/v06/i04/ ) and Badsberg, J.H.: Xlisp+CoCo, Aalborg University, 1996.
enterModel
, makeModel
,
showModel
, makeCurrent
,
returnModel
, returnModelNumber
,
partialAssociations
, optionsCoCo
,
showTest
,
returnTest
, returnDeviance
,
backward
and forward
.
library(MASS); data(anorexia); library(CoCoCg); # showDeviance(model.1 ="[[:Prewt:Postwt][:Treat:Postwt]];", # model.2 ="*;", data = anorexia); # endCoCo(); # CoCoObject <- makeCoCoCg(); enterDataFrame(anorexia, object = CoCoObject); fullModel <- makeModel(enterModel("*", object = CoCoObject)); # fullGraph <- dG(fullModel, title = "Full"); backward(recursive = TRUE, object = CoCoObject); # lastModel <- makeModel("last", object = CoCoObject); # backwardGraph <- dG(lastModel, title = "Last"); showModel("all", object = CoCoObject); makeCurrent(object = CoCoObject); # These two calls does not produce the same deviance, # since the second correctly tests the inhomogene model: showTest(object = CoCoObject); showDeviance(object = CoCoObject); # The default models for showTest and showDeviance are "current" and "base", # but the default models for returnTest and returnDeviance was/are FALSE, # which both will refer to "base" when no models are set: returnTest(object = CoCoObject); # model.1 was/is by default "base"!!! returnTest(FALSE, FALSE, object = CoCoObject); returnTest("base", "base", object = CoCoObject); returnTest("current", object = CoCoObject); returnTest("current", "base", object = CoCoObject); returnTest(model.2 = 1, object = CoCoObject); returnTest("base", "current", object = CoCoObject); returnDeviance(object = CoCoObject); # model.1 was/is by default "base"!!! returnDeviance(FALSE, FALSE, object = CoCoObject); returnDeviance("base", "base", object = CoCoObject); returnDeviance("current", object = CoCoObject); returnDeviance("current", "base", object = CoCoObject); returnDeviance(model.2 = 1, object = CoCoObject); returnDeviance("base", "current", object = CoCoObject); showTest(3, 1, object = CoCoObject); showDeviance(3, 1, object = CoCoObject); returnTest(3, 1, object = CoCoObject); returnDeviance(3, 1, object = CoCoObject); showTest(1, 1, object = CoCoObject); showDeviance(1, 1, object = CoCoObject); returnTest(1, 1, object = CoCoObject); returnDeviance(1, 1, object = CoCoObject); # inhoModel <- makeModel(enterModel("[[:Treat]] / [[:Prewt][:Treat:Postwt]] / [[:Treat:Prewt:Postwt]]", object = CoCoObject)); # inhoModel <- makeModel(enterModel("[[:Treat]] / [[:Prewt][:Treat:Postwt]] / [:Prewt^][:Treat:Prewt:Postwt][:Treat:Postwt^]]", object = CoCoObject)); # inhoModel <- makeModel(enterModel("[[:Treat]] / [[:Prewt][:Treat:Postwt]] / [:Prewt^][:Prewt:Postwt][:Treat:Postwt^]]", object = CoCoObject)); inhoModel <- makeModel(enterModel("[[:Prewt:Postwt][:Treat:Postwt]]", object = CoCoObject)); showModel("all", object = CoCoObject); makeCurrent(object = CoCoObject); showTest(object = CoCoObject); showDeviance(object = CoCoObject); returnTest("current", "base", object = CoCoObject); returnDeviance("current", "base", object = CoCoObject); exportCoCo("anorexia.xpt"); endCoCo(object = CoCoObject);