<\/p>\n
Als we willen testen of een enkele variabele normaal verdeeld is of niet, kunnen we een QQ-plot<\/a> maken om de verdeling te visualiseren of een formele statistische test uitvoeren, zoals een Anderson Darling-test<\/a> of een Jarque-Bera-test<\/a> .<\/span><\/p>\n Als we echter willen testen of meerdere<\/em> variabelen normaal verdeeld zijn als groep, moeten we een multivariate normaliteitstest<\/strong> uitvoeren.<\/span><\/p>\n In deze zelfstudie wordt uitgelegd hoe u de volgende multivariate normaliteitstests uitvoert voor een bepaalde gegevensset in R:<\/span><\/p>\n Gerelateerd:<\/strong> Als we uitschieters in een multivariate context willen identificeren, kunnen we Mahalanobis-afstand<\/a> gebruiken.<\/span><\/p>\n De Mardia-test<\/strong> bepaalt of een groep variabelen al dan niet een multivariate normale verdeling volgt. De nul- en alternatieve hypothesen van de test zijn als volgt:<\/span><\/p>\n H 0<\/sub> (nul): De variabelen volgen een multivariate normale verdeling.<\/span><\/p>\n H a<\/sub> (alternatief): De variabelen volgen geen<\/em> multivariate normale verdeling.<\/span><\/p>\n De volgende code laat zien hoe u deze test in R kunt uitvoeren met behulp van het QuantPsyc-<\/strong> pakket:<\/span><\/p>\n De functie mult.norm()<\/strong> test op multivariate normaliteit in zowel scheefheid als kurtosis van de dataset. Omdat beide p-waarden niet minder dan 0,05 zijn, slagen we er niet in de nulhypothese van de test te verwerpen. We hebben geen bewijs dat de drie variabelen in onze dataset geen multivariate verdeling volgen.<\/span><\/p>\n Een energietest<\/strong> is een andere statistische test die bepaalt of<\/strong> een groep variabelen al dan niet een multivariate normale verdeling volgt. De nul- en alternatieve hypothesen van de test zijn als volgt:<\/span><\/p>\n H 0<\/sub> (nul): De variabelen volgen een multivariate normale verdeling.<\/span><\/p>\n H a<\/sub> (alternatief): De variabelen volgen geen<\/em> multivariate normale verdeling.<\/span><\/p>\n De volgende code laat zien hoe u deze test in R kunt uitvoeren met behulp van het energiepakket<\/b> :<\/span><\/p>\n De p-waarde van de test is 0,31<\/strong> . Omdat dit cijfer niet minder dan 0,05 bedraagt, slagen we er niet in de nulhypothese van de test te verwerpen. We hebben geen bewijs dat de drie variabelen in onze dataset geen multivariate verdeling volgen.<\/span><\/p>\n Opmerking:<\/strong> het argument R=100 specificeert 100 boost-toegewezen replica’s die moeten worden gebruikt bij het uitvoeren van de test. Voor gegevenssets met kleinere steekproeven kunt u dit aantal verhogen om een betrouwbaardere schatting van de teststatistiek te verkrijgen.<\/em><\/span><\/p>\n Hoe u een QQ-plot in R maakt en interpreteert<\/a> Als we willen testen of een enkele variabele normaal verdeeld is of niet, kunnen we een QQ-plot maken om de verdeling te visualiseren of een formele statistische test uitvoeren, zoals een Anderson Darling-test of een Jarque-Bera-test . Als we echter willen testen of meerdere variabelen normaal verdeeld zijn als groep, moeten we een multivariate normaliteitstest […]<\/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-1068","post","type-post","status-publish","format-standard","hentry","category-gids"],"yoast_head":"\n\n
Voorbeeld: Mardia’s test in R<\/strong><\/span><\/h3>\n
library<\/span> (QuantPsyc)\n\n#create dataset\nset.seed(0)\n\ndata <- data.frame(x1 = rnorm(50),\n x2 = rnorm(50),\n x3 = rnorm(50))<\/span>\n\n#perform Multivariate normality test\nmult.norm(data)$ mult.test\n\n<\/span><\/span><\/span>Beta-hat kappa p-val\nSkewness 1.630474 13.5872843 0.1926626\nKurtosis 13.895364 -0.7130395 0.4758213\n<\/strong><\/pre>\n Voorbeeld: energietest in R<\/strong><\/span><\/h3>\n
library<\/span> (energy)\n\n#create dataset\nset.seed(0)\n\ndata <- data.frame(x1 = rnorm(50),\n x2 = rnorm(50),\n x3 = rnorm(50))<\/span>\n\n#perform Multivariate normality test\nmvnorm.etest(data, R= 100<\/span> )\n\n<\/span><\/span>Energy test of multivariate normality: estimated parameters\n\ndata: x, sample size 50, dimension 3, replicates 100\nE-statistic = 0.90923, p-value = 0.31\n<\/strong><\/pre>\n Aanvullende bronnen<\/strong><\/span><\/h3>\n
Hoe een Anderson-Darling-test uit te voeren in R<\/a>
Een Jarque-Bera-test uitvoeren in R<\/a>
Hoe een Shapiro-Wilk-test uit te voeren in R<\/a><\/p>\n","protected":false},"excerpt":{"rendered":"