Пример расчета и оформления результатов

Рассматривается ряд наблюдений за суммами осадков за ноябрь по метеостанции Новый Порт. Хронологический график ряда наблюдений приведен на рис.1.

Рис.1. Хронологический график сумм осадков за ноябрь по метеостанции Новый Порт.

Прежде всего, был проведен анализ однородности эмпирического распределения осадков по статистическим критериям Диксона и Смирнова - Граббса. Для этой цели ряд наблюдений был ранжирован и полученные значения приведены в табл.4. Эмпирическая обеспеченность ранжированного ряда (Р%) определялась по формуле:

, ( 24 )

где m - порядковый номер членов ряда метеорологической характеристики, расположенного в убывающем порядке; n - общее число членов ряда.

Таблица 4.

Рассчитанные значения эмпирической обеспеченности и ряд ранжированных осадков за ноябрь по метеостанции Новый Порт.

Обеспеченность, P %	Значение ранжированных осадков (мм)	Год
1.370	57.0
2.740	52.0
4.110	46.0
5.479	44.0
6.849	44.0
8.219	40.0
9.589	38.0
10.959	37.0
12.329	37.0
13.699	35.0
15.068	35.0
16.438	33.0
17.808	31.0
19.178	29.0
20.548	29.0
21.918	29.0
23.288	28.0
24.658	26.0
26.027	25.0
27.397	23.0
28.767	22.0
30.137	21.0
31.507	20.0
32.877	20.0
34.247	20.0
35.616	20.0
36.986	20.0
38.356	19.0
39.726	19.0
41.096	19.0
42.466	18.0
43.836	18.0
45.205	17.0
46.575	17.0
47.945	17.0
49.315	17.0
50.685	16.0
52.055	16.0
53.425	15.0
54.795	15.0
56.164	15.0
57.534	13.0
58.904	13.0
60.274	13.0
61.644	13.0
63.014	13.0
64.384	13.0
65.753	13.0
67.123	13.0
68.493	12.0
69.863	12.0
71.233	12.0
72.603	10.0
73.973	10.0
75.342	10.0
76.712	9.00
78.082	9.00
79.452	9.00
80.822	9.00
82.192	8.00
83.562	8.00
84.932	8.00
86.301	7.00
87.671	7.00
89.041	7.00
90.411	7.00
91.781	6.00
93.151	6.00
94.521	4.00
95.890	4.00
97.260	3.00
98.630	1.00

Эмпирическое распределение осадков за ноябрь по метеостанции Hовый Порт приведено на рис.2.

.

Рис.2. Эмпирическое распределение осадков за ноябрь по метеостанции Повый Порт.

По ряду наблюдений сумм осадков за ноябрь на метеостанции Новый Порт за период 1924-2003 гг. были получены следующие значения параметров: среднее ряда = 19.2 мм, Cv=0.6424, Cs = 1.0804, Cs/Cv = 1.6818, r(1)= 0.3667. При объеме выборки 72 года (для автокорреляции n=71 год) и числе степеней свободы 69 (n-2) автокорреляция, как следует из формулы (18) и табл.1, является статистически значимой величиной. По формулам (1)-(12) были вычислены расчетные значения статистик Диксона и Смирнова-Граббса, которые приведены в табл.5. Критические значения были получены из таблиц приложения при имеющемся объеме ряда наблюдений (72 года), вычисленным значения коэффициентов асимметрии и автокорреляции и принятом уровне значимости 5%. Полученные критические значения статистик критериев также приведены в табл.5.

Таблица 5

Результаты оценки однородности эмпирического распределения сумм осадков за ноябрь на метеостанции Новый Порт по критериям Диксона и Смирнова-Граббса.

Экстремум	Критерий	Расчетн.знач.	Критич.знач.	Уровень значимости расч.	Вывод
max	Диксон 1	0.0893	0.3268	11.0000	однороден
max	Диксон 2	0.0926	0.3340	11.0000	однороден
max	Диксон 3	0.2037	0.4102	11.0000	однороден
max	Диксон 4	0.2075	0.4068	11.0000	однороден
max	Диксон 5	0.1964	0.4041	11.0000	однороден
min	Диксон 1	0.0357	0.0629	11.0000	однороден
min	Диксон 2	0.0392	0.0729	11.0000	однороден
min	Диксон 3	0.0588	0.0913	11.0000	Однороден
min	Диксон 4	0.0667	0.1006	11.0000	Однороден
min	Диксон 5	0.0536	0.0824	11.0000	Однороден
max	Cмирнов-Граббс	3.0692	4.3010	11.0000	Однороден
min	Cмирнов-Граббс	1.4754	1.8608	11.0000	Однороден

