1. Wczytuje dane z pliku
dane := READPRN("15.txt" )
T
0 1 2 3 4 5 6 7 8 9
dane =
0 360.9 422.9 283.8 318 283.4 580.9 572.1 414.1 413.9 315.9
2. Porzadkujemy dane
d := sort(dane)
T
0 1 2 3 4 5 6 7 8 9
d =
0 251.6 283.4 283.8 315.9 318 321.2 341.8 360.9 385.5 413
3. Ile razy wystapila kazda dana
w := rows(d) w = 30
k := cols(d) k = 1
j := 0,1.. w - 1
Mj := 1
ile_razy_wystj := 1
T
0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19
M =
0 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
Teraz czestosci wystepowania dodajemy do siebie
j
czest_wystj := ile_razy_wystj
"
j = 0
T
0 1 2 3 4 5 6 7 8 9
czest_wyst =
0 1 2 3 4 5 6 7 8 9 10
Teraz zapisujemy w osobnym pliku czest_wyst
WRITEPRN("C_W.dat") := czest_wyst
WRITEPRN("D_ZRED.dat" ) := d
0
Test serii:
0 360.9
dane := READPRN("15.txt." )
1 422.9
2 283.8
Liczba pomiarów:
3 318
4 283.4
5 580.9
n := rows(dane) n = 30
6 572.1
dane =
7 414.1
Zakres iteracji
8 413.9
9 315.9
i := 0,1.. n - 1
10 434.3
Mediana
11 423.2
12 468.4
med := median(dane) med = 423.75
13 617.6
Serie czyli ilosc zer i jedynek 14 421.5
15 656.7
yi := if danei > med, 1, 0
( )
y1 := yi y1 = 15
"
i
y0 := n - y1 y0 = 15
1
yi 0.5
0
0 5 10 15 20 25 30
i
Liczba serii
1
yi
0.5
yi
0
0 5 10 15 20 25 30
i
n-1
serie := 1 + if yi-1 `" yi, 1,0 serie = 12
( )
"
i = 1
Wzor ten musi miec postac 1+..., gdyz w rzeczywistosci sumujemy nie serie tylko
zmiany wartosci z 1 na 0 oraz z 0 na 1.
Liczba tych zmian to wlasnie liczba serii
Wniosek
2<=12<=15+15 - warunek ogólny
12<=!11(U0,05)
Hipoteze o niezaleznosci wynikow (na poziomie istotnosci 0,05) przyjmujemy , poniewaz
liczba serii jest wieksza od wartosci krytycznej liczby serii pomiarow
Test bledow grubych
0
0 360.9
1 422.9
dane := READPRN("15.txt." )
2 283.8
3 318
Zakres iteracji
4 283.4
5 580.9
i := 0,1.. n - 1
6 572.1
dane =
7 414.1
Wartosc minimalna
8 413.9
9 315.9
w_min := min(dane) w_min = 251.6
10 434.3
Wartosc max
11 423.2
12 468.4
w_max := max(dane) w_max = 686.6
13 617.6
Wartosc srednia
14 421.5
w_sred := mean(dane) w_sred = 446.34
15 656.7
Odchylenie srednie skorygowane
odch_sred := Stdev(dane) odch_sred = 117.729
Dane nieuporzadkowane
800
600
danei
w_sred
400
200
0 5 10 15 20 25 30
i
Dane uporzadkowane
dane_up := sort(dane)
800
600
dane_upi
w_sred
400
200
0 5 10 15 20 25 30
i
Obliczanie statystyk B4
(w_max - w_sred)
B4_plus := B4_plus = 2.041
odch_sred
(w_sred - w_min)
B4_minus := B4_minus = 1.654
odch_sred
Wartosc krytyczna (str. 103)
dla poziomu istotnosci 0.05 wartosc B4_plus i B4_minus nie sa obarczone bledem grubym
poniewaz nie jest spelniona nierownosc
Bi<=bi(a,n)