Create factorMarginalization.R