Hoe een t-verdeling in r te plotten

Om de kansdichtheidsfunctie voor de at-verdeling in R te plotten, kunnen we de volgende functies gebruiken:

• dt(x, df) om de kansdichtheidsfunctie te creëren
• curve(functie, van = NULL, tot = NULL) om de kansdichtheidsfunctie uit te zetten

Om de waarschijnlijkheidsdichtheidsfunctie te plotten, moeten we df (vrijheidsgraden) specificeren in de dt() functie, evenals de from en to waarden in de curve() functie.

De volgende code illustreert bijvoorbeeld hoe u een kansdichtheidsfunctie plot voor een verdeling met 10 vrijheidsgraden waarbij de x-as van de grafiek tussen -4 en 4 ligt:

 curve(dt(x, df=10), from=-4, to=4)

Net als bij de normale verdeling is de t-verdeling symmetrisch rond een gemiddelde van 0.

We kunnen een titel toevoegen, het label op de Y-as wijzigen, de rijbreedte vergroten en zelfs de rijkleur wijzigen om de plot esthetisch aantrekkelijker te maken:

 curve(dt(x, df=10), from=-4, to=4, 
    main = 't Distribution (df = 10)', #add title
    ylab = 'Density', #change y-axis label
    lwd = 2, #increase line width to 2
    col = 'steelblue') #change line color to steelblue 

We kunnen ook meerdere curven aan de grafiek toevoegen om t-verdelingen met verschillende vrijheidsgraden te vergelijken. Met de volgende code worden bijvoorbeeld t-verdelingsplots gemaakt met df = 6, df = 10 en df = 30:

 curve(dt(x, df=6), from=-4, to=4, col='blue') 
curve(dt(x, df=10), from=-4, to=4, col='red', add=TRUE)
curve(dt(x, df=30), from=-4, to=4, col='green', add=TRUE)

We kunnen een legenda aan de plot toevoegen met behulp van de legend()- functie, die de volgende syntaxis heeft:

legenda(x, y=NULL, legenda, fill, col, bg, lty, cex)

Goud:

• x, y: de x- en y-coördinaten die worden gebruikt om de legenda te positioneren
• legenda: de tekst die in de legenda moet worden geplaatst
• vulling: vulkleur binnen de legenda
• col: de lijst met kleuren die moeten worden gebruikt voor de lijnen in de legenda
• bg: de achtergrondkleur van de legenda
• lty: lijnstijl
• cex: grootte van de tekst in de legenda

In ons voorbeeld gebruiken we de volgende syntaxis om een legenda te maken:

 #create density plots
curve(dt(x, df=6), from=-4, to=4, col='blue') 
curve(dt(x, df=10), from=-4, to=4, col='red', add=TRUE)
curve(dt(x, df=30), from=-4, to=4, col='green', add=TRUE)

#add legend
legend(-4, .3, legend=c("df=6", "df=10", "df=30"),
       col=c("blue", "red", "green"), lty=1, cex=1.2)