## Graphical Gaussian model
data(carcass)
cm1<-cmod(~.^., carcass)
## Stepwise selection based on BIC
cm2<-backward(cm1,k=log(nrow(carcass)))
## Stepwise selection with fixed edges
cm3<-backward(cm1,k=log(nrow(carcass)),
fixinMAT=matrix(c("LeanMeat","Meat11","Meat12","Meat13","LeanMeat","Fat11","Fat12","Fat13"),
ncol=2))
Run the code above in your browser using DataLab