Я запускаю модерируемую регрессию в R. Я хочу проверить модерацию с помощью plotSlopes, но продолжаю сталкиваться с ошибками независимо от того, сколько проблем я снимаю. Да, это для Uni-назначения. Я не прошу толкований, я просто застрял на этом шаге, и мой гугл-фу в лучшем случае умеренный. Извините, если код выглядит неряшливо, но вот мои команды:
data <- data.frame(ebi.ges.t.z, pq.z, c2.s)
## n = 56
## relevant variables for moderated regression:
predictor <- data$ebi.ges.t.z
criterion <- data$c2.s
moderator <- data$pq.z
## moderated regression model using variables:
model <- lm(data$c2.s ~ data$ebi.ges.t.z * data$pq.z )
summary(model)
## analyzing the interaction between moderator and predictor:
slopes <- plotSlopes(model, plotx = "data$ebi.ges.t.z" , modx = "data$pq.z")
После некоторых настроек я получаю сообщение об ошибке: use modx argument to give name of the variable on the x axis или: undefined columns selected (когда не используется "data $ ..." "-Часть)
Что я и сделал. Я использовал modx = pq.z в качестве переменной модератора (по центру), ebi.ges.tz - центрированный предиктор. Я пробовал разные способы набора текста. Использование функции attach для присоединения набора данных и исключения «данных $ ...» из кода. Я пытался без attach -функции. Я пытался положить его в (), (""). Я просто конвертировал переменную, используя pq.z<-as.numeric(data$pq.z), и снова запустил модель. Я добавил plotx = data $ ebi.ges.tz во всех разновидностях (с (), "" и (""), с присоединением (данными), без, с данными $, без), я пробовал без команды plotx , Это всегда дает мне одну и ту же ошибку. Я не оказываю существенного влияния на коэффициенты, и сама модель не велика, но я получаю цифры. Переменная pq.z-moderator влияет на критерий, самый большой из протестированных мною. Коэффициент взаимодействия между предсказателем и модератором намного меньше. R не нравится это? Разве это не должно игнорировать такие неудобства и просто заставлять цифры давать мне что-то, чтобы я мог эмпирически сообщить, что это не работает? Мне нужен сюжет, чтобы я мог убедиться в отсутствии существенного взаимодействия и интерпретировать соответственно. Есть ли другая альтернатива анализу и составлению графика взаимодействия? Мой учебник дает только функцию plotSlopes.
Спасибо всем.
Вот значения переменных, которые я использовал (наблюдения 1-10)
### Values of the variables
ebi.ges.t.z <- [1] -4.3571429 13.6428571 -1.3571429 13.6428571 11.6428571 -4.3571429 0.6428571
[8] -3.3571429 -0.3571429 -2.3571429
pq.z <- [1] 0.42467286 0.42467286 -0.57404509 0.37467286 0.75928824 -0.52276304 -0.70225022
[8] -0.70225022 -0.70225022 -0.70225022
c2.s <- [,1]
[1,] -0.5661504112
[2,] 1.4740801521
[3,] 0.0036676094
[4,] -0.2867822459
[5,] 4.8415981390
[6,] -0.3031566664
[7,] 0.4700211885
[8,] -0.4694789013
[9,] -0.0396110456
[10,] -0.5479977102