주기 함수의 플롯 밀도

Question

나는 여러 사건이 일어난 시간 (초 자정)을 나타내는 벡터를 얻었으며 그날의 밀도를 그래프로 나타내려고한다. 다음과 같은 한 가지 방법이 있습니다.

rs <- 60*60*24*c(rbeta(5000, 2, 5), runif(10000, 0, 1)) 
den <- density(rs, cut=0) 
plot(den, ylim=range(0,den$y)) 

순환 점 기능이므로 끝점 밀도가 잘못 표시됩니다. 당신이 연속 3 개 기간 플롯 경우에 당신은 중간 기간의 진정한 밀도를 참조하십시오

den <- density(c(rs, rs+60*60*24, rs+2*60*60*24), cut=0) 
plot(den, ylim=range(0,den$y)) 

내 질문에 원본 데이터로부터의 중간 덩어리의 밀도를 얻을 수있는 몇 가지 [더] 방법은 세 겹으로하지 않고, 거기 여부 내가 한 것처럼 관측 수. 물론 종점 근처에 관측치가없는 경우에는 기간의 길이를 제공해야합니다.

Answer 1

커브가 줄을 따라야한다고하는 증거는 반복되는 스플라인 맞춤의 일부가 매력적이라고 ​​생각하지 않습니다. 시간 경계에서 구분을 지정하는 동일한 오브젝트에서 hist()의 결과를 조사해야합니다. logspline 기능은 HTE 데이터에 지정된 범위와 밀도 추정의 계산을 할 수 있습니다 : 제대로 끝 - 중 - 일의 밀도가 일의 첫 번째보다 낮다는 사실을 포착하기 때문에

hist(c(c(rs, rs+60*60*24, rs+2*60*60*24), breaks= 24*3) 
require(logspline) 
?logspline 
fit <- logspline(c(rs), lbound=0, ubound=60*60*24) 
plot(fit) 

더 잘 맞는있는 그 트리플 데이의 밀도 그림으로 실제로 제대로 캡쳐되지는 않습니다.