1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45
|
invisible(sapply(c(
"coin", "ggplot2", "ggthemes", "latex2exp", "magrittr", "SuppDists", "pgirmess"),
function(z) library(z, character.only = TRUE)
))
biomasse <- read.delim("biomasse.tsv")
biomaSansValManq <- sapply(biomasse, function(z) {
tmp <- na.exclude(z)
attr(tmp, "na.action") <- NULL
return(tmp)
})
biomaSansValManq$relevé <- unlist(biomaSansValManq)
biomaSansValManq$effectif <- sapply(biomasse, function(z) na.exclude(z) %>% length)
biomaSansValManq$station <- rep(colnames(biomasse), biomaSansValManq$effectif) %>% factor
###stats
sorties$stats <- kruskal.test(x = biomasse, na.action = "na.exclude")
print(sorties$stats)
print(sorties$stats$statistic)
print(sorties$stats$p.value)
print(sorties$stats$parameter)
sorties$alpha <- 1e-3
sorties$décision <- ifelse(
1L - pchisq(sorties$stats$statistic, sorties$stats$parameter) < sorties$alpha,
paste0("Rejet de H0 avec alpha = ", sorties$alpha, ".\n"),
"Non-rejet de H0.\n"
)
cat(sorties$décision)
###coin
sorties$coin <- kruskal_test(
formula = biomaSansValManq$relevé ~ biomaSansValManq$station,
distribution = approximate(B = 1000L)
)
show(sorties$coin)
statistic(sorties$coin, type = "linear")
statistic(sorties$coin, type = "standardized")
statistic(sorties$coin, type = "test")
statistic(sorties$coin)
pvalue(sorties$coin)
expectation(sorties$coin) |
Partager