(numpy.random.standard\)如何确定自由度？

Question

请告诉我如何确定自由度

1。我使用下面的代码绘制一个带有直方图的student-t分布图

# distribution fitting
param = t.fit(log_rets) 
# now get theoretical values in our interval      
pdf_t = stats.t.pdf(lnspc, param[0], param[1], param[2])
# plot it      
plt.plot(lnspc, pdf_t, label="student t") 
plt.show()

2。我发现参数[3]如下所示

^{pr2}$

3。我会用student-t生成一个随机数，但不知道如何确定自由度

retvec = []
for i in range (n_obs):
    retvec.append(np.random.standard_t(**<degrees of freedom>**, size=n_assets))

请随时发表评论。在

你好，纳尔逊。在

Answer 1

这是一个老问题，但是对于任何阅读它的人来说，params[0]这里是自由度（例如，尾部大小），因此您可以使用scipy随机变量生成随机值：

roll = stats.t.rvs(df=params[0], loc=params[1], scale=params[2], size=n_assets)