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.
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