R kode som illustrerer korrelasjon (avsnitt 12.5)

 

 

# Vi vil illustrere den empiriske korrelasjonskoeffisienten ved ? se p? data der "peak expiratory flow" for 10 pasienter er m?lt med to ulike peak flow meter, som vi her for enkelhets skyld vil kalle apparat 1 og apparat 2. Vi er interessert i ? studere samvariasjonen i m?leresultatene med de to apparatene.

 

# Vi leser inn pef m?lingene for pasientene med de to apparatene:

app1=c(494, 395, 516, 434, 476, 557, 413, 442, 650, 433)

app2=c(512, 430, 520, 428, 500, 600, 364, 380, 658, 445)

 

# Merk at observasjonene h?rer sammen to og to:

cbind(app1,app2)

 

# Vi lager spredningsplott av dataene:

plot(app1,app2)

 

# Vi bestemmer den empiriske korrelasjonskoeffisienten:

cor(app1,app2)

 

# Merk at den empiriske korrelasjonen er symmetrisk, s? vi f?r samme resultat ved kommandoen:

cor(app2,app1)

 

# Test for nullhypotesen H₀: r = 0 mot alternativet at r ikke er lik null (se side 654 i l?reboka):

cor.test(app2,app1)

 

# Kommandoen gir ogs? et konfidensintervall for r (se sidene 655-657 i l?reboka), men det er ikke pensum.