# NOT RUN {
## DAGs
dagNEL <- dag(~ a:b:c + c:d:e, result="graphNEL")
dagMAT <- dag(~ a:b:c + c:d:e, result="matrix")
dagMATS <- dag(~ a:b:c + c:d:e, result="dgCMatrix")
## Undirected graphs
ugNEL <- ug(~a:b:c + c:d:e, result="graphNEL")
ugMAT <- ug(~a:b:c + c:d:e, result="matrix")
ugMATS <- ug(~a:b:c + c:d:e, result="dgCMatrix")
## Is it a DAG?
is.DAG(dagNEL)
is.DAG(dagMAT)
is.DAG(dagMATS)
is.DAG(ugNEL)
is.DAG(ugMAT)
is.DAG(ugMATS)
## Is it an undirected graph
is.UG(dagNEL)
is.UG(dagMAT)
is.UG(dagMATS)
is.UG(ugNEL)
is.UG(ugMAT)
is.UG(ugMATS)
## Is it a triangulated (i.e. chordal) undirected graph
is.TUG(dagNEL)
is.TUG(dagMAT)
is.TUG(dagMATS)
is.TUG(ugNEL)
is.TUG(ugMAT)
is.TUG(ugMATS)
## Example where the graph is not triangulated
ug2NEL <- ug(~ a:b + b:c + c:d + d:a, result="graphNEL")
ug2MAT <- ug(~ a:b + b:c + c:d + d:a, result="matrix")
ug2MATS <- ug(~ a:b + b:c + c:d + d:a, result="dgCMatrix")
is.TUG(ug2NEL)
is.TUG(ug2MAT)
is.TUG(ug2MATS)
## Bidirected graphs
graph::edgemode(ugNEL)
graph::edgemode(ugNEL) <- "directed"
graph::edgemode(ugNEL)
is.DAG(ugNEL)
is.UG(ugNEL)
# }
Run the code above in your browser using DataLab