Добавление условий взаимодействия для отображения окончательных значений на лесных участках? (sjPlot и ggeffects) - PullRequest
0 голосов
/ 02 августа 2020

Итак, у меня есть два графика ниже: один - это коэффициенты для значений преступлений против собственности (1-й уровень), а второй, справа, - коэффициенты для насильственных преступлений. Однако коэффициенты справа - это только разница между коэффициентами преступлений против собственности и насильственных преступлений, а не фактические коэффициенты насильственных преступлений. Мне интересно, есть ли способ добавить коэффициенты преступлений против собственности к коэффициентам насильственных преступлений, чтобы я мог получить фактические коэффициенты насильственных преступлений.

Вот моя модель:

lmer(log(COUNT + 1) ~ CRIME*( cent.log.pop + cent.log.pop.dens + per_capita_income*black + per_capita_income*white + no.grad.hs  + prop5.17.pov + ma.plus + hs + ba + median_gross_rent+cent.EXP_STUDENT +unemployment_rate + asian +  diff.dem) + officers  + (year|PLACE_ID) + (1|COUNTY_ID) + (1|STATE), control = lmerControl(optimizer= "nloptwrap", calc.derivs = FALSE), na.action = 'na.omit', REML = FALSE, city.v.p = data)

введите описание изображения здесь

plot_model(model, rm.terms = c("cent.log.pop", "CRIMECRIME_VIOLENT:cent.log.pop", "CRIMECRIME_VIOLENT:cent.log.pop.dens", "CRIMECRIME_VIOLENT:per_capita_income", "CRIMECRIME_VIOLENT:black", "CRIMECRIME_VIOLENT:white", "CRIMECRIME_VIOLENT:no.grad.hs", "CRIMECRIME_VIOLENT:prop5.17.pov", "CRIMECRIME_VIOLENT:ma.plus", "CRIMECRIME_VIOLENT:hs", "CRIMECRIME_VIOLENT:ba", "CRIMECRIME_VIOLENT:median_gross_rent", "CRIMECRIME_VIOLENT:cent.EXP_STUDENT", "CRIMECRIME_VIOLENT:unemployment_rate", "CRIMECRIME_VIOLENT:asian", "CRIMECRIME_VIOLENT:diff.dem", "CRIMECRIME_VIOLENT:per_capita_income:black", "CRIMECRIME_VIOLENT:per_capita_income:white", "CRIMECRIME_VIOLENT"), group.terms = c(1,1,2,3,3,1,2,1,3,1,1,3,1,2,3,3,3), axis.labels = labels, title = "9a: Effect Size:  Property Crime", dot.size = .5,line.size = .2)...
...