Here is the R output for the multiple regression model using anxiety level and age group as predictors.

Call: Im(formula = satisfaction ~ anxiety + AgeGroup, data = patients) Residuals : Min 1Q Median 30 Max -20.525 -7.600 1.469 7.496 19.975 Coefficients : Estimate Std. Error t value Pr(>Itl) (Intercept) 129 .041 12.563 10.272 5. 01e-13 *** anxiety -25 . 007 5. 700 -4.387 7.57e-05 * * * AgeGroupage35-49 -15.549 3.457 -4. 498 5.33e-05 *** AgeGroupage50-64 -24.441 5.305 -4. 607 3.77e-05 * ** Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '. ' 0.1 ' ' 1 Residual standard error: 10.48 on 42 degrees of freedom Multiple R-squared: 0.6551, Adjusted R-squared: 0.6304 F-statistic: 26.59 on 3 and 42 DF, p-value: 8.457e-10