R kode som illustrerer t test og konfidensintervall for to utvalg (avsnitt 10.2)

 

 

# Vi vil illustrere teorien i avsnitt 10.2 i l�reboka ved � se p� data om "peak expiratory flow" (PEF) for et utvalg av 10 menn og 11 kvinner.

# PEF m�les med et peak flow meter og angir den kraften en kan bl�se luften ut av lungene med (m�leenheten er liter per minutt).

 

# Vi leser inn pef m�lingene:

pef.mann=c(603, 628, 610, 639, 579, 637, 653, 612, 623, 550)

pef.kvinne=c(455,456,523,458, 460, 435, 489, 399, 525, 464, 450)

 

# Vi bruker kommandoen t.test til � teste om kvinner og menn har samme forventede PEF verdi og til � lage et konfidensintervall for differansen i forventede PEF verdier.

 

# Vi bruker f�rst de klassiske metodene som forutsetter at variansen er den samme for kvinner og menn (jf. sidene 504-505 i l�reboka):

t.test(pef.mann,pef.kvinne,var.equal=T)

 

# Vi bruker s� de modifiserte metodene som ikke krever at variansen er den samme for kvinner og menn (jf. sidene 499-503 i l�reboka):

t.test(pef.mann,pef.kvinne,var.equal=F)