gss_ TukeyHSD(data . aov1) Tukey multiple comparisons of means 95% family-wise confidence level Fit: aov(formula = gss123$numpets ~ gss123$marital, data = gss123) $ gss123$marital diff Lwr upr p adj Widowed-Married -0. 57086337 -1. 4973297 0. 3556029 0. 4445954 Divorced-Married -0. 30627127 -1. 0140053 0.4014628 0. 7615689 Separated-Married -0. 09306504 -1.3612677 1. 1751376 0. 9996401 Never married-Married -0. 78391514 -1. 3680985 -0. 1997318 0.0023877 Divorced-Widowed 0. 26459210 -0. 7746596 1. 3038438 0.9574647 Separated-Widowed 0.47779833 -1. 0012198 1.9568164 0. 9033503 Never married-Widowed -0. 21305176 -1. 1724405 0. 7463370 0. 9740534 Separated-Divorced 0. 21320623 -1. 1395843 1.5659968 0. 9928421 Never married-Divorced -0.47764387 -1. 2279601 0. 2726724 0.4102036 Never married-Separated -0.69085010 -1.9832993 0.6015991 0.5885337 marital numpes.png 5,900 APR 30 esc C C Search oGo to file/ function B - Addins . Source Console Terminal Jobs ~ 1 7 > data. aov1 F) gss123$marital 130 32 . 47 3. 587 0. 00648 * * Residuals 1141 10327 9 .05 Signif. codes: 0 '***' 01001 '**' 0.01 '*' 0.05 . . ' 0.1 . ' 1 63668 observations deleted due to missingness > by(gss123$numpets, gss123$marital, mean, na. rm=T) gss123$marital : Married [1] 2. 049587 gss123$marital: Widowed [1] 1.478723 gss123$marital: Divorced [1] 1. 743316 gss123$marital: Separated [1] 1.956522 gss123$marital: Never married [1] 1.265672 > TukeyHSD(data. aov1) Tukey multiple comparisons of means 95% family-wise confidence level Fit: aov(formula = gss123$numpets ~ gss123$marital, data = gss123) 20220430_2304....jpg marital numpes.png 5,902 MAY 1 MacBook esc C G Search or type URL