Hier volgen enkele voorbeelden van responsvariabelen die discrete telresultaten vertegenwoordigen:<\/span><\/p>\n Als de variantie ongeveer gelijk is aan het gemiddelde, past een Poisson-regressiemodel over het algemeen goed bij een dataset.<\/span><\/p>\n Als de variantie echter aanzienlijk groter is dan het gemiddelde, kan een negatief binomiaal regressiemodel over het algemeen beter bij de gegevens passen.<\/span><\/p>\n Er zijn twee technieken die we kunnen gebruiken om te bepalen of Poisson-regressie of negatieve binomiale regressie geschikter is voor een bepaalde dataset:<\/span><\/p>\n 1. Restpercelen<\/strong><\/span><\/p>\n We kunnen een grafiek maken van de gestandaardiseerde residuen tegen de voorspelde waarden uit een regressiemodel.<\/span><\/p>\n Als de meerderheid van de gestandaardiseerde residuen tussen -2 en 2 ligt, is een Poisson-regressiemodel waarschijnlijk geschikt.<\/span><\/p>\n Als veel residuen echter buiten dit bereik vallen, zal een negatief binomiaal regressiemodel waarschijnlijk een betere fit opleveren.<\/span><\/p>\n 2. Test van de waarschijnlijkheidsratio<\/strong><\/span><\/p>\n We kunnen een Poisson-regressiemodel en een negatief binomiaal regressiemodel in dezelfde dataset passen en vervolgens een waarschijnlijkheidsratiotest uitvoeren.<\/span><\/p>\n Als de p-waarde van de test onder een bepaald significantieniveau ligt (bijvoorbeeld 0,05), kunnen we concluderen dat het negatief binomiale regressiemodel een significant betere fit oplevert.<\/span><\/p>\n Het volgende voorbeeld laat zien hoe u deze twee technieken in R kunt gebruiken om te bepalen of het beter is om een Poisson-regressie- of een negatief binomiaal regressiemodel te gebruiken voor een bepaalde gegevensset.<\/span><\/p>\n Stel dat we willen weten hoeveel beurzen een honkbalspeler op een middelbare school in een bepaalde provincie ontvangt op basis van zijn schoolafdeling („A“, „B“ of „C“) en zijn schoolcijfer. toelatingsexamen voor de universiteit (gemeten van 0 tot 100). ).<\/span><\/p>\n Gebruik de volgende stappen om te bepalen of een negatief binomiaal regressiemodel of een Poisson-regressiemodel beter bij de gegevens past.<\/span><\/p>\n Stap 1: Cre\u00eber de gegevens<\/strong><\/span><\/p>\n Met de volgende code wordt de dataset gemaakt waarmee we gaan werken, die gegevens over 1000 honkbalspelers bevat:<\/span><\/p>\n Stap 2: Pas een Poisson-regressiemodel en een negatief binomiaal regressiemodel aan<\/strong><\/span><\/p>\n De volgende code laat zien hoe u zowel een Poisson-regressiemodel als een negatief binomiaal regressiemodel aan de gegevens kunt aanpassen:<\/span><\/p>\n Stap 3: Maak resterende plots<\/strong><\/span><\/p>\n De volgende code laat zien hoe u restplots voor beide modellen kunt maken.<\/span> <\/p>\n Uit de grafieken kunnen we zien dat de residuen meer verspreid zijn voor het Poisson-regressiemodel (merk op dat sommige residuen verder reiken dan 3) vergeleken met het negatief binomiale regressiemodel.<\/span><\/p>\n Dit is een teken dat een negatief binomiaal regressiemodel waarschijnlijk geschikter is, aangezien de residuen van dit model kleiner zijn.<\/span><\/p>\n Stap 4: Voer een waarschijnlijkheidsratiotest uit<\/strong><\/span><\/p>\n Ten slotte kunnen we een waarschijnlijkheidsratiotest uitvoeren om te bepalen of er een statistisch significant verschil is in de pasvorm van de twee regressiemodellen:<\/span><\/p>\n De p-waarde van de test blijkt 3,508072e-29<\/strong> te zijn, wat aanzienlijk minder is dan 0,05.<\/span><\/p>\n We zouden dus kunnen concluderen dat het negatieve binomiale regressiemodel een significant betere fit met de gegevens oplevert vergeleken met het Poisson-regressiemodel.<\/span><\/p>\n Een inleiding tot de negatieve binominale verdeling<\/a> Negatieve binomiale regressie en Poisson-regressie zijn twee soorten regressiemodellen die moeten worden gebruikt wanneer deresponsvariabele wordt weergegeven door discrete telresultaten. Hier volgen enkele voorbeelden van responsvariabelen die discrete telresultaten vertegenwoordigen: Het aantal studenten dat een bepaalde opleiding afrondt Het aantal verkeersongevallen op een bepaald kruispunt Het aantal deelnemers dat een marathon voltooit Het aantal retourzendingen […]<\/p>\n","protected":false},"author":1,"featured_media":0,"comment_status":"open","ping_status":"open","sticky":false,"template":"","format":"standard","meta":{"footnotes":""},"categories":[11],"tags":[],"class_list":["post-1534","post","type-post","status-publish","format-standard","hentry","category-gids"],"yoast_head":"\n\n
Voorbeeld: negatieve binomiale regressie versus Poisson-regressie<\/strong><\/span><\/h3>\n
#make this example reproducible\n<\/span>set. seeds<\/span> (1)\n\n#create dataset\n<\/span>data <- data. frame<\/span> (offers = c(rep(0, 700), rep(1, 100), rep(2, 100),\n rep(3, 70), rep(4, 30)),\n division = sample(c(' A<\/span> ', ' B<\/span> ', ' C<\/span> '), 100, replace = TRUE<\/span> ),\n exam = c(runif(700, 60, 90), runif(100, 65, 95),\n runif(200, 75, 95)))\n\n#view first six rows of dataset\n<\/span>head(data)\n\n offers division exam\n1 0 A 66.22635\n2 0 C 66.85974\n3 0 A 77.87136\n4 0 B 77.24617\n5 0 A 62.31193\n6 0 C 61.06622<\/strong><\/pre>\n #fit Poisson regression model\n<\/span>p_model <- glm(offers ~ division + exam, family = ' fish<\/span> ', data = data)\n\n#fit negative binomial regression model\n<\/span>library<\/span> (MASS)\n\nnb_model <- glm. nb<\/span> (offers ~ division + exam, data = data)\n<\/strong><\/pre>\n #Residual plot for Poisson regression\n<\/span>p_res <- resid<\/span> (p_model)\nplot(fitted(p_model), p_res, col=' steelblue<\/span> ', pch=16,\n xlab=' Predicted Offers<\/span> ', ylab=' Standardized Residuals<\/span> ', main=' Poisson<\/span> ')\nabline(0,0)\n\n#Residual plot for negative binomial regression<\/span>\nnb_res <- resid<\/span> (nb_model)\nplot(fitted(nb_model), nb_res, col=' steelblue<\/span> ', pch=16,\n xlab=' Predicted Offers<\/span> ', ylab=' Standardized Residuals<\/span> ', main=' Negative Binomial<\/span> ')\nabline(0,0)<\/strong> <\/pre>\n![\"Negatieve](\"https:\/\/statorials.org\/wp-content\/uploads\/2023\/08\/negatifbinom1.png\")
<\/p>\n pchisq(2 * ( logLik<\/span> (nb_model) - logLik<\/span> (p_model)), df = 1, lower. tail<\/span> = FALSE<\/span> )\n\n'log Lik.' 3.508072e-29 (df=5)\n<\/strong><\/pre>\n Aanvullende bronnen<\/strong><\/span><\/h3>\n
Een inleiding tot de Poisson-verdeling<\/a><\/p>\n","protected":false},"excerpt":{"rendered":"