R에서 t 분포를 그리는 방법

R의 at 분포에 대한 확률 밀도 함수를 플롯하려면 다음 함수를 사용할 수 있습니다.

• 확률밀도함수를 생성하기 위한 dt(x, df)
• curve(function, from = NULL, to = NULL)를 사용하여 확률 밀도 함수를 플롯합니다.

확률 밀도 함수를 플롯하려면 dt() 함수에 df (자유도)를 지정하고 curve() 함수에 from 및 to 값을 지정해야 합니다.

예를 들어, 다음 코드는 도표의 x축이 -4와 4 사이인 자유도가 10인 분포에 대한 확률 밀도 함수를 도표화하는 방법을 보여줍니다.

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

정규 분포 와 유사하게 t 분포는 평균 0을 기준으로 대칭입니다.

제목을 추가하고, Y축 레이블을 변경하고, 행 너비를 늘리고, 행 색상을 변경하여 플롯을 더욱 보기 좋게 만들 수 있습니다.

 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 

또한 그래프에 여러 곡선을 추가하여 다양한 자유도로 t 분포를 비교할 수도 있습니다. 예를 들어, 다음 코드는 df = 6, df = 10, df = 30인 t 분포도를 만듭니다.

 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() 함수를 사용하여 플롯에 범례를 추가할 수 있습니다.

범례(x, y=NULL, 범례, 채우기, col, bg, lty, cex)

금:

• x, y: 범례 위치를 지정하는 데 사용되는 x 및 y 좌표
• legend: 범례에 넣을 텍스트
• 채우기: 범례 내부의 채우기 색상
• col: 범례 내부 선에 사용할 색상 목록
• bg: 범례의 배경색
• lty: 선 스타일
• cex: 범례의 텍스트 크기

이 예에서는 다음 구문을 사용하여 범례를 만듭니다.

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