Een meer formele manier om naar de normaliteit te kijken, is door te testen of de kurtosis en scheefheid significant verschillen van nul.
Om dit te doen, hebben we het volgende nodig:
kurtosis.test <- function (x) {m4 <- sum ((x-mean (x)) ^ 4) / lengte (x) s4 <- var (x) ^ 2kurt <- (m4 / s4) - 3sek <- sqrt (24 / lengte (x)) totest <- kurt / sekpvalue <- pt (totest, (length (x) -1)) pvalue}
voor kurtosis, en:
skew.test <- functie (x) {m3 <- som ( (x-gemiddelde (x)) ^ 3) / lengte (x) s3 <- sqrt (var (x)) ^ 3skew <- m3 / s3ses <- sqrt (6 / lengte (x)) totest <- scheef / sespt (totest, (length (x) -1)) pval <- pt (totest, (length (x) -1)) pval}
voor Skewness.
Beide tests zijn eenzijdig, dus u moet de p-waarde met 2 vermenigvuldigen om tweezijdig te worden. Als uw p-waarde groter wordt dan één, moet u 1-kurtosis.test () gebruiken in plaats van kurtosis.test.
Als je nog andere vragen hebt, kun je me mailen op j.bredman@gmail.com