У меня есть перекрестная классификация модели, как показано ниже.Я хотел бы представить случайное влияние некоммерческой деятельности на отраслевом уровне на гусеничный сюжет.Какие шаги?Спасибо!
M <- lmer(salary ~ 1+education+male+nonprofit+(1+nonprofit|ind)+(1+nonprofit|occ), data, reml=TRUE)