BICcovar: BIC with covariates

Description Usage Arguments Value Examples

Description

Calculates BIC for the given network and data frame.

Usage

1
BICcovar(nw, dataset)

Arguments

nw

mhnet network

dataset

data frame with attributes given as 'covs' attribute

Value

real number

Examples

1
2
3
4
5
6
7
8
9
data(Gm17403)
data <- Gm17403[,5:7]
covs <- list(Sex = cbind(as.numeric(Gm17403$sex=="M")),
             Diet = cbind(as.numeric(Gm17403$diet=="HF")),
             Batch = cbind(as.numeric(Gm17403$batch==1)))
attr(data, "covs") <- covs

nw.start <- empty.network(c(names(attr(data, "covs")), names(data)))
BICcovar(nw.start, data)

simecek/mhnet documentation built on May 29, 2019, 10:01 p.m.