Получение среднего балла для каждой группы из линейного регрессионного выхода

Question

Я управляю линейной регрессией, прогнозирующей удовлетворенность жизнью по полу, расе и ее взаимодействию.

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

Вот выход:

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 

Я обязан предоставить средний балл удовлетворенности жизнью для (1) каждой половой группы, а также (2) каждой расовой группы (4 в итоге).

Итак, как я могу это сделать, используя R? Я знаю, что могу просто агрегировать данные, но есть намек на то, что я могу использовать некоторые коэффициенты для определения среднего уровня удовлетворенности как для пола, так и для групп расы.

спасибо, что заранее.

Answer 1

Один быстрый способ сделать это:

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