Строки регрессии участка с взаимодействием в R

Question

Я пытаюсь получить график диаграммы рассеяния, показывающий 2 линии регрессии, чтобы эффекты взаимодействия были ясными. Мой код выглядит следующим образом:

lm <- lm(Cage_Int ~ Sex + AVP_mvBNST_ros_caud + Sex:AVP_mvBNST_ros_caud, data=data) 

plot(data$Cage_Int, data$AVP_mvBNST_ros_caud, xlab="V1aR in BNSTmv", 
    ylab="Interaction time", col=colors, pch=20) 
curve(cbind(x,1,x,1*x) %*% coef(lm), add=TRUE, col="red") 
curve(cbind(1,0,x,0*x) %*% coef(lm), add=TRUE, col="blue") 

это дает мне право диаграммы рассеяния график, но линии регрессии не показывают, что это не то, что я хочу (я хочу обе линии регрессии, чтобы показать, один для секса = 0 и один для секса = 1)

Большое спасибо!

Answer 1

Запустите регрессию с помощью sex=0 и sex=1 и создайте две отдельные линейные модели, используя lm(). После того, как вы определили две модели model1 = lm(<insert equation>) и model2 = lm(<insert second equaiton>), вы можете просто объявление линию с использованием abline() или, в данном случае:

abline(model1) abline(model2)

после функции сюжета и что следует сделать трюк. Вы можете добавить цвета и различные линейные эффекты, используя abline().