textmodel_wordfish(data_dfm_LBGexample, dir = c(1,5))
## Not run: ------------------------------------
# ie2010dfm <- dfm(data_corpus_irishbudget2010, verbose = FALSE)
# (wfm1 <- textmodel_wordfish(ie2010dfm, dir = c(6,5)))
# (wfm2a <- textmodel_wordfish(ie2010dfm, dir = c(6,5),
# dispersion = "quasipoisson", dispersion_floor = 0))
# (wfm2b <- textmodel_wordfish(ie2010dfm, dir = c(6,5),
# dispersion = "quasipoisson", dispersion_floor = .5))
# plot(wfm2a@phi, wfm2b@phi, xlab = "Min underdispersion = 0", ylab = "Min underdispersion = .5",
# xlim = c(0, 1.0), ylim = c(0, 1.0))
# plot(wfm2a@phi, wfm2b@phi, xlab = "Min underdispersion = 0", ylab = "Min underdispersion = .5",
# xlim = c(0, 1.0), ylim = c(0, 1.0), type = "n")
# underdispersedTerms <- sample(which(wfm2a@phi < 1.0), 5)
# which(featnames(ie2010dfm) %in% names(topfeatures(ie2010dfm, 20)))
# text(wfm2a@phi, wfm2b@phi, wfm2a@features,
# cex = .8, xlim = c(0, 1.0), ylim = c(0, 1.0), col = "grey90")
# text(wfm2a@phi[underdispersedTerms], wfm2b@phi[underdispersedTerms],
# wfm2a@features[underdispersedTerms],
# cex = .8, xlim = c(0, 1.0), ylim = c(0, 1.0), col = "black")
# if (require(austin)) {
# wfmodelAustin <- austin::wordfish(quanteda::as.wfm(ie2010dfm), dir = c(6,5))
# cor(wfm1@theta, wfmodelAustin$theta)
# }
## ---------------------------------------------
Run the code above in your browser using DataLab