Scipy - Stats - 概率分布å‚æ•°çš„å«ä¹‰

时间:2013-07-23 15:19:38

标签: python numpy scipy

Scipy docs将指数使用的分布形å¼è¡¨ç¤ºä¸ºï¼š

expon.pdf(x) = lambda * exp(- lambda*x)

然而,拟åˆå‡½æ•°éœ€è¦ï¼š

fit(data, loc=0, scale=1)

rvs函数需è¦ï¼š

rvs(loc=0, scale=1, size=1)

问题1: 为什么无关的ä½ç½®å˜é‡ï¼Ÿæˆ‘知é“指数åªæ˜¯æ›´ä¸€èˆ¬åˆ†å¸ƒï¼ˆä¼½çŽ›ï¼‰çš„特定形å¼ï¼Œä½†ä¸ºä»€ä¹ˆè¦åŒ…括ä¸éœ€è¦çš„ä¿¡æ¯å‘¢ï¼Ÿç”šè‡³gamma也没有ä½ç½®å‚数。

问题2: æ‹Ÿåˆï¼ˆ...)的输出与输入å˜é‡çš„顺åºç›¸åŒã€‚我的æ„æ€æ˜¯ 如果我这样åšï¼š

t = fit([....]) , t will have the form t[0], t[1]

我应该将t [0]解释为形状,将t 1解释为比例。

这适用于所有å‘行版å—?

伽玛如何:

fit(data, a, loc=0, scale=1)

1 个答案:

答案 0 :(得分:4)

  1. æ¯ä¸ªå•å˜é‡æ¦‚率分布,无论其通常的公å¼å¦‚何,都å¯ä»¥æ‰©å±•åˆ°åŒ…括ä½ç½®å’Œæ¯”例å‚数。有时,这需è¦å°†åˆ†å¸ƒçš„支æŒä»Žæ­£/éžè´Ÿå®žæ•°æ‰©å±•åˆ°æ•´ä¸ªå®žæ•°è¡Œï¼Œå½“低于loc值时,仅将PDF值设置为0。 scipy.stats执行此æ“作,将locå’Œscale的所有处ç†æ–¹å¼ç§»è‡³æ‰€æœ‰å‘行版共享的公共方法。有人建议删除它,并使gamma loc之类的分布更少éµå¾ªå®ƒä»¬çš„规范表述。然而,事实è¯æ˜Žæœ‰äº›äººå®žé™…上使用éžé›¶locå‚数的“移ä½ä¼½é©¬â€åˆ†å¸ƒæ¥æ¨¡æ‹Ÿå¤ªé˜³é»‘å­çš„大å°ï¼Œå¦‚果我没记错的è¯ï¼Œscipy.stats的当å‰è¡Œä¸ºå¯¹ä»–们æ¥è¯´æ˜¯å®Œç¾Žçš„。所以我们ä¿ç•™å®ƒã€‚

  2. 如果有fit()形状å‚数,(shape0, shape1, ..., shapeN, loc, scale)方法的输出是Nå½¢å¼çš„元组。对于没有形状å‚æ•°çš„æ­£æ€åˆ†å¸ƒï¼Œå®ƒå°†ä»…返回(loc, scale)。对于具有一个的伽玛分布,它将返回(shape, loc, scale)。多个形状å‚数的顺åºä¸Žåˆ†é…中的æ¯ä¸ªå…¶ä»–方法的顺åºç›¸åŒã€‚ This holds for all distributions.