showDeviance {CoCoRaw}R Documentation

Show the difference between likelihoods, dimension, etc.

Description

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.

Usage

showDeviance(model.1 = "current", model.2 = "base", data = NULL, 
             object = .object.of.models(model.1, model.2, data = data, ...), ...)

Arguments

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.

Value

The value TRUE is returned. If an error occurs then a warning is printed.

Author(s)

Jens Henrik Badsberg

References

Badsberg, J.H.: A guide to CoCo, JSS, 2001 ( http://www.jstatsoft.org/v06/i04/ ) and Badsberg, J.H.: Xlisp+CoCo, Aalborg University, 1996.

See Also

enterModel, makeModel, showModel, makeCurrent, returnModel, returnModelNumber, partialAssociations, optionsCoCo, showTest, returnTest, returnDeviance, backward and forward.

Examples


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 <- dynamic.Graph(fullModel, title = "Full");

backward(recursive = TRUE, object = CoCoObject);

# lastModel <- makeModel("last", object = CoCoObject);
# backwardGraph <- dynamic.Graph(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);

[Package Contents]