## ----include=FALSE, cache=FALSE--------------------------------------------
library(CoGAPS)

## --------------------------------------------------------------------------
data(SimpSim)
simulationName <- createGWCoGAPSSets(SimpSim.D, SimpSim.S, nSets=2, "example")

## --------------------------------------------------------------------------
result <- GWCoGAPS(simulationName, nFactor=3, nEquil=500, nSample=500,
    nOutput=250)

## --------------------------------------------------------------------------
plotGAPS(result$Amean, result$Pmean)