d <- read.csv("cells.csv")
cell <- t.test(value ~ factor(group, levels = c("treated", "control")), data = d)
m <- aggregate(value ~ group + mouse, data = d, FUN = mean)
mouse <- t.test(value ~ factor(group, levels = c("treated", "control")), data = m)
cat(sprintf("cells as n:  n = %d + %d, t = %.2f, df = %.1f, p = %.3f\n", sum(d$group == "control"), sum(d$group == "treated"), cell$statistic, cell$parameter, cell$p.value))
cat(sprintf("mouse means: n = %d + %d, t = %.2f, df = %.2f, p = %.2f\n", sum(m$group == "control"), sum(m$group == "treated"), mouse$statistic, mouse$parameter, mouse$p.value))