> pairwise <- pairs(means) > pairwise contrast estimate SE df t.ratio p.value A - B 0.0747 0.190 3317 0.392 0.9795 A - C -0.0476 0.190 3317 -0.250 0.9945 A - D 0.0215 0.191 3317 0.113 0.9995 B - C -0.1222 0.190 3317 -0.643 0.9182 B - D -0.0532 0.191 3317 -0.279 0.9924 C - D 0.0691 0.190 3317 0.363 0.9836 P value adjustment: tukey method for comparing a family of 4 estimates please create these results in a Latex table

Pairwise t-test

$`A - B`
	Welch Two Sample t-test

data:  mean_age by group
t = 0.3924, df = 3316.7, p-value = 0.97952
alternative hypothesis: true difference in means is not equal to 0
95 percent confidence interval:
 -0.3199497  0.4645497
sample estimates:
mean in group A mean in group B 
         53.065          52.991 


$`A - C`
	Welch Two Sample t-test

data:  mean_age by group
t = -0.25048, df = 3316.7, p-value = 0.9945
alternative hypothesis: true difference in means is not equal to 0
95 percent confidence interval:
 -0.3165497  0.2445497
sample estimates:
mean in group A mean in group C 
         53.065          53.123 


$`A - D`
	Welch Two Sample t-test

data