Comment exécuter un test t sur un sous-ensemble de données

Question

D'abord, cet ensemble de données est-il bien rangé pour un test t?

https://i.stack.imgur.com/tMK6R.png

Deuxièmement, je suis en train de faire un test t échantillon pour comparer les moyens à temps 3 d'un traitement et b pour « résultat 1 ». Comment ferais-je cela?

données Exemple:

df <- structure(list(code = c(100, 100, 100, 101, 101, 101, 102, 102, 
     102, 103, 103, 103), treatment = c("a", "a", "a", "b", "b", "b", 
     "a", "a", "a", "b", "b", "b"), sex = c("f", "f", "f", "m", "m", 
     "m", "f", "f", "f", "f", "f", "f"), time = c(1, 2, 3, 1, 2, 3, 
     1, 2, 3, 1, 2, 3), `outcome 1` = c(21, 23, 33, 44, 45, 47, 22, 
     34, 22, 55, 45, 56), `outcome 2` = c(21, 32, 33, 33, 44, 45, 
     22, 57, 98, 65, 42, 42), `outcome 3` = c(62, 84, 63, 51, 45, 
     74, 85, 34, 96, 86, 45, 47)), .Names = c("code", "treatment", 
     "sex", "time", "outcome 1", "outcome 2", "outcome 3"), 
     class = c("tbl_df", "tbl", "data.frame"), row.names = c(NA, -12L)) 

Answer 1

D'abord, vous devrez définir les sous-ensembles que vous voulez tester, vous pouvez exécuter le test t. Vous n'avez pas nécessairement à stocker les sous-ensembles dans les variables comme je l'ai fait, mais cela rend la sortie du test t plus claire.

Normalement, les questions t-test, je vous recommande de l'aide fournie par t.test, mais étant donné que cela implique plus complexe subsetting, j'ai compris comment faire ici:

var_a <- df$`outcome 1`[df$treatment=="a" & df$time==3] 
var_b <- df$`outcome 1`[df$treatment=="b" & df$time==3] 

t.test(var_a,var_b) 

sortie :

Welch Two Sample t-test 

data: var_a and var_b 
t = -3.3773, df = 1.9245, p-value = 0.08182 
alternative hypothesis: true difference in means is not equal to 0 
95 percent confidence interval: 
-55.754265 7.754265 
sample estimates: 
mean of x mean of y 
    27.5  51.5