Melt <- melt(t(model$phi[, , k1] - model$phi[, , k2]))
gDiff <- ggplot(data = Melt, aes(x = Var1, y = Var2, fill = value))
+ geom_tile()
Melt <- melt(t(model$phi[, , k1] - model$phi[, , k2]))
gDiff <- ggplot(data = Melt, aes(x = Var1, y = Var2, fill = value))
+ geom_tile()