khoảng tin cậy của sinh viên t trong python

Question

Tôi quan tâm đến việc sử dụng python để tính toán khoảng tin cậy từ sinh viên t.

Tôi đang sử dụng các chức năng StudentTCI() trong Mathematica và bây giờ cần phải mã hóa các chức năng tương tự trong python http://reference.wolfram.com/mathematica/HypothesisTesting/ref/StudentTCI.html

Tôi không khá chắc chắn làm thế nào để xây dựng chức năng này bản thân mình, nhưng trước khi tôi bắt tay vào đó, là thế này chức năng trong python một nơi nào đó? Như numpy? (Tôi đã không sử dụng numpy và cố vấn của tôi khuyên không sử dụng gumpy nếu có thể).

Cách dễ nhất để giải quyết vấn đề này là gì? Tôi có thể sao chép mã nguồn từ StudentTCI() trong numpy (nếu nó tồn tại) vào mã của tôi như là một định nghĩa chức năng?

chỉnh sửa: Tôi sẽ cần phải xây dựng TCI Sinh viên bằng mã python (nếu có thể). Cài đặt scipy đã biến thành một ngõ cụt. Tôi có cùng một vấn đề mà mọi người khác đang gặp phải, và không có cách nào tôi có thể yêu cầu Scipy cho mã tôi phân phối nếu mất quá nhiều thời gian để thiết lập.

Bất kỳ ai biết cách xem mã nguồn cho thuật toán trong phiên bản scipy? Tôi đang suy nghĩ tôi sẽ tái cấu trúc nó thành một định nghĩa python.

Answer 1

Tôi đoán bạn có thể sử dụng scipy.stats.t và nó interval phương pháp:

In [1]: from scipy.stats import t 
In [2]: t.interval(0.95, 10, loc=1, scale=2) # 95% confidence interval 
Out[2]: (-3.4562777039298762, 5.4562777039298762) 
In [3]: t.interval(0.99, 10, loc=1, scale=2) # 99% confidence interval 
Out[3]: (-5.338545334351676, 7.338545334351676) 

Chắc chắn, bạn có thể thực hiện chức năng của riêng bạn nếu bạn muốn. Chúng ta hãy làm cho nó trông giống như trong Mathematica:

from scipy.stats import t 


def StudentTCI(loc, scale, df, alpha=0.95): 
    return t.interval(alpha, df, loc, scale) 

print StudentTCI(1, 2, 10) 
print StudentTCI(1, 2, 10, 0.99) 

Kết quả:

(-3.4562777039298762, 5.4562777039298762) 
(-5.338545334351676, 7.338545334351676)