Блог

Я пытаюсь сгенерировать графики условий взаимодействия на основе линейной модели

В категориальных переменных — предикторах есть некоторые недостатки — например, не каждый «region_new» присутствует для каждого «состояния», поэтому будут некоторые термины взаимодействия NA, например, из сводки модели:

 stateACT:region_newRegional            NA         NA      NA       NA 
stateNT:region_newCities               NA         NA      NA       NA       
stateACT:region_newRemote              NA         NA      NA       NA    
stateTAS:region_newRemote              NA         NA      NA       NA    
stateVIC:region_newRemote              NA         NA      NA       NA    
stateACT:region_newVery remote         NA         NA      NA       NA    
stateTAS:region_newVery remote         NA         NA      NA       NA    
stateVIC:region_newVery remote         NA         NA      NA       NA    
 

Это связано с тем, что в штате «ЗАКОН» нет данных по «Региональным» районам, в штате «NT» нет «Городов» и т. Д

Тем не менее, когда я использую функцию sjPlot::plot_model(…, type = ‘int’), как показано ниже, она отображает прогнозируемые значения для взаимодействия stateACT:region_newRegional , даже если оно должно быть пустым/отсутствующим (среди других терминов взаимодействия, которые также должны быть пустыми/отсутствующими)

 writing_lm <- lm(writing ~ (year   state   region_new   grade)^2, data = scores_df)

summary(writing_lm)


plot_model(writing_lm, type = "int",colors = "Set1")
 

График, созданный с предсказанными условиями взаимодействия, которые должны отсутствовать

Есть ли способ, которым я могу сказать sjPlot::plot_model , чтобы игнорировать/не отображать прогнозируемые значения для этих терминов там, где они должны отсутствовать? Я просмотрел документацию по упаковке от CRAN и искал здесь другие вопросы, но пока безрезультатно.