Distribuzioni continue
Distribuzioni continue
Uniforme
Esempio:
layout(matrix(c(1,2,3,4), 2, 2))
plot(seq(0,1, by=0.01), dunif(seq(0,1, by=0.01), min=0, max=1), type = "s", ylab="dunif",
xlab="La probabilità di ciascun esito")
plot(seq(0,1, by=0.01), punif(seq(0,1, by=0.01), min=0, max=1), type = "s", ylab="punif",
xlab="La probabilità di un esito minore o uguale")
plot(seq(0, 1, by=0.2), qunif(seq(0, 1, by=0.2), min=0, max=1), type = "s",
ylim = c(0,1), ylab="qunif", xlab="Suddivisione per percentili")
plot(table(floor(runif(100, min=0, max = 1)*100)/100), type="h",
ylim = c(0, 1), ylab="runif", xlab="La distribuzione di 100 lanci casuali")
Normale
Esempio: l’altezza della popolazione umana
layout(matrix(c(1,2,3,4), 2, 2))
plot(115:235, dnorm(115:235, mean = 175, sd=10), type = "s", ylab="dnorm",
xlab="La probabilità di ciascun esito")
plot(115:235, pnorm(115:235, mean = 175, sd=10), type = "s", ylab="pnorm",
xlab="La probabilità di un esito minore o uguale")
plot(seq(0, 1, by=0.02), qnorm(seq(0, 1, by=0.02), mean = 175, sd=10),
type = "s", ylab="qnorm", xlab="Suddivisione per percentili")
plot(table(floor(rnorm(1000, mean = 175, sd=10))), type="h",
ylab="rnorm", xlab="La distribuzione di 100 lanci casuali")
Gamma
Esempio:
layout(matrix(c(1,2,3,4), 2, 2))
plot(seq(1,5, by = 0.02), dgamma(seq(1,5, by = 0.02), shape = 5, rate = 2), type = "s",
ylab="dgamma", xlab="La probabilità di ciascun esito")
lines(seq(1,5, by = 0.02), dgamma(seq(1,5, by = 0.02), shape = 6, rate = 2.2))
lines(seq(1,5, by = 0.02), dgamma(seq(1,5, by = 0.02), shape = 7, rate = 2.4))
plot(seq(1,5, by = 0.02), pgamma(seq(1,5, by = 0.02), 5, rate = 2), type = "s", ylab="pgamma",
xlab="La probabilità di un esito minore o uguale")
plot(seq(0,1, by = 0.02), qgamma(seq(0,1, by = 0.02), 5, rate = 2), type = "s", ylab="qgamma", xlab="Suddivisione per percentili")
plot(table(floor(rgamma(seq(1,5, by = 0.02), 5, rate = 2))), type="h",
ylab="rgamma", xlab="La distribuzione di 100 lanci casuali")
Esponenziale
Esempio:
layout(matrix(c(1,2,3,4), 2, 2))
plot(seq(1,5, by = 0.02), dexp(seq(1,5, by = 0.02), rate = 2), type = "s",
ylab="dexp", xlab="La probabilità di ciascun esito")
lines(seq(1,5, by = 0.02), dexp(seq(1,5, by = 0.02), rate = 2.2))
lines(seq(1,5, by = 0.02), dexp(seq(1,5, by = 0.02), rate = 2.4))
plot(seq(1,5, by = 0.02), pexp(seq(1,5, by = 0.02), rate = 2), type = "s", ylab="pexp",
xlab="La probabilità di un esito minore o uguale")
plot(seq(0,1, by = 0.02), qexp(seq(0,1, by = 0.02), rate = 2), type = "s", ylab="qexp", xlab="Suddivisione per percentili")
plot(table(floor(rexp(seq(1,5, by = 0.02), rate = 0.1))), type="h",
ylab="rexp", xlab="La distribuzione di 100 lanci casuali")