R中概率密度(pdf)的提取

i 表示y的每个预测值(连同密度),导出 model_i t .

分母是上述乘积的积分。我的问题是:如何估计密度。要得到变量的密度,需要一些数据点。到目前为止我有这个:

y<-c(-0.00604,-0.00180,0.00292,-0.0148)
forecastsy_model1<-c(-0.0183,0.00685) # respectively time t=1 and t=2 of the forecasts
forecastsy_model2<-c(-0.0163,0.00931) # similarly
all.y.1<-c(y,forecasty_model1) #together in one vector 
all.y.2<-c(y,forecasty_model2) #same

x1 暂时的 t=1 t=6 为了做产品。我已经考虑过用这个来计算密度:

dy1<-density(all.y.1)
which(dy1$x==0.00685)
integer(0) #length(dy1$x) : 512

具有 dy1$x 包含 n 根据文件估计密度的点的坐标。不应该 为6,或至少包含 y 我提供的?什么是提取密度(pdf)的正确方法 是的 ?