我正在写一些模拟,我发现我需要对正态分布的尾部进行过采样,这样才能得到足够的样本,对于一个特定的变量来说,样本值较低。还有什么比这个更好的吗?
from scipy.stats import norm, uniform
tail_high = .01
n_samples = 1000
tail_rvs = norm.ppf(uniform.rvs(0, tail_high, n_samples))发布于 2015-03-05 01:42:11
假设您确实需要从正态分布中取样,您可能可以使用DIY method或http://en.m.wikipedia.org/wiki/Box-Muller_transform。
像目前在for https://github.com/scipy/scipy/issues/2477中实现的一样,截断规范也有一个悬而未决的问题。原始票据提供了几个替代实现的链接。
https://stackoverflow.com/questions/28857583
复制相似问题