В табл.5 приведены названия критериев, расчетные значения статистик, их критические значения и уровень значимости, соответствующий расчетным значениям, который также определяется интерполяцией из таблиц критических значений Приложения. В данном случае при проверке экстремумов уровень значимости составляет более 10% и поэтому в таблице приведено значение 11%, свидетельствующее о том, что гипотеза об однородности принимается с высокой степенью достоверности. В остальных случаях, если уровень значимости, соответствующий расчетному значению статистки, находится в интервале от 1% до 10% - приводится его численное значение. Это позволяет дифференцировать вывод о принятии или отклонении нулевой гипотезы об однородности. Так, если уровень значимости близок к 5% (находится в диапазоне от 4% до 5%), то гипотеза принимается, но с учетом этого обстоятельства и вывод помещается в скобки: (+).

Для оценки стационарности дисперсий и средних значений многолетний ряд наблюдений вначале был разбит на две равные части или подвыборки: первая с 1924 по 1964 гг., вторая с 1965 по 2003 гг. Для каждой подвыборки были рассчитаны средние значения и дисперсии, которые соответственно равны: Ycp_I.=13.7 мм, σ ²_I. =115.98 мм² и Ycp_II.=24.6 мм, σ ²_II. =130.75 мм². По полученным дисперсиям и средним значениям были определены расчетные значения статистик Фишера и Стьюдента (формулы (19) и (22)), которые приведены в табл.6. На основе вычисленных значений новых степеней свободы (формулы (20) и (21), таблицы 2 и 3) и таблиц приложения были определены критические значения статистик Фишера и Стьюдента, которые также даны в табл.6. Из результатов табл.6 следует, что расчетная статистика критерия Фишера меньше критической, поэтому вывод о стационарности дисперсий принимается. Вместе с тем, расчетное значение статистики Стьюдента значительно превышает его критическое значение и гипотеза о стационарности средних значений двух частей временного ряда осадков отклоняется с высокой степенью достоверности.

Таблица 6.

Оценка стационарности средний значений и дисперсий ряда наблюдений сумм осадков за ноябрь на метеостанции Новый Порт.

Критерий	Расчетное зн.	Критич.зн.	Уровень значимости расч.	Вывод
Критерий Фишера	1.1274	1.9141	11.0000	однороден
Критерий Стъюдента	4.1116	2.8580	0.9000	неоднороден

Вторым способом временной ряд осадков, приведенный на рис.1, был разделен по дате разрыва данных наблюдений, которые закончились в 1955 г. и были возобновлены в 1961 г. Поэтому ряд наблюдений был разбит на следующие две части: 1924 -1955 гг. и 1961 – 2003 г. Для каждой подвыборки были рассчитаны средние значения и дисперсии, которые соответственно равны: Ycp_I.=10.7 мм, σ ²_I. =25.6 мм² и Ycp_II.=26.0 мм, σ ²_II. =149.7 мм². Как следует из сравнения параметров, полученных первым и вторым способами, во втором случае как существенно отличаются средние значения (в 2.6 раза), так уже и дисперсии (в 5.8 раза). Полученные расчетные и критические значения статистик критериев Фишера и Стьюдента при втором способе разбиения ряда приведены в Табл.7 и свидетельствуют о существенной нестационарности как средних значений, так и дисперсий.

Таблица 7.

Оценка стационарности средний значений и дисперсий ряда наблюдений сумм осадков за ноябрь на метеостанции Новый Порт при разбиении ряда по датам производства наблюдений.

Критерий	Расчетное зн.	Критич.зн.	Уровень значимости расч.	Вывод
Критерий Фишера	5.9457	1.9494	0.9000	неоднороден
Критерий Стъюдента	7.2174	2.8727	0.9000	неоднороден

Оценка стационарности осадков были проведена и для остальных месяцев года и результаты анализа помимо индивидуальных таблиц расчетов для каждого месяца (Табл.5- 7) были представлены в виде одной обобщенной таблицы 8. В таблице 8 знаком «+» представлен вывод о принятии гипотезы однородности и стационарности, а знаком «-» - об отклонении. Если расчетное значение статистики критерия близко к критическому, но превышает его, то в таблице также приводится уровень значимости, соответствующий расчетному значению критерия. В этом случае вывод о принятии или отклонении гипотезы однородности и стационарности является сомнительным и заключается в скобки.

Таблица 8.

Результаты оценки стационарности месячных сумм осадков по критериям

Диксона (Д) , Фишера (Ф) и Стьюдента (Ст) на метеостанции Новый Порт

Д	Ф	Ст	Д	Ф	Ст	Д	Ф	Ст
январь	февраль	Март
+	+	-	+	+	+	+	-	-
апрель	Май	Июнь
+	(-)3.0	-	+	-	+	+	+	-
июль	август	Сентябрь
+	+	+	+	+	+	+	-	+
октябрь	Ноябрь	Декабрь
+	+	-	+	+	-	+	(+)4.4	-

Анализ табл.8 позволяет сделать выводы, что все распределения сумм осадков по месяцам являются однородными. Установлена нестационарность дисперсий и особенно средних в холодные месяцы года. Полученная нестационарность вполне может быть обусловлена сменой приборов для регистрации осадков как раз в 1960е годы, причем предшествующий осадкомер не позволял нейтрализовать ветровое выдувание осадков, что особенно характерно для твердых осадков и для станций, находятся на морском побережье в зоне больших ветров.