R에서 로그 정규 분포를 그리는 방법

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

• dlnorm(x, 평균 로그 = 0, sdlog = 1)을 사용 하여 확률 밀도 함수를 생성합니다.
• curve(function, from = NULL, to = NULL)를 사용하여 확률 밀도 함수를 플로팅합니다.

예를 들어, 다음 코드는 플롯의 x축이 0에서 10까지인 평균 = 0이고 표준 편차 = 1(로그 척도)인 로그 정규 분포에 대한 확률 밀도 함수를 플롯하는 방법을 보여줍니다.

 curve(dlnorm(x, meanlog=0, sdlog=1), from=0, to=10)

기본적으로 평균 로그 = 0 및 sdlog = 1입니다. 이는 dlnorm() 함수에서 이러한 매개변수를 지정하지 않고도 정확히 동일한 플롯을 생성할 수 있음을 의미합니다.

 curve(dlnorm(x), from=0, to=10)

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

 curve(dlnorm(x), from=0, to=10, 
    main = 'Log Normal Distribution', #add title
    ylab = 'Density', #change y-axis label
    lwd = 2, #increase line width to 2
    col = 'steelblue') #change line color to steelblue 

또한 로그 정규 분포를 다양한 표준 편차와 비교하기 위해 그래프에 여러 곡선을 추가할 수도 있습니다. 예를 들어, 다음 코드는 sdlog = 0.3, sdlog = 0.5, sdlog = 1인 정규 분포 도표를 만듭니다.

 curve(dlnorm(x, meanlog=0, sdlog=.3), from=0, to=10, col='blue')
curve(dlnorm(x, meanlog=0, sdlog=.5), from=0, to=10, col='red', add=TRUE)
curve(dlnorm(x, meanlog=0, sdlog=1), from=0, to=10, col='purple', add=TRUE)

다음 구문을 사용하는 legend() 함수를 사용하여 플롯에 범례를 추가할 수 있습니다.

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

금:

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

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

 #create density plots
curve(dlnorm(x, meanlog=0, sdlog=.3), from=0, to=10, col='blue')
curve(dlnorm(x, meanlog=0, sdlog=.5), from=0, to=10, col='red', add=TRUE)
curve(dlnorm(x, meanlog=0, sdlog=1), from=0, to=10, col='purple', add=TRUE)

#add legend
legend(6, 1.2, legend=c("sdlog=.3", "sdlog=.5", "sdlog=1"),
       col=c("blue", "red", "purple"), lty=1, cex=1.2)