lm
step
\[ \large y_{i} = \beta_0 + \beta_1 x_{1i} + \beta_2 x_{2i} + \beta_3 x_{1i}x_{2i} + \epsilon_{i} \]
\[ \large fm \leftarrow lm(SBP \sim BMI + DM + BMI:DM, \space data=BP) \]
\[ \large anova(fm) \]
\[ \large summary(fm) \]