Jak wykreślić rozkład t w r

Aby wykreślić funkcję gęstości prawdopodobieństwa dla rozkładu at w R, możemy skorzystać z następujących funkcji:

• dt(x, df) , aby utworzyć funkcję gęstości prawdopodobieństwa
• krzywa(funkcja, od = NULL, do = NULL), aby wykreślić funkcję gęstości prawdopodobieństwa

Aby wykreślić funkcję gęstości prawdopodobieństwa, musimy określić df (stopnie swobody) w funkcji dt() oraz wartości od i do w funkcji curve() .

Na przykład poniższy kod ilustruje sposób wykreślenia funkcji gęstości prawdopodobieństwa dla rozkładu o 10 stopniach swobody, gdzie oś x wykresu mieści się w przedziale od -4 do 4:

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

Podobnie do rozkładu normalnego , rozkład t jest symetryczny wokół średniej wynoszącej 0.

Możemy dodać tytuł, zmienić etykietę osi Y, zwiększyć szerokość wiersza, a nawet zmienić kolor wiersza, aby wykres był bardziej estetyczny:

 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 

Możemy także dodać do wykresu wiele krzywych, aby porównać rozkłady t o różnych stopniach swobody. Na przykład poniższy kod tworzy t wykresów rozkładu z df = 6, df = 10 i 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)

Legendę do wykresu możemy dodać za pomocą funkcji legend() , która przyjmuje następującą składnię:

legenda(x, y=NULL, legenda, wypełnienie, col, bg, lty, cex)

Złoto:

• x, y: współrzędne x i y użyte do umiejscowienia legendy
• legenda: tekst, który należy umieścić w legendzie
• fill: kolor wypełnienia legendy
• col: lista kolorów linii wewnątrz legendy
• bg: kolor tła legendy
• lty: styl linii
• cex: rozmiar tekstu w legendzie

W naszym przykładzie użyjemy następującej składni do utworzenia legendy:

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