GGMselect (version 0.1-0)

GGMselect-package: Gaussian Graphs Models selection

Description

GGMselect is a package dedicated to graph estimation in Gaussian Graphical Models. The main functions return the adjacency matrix of an undirected graph estimated from a data matrix.

This package is developed in the Applied Mathematics and Informatics (http://www.jouy.inra.fr/mia_eng/) Lab of INRA - Jouy-en-Josas, France.

To cite GGMselect, please use citation("GGMselect").

Details

ll{ Package: GGMselect Version: 0.1-0 Date: 2009-07-06 Depends: mvtnorm, lars, mvbutils (>= 2.0), gtools Suggests: network, glasso License: GPL (>= 3) URL: http://www.inra.fr/miaj/public/logiciels/GGMselect }

References

Please use citation("GGMselect").

See Also

selectFast, selectQE, selectMyFam,convertGraph, simulateGraph, penalty

Examples

p=30
n=30
# simulate graph
eta=0.11
Gr <- simulateGraph(p,eta)
# simulate data
X <- rmvnorm(n, mean=rep(0,p), sigma=Gr$C)
# estimate graph
GRest <- selectFast(X)

# plot result
library(network)
par(mfrow=c(1,2))
gV <- network(Gr$G)
plot(gV,jitter=TRUE, usearrows = FALSE, label=1:p,displaylabels=TRUE)
g <- network(GRest$EW$G)
plot(g, jitter=TRUE, usearrows = FALSE,  label=1:p,displaylabels=TRUE)