250x250
Link
๋‚˜์˜ GitHub Contribution ๊ทธ๋ž˜ํ”„
Loading data ...
Notice
Recent Posts
Recent Comments
๊ด€๋ฆฌ ๋ฉ”๋‰ด

Data Science LAB

[Python] ๋“ฑ๋ถ„์‚ฐ๊ฒ€์ • (์ •๊ทœ๋ถ„ํฌ์˜ ๋ชจ๋ถ„์‚ฐ์— ๋Œ€ํ•œ ๊ฒ€์ •) ๋ณธ๋ฌธ

๐Ÿ›  Machine Learning/๊ธฐ์ดˆ ํ†ต๊ณ„

[Python] ๋“ฑ๋ถ„์‚ฐ๊ฒ€์ • (์ •๊ทœ๋ถ„ํฌ์˜ ๋ชจ๋ถ„์‚ฐ์— ๋Œ€ํ•œ ๊ฒ€์ •)

ใ…… ใ…œ ใ…” ใ…‡ 2022. 9. 20. 05:02
728x90

๋“ฑ๋ถ„์‚ฐ๊ฒ€์ •

1. ์–‘์ธก๊ฒ€์ • 
- ๊ท€๋ฌด๊ฐ€์„ค : ํ‘œ๋ณธ์˜ ๋ถ„์‚ฐ์€ ๊ฐ™๋‹ค.
- ๋Œ€๋ฆฝ๊ฐ€์„ค : ํ‘œ๋ณธ์˜ ๋ถ„์‚ฐ์€ ๊ฐ™์ง€ ์•Š๋‹ค.

2. ๋‹จ์ธก ๊ฒ€์ •
- ๊ท€๋ฌด๊ฐ€์„ค : ํ‘œ๋ณธ์˜ ๋ถ„์‚ฐ์€ ๊ฐ™๋‹ค.
- ๋Œ€๋ฆฝ๊ฐ€์„ค : ํ‘œ๋ณธ์˜ ๋ถ„์‚ฐ์€ ์ž‘๋‹ค. or ํฌ๋‹ค

 

 

ํŒŒ์ด์ฌ์œผ๋กœ ํ•จ์ˆ˜ ๊ตฌํ˜„

from scipy.stats import chi2

def var_test(x, va0, direction = 'two-tailed', alpha = 0.05):
    n = len(x)
    Q = (n-1) * np.var(x) / va0
    print('direction', direction)
    
    if direction == 'lower':
        q = chi2.ppf(alpha, n-1)
        
        if Q <= q:
            return 'H_0 rejected'
        else:
            return 'H_0 not rejected'
        
    elif direction == 'upper':
        q = chi2.ppf(1-alpha, n-1)
        if Q >= q:
            return 'H_0 rejected'
        else:
            return 'H_0 not rejected'
        
    else:
        q1 = chi2.ppf(alpha/2, n-1)
        q2 = chi2.ppf(1 - (alpha/2), n-1)
        
        print('Q',Q)
        print('q1',q1)
        print('q2',q2)
        
        if Q <= q1 or Q >= q2:
            return 'H_0 rejected'
        else:
            return 'H_0 not rejected'
        
var_test(df2['content'], va0 = 1.3)

 

728x90
Comments