Basisvaardigheden in R: Grafieken: de basis
Histogram
Hieronder volgt een voorbeeld van een eenvoudig histogram gemaakt met de functie hist(). Opnieuw gebruiken we de luchtkwaliteit-dataset.
Ozon voorbeeld van een histogram
R script
with(airquality,
+ hist(Ozone,
+ xlab="Gemiddeld ozonniveau (ppb)",
+ ylab = "Frequentie",
+ main="Histogram van ozonniveaus"
+ )
+ )
Dit voorbeeld geeft een eenvoudig histogram weer dat de frequentie toont waarmee bepaalde ozonwaarden voorkomen in de dataset.
Histogram
Het volgende voorbeeld gaat over het maken van een histogram van een reeks willekeurig gegenereerde getallen voor een Poissonverdeling met dichtheidsfunctie \[p(k)=\frac{\lambda^{k}e^{-\lambda}}{k!}\] voor \(k=0,1,2,3\ldots\) en een vaste parameterwaarde \(\lambda\).
Histogram van data uit een Poissonverdeling
> set.seed(123) > rpois(n = 10,lambda = 5.5) # een voorbeeldrijtje van willekerig getrokken getallen [1] 4 7 5 8 9 2 5 8 6 5
> data <- rpois(n = 1000, lambda = 5.5)
> hist(data, breaks = 11, ylab = "frequentie", main = "histogram van data")
# suggestie van 10 bins, maar in werkelijkheid 15 gebruikt
> forced_breaks <- seq(min(data), max(data), length.out = 11) # 10 bins geforceerd
> forced_breaks
[1] 0.0 1.5 3.0 4.5 6.0 7.5 9.0 10.5 12.0 13.5 15.0
> hist(data, breaks = forced_breaks, ylab = "frequentie", main = "histogram van data")
> # dichtheidshistogram + grafiek van de dichtheidsfunctie van de Poisson distributie
> hist(data, breaks = 16, freq = FALSE, ylab = "dichtheid", main = "dichtheidshistogram")
> curve(5.5^x*exp(-5.5)/factorial(x), from = 0, to = 15, col = "blue", lwd = 3, add = TRUE)
Ontgrendel volledige toegang
