1 / 30

Help! Statistiek!

Doorlopende serie laagdrempelige lezingen, voor iedereen vrij toegankelijk. Help! Statistiek!. Doel: Informeren over statistiek in klinisch onderzoek. Tijd: Derde woensdag in de maand, 12-13 uur 21 januari : Poisson regressie 18 februari : Graven naar causaliteit

gudrun
Download Presentation

Help! Statistiek!

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript


  1. Doorlopende serie laagdrempelige lezingen, voor iedereen vrij toegankelijk. Help! Statistiek! Doel: Informeren over statistiek in klinisch onderzoek. Tijd: Derde woensdag in de maand, 12-13 uur 21 januari : Poisson regressie 18 februari : Graven naar causaliteit 18 maart : Betrouwbaarheidsintervallen Sprekers: Vaclav Fidler, Hans Burgerhof, Wendy Post, Sacha la Bastide DG Epidemiologie

  2. Overzicht • Welke soort onderzoeksvragen • Introductie van voorbeeld • Waarom geen gewone lineaire regressie? • Wat is het Poisson regressiemodel • De Poisson verdeling • Specificatie van het regressiemodel • Interpretatie van parameters • Schatten en toetsen • Model fit • Hoe in SPSS • Referenties

  3. Onderzoeksvragen Mogelijke vragen: Zijn er in Nederland economische determinanten die het aantal kinderen voorspellen? Welke variabelen voorspellen het aantal blessures in een bepaald sport seizoen? Welke variabelen bepalen het aantal nieuwe tumorgevallen in een bepaald gebied in een jaar? Algemeen: Men is geïnteresseerd in relatie tussen Y: aantal events (afhankelijke variabele); niet negatieve gehele getallen!!!! X1 tm Xk : k continue en/of categoriale variabelen

  4. Onderzoeksvragen Regressie modellen Altijd relatie tussen afhankelijke variabele Y en onafhankelijke variabelen Lineaire regressie: Y is continu Logistische regressie: Y is dichotoom Poisson regressie: Y is een aantal niet negatief gehele getallen Schat gemiddelde van Y als functie van predictoren

  5. data Gegevens (gemanipuleerde data van Michel Brink): 50 topsporters in leeftijd 15-19 jaar Afhankelijke variabele: aantal blessures in een seizoen Achtergrond variabelen: Geslacht: 23 vrouwen (“0”), 27 mannen (“1”) Vetpercentage: gemiddeld 8; sd = 2; range 4-13 Onderzoeksvraag: Wordt het aantal blessures bepaald door geslacht en vetpercentage?

  6. data

  7. Poisson verdeling Bij niet vaak voorkomende gebeurtenissen is de Poisson verdeling geschikt vanwege 1. geen negatieve getallen 2. Positief verwachte waarde (gemiddelde) e -µ µy Pr(Y= y) = -------------, voor µ > 0 y! Gemiddelde (verwachting) = µ Speciale eigenschap: gemiddeld aantal, µ = variantie van aantal= var(Y)

  8. Poisson verdeling Simulatie van 100 waarnemingen Uit Poisson verdeling met verschillende verwachtingen mu

  9. Poisson verdeling Poisson verdeling: Gerelateerd aan binomiale verdeling bin(n,p) Voor n heel groot, en p heel klein: Verwachting = np, variantie np(1-p)  np Wanneer voldoet de Poissonverdeling? De kans op het optreden van tenminste één gebeurtenis in een zeker tijdsinterval is proportioneel aan de lengte van dat tijdsinterval De kans op meer dan 1 gebeurtenis in een heel klein tijdsinterval is verwaarloosbaar klein De aantallen gebeurtenissen in verschillende tijdsintervallen zijn onderling onafhankelijk

  10. Poisson regressie model Onze data over blessures: steekproef 50 personen Het aantal blessures voor persoon i, Yi heeft een Poisson verdeling met parameter µi Notatie: Yi ~ Poisson(µi) , Let op µi = verwachting (gemiddelde) = var(Yi) > 0

  11. Poisson regressie model Regressiemodel: Schat gemiddelde als functie van predictoren Lineaire model: µi = β0 + β1*geslachti + β2*vetpercentagei Probleem: negatieve schattingen zijn mogelijk Mogelijke oplossing: Schat log µi als functie van predictoren: Dus, Log(µi)= β0 + β1*geslachti + β2*vetpercentagei

  12. Vergelijking logistisch model Modelleren van gemiddelde in logistisch model: Een proportie  dat ligt tussen 0 en 1 Vandaar: log( /1- ) = β0 + β1*x1 + β2*x2 Zowel logistische regressie als Poisson regressie behoren tot de klasse van generalized linear models Gebruik van verschillende linkfuncties Logistisch model: log( /1- ) Poisson model: log(µi)

  13. Poisson regressie model Poisson regressiemodel: Yi ~ Poisson(µi), waarbij µi= exp(β0 + β1*geslachti + β2*vetpercentagei) Dus, Yi ~ Poisson(exp(β0 + β1*geslachti + β2*vetpercentagei)) Interpretatie??

  14. Poisson regressie model: voorbeeld Onze data van 50 sporters met de vraag : Zijn geslacht en vetpercentage van invloed op aantal blessures? Resultaten van Poisson analyse in R (maximum likelihood)

  15. Poisson regressie model: voorbeeld Resultaten van Poisson analyse in R: µ = exp(β0 + β1*geslacht + β2*vetpercentage) µdak = exp(-17 + 0.42*man + 0.11*vet) = exp(-0.17) * exp(0.42*man) * exp(0.11*vet) = 0.84 * 1.52man * 1.12vet Vrouw: 0.84 * 1.12vet Man: 0.84 * 1.52* 1.12vet

  16. Poisson regressie modelscatterplot + schattingen

  17. Poisson regressie model In de niet gemanipuleerde data: Elke sporter heeft een verschillende follow-up periode Oplossing voor dat probleem: ti = lengte follow-up voor persoon i: Schat log (µi/ti) als functie van predictoren. log(µi/ti)= β0 + β1*geslachti + β2*vetpercentagei log(µi) – log(ti)= β0 + β1*geslachti + β2*vetpercentagei log(µi) = log(ti) + β0 + β1*geslachti + β2*vetpercentagei log(ti) = offset gemiddelde proportioneel met follow-up Verdubbeling van follow-up betekent verdubbeling van aantal (mits andere predictoren zelfde blijven)

  18. Poisson regressie modeltoetsen Wald test voor elke predictor De interactieterm was niet significant: p-value 0.08

  19. Poisson regressie modeltoetsen Likelihood ratio test: gebaseerd op likelihood Vergelijken van geneste modellen: Verschil ~ 2 verdeeld

  20. Poisson regressie modelgoodness of fit Goodness of fit: Voorspelde waarden (ydak) vergelijken met geobserveerde waarnemingen (y): (yi – ydaki) ei = residu(gestandaardiseerd)i = --------------- , (ydaki) Onder Poisson model: Gemiddeld 0 en variantie 1:

  21. Poisson regressie modelmodelfit

  22. Overdispersion Bij Poisson verdeling: Gemiddelde = variantie! In veel gevallen: variantie > gemiddelde (bijvoorbeeld bij veel nullen) Dit heet overdispersion: Kan je checken en toetsen gestandaardiseerd residu is basis voor toetsing Bij overdispersion: in R: met optie quasipoisson werken

  23. R-syntax model1 = glm(formula = n_blessures ~ geslacht + vet, family =poisson) summary(model1) model2 = glm(formula = n_blessures ~ geslacht + vet, family = quasipoisson) summary(model3) In onze data: geen overdispersion.

  24. Poisson regressie modelin SPSS

  25. Poisson regressie modelSPSS

  26. Poisson regressie modelSPSS

  27. Poisson regressie modelSPSS

  28. Poisson regressie modelSPSS

  29. Referenties • Matthews D.E. & Farewell V.T. Using and Understanding Medical Statistics (hoofdstuk 12) • Gelman A. & Hill J. • Data analysis Using Regression and multilevel/hierarchical Models (hoofdstuk 6) • Mc Cullagh P. & Nelder J.A. • Generalized linear models (hoofdstuk 6)

  30. Volgende keer 18 februari: Graven naar causaliteit Zaal: 16

More Related