Föreläsningar (i efterhand)

Kursöversikt
Statistikens roll
Dos-respons
Design-issues
Farmakokinetik
Matematisk modellering med statistiska tillämpningar
– en projektkurs
Anders Källén
Matematikcentrum LTH
8 april 2015
Anders Källén
Matematisk modellering med statistiska tillämpningar – en projektkurs
Kursöversikt
Statistikens roll
Dos-respons
Design-issues
Farmakokinetik
Kursöversikt
(A) Om vad statistiska resultat betyder och inte betyder
1. Hur används statistik inom medicinsk vetenskap? – med
speciellt betoning på läkemedelsindustring
2. Vad betyder resultaten egentligen? Vad mäter ett p-värde?
Statistisk signifikans vs. medicinsk signifikans
3. Punktskattning, konfidens och verklighet – hur hänger de ihop?
4. Multiplicitetsproblemet – varför man ofta behöver förspecificera
vad man ska göra
5. Publikationsbias – varför även icke-konklusiva resultat borde
publiceras
6. Hur stora ska studier göras?
(B) Om svårigheterna med att välja dos för ett läkemedel
1. Hur kan dose-respons modelleras? Hur ska man välja dos om
det inte finns någon slumpvariation?
2. Vad är praktiskt och rimligt att göra?
Anders Källén
Matematisk modellering med statistiska tillämpningar – en projektkurs
Kursöversikt
Statistikens roll
Dos-respons
Design-issues
Farmakokinetik
Projekt nr 1 – ett mini-projekt
1. Handlar både om att förstå vad statistiska resultat betyder och hur de
statistiskt oskolade kan tänkas misstolka dem
2. Ni ska få ett uppdrag att tolka ett studieresultat för en dos-respons-studie
för att avgöra hur mycket information den ger.
3. Dessutom ska ni genomskåda hur två andra instanser har tolkat resultaten
4. Slutligen ska ni identifiera vilken mer information som vore önskvärd. I det
ingår att föreslå en ytterligare studie för att få kompletterande information.
5. Ni ska arbeta två och två och lämna in en gemensam rapport som ska
vara förklarande för en i statistik oskolad person.
Anders Källén
Matematisk modellering med statistiska tillämpningar – en projektkurs
Kursöversikt
Statistikens roll
Dos-respons
Design-issues
Farmakokinetik
(C) Även när man räknat rätt kan resultatet vara fel
1. Varför ska studier vara randomiserade (och blinda)? –
behovet av en “concurrent control”
2. Etiska aspekter
3. Föreläsning(ar) med många exempel
(D) Lite farmakokinetik – behövs för det stora projektet
1. Absorption, fördelning och elimination av ett läkemedel –
farmakokinetikernas sätt att beskriva detta
2. Några tips för MATLAB-användare
Anders Källén
Matematisk modellering med statistiska tillämpningar – en projektkurs
Kursöversikt
Statistikens roll
Dos-respons
Design-issues
Farmakokinetik
Projekt nr 2 – huvudprojektet
1. Ni ska arbeta i grupper om fyra. Projektet innehåller både teoretiska och
praktiska (MATLAB) utmaningar.
2. Det finns ingen entydig lösning på projektet!
3. Slutlig rapport inlämnad till kurschef (per epost) senast 13/5.
Presentation av den ska ske i större grupp den 26/5.
4. Lägesrapport till kurschef (per epost) senast 20/4.
5. Projektet handlar om att analysera plasmakoncentrationsdata på ett
(fiktivt) läkemedel i syfte att identifiera hur det bör administreras. Och hur
man ska hantera det när biverkningar uppträder.
6. Rapporten ska skrivas som en vetenskaplig rapport, och för insikter om
vad det innebär (struktur, vetenskaplig ärlighet etc) finns en obligatorisk
gästföreläsning den 22/4 på det temat!
Anders Källén
Matematisk modellering med statistiska tillämpningar – en projektkurs
Kursöversikt
Statistikens roll
Dos-respons
Design-issues
Farmakokinetik
Inledning
I
Många beslut fattas på basen av data som har genomgått en statistisk
analys
I
Beslutsfattare ofta oskolade i statistikens metoder och filosofi – vill ha
enkla svar
I
Det är därför viktigt att den som gjort analyserna “kan prata med bönder
på bönders språk, och de lärde på latin”.
I
Detta innefattar en förståelse för de vanligaste misstolkningarna av
statistik
I
Vilket är vad denna första föreläsning handlar om..
Vi ska använda medicinsk statistik, vilket bl.a. innefattar läkemedelsindustrins
användning, som exempel eftersom den på ett ganska renodlat sätt kan både
definiera problemen och illustrera dem.
Anders Källén
Matematisk modellering med statistiska tillämpningar – en projektkurs
Kursöversikt
Statistikens roll
Dos-respons
Design-issues
Farmakokinetik
Från åderlåtning till EBM
I
Varför har läkekonsten kunna använda ineffektiva metoder under många
hundra år utan att genomskådas?
Anders Källén
Matematisk modellering med statistiska tillämpningar – en projektkurs
Kursöversikt
Statistikens roll
Dos-respons
Design-issues
Farmakokinetik
Från åderlåtning till EBM
I
Varför har läkekonsten kunna använda ineffektiva metoder under många
hundra år utan att genomskådas?
I
Första svaret: placebo-effekten – ibland sker spontant tillfrisknande
I
Omvänt, de flesta (bra) terapierna fallerar ibland
I
Med andra ord: det är stor slumpvariation i effekten av en viss behandling
I
Och inte ens läkarna hade tillräcklig erfarenhet för att se signal genom
bruset
Anders Källén
Matematisk modellering med statistiska tillämpningar – en projektkurs
Kursöversikt
Statistikens roll
Dos-respons
Design-issues
Farmakokinetik
Från åderlåtning till EBM
I
Varför har läkekonsten kunna använda ineffektiva metoder under många
hundra år utan att genomskådas?
I
Första svaret: placebo-effekten – ibland sker spontant tillfrisknande
I
Omvänt, de flesta (bra) terapierna fallerar ibland
I
Med andra ord: det är stor slumpvariation i effekten av en viss behandling
I
Och inte ens läkarna hade tillräcklig erfarenhet för att se signal genom
bruset
I
Leta liten signal i mycket brus ⇒ stora datamaterial + statistik
I
Dessutom måste data samlas in på ett korrekt sätt – annars blir det fel
I
EBM = Evidence Based Medicine: endast behandlingar som har bevisat
sitt värde i riktiga kliniska studier får användas
Anders Källén
Matematisk modellering med statistiska tillämpningar – en projektkurs
Kursöversikt
Statistikens roll
Dos-respons
Design-issues
Farmakokinetik
Vad innebär EBM?
I
Först måste vi bestämma oss för vad vi vill visa!
I
Sedan måste vi skriva en detaljerad studieplan som talar om både hur
studien ska genomföras i detalj, och hur vi tänker använda de data som
den genererar för att dra de slutsatser vi specificerat.
I
Nyckelordet här är förspecificera – fundera på varför?
Anders Källén
Matematisk modellering med statistiska tillämpningar – en projektkurs
Kursöversikt
Statistikens roll
Dos-respons
Design-issues
Farmakokinetik
Vad innebär EBM?
I
Först måste vi bestämma oss för vad vi vill visa!
I
Sedan måste vi skriva en detaljerad studieplan som talar om både hur
studien ska genomföras i detalj, och hur vi tänker använda de data som
den genererar för att dra de slutsatser vi specificerat.
I
Nyckelordet här är förspecificera – fundera på varför?
De statistiska analyserna beräknar sedan ett tvåsidigt p-värde till
huvudhypotesen med enkel logik:
I
1. om p < 5% så finns en effekt
2. om p > 5% har vi inte lyckats påvisa en effekt
I
Att testet ska vara tvåsidigt är bestämt av FDA och betyder i praktiken
att signifikansnivån är 2.5%
Anders Källén
Matematisk modellering med statistiska tillämpningar – en projektkurs
Kursöversikt
Statistikens roll
Dos-respons
Design-issues
Farmakokinetik
Vad innebär EBM?
I
Först måste vi bestämma oss för vad vi vill visa!
I
Sedan måste vi skriva en detaljerad studieplan som talar om både hur
studien ska genomföras i detalj, och hur vi tänker använda de data som
den genererar för att dra de slutsatser vi specificerat.
I
Nyckelordet här är förspecificera – fundera på varför?
De statistiska analyserna beräknar sedan ett tvåsidigt p-värde till
huvudhypotesen med enkel logik:
I
1. om p < 5% så finns en effekt
2. om p > 5% har vi inte lyckats påvisa en effekt
I
Att testet ska vara tvåsidigt är bestämt av FDA och betyder i praktiken
att signifikansnivån är 2.5%
I
Vad var nu p-värden för något?
Anders Källén
Matematisk modellering med statistiska tillämpningar – en projektkurs
Kursöversikt
Statistikens roll
Dos-respons
Design-issues
Farmakokinetik
Ett statistiskt tests anatomi
I
Huvudhypotesen definierar en parameter som man vill skatta: ofta en
behandlingseffekt.
Anders Källén
Matematisk modellering med statistiska tillämpningar – en projektkurs
Kursöversikt
Statistikens roll
Dos-respons
Design-issues
Farmakokinetik
Ett statistiskt tests anatomi
I
Huvudhypotesen definierar en parameter som man vill skatta: ofta en
behandlingseffekt.
I
Denna parameter anses ha ett visst värde θ i den population man studerar
(t.ex. genomsnittlig behandlingseffekt)
Anders Källén
Matematisk modellering med statistiska tillämpningar – en projektkurs
Kursöversikt
Statistikens roll
Dos-respons
Design-issues
Farmakokinetik
Ett statistiskt tests anatomi
I
Huvudhypotesen definierar en parameter som man vill skatta: ofta en
behandlingseffekt.
I
Denna parameter anses ha ett visst värde θ i den population man studerar
(t.ex. genomsnittlig behandlingseffekt)
I
Huvudhypotesen innebär att vi specificerar ett värde θ0 (t.ex. att
behandingseffekten = 0).
Anders Källén
Matematisk modellering med statistiska tillämpningar – en projektkurs
Kursöversikt
Statistikens roll
Dos-respons
Design-issues
Farmakokinetik
Ett statistiskt tests anatomi
I
Huvudhypotesen definierar en parameter som man vill skatta: ofta en
behandlingseffekt.
I
Denna parameter anses ha ett visst värde θ i den population man studerar
(t.ex. genomsnittlig behandlingseffekt)
I
Huvudhypotesen innebär att vi specificerar ett värde θ0 (t.ex. att
behandingseffekten = 0).
Ur insamlade data x1 , . . . , xn kan man sedan beräkna en teststorhet som
har en kumulativ fördelningsfunktion F (x, θ) som beror på värdet av θ.
I
1. x → F (x, θ), fixt θ, är sannolikhetsteori
2. θ → F (x, θ), fixt x (observerat värdet) är statistik.
Anders Källén
Matematisk modellering med statistiska tillämpningar – en projektkurs
Kursöversikt
Statistikens roll
Dos-respons
Design-issues
Farmakokinetik
Ett statistiskt tests anatomi
I
Huvudhypotesen definierar en parameter som man vill skatta: ofta en
behandlingseffekt.
I
Denna parameter anses ha ett visst värde θ i den population man studerar
(t.ex. genomsnittlig behandlingseffekt)
I
Huvudhypotesen innebär att vi specificerar ett värde θ0 (t.ex. att
behandingseffekten = 0).
Ur insamlade data x1 , . . . , xn kan man sedan beräkna en teststorhet som
har en kumulativ fördelningsfunktion F (x, θ) som beror på värdet av θ.
I
1. x → F (x, θ), fixt θ, är sannolikhetsteori
2. θ → F (x, θ), fixt x (observerat värdet) är statistik.
I
Ett ensidigt p-värde (för mothypotesen θ > θ0 ) beräknas som
p = P(T ≥ x|θ0 ) = 1 − F (x, θ0 )
I
Motsvarande 2-sidiga p-värde fås genom att dubbla detta
Anders Källén
Matematisk modellering med statistiska tillämpningar – en projektkurs
Kursöversikt
Statistikens roll
Dos-respons
Design-issues
Farmakokinetik
Ett statistiskt tests anatomi
I
Huvudhypotesen definierar en parameter som man vill skatta: ofta en
behandlingseffekt.
I
Denna parameter anses ha ett visst värde θ i den population man studerar
(t.ex. genomsnittlig behandlingseffekt)
I
Huvudhypotesen innebär att vi specificerar ett värde θ0 (t.ex. att
behandingseffekten = 0).
Ur insamlade data x1 , . . . , xn kan man sedan beräkna en teststorhet som
har en kumulativ fördelningsfunktion F (x, θ) som beror på värdet av θ.
I
1. x → F (x, θ), fixt θ, är sannolikhetsteori
2. θ → F (x, θ), fixt x (observerat värdet) är statistik.
I
Ett ensidigt p-värde (för mothypotesen θ > θ0 ) beräknas som
p = P(T ≥ x|θ0 ) = 1 − F (x, θ0 )
I
Motsvarande 2-sidiga p-värde fås genom att dubbla detta
I
Innebörden: det är så (o)sannolikt att vi råkar påvisa en behandlingseffekt
även om det inte finns någon på riktigt.
Anders Källén
Matematisk modellering med statistiska tillämpningar – en projektkurs
Kursöversikt
Statistikens roll
Dos-respons
Design-issues
Farmakokinetik
I
I
Istället för att beräkna p-värden kan man fråga sig vad man lärt sig om
parametern θ.
Detta definierar konfidensfunktionen
C (θ) = 1 − F (x, θ)
Från vilken vi kan beräkna konfidensintervall
1
0.9
Confidence function C(θ)
I
(ensidig).
0.8
0.7
0.6
0.5
0.4
0.3
0.2
0.1
0
θ1 θ0
θ2
θ∗
Population parameter θ
Anders Källén
Matematisk modellering med statistiska tillämpningar – en projektkurs
Kursöversikt
Statistikens roll
Dos-respons
Design-issues
Farmakokinetik
Illustration - ett vanligt t-test
1. Vi har två stickprov x1 , . . . , xm ∈ N(m, σ) och y1 , . . . , yn ∈ N(m + θ, σ).
2. Vi beräknar medelvärden och stickprovsvarians vars fördelningar vi känner
3. Vår stickprovsfunktion blir ȳ − x̄ (s 2 är en “nuisance parameter”).
4. Statistiken lär oss att
ȳ − x̄ − θ
q
∈ t(n + m − 2) ≈ N(0, 1).
1
s m
+ n1
5. Konfidensfunktionen ges (approximativt) av
ȳ − x̄ − θ
θ − (ȳ − x̄)
C (θ) = 1 − Φ( q
) = Φ( q
)
1
1
1
s m+n
s m
+ n1
6. Egentligen ska här ingå en icke-central t-fördelning..
Anders Källén
Matematisk modellering med statistiska tillämpningar – en projektkurs
Kursöversikt
Statistikens roll
Dos-respons
Design-issues
Farmakokinetik
Några vanliga missförstånd
1. Notera att funktionen F (x, θ) beror av stickprovsstorleken n
Anders Källén
Matematisk modellering med statistiska tillämpningar – en projektkurs
Kursöversikt
Statistikens roll
Dos-respons
Design-issues
Farmakokinetik
Några vanliga missförstånd
1. Notera att funktionen F (x, θ) beror av stickprovsstorleken n
2. Ett p-värde kan vara litet både därför att effekten är stor och därför att
stickprovet är stort..
Anders Källén
Matematisk modellering med statistiska tillämpningar – en projektkurs
Kursöversikt
Statistikens roll
Dos-respons
Design-issues
Farmakokinetik
Några vanliga missförstånd
1. Notera att funktionen F (x, θ) beror av stickprovsstorleken n
2. Ett p-värde kan vara litet både därför att effekten är stor och därför att
stickprovet är stort..
3. så ett litet p-värde garanterar inte en stor effekt
Anders Källén
Matematisk modellering med statistiska tillämpningar – en projektkurs
Kursöversikt
Statistikens roll
Dos-respons
Design-issues
Farmakokinetik
Några vanliga missförstånd
1. Notera att funktionen F (x, θ) beror av stickprovsstorleken n
2. Ett p-värde kan vara litet både därför att effekten är stor och därför att
stickprovet är stort..
3. så ett litet p-värde garanterar inte en stor effekt
4. p-värden mäter ö.h.t. inte hur stor en effekt är utan hur säkra vi är på att
det finns någon effekt över huvud taget
5. Det du vet om effektstorleken beskrivs av konfidensintervallet.
Anders Källén
Matematisk modellering med statistiska tillämpningar – en projektkurs
Kursöversikt
Statistikens roll
Dos-respons
Design-issues
Farmakokinetik
Några vanliga missförstånd
1. Notera att funktionen F (x, θ) beror av stickprovsstorleken n
2. Ett p-värde kan vara litet både därför att effekten är stor och därför att
stickprovet är stort..
3. så ett litet p-värde garanterar inte en stor effekt
4. p-värden mäter ö.h.t. inte hur stor en effekt är utan hur säkra vi är på att
det finns någon effekt över huvud taget
5. Det du vet om effektstorleken beskrivs av konfidensintervallet.
6. Omvänt, om p-värdet inte är signifikant, betyder det inte att läkemedlet
inte har effekt bara att du misslyckats med att visa att det har effekt.
7. Du kan alltså inte säga p > 5% ⇒ läkemedlet fungerar inte!
I
Varför? Vilken är den alternativa förklaringen?
Anders Källén
Matematisk modellering med statistiska tillämpningar – en projektkurs
Kursöversikt
Statistikens roll
Dos-respons
Design-issues
Farmakokinetik
Varför måste vi förspecificera?
Betrakta följande problem.
En kvinna har våldtagits och man har hittat främmande DNA hos
henne.
1. Vittnen pekar ut en viss man, och när han DNA-undersöks hittas
en match.
2. Det hittade DNAt jämförs med en DNA-databas omfattande
20000 män. En matchning hittas.
Antag att om man slumpmässigt jämför två DNA profiler så får man
en match en gång på 10000.
Hur stort är DNAts bevisvärde i de två situationerna? När är det rätt att dra
mannen inför rätta?
Anders Källén
Matematisk modellering med statistiska tillämpningar – en projektkurs
Kursöversikt
Statistikens roll
Dos-respons
Design-issues
Farmakokinetik
Varför måste vi förspecificera?
Betrakta följande problem.
En kvinna har våldtagits och man har hittat främmande DNA hos
henne.
1. Vittnen pekar ut en viss man, och när han DNA-undersöks hittas
en match.
2. Det hittade DNAt jämförs med en DNA-databas omfattande
20000 män. En matchning hittas.
Antag att om man slumpmässigt jämför två DNA profiler så får man
en match en gång på 10000.
Hur stort är DNAts bevisvärde i de två situationerna? När är det rätt att dra
mannen inför rätta?
Slutsats: p-värdets historia är fundamental för dess tolkning!
Anders Källén
Matematisk modellering med statistiska tillämpningar – en projektkurs
Kursöversikt
Statistikens roll
Dos-respons
Design-issues
Farmakokinetik
Dags att larma?
På en viss arbetsplats har man noterat att ovanligt många av deras
anställda har fått barn med allvarliga missbildningar. Man samlade
därför in data för de föregående 5 åren och fann att under den
perioden hade 50 kvinnliga medarbetare blivit gravida och fött barn
och av dem hade så många som 5 fått allvarliga missbildningar. Så
allvarliga missbildningar förekommer normalt i 2% av alla barn.
För en Bin(50, 0.02)-fördelning gäller att
Utfall
Sannolikhet
0
0.364
1
0.372
Anders Källén
2
0.186
3
0.0607
4
0.015
≥5
0.00321
Matematisk modellering med statistiska tillämpningar – en projektkurs
Kursöversikt
Statistikens roll
Dos-respons
Design-issues
Farmakokinetik
Det tragiska fallet med Sally Clark
I
1999 dömdes Sally Clark för mord på sina två barn (11 resp 8 veckor).
Försvaret hävdade att det var två fall av SIDS.
I
Domen baserades huvudsakligen på att expertvittnet Roy Meadow, en
barnläkare, hävdade att sannolikheten att två barn i samma familj skulle
dö av SIDS är ungefär 1 : 73 · 106 , varför det måste vara mord
I
Är detta ett relevant sätt att använda p-värden?
Anders Källén
Matematisk modellering med statistiska tillämpningar – en projektkurs
Kursöversikt
Statistikens roll
Dos-respons
Design-issues
Farmakokinetik
Det tragiska fallet med Sally Clark
I
1999 dömdes Sally Clark för mord på sina två barn (11 resp 8 veckor).
Försvaret hävdade att det var två fall av SIDS.
I
Domen baserades huvudsakligen på att expertvittnet Roy Meadow, en
barnläkare, hävdade att sannolikheten att två barn i samma familj skulle
dö av SIDS är ungefär 1 : 73 · 106 , varför det måste vara mord
I
Är detta ett relevant sätt att använda p-värden?
I
Vad hände sedan? En högre domstol upphävde domen, men Sally hämtade
sig aldrig utan dog av alkoholism i sitt hem 8 år senare, 42 år gammal.
Anders Källén
Matematisk modellering med statistiska tillämpningar – en projektkurs
Kursöversikt
Statistikens roll
Dos-respons
Design-issues
Farmakokinetik
Publikationsbias
Om tidskrifter endast tar in resultat som är statistiskt signifikanta, vilken blir
effekten?
105
Precision (No of patients)
ISIS-4
ISIS-2
GISSI-1
104
103
102
101 −1
10
100
Odds Ratio
101
Figuren illustrerar skillnaden mellan vad man sett i flera små studier och vad
man såg en stor studie.
Anders Källén
Matematisk modellering med statistiska tillämpningar – en projektkurs
Kursöversikt
Statistikens roll
Dos-respons
Design-issues
Farmakokinetik
Typfel i statistiken och annars
När man testar om det finns en signal kan olika saker hända:
Signal
Test
Ja
Nej
Finns
Korrekt
Miss (Typ II-fel)
Saknas
Falskt alarm (Typ I-fel)
Korrekt
Denna tabell är fundamental i en mängd olika situationer, t.ex.
I
Vid medicinsk diagnostik
I
När man ska fälla en misstänkt brottsling
I
När man ska tolka statistik
Anders Källén
Matematisk modellering med statistiska tillämpningar – en projektkurs
Kursöversikt
Statistikens roll
Dos-respons
Design-issues
Farmakokinetik
Typfel i statistiken och annars
När man testar om det finns en signal kan olika saker hända:
Signal
Test
Ja
Nej
Finns
Korrekt
Miss (Typ II-fel)
Saknas
Falskt alarm (Typ I-fel)
Korrekt
Denna tabell är fundamental i en mängd olika situationer, t.ex.
I
Vid medicinsk diagnostik
I
När man ska fälla en misstänkt brottsling
I
När man ska tolka statistik
Den romerska lagens grundpelare: du är oskyldig tills bevisad skyldig!
vilket betyder att ett läkemedel är verkningslöst tills bevisat effektivt.
Anders Källén
Matematisk modellering med statistiska tillämpningar – en projektkurs
Kursöversikt
Statistikens roll
Dos-respons
Design-issues
Farmakokinetik
Typ I och Typ II fel i statistik
1. Antag att vi har en teststorhet X som har en fördelningsfunktion F (x, θ)
(stora värden på X är bra)
Anders Källén
Matematisk modellering med statistiska tillämpningar – en projektkurs
Kursöversikt
Statistikens roll
Dos-respons
Design-issues
Farmakokinetik
Typ I och Typ II fel i statistik
1. Antag att vi har en teststorhet X som har en fördelningsfunktion F (x, θ)
(stora värden på X är bra)
2. Om nollhypotesen är θ = 0 och observerat värde på X är x ∗ så ges Typ
I-felet av p = 1 − F (x ∗ , 0).
Anders Källén
Matematisk modellering med statistiska tillämpningar – en projektkurs
Kursöversikt
Statistikens roll
Dos-respons
Design-issues
Farmakokinetik
Typ I och Typ II fel i statistik
1. Antag att vi har en teststorhet X som har en fördelningsfunktion F (x, θ)
(stora värden på X är bra)
2. Om nollhypotesen är θ = 0 och observerat värde på X är x ∗ så ges Typ
I-felet av p = 1 − F (x ∗ , 0).
3. Fixera α och låt xc lösa ekvationen F (xc , 0) = 1 − α.
Anders Källén
Matematisk modellering med statistiska tillämpningar – en projektkurs
Kursöversikt
Statistikens roll
Dos-respons
Design-issues
Farmakokinetik
Typ I och Typ II fel i statistik
1. Antag att vi har en teststorhet X som har en fördelningsfunktion F (x, θ)
(stora värden på X är bra)
2. Om nollhypotesen är θ = 0 och observerat värde på X är x ∗ så ges Typ
I-felet av p = 1 − F (x ∗ , 0).
3. Fixera α och låt xc lösa ekvationen F (xc , 0) = 1 − α.
4. Då är Typ I-felet < α precis då x ∗ > xc .
Anders Källén
Matematisk modellering med statistiska tillämpningar – en projektkurs
Kursöversikt
Statistikens roll
Dos-respons
Design-issues
Farmakokinetik
Typ I och Typ II fel i statistik
1. Antag att vi har en teststorhet X som har en fördelningsfunktion F (x, θ)
(stora värden på X är bra)
2. Om nollhypotesen är θ = 0 och observerat värde på X är x ∗ så ges Typ
I-felet av p = 1 − F (x ∗ , 0).
3. Fixera α och låt xc lösa ekvationen F (xc , 0) = 1 − α.
4. Då är Typ I-felet < α precis då x ∗ > xc .
5. Type II-felet borde använda F (x, θs ) där θs är det sanna värdet.
Anders Källén
Matematisk modellering med statistiska tillämpningar – en projektkurs
Kursöversikt
Statistikens roll
Dos-respons
Design-issues
Farmakokinetik
Typ I och Typ II fel i statistik
1. Antag att vi har en teststorhet X som har en fördelningsfunktion F (x, θ)
(stora värden på X är bra)
2. Om nollhypotesen är θ = 0 och observerat värde på X är x ∗ så ges Typ
I-felet av p = 1 − F (x ∗ , 0).
3. Fixera α och låt xc lösa ekvationen F (xc , 0) = 1 − α.
4. Då är Typ I-felet < α precis då x ∗ > xc .
5. Type II-felet borde använda F (x, θs ) där θs är det sanna värdet.
6. Men θs är okänt, så man betraktar styrkefunktionen
β(θ) = 1 − F (xc , θ)
(β(0) = α) . . .
7. ... vilket blir sannolikheten att vi förkastar hypotesen θ = 0 när det sanna
värdet är θ.
Anders Källén
Matematisk modellering med statistiska tillämpningar – en projektkurs
Kursöversikt
Statistikens roll
Dos-respons
Design-issues
Farmakokinetik
1. Antag normalfördelning, så att fördelningsfunktionen är
F (x − θ) = Φ( x−θ
)
σ
2. Då gäller att F (0) = 0.5 och att β(θ) = 1 − F (xc − θ) = F (θ − xc ).
3. Om vi sätter θc = xc så gäller alltså att
β(θc ) = 0.5,
β(θ) = F (θ − θc ).
β(θ) = Probability for success
1
0.9
0.8
0.7
0.6
0.5
0.4
0.3
0.2
0.1
0
θc
θ0
Population parameter θ
Anders Källén
2θc
Matematisk modellering med statistiska tillämpningar – en projektkurs
Kursöversikt
Statistikens roll
Dos-respons
Design-issues
Farmakokinetik
Hur stor ska vi göra en studie?
1. Målet är att göra en så liten studie som möjligt som har realistiska
möjligheter att påvisa en effekt
Anders Källén
Matematisk modellering med statistiska tillämpningar – en projektkurs
Kursöversikt
Statistikens roll
Dos-respons
Design-issues
Farmakokinetik
Hur stor ska vi göra en studie?
1. Målet är att göra en så liten studie som möjligt som har realistiska
möjligheter att påvisa en effekt
2. Det kräver att vi specificerar vilka Typ I- och Typ II-fel vi kan
acceptera/vill ha
Anders Källén
Matematisk modellering med statistiska tillämpningar – en projektkurs
Kursöversikt
Statistikens roll
Dos-respons
Design-issues
Farmakokinetik
Hur stor ska vi göra en studie?
1. Målet är att göra en så liten studie som möjligt som har realistiska
möjligheter att påvisa en effekt
2. Det kräver att vi specificerar vilka Typ I- och Typ II-fel vi kan
acceptera/vill ha
3. Konvention för Typ I: signifikansnivå 5% tvåsidigt ⇒ α = 0.025. Allmänt
accepterat!
Anders Källén
Matematisk modellering med statistiska tillämpningar – en projektkurs
Kursöversikt
Statistikens roll
Dos-respons
Design-issues
Farmakokinetik
Hur stor ska vi göra en studie?
1. Målet är att göra en så liten studie som möjligt som har realistiska
möjligheter att påvisa en effekt
2. Det kräver att vi specificerar vilka Typ I- och Typ II-fel vi kan
acceptera/vill ha
3. Konvention för Typ I: signifikansnivå 5% tvåsidigt ⇒ α = 0.025. Allmänt
accepterat!
4. För Typ II används begreppet styrka (power). Beror på studien men
nästan alltid ≥ 80% power (β ≤ 0.2).
Anders Källén
Matematisk modellering med statistiska tillämpningar – en projektkurs
Kursöversikt
Statistikens roll
Dos-respons
Design-issues
Farmakokinetik
Hur stor ska vi göra en studie?
1. Målet är att göra en så liten studie som möjligt som har realistiska
möjligheter att påvisa en effekt
2. Det kräver att vi specificerar vilka Typ I- och Typ II-fel vi kan
acceptera/vill ha
3. Konvention för Typ I: signifikansnivå 5% tvåsidigt ⇒ α = 0.025. Allmänt
accepterat!
4. För Typ II används begreppet styrka (power). Beror på studien men
nästan alltid ≥ 80% power (β ≤ 0.2).
5. Dessutom måste man bestämma sig för vilken effekt vi inte vill missa och
6. hitta bra skattningar på “nuisance parameters” (typiskt standardavvikelse)
Anders Källén
Matematisk modellering med statistiska tillämpningar – en projektkurs
Kursöversikt
Statistikens roll
Dos-respons
Design-issues
Farmakokinetik
Hur stor ska vi göra en studie?
1. Målet är att göra en så liten studie som möjligt som har realistiska
möjligheter att påvisa en effekt
2. Det kräver att vi specificerar vilka Typ I- och Typ II-fel vi kan
acceptera/vill ha
3. Konvention för Typ I: signifikansnivå 5% tvåsidigt ⇒ α = 0.025. Allmänt
accepterat!
4. För Typ II används begreppet styrka (power). Beror på studien men
nästan alltid ≥ 80% power (β ≤ 0.2).
5. Dessutom måste man bestämma sig för vilken effekt vi inte vill missa och
6. hitta bra skattningar på “nuisance parameters” (typiskt standardavvikelse)
7. Notera att effekt vi inte vill missa är inte ett statistiskt problem att
bestämma.
Anders Källén
Matematisk modellering med statistiska tillämpningar – en projektkurs
Kursöversikt
Statistikens roll
Dos-respons
Design-issues
Farmakokinetik
Översättning till matematik
1. Man approximerar oftast med normalfördelning med känt σ som beror av hur
många patienter vi har i studien:
σ
θ̂ ∈ N(θ, √ ).
n
2. För test θ ≤ 0 får vi styrkefunktionen
√
n(θ − θc )
σ
och vi bestämmer n genom att lösa en ekvation på formen β(θ0 ) = 0.8 där θ0 är
skillnaden vi inte vill missa.
3. Från detta får man en allmän formel, som i fallet av tvåsidigt test på 5% nivå
och styrka 80% blir
7.85σ 2
n=
.
θ02
β(θ) = Φ
4. Med två grupper med storlekar rn, (1 − r )n har vi
σ0
σ= p
r (1 − r )
där σ0 är standardavvikelsen i en grupp. Speciellt ser vi att styrkan är som mest
om vi fördelar patienterna i lika stora grupper.
Anders Källén
Matematisk modellering med statistiska tillämpningar – en projektkurs
Kursöversikt
Statistikens roll
Dos-respons
Design-issues
Farmakokinetik
Hur bestämmer man vilken dos ett nytt läkemedel ska ges i?
1. Det är en sak att avgöra om en viss dos av ett läkemedel är en effektiv
behandling . . .
2. . . . men en annan att avgöra vilken dos som bör användas på marknaden
3. En given dos kan vara för hög (för mycket biverkningar) eller för låg (bara
liten effekt)
4. Under läkemedelsutvecklingprocessen måste man därför göra (minst) en
dos-respons-studie
5. Man kan endast studera ett relativt litet antal doser, och måste välja dessa
så att man kan få ut maximal information
6. Men hur ska man välja dos från den?
7. Analys? Statistiker jämför ofta doserna parvis – är det vettigt?
8. Detta handlar mini-projektet om, och här kommer lite bakgrund
Anders Källén
Matematisk modellering med statistiska tillämpningar – en projektkurs
Kursöversikt
Statistikens roll
Dos-respons
Design-issues
Farmakokinetik
Hur ser en normal dos-respons kurva ut?
Följande antaganden torde vara okontroversiella:
1. Effekten mäts positiv och placebo har noll effekt (skilj på effekt och
mätvärde!) E (0) = 0
2. Högre dos bör ge högre effekt (om vi kan bortse från biverkningar) (E (D)
växande)
3. Det finns en maximal effekt som kan uppnås med höga doser
(limD→∞ E (D) = Emax )
Anders Källén
Matematisk modellering med statistiska tillämpningar – en projektkurs
Kursöversikt
Statistikens roll
Dos-respons
Design-issues
Farmakokinetik
Hur ser en normal dos-respons kurva ut?
Följande antaganden torde vara okontroversiella:
1. Effekten mäts positiv och placebo har noll effekt (skilj på effekt och
mätvärde!) E (0) = 0
2. Högre dos bör ge högre effekt (om vi kan bortse från biverkningar) (E (D)
växande)
3. Det finns en maximal effekt som kan uppnås med höga doser
(limD→∞ E (D) = Emax )
Den kliniska traditionen ger att om en tablett ger otillräcklig effect, ta två.
Matematiskt (empirisk lag): inom ett visst dosintervall gäller att
E (D) ≈ α + β log D.
Detta gäller i någon omgivning av ED50 , som är den dos som ger 50% av
maximal effekt.
Anders Källén
Matematisk modellering med statistiska tillämpningar – en projektkurs
Kursöversikt
Statistikens roll
Dos-respons
Design-issues
Farmakokinetik
1
0.9
Normalized response
0.8
0.7
0.6
0.5
0.4
0.3
0.2
0.1
0 0
10
101
102
Dose (on log-scale)
103
En sådan funktion kan skrivas
E (D) = Emax Ψ(a + b ln D)
där Ψ(x) är en CDF med Ψ(0) = 1/2 (bekvämlighetsval). Vanliga val:
Ψ(x) = Φ(x),
Anders Källén
Ψ(x) =
1
.
1 + e −x
Matematisk modellering med statistiska tillämpningar – en projektkurs
Kursöversikt
Statistikens roll
Dos-respons
Design-issues
Farmakokinetik
Ett exempel från toxikologi
För en viss substans har vi följande data (varje grupp innehöll 10 råttor):
dose:
deaths:
5
1
25
3
100
8
800
10
Vi vill skatta LD50 (den dos som dödar 50% av råttorna). Låt p(D) vara sannolikheten
att en råtta dör när den får dosen D.
1. Antag först att p(D) = Ψ(a + b ln D), där Ψ(x) är logistisk CDF. Ger svart kurva
och skattningen 37.6
2. För att ta hänsyn till bakgrundsmortalitet inför vi en parameter till:
p(D) = π + (1 − π)Ψ(a + b ln D).
Ger grå kurva och skattningen 53.6.
1
0.9
Percent mice dead
0.8
0.7
0.6
0.5
0.4
0.3
0.2
Background mortality
0.1
0 0
10
101
Anders Källén
102
Dose of drug
103
Matematisk modellering med statistiska tillämpningar – en projektkurs
Kursöversikt
Statistikens roll
Dos-respons
Design-issues
Farmakokinetik
Vad måste man ta hänsyn till när man bestämmer dos?
1. Höga doser är alltid en risk, så dos bör ligga på log-linjär del av
dos-responskurvan
2. Det finns alltid flera effekter att ta ställning till: positiva effekter och
biverkningar
3. Att hitta en kompromiss kan vara svårt och ofta blir farmaceutiska
faktorer och ekonomi viktiga spelare
4. Under alla omständigheter återstår det att sedan i ett Fas III-program visa
att de(n) fastställd(a) dosen är tillräckligt effektiv och säker.
Nu är det dags för mini-projektet!
Anders Källén
Matematisk modellering med statistiska tillämpningar – en projektkurs
Kursöversikt
Statistikens roll
Dos-respons
Design-issues
Farmakokinetik
Vilka slutsatser får man dra från en statistisk analys?
I
Även när den statistiska analysen är korrekt genomförd är det inte klart
vad man får dra för slutsatser
I
För det spelar designen av studien en viktig roll. (Ursäkta svengelskan)
I
En analys talar om att två grupper (nästan säkert) skiljer sig åt. Det man
vill uttala sig om är att det är något som orsakar denna skillnad.
I
När får man lov att göra det?
I
Nedan följer ett antal exempel som delvis är historiska milstolpar mot det
vi idag kallar EBM.
Anders Källén
Matematisk modellering med statistiska tillämpningar – en projektkurs
Kursöversikt
Statistikens roll
Dos-respons
Design-issues
Farmakokinetik
Att göra en gallup
Varför är Gallup-undesökningar ett begrepp och t.o.m. ett verb?
I
1936 års presidentval i USA stod mellan Alfred Landon (R) och Franklin
D. Roosevelt (D)
I
Literary Digest tillfrågade 2 miljoner röstare per post. Resultat: Landon får
en jordskredsseger (370 mot 161)
Anders Källén
Matematisk modellering med statistiska tillämpningar – en projektkurs
Kursöversikt
Statistikens roll
Dos-respons
Design-issues
Farmakokinetik
Att göra en gallup
Varför är Gallup-undesökningar ett begrepp och t.o.m. ett verb?
I
1936 års presidentval i USA stod mellan Alfred Landon (R) och Franklin
D. Roosevelt (D)
I
Literary Digest tillfrågade 2 miljoner röstare per post. Resultat: Landon får
en jordskredsseger (370 mot 161)
I
Det faktiska utfallet var 523 mot 8 för Roosevelt!
I
Varför så fel?
Anders Källén
Matematisk modellering med statistiska tillämpningar – en projektkurs
Kursöversikt
Statistikens roll
Dos-respons
Design-issues
Farmakokinetik
Att göra en gallup
Varför är Gallup-undesökningar ett begrepp och t.o.m. ett verb?
I
1936 års presidentval i USA stod mellan Alfred Landon (R) och Franklin
D. Roosevelt (D)
I
Literary Digest tillfrågade 2 miljoner röstare per post. Resultat: Landon får
en jordskredsseger (370 mot 161)
I
Det faktiska utfallet var 523 mot 8 för Roosevelt!
I
Varför så fel? Selektionsbias: urval baserat på läsare, bilägare och folk med
telefon.
I
George Gallup’s organisation undersökte endast 50000 och lyckades få en
korrekt prediktion!
I
Man kan inte kompensera bias med stora antal!!
Anders Källén
Matematisk modellering med statistiska tillämpningar – en projektkurs
Kursöversikt
Statistikens roll
Dos-respons
Design-issues
Farmakokinetik
Dr Snow och kolerans orsaker
I
John Snow räknas ibland som epidemiologins fader (vattenpumpen på
Broad Street)
I
Frågan var vid den tiden: hur spreds kolera? Dålig luft??
I
Snow jämförde antal fall i hushåll som fick sitt vatten från två olika ställen
på Themsen. Det ena nära ett avlopp, det andra rent.
Avloppsvatten Rent vatten
Kolera fall
4 093
461
Totalantal
266 516
173 749
Relativ risk: 5.8 med 95% CI (5.3, 6.4)
I
Vad bevisar detta?
Anders Källén
Matematisk modellering med statistiska tillämpningar – en projektkurs
Kursöversikt
Statistikens roll
Dos-respons
Design-issues
Farmakokinetik
Dr Snow och kolerans orsaker
I
John Snow räknas ibland som epidemiologins fader (vattenpumpen på
Broad Street)
I
Frågan var vid den tiden: hur spreds kolera? Dålig luft??
I
Snow jämförde antal fall i hushåll som fick sitt vatten från två olika ställen
på Themsen. Det ena nära ett avlopp, det andra rent.
Avloppsvatten Rent vatten
Kolera fall
4 093
461
Totalantal
266 516
173 749
Relativ risk: 5.8 med 95% CI (5.3, 6.4)
I
Vad bevisar detta?
I
Inte randomiserad, men alla sociala klasser och regioner fick vatten från
båda bolagen. Kvar finns därför endast förklaringen förorenat vatten.
Anders Källén
Matematisk modellering med statistiska tillämpningar – en projektkurs
Kursöversikt
Statistikens roll
Dos-respons
Design-issues
Farmakokinetik
Leder rökning till lungcancer?
I
Folk började dö i lungcancer i början av 1900-talet. Varför?
I
Richard Doll och Bradford Hill bestämde sig för att undersöka relationen
mellan rökning och lungcancer
I
De fick ihop följande tabell (1950)
Fall
Kontroller
Rökare
647
622
icke-rökare
2
27
Odds ratio 14 med 95% CI (4.2,87.3)
I
Studien är en fall-kontroll-studie på sjukhus i 4 sjukhus, fall=lungcancer,
kontroll=ej lungcancer
I
Kan vi dra slutsatsen att rökning orsakar lungcancer? (Fisher accepterade
aldrig sambandet)
Anders Källén
Matematisk modellering med statistiska tillämpningar – en projektkurs
Kursöversikt
Statistikens roll
Dos-respons
Design-issues
Farmakokinetik
Simpsons paradox
Två metoder att operera njursten jämfördes. Vi rak jämförelse lyckades
behandling A (öppen kirurgi) i 78% av fallen, medan den mindre invasiva
metoden B lyckades i 83% av fallen. Betyder det att B är en bättre metod?
Anders Källén
Matematisk modellering med statistiska tillämpningar – en projektkurs
Kursöversikt
Statistikens roll
Dos-respons
Design-issues
Farmakokinetik
Simpsons paradox
Två metoder att operera njursten jämfördes. Vi rak jämförelse lyckades
behandling A (öppen kirurgi) i 78% av fallen, medan den mindre invasiva
metoden B lyckades i 83% av fallen. Betyder det att B är en bättre metod?
Om man delar upp resultatet efter storlek på
stenarna får man följande tabell
Success
Failure
Odds
small stones
A
B
81
234
6
36
13.5
6.5
large stones
A
B
192
55
71
25
2.7
2.2
Total
A
B
273
289
77
61
3.5
4.7
Här ser vi att A är bättre både på små och
stora stenar! Hur är detta möjligt?
Anders Källén
Matematisk modellering med statistiska tillämpningar – en projektkurs
Kursöversikt
Statistikens roll
Dos-respons
Design-issues
Farmakokinetik
Simpsons paradox
Två metoder att operera njursten jämfördes. Vi rak jämförelse lyckades
behandling A (öppen kirurgi) i 78% av fallen, medan den mindre invasiva
metoden B lyckades i 83% av fallen. Betyder det att B är en bättre metod?
Om man delar upp resultatet efter storlek på
stenarna får man följande tabell
large stones
A
B
192
55
71
25
2.7
2.2
Total
A
B
273
289
77
61
3.5
4.7
Number of successes
Success
Failure
Odds
small stones
A
B
81
234
6
36
13.5
6.5
B1 + B2
Anders Källén
A2
200
100
A1
B2
0
Här ser vi att A är bättre både på små och
stora stenar! Hur är detta möjligt?
A1 + A2
B1
0
10
20
30
40
50
Number of failures
60
70
Matematisk modellering med statistiska tillämpningar – en projektkurs
Kursöversikt
Statistikens roll
Dos-respons
Design-issues
Farmakokinetik
Simpsons paradox
Två metoder att operera njursten jämfördes. Vi rak jämförelse lyckades
behandling A (öppen kirurgi) i 78% av fallen, medan den mindre invasiva
metoden B lyckades i 83% av fallen. Betyder det att B är en bättre metod?
Om man delar upp resultatet efter storlek på
stenarna får man följande tabell
large stones
A
B
192
55
71
25
2.7
2.2
Total
A
B
273
289
77
61
3.5
4.7
Number of successes
Success
Failure
Odds
small stones
A
B
81
234
6
36
13.5
6.5
B1 + B2
A2
200
100
A1
B2
0
Här ser vi att A är bättre både på små och
stora stenar! Hur är detta möjligt?
A1 + A2
B1
0
10
20
30
40
50
Number of failures
60
70
Vilken analys ska vi tro på?
Anders Källén
Matematisk modellering med statistiska tillämpningar – en projektkurs
Kursöversikt
Statistikens roll
Dos-respons
Design-issues
Farmakokinetik
James Lind och skörbjugg
I
Den skottske läkaren James Lind betraktas gärna som den kliniska
studiens fader
I
Han var Surgeon’s mate på HMS Salisbury när han i maj 1747 ville
undersöka vad som kunde påverka skörbjugg
I
Skörbjugg var en sjukdom som var ett stort hinder för längre resor över de
stora haven
I
Lind to 12 sjömän med likartade stadier av skörbjugg och testade 6 olika
behandlingar, till varje behandling hade han en kontroll
I
De som fick apelsin och citron klarade sig överlägset bäst!
Anders Källén
Matematisk modellering med statistiska tillämpningar – en projektkurs
Kursöversikt
Statistikens roll
Dos-respons
Design-issues
Farmakokinetik
James Lind och skörbjugg
I
Den skottske läkaren James Lind betraktas gärna som den kliniska
studiens fader
I
Han var Surgeon’s mate på HMS Salisbury när han i maj 1747 ville
undersöka vad som kunde påverka skörbjugg
I
Skörbjugg var en sjukdom som var ett stort hinder för längre resor över de
stora haven
I
Lind to 12 sjömän med likartade stadier av skörbjugg och testade 6 olika
behandlingar, till varje behandling hade han en kontroll
I
De som fick apelsin och citron klarade sig överlägset bäst!
I
Idag vet vi att skörbjugg = brist på C-vitamin
I
Konstigt nog dröjde till 1795 innan den brittiska flottan delade ut
citron-juice till sjömännen!
Anders Källén
Matematisk modellering med statistiska tillämpningar – en projektkurs
Kursöversikt
Statistikens roll
Dos-respons
Design-issues
Farmakokinetik
Den första randomiserade studien
I
1946 ville man avgöra om streptomycin fungerade på lung-TB (eftersom
det dödade TB-bakterien in vitro)
I
En grupp behandlades med sängläge i 6 månader(!) enbart, den andra
med sängläge plus streptomycin.
I
Mortalitetsresultaten:
Streptomycin
Endast sängläge
Döda
4
14
Överlevande
51
38
Relativ risk 0.27 med 95% CI (0.08,0.70) (p = 0.006).
Anders Källén
Matematisk modellering med statistiska tillämpningar – en projektkurs
Kursöversikt
Statistikens roll
Dos-respons
Design-issues
Farmakokinetik
Den första randomiserade studien
I
1946 ville man avgöra om streptomycin fungerade på lung-TB (eftersom
det dödade TB-bakterien in vitro)
I
En grupp behandlades med sängläge i 6 månader(!) enbart, den andra
med sängläge plus streptomycin.
I
Mortalitetsresultaten:
Streptomycin
Endast sängläge
Döda
4
14
Överlevande
51
38
Relativ risk 0.27 med 95% CI (0.08,0.70) (p = 0.006).
I
Det banbrytande med studien var emellertid det sätt man valde vem som
fick vilken behandling (nästa sida)
I
Men beslutet var inte lätt - patienter ska alltid får den bästa behandlingen.
Det som räddade studien var att det var ont om streptomycin.
Anders Källén
Matematisk modellering med statistiska tillämpningar – en projektkurs
Kursöversikt
Statistikens roll
Dos-respons
Design-issues
Farmakokinetik
Hur man allokerade patienter i streptomycin-studien
“Determination whether a patient would be treated by streptomycin and
bed-rest (S case) or by bed-rest alone (C case) was made by reference to a
statistical series based on random sampling numbers drawn up for each sex
at each center by Professor Bradford Hill: the details of the series were
unknown to any of the investigators or to the coordinator and were
contained in a set of sealed envelopes, each bearing on the outside only the
name of the hospital and a number. After acceptance of a patient by the
panel, and before admission to the streptomycin centre, the appropriate
numbered envelope was opened at the central office; the card inside told if
the patient was to be an S or a C case, and this information was then given
to the medical officer of the centre. Patients were not told before admission
that they were to get special treatment”
Anders Källén
Matematisk modellering med statistiska tillämpningar – en projektkurs
Kursöversikt
Statistikens roll
Dos-respons
Design-issues
Farmakokinetik
Hormonsersättning och hjärtsjukdom
Hormonersättning (HRT) ges till (vissa) kvinnor runt menopausen i syfte att
reducera obehag.
I
Under 1980-talet kom ett antal små observationsstudier som antydde en
positiv effekt av HRT på hjärtsjukdom. En meta-analys av 31 studier
antydde en 50% reduktion
I
Under 1991-92 konstruerades en dubbel-blind, placebo-kontrollerad klinisk
studie i syfte att få konklusiva resultat
I
Total 16 608 kvinnor randomiserades jämnt mellan två grupper och
studien planerades att pågå i 8.5 år
Anders Källén
Matematisk modellering med statistiska tillämpningar – en projektkurs
Kursöversikt
Statistikens roll
Dos-respons
Design-issues
Farmakokinetik
Hormonsersättning och hjärtsjukdom
Hormonersättning (HRT) ges till (vissa) kvinnor runt menopausen i syfte att
reducera obehag.
I
Under 1980-talet kom ett antal små observationsstudier som antydde en
positiv effekt av HRT på hjärtsjukdom. En meta-analys av 31 studier
antydde en 50% reduktion
I
Under 1991-92 konstruerades en dubbel-blind, placebo-kontrollerad klinisk
studie i syfte att få konklusiva resultat
I
Total 16 608 kvinnor randomiserades jämnt mellan två grupper och
studien planerades att pågå i 8.5 år
I
Avbröts efter 5.2 år p.g.a. en interims-analys p.g.a indikationer på en ökad
risk, skattad till 29%
I
Varför?
Anders Källén
Matematisk modellering med statistiska tillämpningar – en projektkurs
Kursöversikt
Statistikens roll
Dos-respons
Design-issues
Farmakokinetik
Hormonsersättning och hjärtsjukdom
Hormonersättning (HRT) ges till (vissa) kvinnor runt menopausen i syfte att
reducera obehag.
I
Under 1980-talet kom ett antal små observationsstudier som antydde en
positiv effekt av HRT på hjärtsjukdom. En meta-analys av 31 studier
antydde en 50% reduktion
I
Under 1991-92 konstruerades en dubbel-blind, placebo-kontrollerad klinisk
studie i syfte att få konklusiva resultat
I
Total 16 608 kvinnor randomiserades jämnt mellan två grupper och
studien planerades att pågå i 8.5 år
I
Avbröts efter 5.2 år p.g.a. en interims-analys p.g.a indikationer på en ökad
risk, skattad till 29%
I
Varför?
I
Kvinnor som bad om HRT är från högre socialgrupp, vilka ofta har en
bättre vanor med diet och motion. Det var dessas effekt som slog igenom i
observationsstudierna.
Anders Källén
Matematisk modellering med statistiska tillämpningar – en projektkurs
Kursöversikt
Statistikens roll
Dos-respons
Design-issues
Farmakokinetik
Hjälper det mot spädbarnsdödlighet att modern röker?
Följande tabell visar en relation mellan dödfödda barn och
moderns rökvanor
Infant
Stillborn
Liveborn
P(Dead) (%)
Mother
Non-smoker
smoker
40
6
10665
1277
0.374
0.468
Så rökande mödrar har fler dödfödda barn. Många data pekar
på detta.
Anders Källén
Matematisk modellering med statistiska tillämpningar – en projektkurs
Kursöversikt
Statistikens roll
Dos-respons
Design-issues
Farmakokinetik
Hjälper det mot spädbarnsdödlighet att modern röker?
Följande tabell visar en relation mellan dödfödda barn och
moderns rökvanor
Infant
Stillborn
Liveborn
P(Dead) (%)
Mother
Non-smoker
smoker
40
6
10665
1277
0.374
0.468
Så rökande mödrar har fler dödfödda barn. Många data pekar
på detta.
Men om vi tittar på barn som väger < 2.5kg vid födseln har vi
följande tabell
Infant
Stillborn
Liveborn
P(Dead) (%)
Mother
Non-smoker
smoker
19
1
371
73
4.87
1.35
Betyder det att små barn har ökad överlevnadschans om
modern röker?
Anders Källén
Matematisk modellering med statistiska tillämpningar – en projektkurs
Kursöversikt
Statistikens roll
Dos-respons
Design-issues
Farmakokinetik
Hjälper det mot spädbarnsdödlighet att modern röker?
Följande tabell visar en relation mellan dödfödda barn och
moderns rökvanor
Infant
Stillborn
Liveborn
P(Dead) (%)
Mother
Non-smoker
smoker
40
6
10665
1277
0.374
0.468
Så rökande mödrar har fler dödfödda barn. Många data pekar
på detta.
Men om vi tittar på barn som väger < 2.5kg vid födseln har vi
följande tabell
Infant
Stillborn
Liveborn
P(Dead) (%)
Mother
Non-smoker
smoker
19
1
371
73
4.87
1.35
Will Rogers:
When the Okies
left Oklahoma
and moved to
Californa, they
raised the
average
intelligence level
in both states.
Betyder det att små barn har ökad överlevnadschans om
modern röker?
Anders Källén
Matematisk modellering med statistiska tillämpningar – en projektkurs
Kursöversikt
Statistikens roll
Dos-respons
Design-issues
Farmakokinetik
Ger ultraljudsundersökningar hjärnskador?
I
I Sverige görs rutinmässigt ULjus kring vecka 17 under en graviditet för
att säkra god tillväxt och utesluta allvarliga missbildningar
I
Kan detta ge hjärnskador?
I
En grupp ville studera om det fanns en ökad risk för vänsterhänthet pga
ULjus.
I
Studien genomfördes 1985-87 vid 19 kliniker i Sverige.
I
4997 kvinnor randomiserads till två grupper av vilka den ena fick rutin
ULus vecka 15, den andra inte.
I
10 år senare fick mödrarna ett frågeformulär om barnens hänthet.
I
Resultat:
Screened Group
Non-screened Group
I
Right-handed
636
648
Non-Right handed
156
134
OR 1.19 (0.92,1.53). Slutsats?
Anders Källén
Matematisk modellering med statistiska tillämpningar – en projektkurs
Kursöversikt
Statistikens roll
Dos-respons
Design-issues
Farmakokinetik
Varför ge upp? Inspektera data noggrannt
I
Vid inspektion av data noterades att i 103 fall hade det faktiskt gjorts en
ULus före vecka 19 i den andra gruppen.
I
Vad göra?
Anders Källén
Matematisk modellering med statistiska tillämpningar – en projektkurs
Kursöversikt
Statistikens roll
Dos-respons
Design-issues
Farmakokinetik
Varför ge upp? Inspektera data noggrannt
I
Vid inspektion av data noterades att i 103 fall hade det faktiskt gjorts en
ULus före vecka 19 i den andra gruppen.
I
Vad göra?
I
Om dessa flyttas över till första gruppen fås följande tabell
Exposed Group
Non-exposed Group
Right-handed
655
629
Non-Right handed
167
123
I
OR 1.3 med 85% CI (1.01,1.69).
I
Betyder detta att vi hittat bevis för att tidig ULus ger hjärnskador?
Anders Källén
Matematisk modellering med statistiska tillämpningar – en projektkurs
Kursöversikt
Statistikens roll
Dos-respons
Design-issues
Farmakokinetik
Vad är farmakokinetik?
I
Farmakokinetik är läran om hur kroppen hanterar ett läkemedel
I
I
I
I
upptag
fördelning i kroppen
elimination
Försöker besvara frågor som
I
I
I
I
Vilken dos ska ges för att erhålla en viss effekt?
Hur länge varar effekten för denna dos?
Hur ofta måste läkemedlet administreras för att effekten ska bibehållas?
Vid vilken dosering uppträder oönskade biverkningar?
Anders Källén
Matematisk modellering med statistiska tillämpningar – en projektkurs
Kursöversikt
Statistikens roll
Dos-respons
Design-issues
Farmakokinetik
Farmakokinetikens grunder
I
Det man mäter är koncentrationen C (t) i blodet (eller plasma, serum)
I
Läkemedlet sprider sig (oftast) även till andra delar av kroppen
I
Den fundamentala massbalansekvationen för läkemedel kroppen är att
M 0 (t) = a(t) − CL(t)C (t)
där M(t) =mängd i kroppen, a(t) = absorptionshastighet och CL(t) =
clearance (ofta konstant, ibland beroende av C (t)).
I
I den empiriska farmakokinetiken används denna till att härled några
kinetiska parametrar som beskriver aspekter av upptag, fördelning (oftast
volymer) och elimination (clearance).
I
Ett (äldre) alternativ är att modellera själva fördelningsprocessen av
läkemedlet i kroppen. Det är det vi ska göra i projektet.
Anders Källén
Matematisk modellering med statistiska tillämpningar – en projektkurs
Kursöversikt
Statistikens roll
Dos-respons
Design-issues
Farmakokinetik
En enkel fördelningsmodell
I
Den enklaste modellen är att tänka på kroppen som en stor behållare i
snabb jämvikt med blodet. Då är M(t) = VC (t) och vi får en
differentialekvation i C (t).
I
Den näst enklaste är att tänka på kroppen som bestående av två behållare
som kommunicerar med varandra. Antag att deras volymer är fixa och att
de är väl blandade. Konstanta hastighetskonstanter.
I
Systemet är linjärt i den meningen att totaldosen ger plasma
koncentrationer som är summan av vad de individuella ger
I
Detta leder till ett system av differentialekvationer
I
Cc0 (t) =
a(t) − kpc Cp (t) − (kcp + kce )Cc (t)
(1)
Cp0 (t) =
kcp Cc (t) − (kpc + kpe )Cp (t)
(2)
Om a(t) = 0 och Cc (0) = C0 , Cp (0) = 0 (bolus dos) ger detta att
Cc (t) = Ae −λt + Be −µt .
Anders Källén
Matematisk modellering med statistiska tillämpningar – en projektkurs
Kursöversikt
Statistikens roll
Dos-respons
Design-issues
Farmakokinetik
I
Absorptionen antas ofta vara enkel: en tablett i magen sönderfaller med en
hastighet som är proportionell mot hur mycket som finns kvar. Leder till
att
a(t) = Fe −ka t .
I
Det lämnas som övning att tänka ut hur det förändrar ekvationen för Cc .
Det behövs för projektet, där den slutliga formen för Cc (t) ska jämföras
med data från 10 personer. Notera att konstanterna skiljer sig åt mellan
individerna (dock inte jättemycket).
I
En viktig del av projektet är att tänka ut vad som händer när man tar nya
doser – hur ändras plasmakoncentrationerna?
I
En annan är att tänka på hur man sammanfattar information från 10
individer till en slutsats.
I
Det finns många sätt att lösa projektets detaljer – så jag förväntar mig
flera olika rapporter
I
Anpassning till data kräver att man använder en optimerings-funktion i
MATLAB.
Anders Källén
Matematisk modellering med statistiska tillämpningar – en projektkurs
Kursöversikt
Statistikens roll
Dos-respons
Design-issues
Farmakokinetik
Anpassning till data i MATLAB
Antag att vi har följande uppsättning data:
t
x
0
0
1.40
1.49
2.79
0.399
4.19
-0.75
5.58
-0.42
6.98
0.32
8.38
0.32
9.77
-0.10
11.17
-0.21
12.57
0
Av någon anledning tror vi oss veta att det finns ett samband på formen
x = Ce −λ1 t sin(λ2 t)(= f (t, C , λ1 , λ2 )).
Men för vilka konstanter?
Enklast optimering är att använda minsta-kvadrat-methoden: minimera
X
Q(C , λ1 , λ2 ) =
(xi − f (ti , C , λ1 , λ2 ))2 .
i
Anders Källén
Matematisk modellering med statistiska tillämpningar – en projektkurs
Kursöversikt
Statistikens roll
Dos-respons
Design-issues
Farmakokinetik
Enkel kod i MATLAB
Vi definierar först en funktionsfil Qfunc.m som beräkna Q:
function mk_norm = Qfunc(constants,t_data,x_data)
t = t_data;
x = x_data;
C = constants(1);
l1=constants(2);
l2=constants(3);
xnew = C*exp(-l1*t).*sin(l2*t);
mk_norm = norm(xnew-x)
1.5
1
0.5
0
−0.5
Sedan blir själva anropet
−1
0
5
10
15
constants=[1 0.1 0.8]; % initial guess
constants = fminsearch(’Qfunc’,constants,[],t,x);
Anders Källén
Matematisk modellering med statistiska tillämpningar – en projektkurs