Структурное моделирование темы (stm) Ошибка в makeTopMatrix (распространенность, данные): Ошибка создания матрицы модели

Question

Я пытаюсь выполнить начальные шаги этого урока stm

https://github.com/dondealban/learning-stm

с этим набором данных, он является частью исходного

http://www.mediafire.com/file/1jk2aoz4ac84jn6/data.csv/file

install.packages("stm")
library(stm)
load("VignetteObjects.RData") 
data <- read.csv("C:/data.csv") 
head(data)
processed <- textProcessor(data$documents, metadata=data)
out <- prepDocuments(processed$documents, processed$vocab, processed$meta)
docs <- out$documents
vocab <- out$vocab
meta <- out$meta


poliblogPrevFit <- stm(out$documents, out$vocab, K=4, prevalence=~rating+s(day), 
                       max.em.its=200, data=out$meta, init.type="Spectral", 
                       seed=8458159)

Но я продолжаю получать ту же ошибку

Error in makeTopMatrix(prevalence, data) : Error creating model matrix.
                 This could be caused by many things including
                 explicit calls to a namespace within the formula.
                 Try a simpler formula.

Пожалуйста, кто-нибудь может запустить его в 64-битной версии MS Windows R-3.5.2 .. Я не мог даже найти подобные ошибкив любом месте ..

Answer 1

Кажется, ваша проблема заключалась в том, что с выбранной выборкой у вас получился фактор-объект с одним уровнем:

> levels(meta$rating)
[1] "Conservative"

Использование такой переменной не имеет никакого смысла, посколькунет вариаций между делами.Если вы используете исходные данные, ваш код работает абсолютно нормально:

data <- read.csv("https://raw.githubusercontent.com/dondealban/learning-stm/master/data/poliblogs2008.csv")

processed <- textProcessor(data$documents, metadata = data)
out <- prepDocuments(processed$documents, processed$vocab, processed$meta)
docs <- out$documents
vocab <- out$vocab
meta <- out$meta

levels(meta$rating)
[1] "Conservative" "Liberal" 

poliblogPrevFit <- stm(docs, vocab, K = 4, prevalence = ~rating+s(day), 
                       max.em.its = 200, data = out$meta, init.type = "Spectral", 
                       seed = 8458159)