Разделите несколько сюжетов на 4 разных сюжета

Question

Этот фрагмент:

library(stm)
gadarian <- gadarian
K<-c(5,10,15)
temp<-textProcessor(documents=gadarian$open.ended.response,metadata=gadarian)
out <- prepDocuments(temp$documents, temp$vocab, temp$meta)
documents <- out$documents
vocab <- out$vocab
meta <- out$meta
set.seed(02138)
K<-c(5,10,15)
df1 <- searchK(documents, vocab, K, prevalence=~treatment + s(pid_rep), data=meta)

создает график с помощью этой команды:

plot(df1)

Этот график содержит 4 понимания сюжета.

Как можно получить их на 4 разных участка с теми же метками, что и на одном объединенном участке?

Answer 1

Вы можете попробовать это:

plot(df1$results$K, df1$results$heldout, type = "p", main = "Held-Out Likelihood", 
    xlab = "Number of Topics (K)", ylab = "Held-Out Likelihood", ask = TRUE)
lines(df1$results$K, df1$results$heldout, lty = 1, col = 1)

plot(df1$results$K, df1$results$residual, type = "p", main = "Residuals", 
    xlab = "Number of Topics (K)", ylab = "Residuals", ask = TRUE)
lines(df1$results$K, df1$results$residual, lty = 1, col = 1)

plot(df1$results$K, df1$results$semcoh, type = "p", main = "Semantic Coherence", 
      xlab = "Number of Topics (K)", ylab = "Semantic Coherence", ask = TRUE)
lines(df1$results$K, df1$results$semcoh, lty = 1, col = 1)

plot(df1$results$K, df1$results$lbound, type = "p", main = "Lower Bound", 
     xlab = "Number of Topics (K)", ylab = "Lower Bound", ask = TRUE)
lines(df1$results$K, df1$results$lbound, lty = 1, col = 1)

Answer 2

Вы можете использовать что-то вроде

plot(df1$results$K,df1$results$heldout, type = "b", 
     xlab = "Number of Topics (K)", ylab = "Held-out Likelihood",
     main="Held-out Likelihood")

plot(df1$results$K,df1$results$semcoh, type = "b", 
     xlab = "Number of Topics (K)", ylab = "Semantic Coherence",
     main="Semantic Coherence")

plot(df1$results$K,df1$results$residual, type = "b", 
     xlab = "Number of Topics (K)", ylab = "Residual",
     main="Residual")

plot(df1$results$K,df1$results$lbound, type = "b", 
     xlab = "Number of Topics (K)", ylab = "Lower Bound",
     main="Lower Bound")