БИОТЕХНИЧЕСКИЙ ЗАКОН И АДЕКВАТНОСТЬ ГОТОВОЙ МОДЕЛИ

По остаткам размах вариации (изменчивости) R по совокупности остатков равен
. (1)
Среднее линейное отклонение составляет . (2)
Дисперсия остатков с учетом малой выборки, когда ,будет равной
. (3)
При впереди правой части уравнения (3) принимается отношение 1 / n . Среднее квадратичное отклонение остатков равно
. (4)
Коэффициент изменчивости (вариации) , дающий относительную оценку вариабельности остатков, определяется из выражения
. (5)
Нормированное отклонение наблюдения определяется по формуле
. (6)
Принимается значение доверительного интервала -t - сигмальная единица. Например, пусть класс надежности расчетов принят B, то есть расчеты необходимы с малым риском. Поэтому [t]= 1,96. Для решения вопроса об исключения наблюдения и повторном поиске значений параметров модели по рекомендациям [2, с.260] необходимо выяснить сущность отклонения (так называемые «особые» точки, которые нами рассмотрены отдельно). Формально, если выполнено соотношение
, (7)
то такое наблюдение оставляется в статистической совокупности. При невыполнении условия (7) наблюдение исключается и поиск проводится по новому сеансу на ПЭВМ.
Асимметрия (скошенность) A S распределения остатков вычисляется по формуле
.(8)
Если A S > 0 , то наблюдается правосторонняя асимметрия (вершина распределения расположена слева среднего значения), при условии A S < 0 - левосторонняя асимметрия. Чаще всего асимметрия выражает нeнopмaльность распределения или наличие нескольких, сдвинутых вместе, распределений.
Эксцесс (островершинность) E будет равен
.(9)
При островершинности E > 0 , то эксцесс положительный. Если E < 0 , то распределение остатков плосковершинное. По значениям t и σ определяются параметры функции кривой распределения. Уравнение нормированной (стандартной) плотности нормального распределения будет иметь вид [3, с.41]
. (10)
Другие виды распределения, как (хи - квадрат), Стьюдента ( t - распределение), биномиальное, логарифмическое, логнормальное и другие в технических исследованиях встречаются редко. В исследованиях могут встречаться известные устойчивые законы распределения, а также их комбинации.
Предельная ошибка выборки остатков вычисляется по формуле
,(11)
из которой вычисляется допустимое значение n по формуле
, (12)
где допустимое значение t принимается по известным данным. Полученное значение n сравнивается со значением расчетного числа наблюдений.
Такие сложные расчеты необходимы при аппроксимации. В них нет необходимости при процеДypaх идентификации. Потому при поиске параметров модели можно воспользоваться еще более простыми формулами:
- ошибка среднеквадратичного отклонения ;(13)
- показатель точности наблюдений ;(1)
- доверительные интервалы ; .(15)
Адекватность модели. Для вычислений критериев адекватности необходимо определить:
- дисперсию теоретических значений показателя ,(16)
где
; (17)
- дисперсию фактических значений показателя
; (18)
- дисперсию остатков
. (19)
Значение дисперсии остатков можно также вычислять по правилу сложения дисперсией, поэтому
. (20)
Для измерения тесноты связи принимается коэффициент детерминации B yx , равный [3, с.102]
.(21)
Чем ближе , то тем больше совпадают теоретические и фактические выходные результаты. Мера неопределенности регрессии, то есть влияние «необъясненной» дисперсии , определится из выражения
. (22)
Исправленный коэффициент детерминации B * yx определится из формулы [3, с.108]
. (23)
С вероятностью 1 - α , где α - уровень значимости по табл. 2, коэффициент детерминации проверяется по критерию Фишера [3, с.202]
. (24)
Существенность дисперсионной связи вычисляется по критерию Фишера [25, с.33] по формуле
, (25)
где - табличное значение критерия Фишера, n - число наблюдений, N - число параметров модели. Регрессионная модель считается работоспособной, если [4, с.34]. По критерию Фишера сравниваются несколько альтернативных регрессионных моделей, если таковые возникнут при обсуждении проблемы построения конструкций моделей (способами аппроксимации). Например, лучше обычная или нормализованная формы модели, предпочтительнее арифметические или алгебраические оперативные константы.
Автоматический поиск параметров модели. Для оценки сходимости и адекватности модели автоматического поиска параметров модели целесообразно вычислять статистические показатели, соотнесенные к выходным результатам объекта исследования:
- сумма квадратов отклонений
;(26)
- среднеквадратичное отклонение
; (27)
- ошибка среднеквадратичного отклонения
; (28)
- коэффициент изменчивости
; (29)
- показатель точности
;(30)
- относительное отклонение выходных результатов
; (31)
- среднее относительное отклонение выходных результатов
. (32)
Максимальное отклонение Δ max выбирается из значений (31) по абсолютной величине из множества | Δ i | .
Сравнение моделей по критериям адекватности. Все существующие методы анализа адекватности моделей статистическим выборкам предполагают однородность экспериментов, включенных в одно множество. При этом классическая статистика четко предполагает только наличие закона нормального распределения.
Однако полезно сравнить предложенный метод оценки адекватности по максимальной относительной погрешности с известными методами, в частности, с методом Фишера. Он наиболее распространен в сельскохозяйственной и инженерно-экологической науке.
Средняя относительная погрешность по остаткам от модели. Поэтому наиболее простому известному способу суммируются все отношения остатков ξ i по модулю (без учета знаков, то есть | ξ i | ) к фактическим значениям изучаемого показателя . В результате исключается проблема сравнимости моделей и одновременно «проклятие» размерности показателя . Затем вычисляется результат как сумма ∑ и она делится на число n члeнов (наблюдений) статистической выборки.
Тогда получается сопоставимый для любых статистических моделей, полученных аппроксимацией произвольных формул или идентификацией устойчивых законов и закономерностей критерий адекватности - средняя относительная погрешность (%)
. (33)
Однако любая средняя арифметическая величина относится к нормальному закону распределения, поэтому не применима к неоднородным статистическим выборкам, каждый члeн которых имеет собственную индивидуальность. Поэтому по формуле (33) можно получить одно и то же значение относительной погрешности для нескольких выборок, имеющих различные друг от друга законы распределения. При этом размах может быть сильно различающимся. Поэтому лучше всего определять относительную погрешность по формуле (33), а затем выделить . Тогда получается, что требование адекватности по максимальной относительной погрешности в несколько раз (равное отношению ) строже по сравнению с известным методом.
Средняя квадратичная погрешность. Чтобы исключить влияние знаков, было придумано возвести в квадрат остатки, то есть получить выборку в виде ряда ξ 2 .
Но вначале придется определить среднее арифметическое от всех значений изучаемого показателя, то есть применить формулу для нормального распределения случайных событий
. (34)
Среднеарифметического объекта в природе физически не существует. Поэтому эта физическая величина весьма условная, придуманная математиками для сравнения статистических выборок друг с другом. Однако много существует однородных выборок с члeнами, равными по сущности неотличительных признаков, у которых с различным размахом может оказаться одинаковое значение среднего арифметического.
Таким образом, уже изначально ясно, что метод среднеквадратичного отклонения (погрешности) не примелем для выявления адекватности идентификации закономерностей. Однако до сих пор применяют для выявления адекватности аппроксимированных уравнений, не имеющих содержательного смысла по параметрам и отдельным составляющим, формулу (при этом к методу наименьших квадратов не относится)
. (35)
F-критерий Фишера. Этот показатель адекватности относится к дисперсионному анализу. Он оказывает статистическую значимость уравнения регрессии (в нашем случае статистически устойчивой закономерности, построенной на основе устойчивых законов) в целом. Анализ адекватности модели выполняется сравнением фактического F Ф и табличного (критического) F значений критерия Фишера.
При этом фактическое значение критерия вычисляется, с использованием дисперсий на одну степень свободы, по формуле
,(36)
где n - число члeнов (единиц) статистической совокупности, шт., m - количество объясняющих факторов в уравнении (модели), шт. S Ф - факторная дисперсия (по размерности изучаемого показателя), S- остаточная от модели дисперсия (по размерности показателя).
В формуле (36) обращает внимание тот факт, что по выражению в числителе n - m - 1 количество объясняющих переменных должно соответствовать условию .
Многие примеры идентификации показали, что это условие по Фишеру при использовании устойчивых законов необязательное, поэтому количество одновременно учитываемых факторов может быть и больше числа наблюдений (например. 14 факторов при 10 сложных по количеству влияющих факторов наблюдениях).
При этом было замечено, что конструктивная связь между факторами требует наименьшего количества параметров модели. С переходом на мультиколлениарные связи потребуется большее количество параметров модели, а при аддитивной связи будет наблюдаться наибольшее количестве параметров у искомой статистической закономерности.
Факторная дисперсия, то есть сумма квадратов отклонений фактических значений показателя от их среднеарифметического по формуле (7.25), вычисляется по уравнению
, (37)
где y - расчетное по модели значение изучаемого показателя, - среднее арифметическое у фактических значений показателя. Остаточная дисперсия, то есть сумма квадратов отклонений по остаткам , определяется по выражению
.(38)
Общая вариация (дисперсия) результата определяется по формуле
. (39)
Примеры сравнения. Для сопоставления различных критериев адекватности используем табличные данные по трем гидрометрическим створам и пяти показателям свойств у 18 проб травы [5]. Табличное критическое значение критерия Фишера равно 3,68 для всех пяти двухфакторных моделей при степенях свободы для расчета фактического критерия Фишера k 1 = n - m - 1 = 18 - 2 - 1 = 15 и при этом для двух учтенных факторов k 2 = m = 2 .
Удельная масса сырой травы. В табл. 1 приведены результаты расчетов различных критериев адекватности одной и той же модели. Предлагаемый нами критерий адекватности приведен в первом столбце. Он не требует никаких дополнительных расчетов и прост в применении тем, что нужно просмотреть относительную погрешность по всем строкам и выбрать максимальную из значений погрешности, обозначив как Δ max (%).
Поверхность отклика на рис. 1 сложна и она в среднем хаpaктеризуется статистической моделью вида
. (39)
На площадке № 7 показатель по формуле (39) получил наибольшую погрешность Δ max = 54,24% . Относительную погрешность Δ = 38,21% получила травяная проба на пробной площадке № 15. Средняя относительная погрешность (%) требует небольших вычислений. Вначале определяется сумма погрешностей по модулю (второй столбец табл. 1), а затем вычисляется средняя погрешность.
Рис. 1. Поверхность отклика массы срезанной травы в среднем на трех створах реки Ировка
Таким образом, формула (39) получила по предложенному нами критерию адекватности Δ max для седьмой пробной площадки и при этом второй статистический показатель . Средняя квадратичная погрешность равна . Тогда получается, что максимальная относительная погрешность больше средних величин соответственно в 54,24 / 21,14 ≈ 2,6 и 54,24 / 23,80 ≈ 2,3 раза.
Распределение погрешностей только частично подчиняется нормальному закону Гаусса-Лапласа. Однако биотехнический закон адекватнее. Например, нормальное и ранговое распределения значений из второго столбца табл. 1 определяются формулами:
- закон нормального распределения погрешностей по модулю
;(40)
- распределение по биотехническому закону проф. П.М. Мазуркина
. (41)
Сравнение показывает, что нормальное распределение по закону Гаусса-Лапласа при условии (здесь этот закон приводится напрямую по параметрам) относительно оси ординат биотехнический закон (закон гибели как его частный случай) имеет более высокое значение коэффициента корреляции 0,9962 > 0,9425 (рис. 2).
Таблица 1. Критерии адекватности закономерности (39) массы проб сырой травы
|
№ п/п |
Δ max, % |
, % |
Среднее квадратичное |
F -критерий Фишера F Ф =3,20<3,68 |
||
|
|
, % |
S Ф |
S |
|||
|
1 |
16,91 |
16,91 |
880 |
22141,4 |
48796,8 |
22141,4 |
|
2 |
8,76 |
8,76 |
580 |
2580,6 |
14520,3 |
2580,6 |
|
3 |
30,26 |
30,26 |
420 |
16154,4 |
1354,2 |
16154,4 |
|
4 |
1,45 |
1,45 |
475 |
47,6 |
806,6 |
47,6 |
|
5 |
13,92 |
13,92 |
600 |
6972,3 |
38,4 |
6972,3 |
|
6 |
21,70 |
21,70 |
440 |
9120,3 |
635,0 |
9120,3 |
|
7 |
54,24 |
54,24 |
340 |
34003,4 |
198,8 |
34003,4 |
|
8 |
37,08 |
37,08 |
390 |
20909,2 |
590,5 |
20909,2 |
|
9 |
29,60 |
29,60 |
780 |
53314,8 |
1505,4 |
53314,8 |
|
10 |
29,05 |
29,05 |
630 |
33489,0 |
4006,9 |
33489,0 |
|
11 |
19,04 |
19,04 |
490 |
8704,9 |
12905,0 |
8704,9 |
|
12 |
6,52 |
6,52 |
315 |
506,3 |
20391,8 |
2756,3 |
13 |
11,42 |
11,42 |
565 |
4160,3 |
96,0 |
4160,3 |
|
14 |
10,91 |
10,91 |
460 |
2520,0 |
0,0 |
2520,0 |
|
15 |
38,21 |
38,21 |
380 |
21083,0 |
222,0 |
21083,0 |
|
16 |
25,28 |
25,28 |
595 |
22620,2 |
4316,5 |
22620,2 |
|
17 |
8,00 |
8,00 |
415 |
1102,2 |
3856,4 |
1102,2 |
|
18 |
18,26 |
18,26 |
430 |
6162,3 |
3,2 |
6162,3 |
|
Сумма |
380,61 |
9185 |
265592,0 |
114244,0 |
267842,1 |
|
|
Критерий |
21,15 |
510,3 |
23,80 |
|
|
|
Вроде бы асимметрия небольшая, но она имеет колоссальное креативное значение. Поэтому формула (39) может быть применена, то есть, не отброшена как статистически недостоверное, только из-за того, что её конструкция идентична конструкции моделей других показателей у свойств луговой травы. Однако уравнение с конкретными параметрами модели (39) не может быть использованы в различных технологических или иных расчетах.
|
|
|
|
а) |
б) |
Рис. 2. Графики нормального (а) и биотехнического (б) распределения ошибок
Результаты сравнительных расчетов приведены в табл. 2.
Полужирным шрифтом выделены наилучшие значения того или иного критериев адекватности. Заметим здесь, что эмпирическим путем в эконометрике достигнуто требование к средней погрешности эконометрической модели в 8,0%. Этот уровень вполне совпадает с погрешностью 8,81% табл. 2 для относительной влажности.
Сравнение показало, что наипростейшим критерием адекватности и более точным по сущности идентифицируемых природных и природно-антропогенных явлений и процессов является статистический показатель Δ max - максимальная относительная погрешность (в процентах).
В заключение приведем слова Н. Винера [6, с.76]: «... современный аппарат малых выборок, как только он выходит за рамки простого подсчета своих собственных, специально определенных параметров и превращается в метод положительных статистических выводов для новых случаев, уже не внушает мне доверия, Исключение составляет случай, когда этот аппарат применяется статистиком, который явно знает или хотя бы неявно чувствует основные элементы динамики исследуемой ситуации».
Таблица 2. Критерии адекватности по моделям изменения показателей луговой травы
|
Параметр травы |
Δ max, % |
, % |
|
, % |
F Ф |
R 2 |
R |
|
m , г/м2 |
54.24 |
21.15 |
510.3 |
23.80 |
3.20 |
0.299 |
0.547 |
|
m c, г/м2 |
51.31 |
19.22 |
138.6 |
20.62 |
3.82 |
0.461 |
0.679 |
|
m в0, г/м2 |
59.87 |
22.95 |
364.2 |
26.26 |
3.12 |
0.294 |
0.542 |
|
W , % |
26.21 |
8.81 |
250.1 |
10.32 |
4.59 |
0.379 |
0.616 |
|
, г/(м2 ч) |
46.66 |
18.77 |
1.312 |
21.44 |
1.87 |
0.200 |
0.447 |
Статья опубликована при поддержке гранта 3.2.3/4603 МОН РФ
СПИСОК ЛИТЕРАТУРЫ:
- Мазуркин, П.М. Статистическое моделирование. Эвристико-математический подход / П.М. Мазуркин. - Научное издание. - Йошкар-Ола: МарГТУ, 2001. - 100с.
- Пасхавер, И.С. Общая теория статистики. Для программированного обучения /
И.С. Пасхавер, А.Л. Яблочник. М.: Финансы и статистика, 1983. - 432 с. - Фёрстер, Э. Методы корреляционного и регрессионного анализа: Руководство для экономистов / Э. Фёрстер, Б. Рёнц. - М.: Экономика и статистика, 1983. - 302 с.
- Елисеева, И.И. Логика прикладного статистического анализа / И.И. Елисеева, В.О. Рукавишников. - М.: Финансы и статистика, 1982. - 192 с.
- Мазуркин, П.М. Измерение продуктивности травяного покрова пойменного луга /
П.М. Мазуркин, С.И. Михайлова // Современные наукоемкие технологии: материалы заочной электронной конференции. - 2008. - № 7. - С.91-92. - Винер, Н. Кибернетика или управление и связь в животном и машине / Н. Винер. - 2-е изд. - М.: Наука, 1983. - 344 с.
Статья в формате PDF
119 KB...
11 04 2026 20:43:24
09 04 2026 15:17:29
Статья в формате PDF
113 KB...
08 04 2026 15:46:17
Летом 2012 года был проведен мониторинг расхода воды на малом водотоке. Мерный сосуд был принят в виде ковша емкостью один литр. Все измерения проводились вечером с 17-00 часов. Поэтому текущее время берется целыми сутками. Модель динамики имеет две составляющие: первая составляющая является законом экспоненциального роста, а вторая волновым возмущением с переменными амплитудой и частотой колебания. Показана методика моделирования с процеДypaми: 1) выявление постоянного члeна; 2) по остаткам от постоянного члeна, последовательно усложняя конструкцию, идентифицируется волновая функция; 3) постоянный члeн совмещается с волновой функцией; 4) усложняется конструкция тренда до устойчивого не волнового закона.
...
07 04 2026 7:24:51
Для исследования вариаций параметров живых существ, обитающих в биосфере в разных широтных регионах, в частности экваториальных, построена модель экваториального электроджета, основанная на численном решении дифференциальных уравнений второй степени для потенциала, вызванного прострaнcтвенным зарядом.
...
06 04 2026 17:49:46
В статье отражен анализ работы действующей системы предварительного охлаждения воздуха (СПОВ) воздухоразделительной установки (ВРУ) КААр-30М; выявлены основные проблемы действующей СПОВ; предложены технологические решения, которые способствуют преодолению существующих недостатков и проведен сравнительный анализ действующей и модернизированной систем по основным показателям.
...
05 04 2026 14:16:36
Статья в формате PDF
105 KB...
04 04 2026 19:10:56
Статья в формате PDF
147 KB...
03 04 2026 6:58:18
Статья в формате PDF
106 KB...
02 04 2026 1:14:45
В экспериментах по микроэволюции генетически модифицированных бактерий (ГМО) при непрерывном культивировании показано, что при переходе от одного стационарного состояния к другому в открытой биологической системе скорость производства энтропии должна возрастать, а не уменьшаться, как следует из основных положений неравновесной термодинамики. С точки зрения термодинамики проточные культуры микроорганизмов – хемостат и турбидостат – это открытые термодинамические системы, способные находиться в устойчивых стационарных состояниях. Причем, в соответствии с классификацией М.Эйгена (1973), хемостат соответствует случаю постоянных потоков, а турбидостат – случаю постоянной организации. Несмотря на кажущееся разнообразие микроэволюционных переходов в двух типах открытых систем при их изучении обнаруживаются общие закономерности. Важнейшей из них является возрастание потока использованной популяциями свободной энергии, и, следовательно, возрастание теплорассеяния и скорости производства энтропии. Результаты свидетельствуют о необходимости дальнейшего развития термодинамической теории открытых биологических систем, дальнейшего изучения общих закономерностей биологического развития.
...
01 04 2026 11:14:11
Статья в формате PDF
269 KB...
31 03 2026 0:53:41
Статья в формате PDF
126 KB...
30 03 2026 14:12:45
Статья в формате PDF
112 KB...
29 03 2026 11:52:22
Статья в формате PDF
126 KB...
28 03 2026 4:52:24
Статья в формате PDF
321 KB...
26 03 2026 17:44:48
Статья в формате PDF
105 KB...
25 03 2026 16:54:20
Бесплодие является одним из главным заболеванием коров. По причине непригодности к воспроизводству из стад выбывает более половины животных. По этой причине сельскохозяйственные предприятия терпят существенные убытки. В настоящее время в производстве требуются современные методы лечения, которые отличались бы высокой эффективностью, широким спектром действия, низкозатратностью. Авторы считают, что такой инновационной технологией является лечение нарушения воспроизводительной системы коров и телок с использованием метода криотерапии и озонированными гомеопатическими препаратами.
...
24 03 2026 4:52:59
Статья в формате PDF
112 KB...
23 03 2026 13:31:16
Показана возможность использования электрохимически активированной воды (в виде анолита и католита) для повышения урожайности зерновых и овощных (картофеля) культур и улучшения фитосанитарной ситуации с помощью модуля активации оросительной воды. Наиболее энтомоцидным действием в отношении пшеничного трипса обладал анолит с окислительно-восстановительным потенциалом +600 и +900 мВ. Католит с ОВП – 700 мВ способствовал увеличению всхожести до 96%. Хороший результат в борьбе против колорадского жука давала предпосевная обработка клубней картофеля вначале анолитом, а потом католитом. Заселенность кустов колорадским жуком и проволочником снизилась на 37–83%. Наиболее эффективно в плане оптимизации фитосанитарного состояния посевов сочетание предпосевной обработки семян с последующим опрыскиванием стeблестоя католитом или анолитом. ...
22 03 2026 20:49:15
Статья в формате PDF
108 KB...
21 03 2026 13:55:13
20 03 2026 13:41:32
В статье излагаются положения новой концепции на субстрат миндалевидного комплекса, предлагающей рассматривать эту структуру лимбической системы как ядерно-палеокортикальный компонент мозга.
...
19 03 2026 23:21:28
Статья в формате PDF
133 KB...
18 03 2026 18:26:26
Известные значения констант диссоциации одного из самых распространенных природных флавоноидов – кверцетина – отличаются крайней невоспроизводимостью. Одной из причин этого следует считать легкое окисление кверцетина в процессе титрования кислородом воздуха. Для устранения этого эффекта предложен модифицированный вариант потенциометрического титрования с барботированием инертного газа (азот) через титруемый раствор с добавкой в него неионогенного детергента. Полученное таким способом значение pKaI кверцетина равно 6.62 ± 0.04. Из этого следует принципиально важный вывод: в нейтральной среде (при рН ~ 7) кверцетин и, возможно, другие флавонолы, пpaктически полностью диссоциированы.
...
17 03 2026 9:10:57
Статья в формате PDF
143 KB...
16 03 2026 2:51:16
Статья в формате PDF
253 KB...
15 03 2026 11:48:16
Статья в формате PDF
120 KB...
14 03 2026 18:21:13
Распространённость мастопатии в популяции может достигать более 70 % и не зависит от этнического фенотипа. 92,5 % пациенток, самостоятельно обратившихся по поводу мастопатии, – это городские жители из социальной категории «служащие» со средним специальным и высшим гуманитарным образованием. Сопутствующие заболевания органов пищеварения и урогeнитaльной системы, а также девиантные психологические черты личности достоверно чаще регистрируются у женщин с мастопатией, чем в контроле. Более 70 % женщин отмечают усиление симптомов мастопатии после обострения соматических заболеваний и нервных стрессов, а более 80 % испытывают психологический дискомфорт от направления в онкодиспансер.
Необходимы специализированные маммологические кабинеты при женских консультациях и поликлиниках для квалифицированной диагностики, лечения и психологической коррекции пациенток с доброкачественными заболеваниями молочных желез.
...
13 03 2026 16:37:18
Статья в формате PDF
110 KB...
12 03 2026 20:23:25
Статья в формате PDF
308 KB...
11 03 2026 13:45:57
Статья в формате PDF
134 KB...
09 03 2026 1:39:37
08 03 2026 18:59:50
Статья в формате PDF
110 KB...
07 03 2026 23:29:16
Статья в формате PDF
245 KB...
06 03 2026 17:40:38
С помощью геоинформационной системы были получены точные измеренные значения каждого годичного слоя на всем керне древесины сосны. Данные обработаны в математической среде и получена статистическая формула, которая состоит из 16 составляющих, что позволило дать ориентировочный долгосрочный прогноз.
...
05 03 2026 18:10:45
Статья в формате PDF
307 KB...
04 03 2026 21:32:28
Еще:
Поддержать себя -1 :: Поддержать себя -2 :: Поддержать себя -3 :: Поддержать себя -4 :: Поддержать себя -5 :: Поддержать себя -6 :: Поддержать себя -7 :: Поддержать себя -8 :: Поддержать себя -9 :: Поддержать себя -10 :: Поддержать себя -11 :: Поддержать себя -12 :: Поддержать себя -13 :: Поддержать себя -14 :: Поддержать себя -15 :: Поддержать себя -16 :: Поддержать себя -17 :: Поддержать себя -18 :: Поддержать себя -19 :: Поддержать себя -20 :: Поддержать себя -21 :: Поддержать себя -22 :: Поддержать себя -23 :: Поддержать себя -24 :: Поддержать себя -25 :: Поддержать себя -26 :: Поддержать себя -27 :: Поддержать себя -28 :: Поддержать себя -29 :: Поддержать себя -30 :: Поддержать себя -31 :: Поддержать себя -32 :: Поддержать себя -33 :: Поддержать себя -34 :: Поддержать себя -35 :: Поддержать себя -36 :: Поддержать себя -37 :: Поддержать себя -38 ::