4 T-test

4.1 One-sample T-test

Eksempel: Slagvolumen målt ved to metoder (parrede målinger, uge 1)

I slides uge 1 lavede vi et one-sample T-test på differenserne:

data ms; 
  set ms;
  dif=mf-sv;
run;
proc ttest data=ms;
    var dif;
run;    

Her får man bla gennemsnit incl CI (95% CL Mean), t-teststørrelse og p-værdi.

Et one-sampe t-test på differenserne er det samme som et parret test T-test (se at output er fuldstændigt identisk):

proc ttest data=ms;
    paired mf*sv;
run;

4.2 Two-sample t-test

Eksempel: Rygning og lungefunktion for Vietnam veteraner (rygning og lungegfunktion, uge 2)

Two-sample t-test udføres ved proc ttest. I class-linjen angiver man den variabel, som angiver de to grupper. Responsvariablen sættes i var-linjen.

proc ttest data=viet;
    class smoke01;
      var FEV1;
run;

Jeg har en gammel video, hvor jeg ved et t-test sammenligner BMI for mænd og kvinder i SundBy data (benyttet i opgave 1 uge 1: SundBy). Koden til at udføre t-testet er:

data sundby0;
  infile "http://staff.pubhealth.ku.dk/~sr/BasicStatistics/datasets/sundby.txt" URL firstobs=2;
  input kon v75 v76 v17 v24af;
  vaegt=v75;
  hoejde=v76/100;
  bmi = vaegt/hoejde**2;
  gender=kon;
  drop v75 v76 v17 v24af kon;
run;

proc ttest data=sundby0;
  var bmi;
  class gender;
run;

og output er forklaret i denne video (vises kun i html-versionen af denne e-bog)