从线性回归输出获取每组的平均得分

Question

I run the linear regression predicting life satisfaction by sex, race and its interaction.

lm2 <-lm(nids$satisfaction~nids$male+nids$race+nids$male:nids$race)

Here is an output:

Call:
lm(formula = nids$satisfaction ~ nids$male + nids$race + nids$male:nids$race)

Residuals:
    Min      1Q  Median      3Q     Max 
-6.6613 -1.3366 -0.0485  1.7378  4.9515 

Coefficients:
                    Estimate Std. Error t value Pr(>|t|)    
(Intercept)          4.17751    0.05467  76.410  < 2e-16 ***
nids$male            0.39318    0.08564   4.591 4.45e-06 ***
nids$race            0.87095    0.03421  25.459  < 2e-16 ***
nids$male:nids$race -0.17947    0.05261  -3.411 0.000649 ***
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Residual standard error: 2.358 on 12016 degrees of freedom
Multiple R-squared:  0.07414,   Adjusted R-squared:  0.07391 
F-statistic: 320.7 on 3 and 12016 DF,  p-value: < 2.2e-16

I'm required to provide the mean score of life satisfaction for (1) each sex group as well as for (2) each race group (4 in total).

So, how can I do it using R? I know that I can just aggregate the data but there is a hint that I can use some coefficients to figure out the mean of satisfaction level for both sex and race groups.

Thank you very much in advance.

Answer 1

One quick way of doing it:

malenids <- nids[nids$male == 1, ]
femalenids <- nids[nids$male == 0, ]
lapply(split(malenids, malenids$race), function(x) mean(x$satisfaction))
lapply(split(femalenids, femalenids$race), function(x) mean(x$satisfaction))