R에 대한 커널 밀도 추정을 통한 시뮬레이션

Question

R. Dataset의 커널 분포에 기반한 시뮬레이션을 수행하려고했습니다. - HSAUR2 패키지의 CYG OB1이 사용되었습니다. 데이터 세트를 사용하여 커널 밀도 추정치를 찾기위한 분석을 수행했습니다. 이 커널 밀도를 사용하여 이변 수 값을 시뮬레이션하려고했습니다. 지금까지 사용 된 코드

> CYGOB1d <- bkde2D(CYGOB1, bandwidth = sapply(CYGOB1, dpik)) 
> plot(CYGOB1, xlab = "log surface temperature", ylab = "log light intensity") 
> contour(x = CYGOB1d$x1, y=CYGOB1d$x2, z=CYGOB1d$fhat, add = TRUE) 
> persp(x=CYGOB1d$x1, y=CYGOB1d$x2, z = CYGOB1d$fhat, xlab = "log surface 
    temperature", ylab = "log light intensity", zlab ="density") 

어떻게 커널 밀도를 기반으로 시뮬레이션 (1000 회 실행)을 수행 할 수 있습니까? z이 확률 변수 인 경우

Answer 1

, 당신은

d <- density(z,n=512) 
#n is the number of points used for the curve, should be a power of two (512 is default) 

sample(x = d$x,prob = d$y,size=1000,replace=TRUE) 
#samples from the n values of x, according to prob y (the density) 

이는 계산을 수행하는 방법에 대한 다양한 옵션을 ?density를 참조 할 수 있습니다.

Answer 2

답장을 보내 주셔서 감사합니다. I 위의 코드 및 I는 다음과 얻고 용액을 사용하여 시도 (화상으로서 부착 된) 데이터가 이변 량 데이터이다

simulation results

동안 다음과 같이

 logst logli 
[1,] 4.37 5.23 
[2,] 4.56 5.74 
[3,] 4.26 4.93 
[4,] 4.56 5.74 
[5,] 4.30 5.19 
[6,] 4.46 5.46 
[7,] 3.84 4.65 
[8,] 4.57 5.27 
[9,] 4.26 5.57 
[10,] 4.37 5.12 

나는 시뮬레이션 결과가 내게 비슷한 숫자의 쌍을 줄 것이라고 기대했다.