[도구 R과 3D 그래픽스] 3. 표면을 그려주는 persp3d함수

[도구 R과 3D 그래픽스] 3. 표면을 그려주는 persp3d함수

persp3d함수를 이용하여 이변수 표준정규분포 그래프를 그려봅시다. 

x=seq(-3,3,0.1)

y=seq(-3,3,0.1)

outer 함수를 사용하여 x,y 격자에 해당하는 z값을 얻습니다. outer 함수의 사용 방법을 예로 들면, 아래와 같습니다. 둘의 곱을 계산하는 경우입니다. 

x=c(1,2,3)

y=c(1,2,3)

z=outer(x,y,"*")

> z

     [,1] [,2] [,3]

[1,]    1    2    3

[2,]    2    4    6

[3,]    3    6    9

outer 함수를 우리 예제에 적용합시다. 

x=seq(-3,3,0.1)

y=seq(-3,3,0.1)

f=function(x,y){dnorm(x)*dnorm(y)}

z=outer(x,y,f)

그래프를 그려봅시다

persp3d(x, y, z, col = "yellow",

        xlab = "X", ylab = "Y", zlab = "probability density")

만약 축의 길이에 실제 값이 반영되게 하고 싶으면 aspect="iso" 옵션을 추가하면 됩니다.