Hoewel de interactieterm niet significant is, zullen we deze toch in ons model laten. We zullen het marginale dosis-effect schatten. Dit is het gemiddelde effect over alle vissen heen en is gelijk aan het effect voor een vis met gemiddeld gewicht, m.a.w. \(\beta_d +\beta_{d:g}*\bar{x}_g\)
Eerst kunnen we dit doen op basis van de coëfficiënten die we vinden in de summary:
summary(lmGoudInt)
0.3691+wBar*(-0.5146)
#Call:
#lm(formula = log2minsurv ~ dosis * gewicht, data = poison %>%
# filter(soort == 1))
#
#Residuals:
# Min 1Q Median 3Q Max
#-0.70747 -0.27768 -0.06423 0.29306 0.91394
#
#Coefficients:
# Estimate Std. Error t value Pr(>|t|)
#(Intercept) -1.6571 2.7690 -0.598 0.554
#dosis 0.3691 1.6541 0.223 0.825
#gewicht 2.3961 1.2161 1.970 0.057 .
#dosis:gewicht -0.5146 0.7290 -0.706 0.485
#---
#Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
#
#Residual standard error: 0.3891 on 34 degrees of freedom
#Multiple R-squared: 0.6406, Adjusted R-squared: 0.6088
#F-statistic: 20.2 on 3 and 34 DF, p-value: 1.078e-07
#
#[1] -0.789865