2011-05-03 3 views
1

나는 이산 데이터의 벡터를 가지고 있으며,이 데이터와 관련된 경험적 분포를 시뮬레이트하고 싶습니다. 적합을 수행 한 후 rlogspline 함수로 시뮬레이션했습니다. < -logspline (vector_of_data) 여기서 vector_of_data는 데이터가 연속 배포에서 오는 것입니다. 그 때문에 logpline을 사용했지만,이 벡터를 사용하면이 벡터의 값이 개별적인 특성을 갖기 때문에 Logpline을 사용하여 해당 값에 "적합"을 조정할 수 없습니다 .이산 데이터 벡터로 시뮬레이션하기

기본적으로 내가하고 싶은 것은 관찰 된 데이터의 "적합"을 조정 한 다음 그 값을 시뮬레이트하기 위해 그 적합성을 사용하는 것입니다. 이것이 R에서 할 수 있다고 생각합니까?

도움 주셔서 감사합니다. 나는 sample(x,...,replace=TRUE)가 (교체와 샘플링) 경험적 분포 시뮬레이션해야한다고 생각

+0

당신은 당신이 logspline에서 얻을 시뮬레이션에) (원형 적용 단지 수없는 이유는 무엇입니까? –

+0

그래, 그럴 수는 있어도 문제는 단지 양의 이산 값 벡터 만 있으면 그 값이 음수 인 시뮬레이션으로 끝날 수 있다는 것입니다. 왜냐하면 스플라인이 0에서 계산되는 방식 때문입니다. – natorro

+0

아니요. logpline에 대한 두 번째 및 세 번째 인수를 신중하게 고려하지 않았습니다. –

답변

2

...

+0

네,하지만 여기서 문제는 제가 이미 관찰 한 값을 돌려 줄 것입니다. 이미 관찰되지 않은 "실제"밀도의 값을 갖고 싶습니다. 밀도. – natorro

0

나는 나는 당신이하려고하는,하지만 당신은 예를 들어, quantilerunif 같은 것을 사용할 수 있습니다 정확히 완전히 명확하지 :

obs <- c(125,110,115,100,150)    # original observations 
sim <- quantile(obs, runif(10000))  # simulations 
hist(sim, freq=FALSE) 

sim histogram

관련 문제