Distribuzioni continue

Posted by On Filed under Probabilità & statistica, R

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")