Python'da özel olasılık yoğunluk işlevini tanımlama

Kendi olasılık yoğunluk fonksiyonumu tanımlamak için bazı yerleşik Python paketini (örneğin SciPy) kullanarak (daha önce herhangi bir veri olmadan, sadece ), bu yüzden onunla hesaplamalar yapabilirim (elde etme gibi) sürekli rasgele değişkenin varyansı)? Tabii ki, SymPy veya Sage alabilir, sembolik bir işlev yaratabilir ve işlemleri yapabilirim, ama tüm bu işi kendim yapmak yerine zaten uygulanmış bir paketi kullanıp kullanamayacağımı merak ediyorum.$f(x) = a x + b$

Scvy.stats içinde rv_continuous sınıfını alt sınıflara ayırmanız gerekir

import scipy.stats as st

class my_pdf(st.rv_continuous):
    def _pdf(self,x):
        return 3*x**2  # Normalized over its range, in this case [0,1]

my_cv = my_pdf(a=0, b=1, name='my_pdf')

şimdi my_cv verilen PDF ve aralık ile sürekli rasgele bir değişkendir [0,1]

Not bu örnekte my_pdfve my_cvkeyfi isimler (yani şey olabilirdi) ama vardır _pdfolduğunu değil keyfi; o ve _cdfiçinde yöntemlerdir st.rv_continuousişe subclassing için sırayla üzerine yazılır olması gereken biri.

Sympy.stats sayfasını kontrol etmelisiniz. Rastgele değişkenlerle başa çıkmak için bir arayüz sağlar. Aşağıdaki örnek, Xyoğunluk ile birim aralığında tanımlanan rastgele bir değişken sağlar2x

In [1]: from sympy.stats import *
In [2]: x = Symbol('x')
In [3]: X = ContinuousRV(x, 2*x, Interval(0, 1))

In [4]: P(X>.5) 
Out[4]: 0.750000000000000

In [5]: Var(X) # variance
Out[5]: 1/18

In [6]: E(2*cos(X)+X**2) # complex expressions are ok too
Out[6]: -7/2 + 4⋅cos(1) + 4⋅sin(1)

Eğer ilgileniyorsanız bu soyutlama bazı oldukça karmaşık manipülasyonlarla başa çıkabilir.