R kode som viser beregningene i eksempel 10.12

 

 

# I eksempel 10.12 i l?reboka ser en p? det ber?mte vaksinefors?ket til Jonas Salk.

 

# Salk ville teste ut en vaksine mot poliomyelitt. I et dobbelt blindt fors?k fikk behandlingsgruppen den nye vaksinen, men kontrollgruppen fikk en inaktiv placebo injeksjon. Sp?rsm?let var da hvor store grupper en m?tte ha.

 

# La p₁ v?re sannsynligheten for at et tilfeldig valgt barn i kontrollgruppen f?r poliomyelitt, mens p₂ er den tilsvarende sannsynligheten for et barn i kontrollgruppen.

 

# Vi vil finne hvor store de to gruppene m? v?re for at sannsynligheten for feil av type II skal v?re 10% for en test med signifikansniv? 5% hvis p₁=0.00030 og p₂=0.00015. Ved ? bruke formel (10.7) i l?reboka, finner vi:

z.alfa=qnorm(0.95)

z.beta=qnorm(0.90)

p1=0.00030

q1=1-p1

p2=0.00015

q2=1-p2

d=p1-p2

n=(z.alfa*sqrt((p1+p2)*(q1+q2)/2)+z.beta*sqrt(p1*q1+p2*q2))^2/d^2

 

 

# R har en kommando som gir resultatet direkte:

power.prop.test(p1=0.00030,p2=0.00015,power=0.90,alternative="one.sided")