让学生适应不完整的分布

Question

感谢封闭形式公式（我研究风险中性密度，这个公式的王：RND formula，第8页），我有这种类型的不完整分布：

我的想法是将这种密度与学生一致。

我已经尝试了MASS和fitdistrplus套餐，但却找不到如何执行我的任务。我现在能做的一切就是得到拟合的参数（m = 1702.041，s = 6.608536，df = 15.18036），但是从这里我不知道如何得到我的分布值。

代码示例：

temp = matrix(nrow=1000, ncol=3)
colnames(temp) = c("strikes", "first_density", "mulitply_first_density")
temp = as.data.frame(temp)

# we generate fake data
temp$strikes = seq(1000,2000,length=1000)
temp$first_density = runif(1000,max=0.006, min=1e-10)

# we multiply our first density to generate our sample
temp$mulitply_first_density = temp$first_density*1000000

# we generate our sample
vec = vector()
for (i in 1:nrow(temp))
{
    vec = c(vec, rep(temp$strike[i], temp$mulitply_first_density[i]))
}

# we laod our library
library("MASS")

# we fir our parameters
fitted_parameters = fitdistr(vec, "t")

Answer 1

使用位置和比例参数的t密度函数的公式在文档示例中给出为Map。

mydt