у і параметрів закону розподілу випадкової величини. Для ряду практично важливих випадків цей закон можна визначити з теоретичних міркувань.
3.3 Реалізація на мові програмування Pascal
бутстрап комп'ютерний інтервал статистика
Program Ravn_V_K; Uses CRT;
{Дослідження бутстраповскіх інтервальних оцінок для вибірки з рівномірним розподілом}
Const T5=1.67655089261685; T2_5=2.00957523712924; {t-статистики на 49 ст. св}
{Всі константи розраховані на вибірку об'ємом 50 елементів}
Function Ravn (A, B: real): real;
Begin
Result:=A + (B-A) * Random;
end ; {Ravn}
Var X: array [1..1000] of real; {Вибірка}, S, A_, B_: real; {Параметри вибірки}
N, M, G, B: integer; {Параметри експерименту}
PA: array [1..3] of integer; {Оцінки Класичні }: array [1..3, 0..20] of integer; {Бутстраповскіе оцінки}
Ab: array [1..10000] of real; {Бутстраповскіе статистики}, J, K, L, Q, H: integer; {Допоміжні змінні} _, S_, R: real; {Допоміжні змінні}: Text; {Файл з вхідними, а потім і вихідними даними}
St: String;
BEGIN
{Налаштування моделі} (Dat, Param.txt ); ReSet (Dat); (Dat, N); ReadLn (Dat, A); ReadLn (Dat, S); ReadLn (Dat, M); (Dat, G); ReadLn (Dat, B); Close (Dat); ( бутстрапа рівномірна модель ( + intToStr (N) + + floatToStr (A) + + floatToStr (S) + ) ); ; Randomize;
{Ініціалізація лічильників}
For I:=1 to 3 do begin
PA [I]:=0;
For J:=0 to 20 do PbA [I , J]:=0;
end ; _:=A - S * sqrt (3); B_:=A + S * sqrt (3);
{Прогони моделі}
For I:=1 to M do Begin (I mod 10)=0 then Begin
GoToXY (33,12); Write (I)
end ;
{Генеруємо вибірку}
For J:=1 to N do X [J ]:=Ravn (A_, B_);
{Реалізуємо класику} _:=0;
For J:=1 to N do Begin
X _ +=X [J]; S _ +=Sqr (X [J])
End ; _ /=N; S _:=(S_-Sqr (X _) * N)/(N - 1);
If (A gt; X_-Sqrt (S_/N) * T5) then Inc (PA [1]) ;
If (A lt; X_ + Sqrt (S_/N) * T5) then Inc (PA [2]) ;
If (A gt; X_-Sqrt (S_/N) * T2_5) and (A lt; X_ + Sqrt ( S_/N) * T2_5) then Inc (PA [3]);
{Реалізуємо бутстрап}
For K:=0 to 20 do begin
Q:=K + B; {Обсяг бутстраповской вибірки}
For L:=1 to G do begin
Ab [L]:=0;
for J:=1 to Q do Begin
H:=Random (N) +1; Ab [L] +=X [H];
end ; [L] /=Q;
end ;
{Сортування результатів}
For L:=1 to G - 1 do begin
R:=Ab [L]; H:=L;
for J:=L + 1 to G do R gt ; Ab [J] then begin
R:=Ab [J]; H:=J
end ; [H]:=Ab [L]; Ab [L]:=R
end ;
IF (A gt; Ab [trunc (G * 0.05)]) then inc (PbA [1, K] );
IF (A lt; Ab [trunc (G * 0.95)]) then in...