Функция ФИШЕР
В этой статье описаны синтаксис формулы и использование функции ФИШЕР в Microsoft Excel.
Описание
Возвращает преобразование Фишера для аргумента x. Это преобразование строит функцию, которая имеет нормальное, а не асимметричное распределение. Данная функция используется для проверки гипотез с помощью коэффициента корреляции.
Синтаксис
ФИШЕР(x)
Аргументы функции ФИШЕР описаны ниже.
Замечания
-
Если x не является числом, фишер возвращает #VALUE! значение ошибки #ЗНАЧ!.
-
Если x ≤ -1 или x ≥ 1, фишер возвращает #NUM! значение ошибки #ЗНАЧ!.
-
Уравнение для преобразования Фишера имеет следующий вид:
Пример
Скопируйте образец данных из следующей таблицы и вставьте их в ячейку A1 нового листа Excel. Чтобы отобразить результаты формул, выделите их и нажмите клавишу F2, а затем — клавишу ВВОД. При необходимости измените ширину столбцов, чтобы видеть все данные.
Формула |
Описание |
Результат |
=ФИШЕР(0,75) |
Преобразование Фишера для аргумента 0,75 |
0,9729551 |
понимаем лучше экономику и финансы
Сергей Кикевич Все авторы
В 1911 году американский финансист Ирвинг Фишер в своей книге «Покупательная сила денег» опубликовал результат своих наблюдений за инфляцией, денежной массой и объемом производства в виде простой но очень емкой формулы:
M V = P Q
М – объем денежной массы
V – скорость оборачиваемости денег
P – уровень цен
Q – объём производства
С тех пор эта формула используется повсеместно для анализа финансовой ситуации и в целях формирования монетарной политики.
Использование формулы Фишера
Уравнение может иметь несколько интерпретаций. Одна из них, возможно, наиболее важная:
При увеличении денежной массы (левая сторона уравнения) возможным результатом может стать как рост объема производства и товаров на рынке, так и рост уровня цен (правая сторона уравнения).
В реальной жизни могут происходить оба явления, но в разных пропорциях. И эти пропорции довольно качественно характеризуют свойства экономики. В развитых странах, где ВВП в основном складывается из промышленных и высоко технологичных продуктов и услуг, рынок спокойно «съедает» очередную порцию денег, напечатанную государством, без значительного роста цен. «Лишние» деньги, которые не может осовоить экономика напрямую, уходят в другие виды долгосрочных активов: акции, облигации, взаимные фонды, пенсионные накопления и т.п.
В других странах, зависящих от природных ресурсов с низкой зависимостью ВВП от реального сектора, наблюдается обратное явление. Сколько денег в экономику не вкачивай, на производство это не оказывает ни малейшего влияния.
Уравнение Фишера и ситуация с денежной массой в России
Становится довольно понятно, почему опыт США с их «количественным смягчением» и колоссальным вливанием новых денег в экономику не получается применить в России.
Со времен Кудрина Минфин и ЦБ взяли обязательство контролировать инфляцию доступными для них методами, в т.ч. ограничивая денежную массу. До этого, кажется, в нашей стране никто уравнение Фишера не изучал. Тем не менее даже сейчас регулярно слышатся призывы (в основном от политиков левого толка) вроде «давайте напечатаем денег и оживим экономику». Разобравшись в смысле уравнения Фишера, становится очевидным, почему простое печатание денег никогда не приводит к желаемым эффектам.
Репетитор оценщика — Формула Фишера. Номинальная и реальная ставки
ФОРМУЛА ФИШЕРА. Перевод номинальной ставки в реальную и наоборот.
В процессе оценки необходимо учитывать, что номинальные и реальные (то есть, включающие и не включающие инфляционный компонент) безрисковые ставки.
Номинальная ставка процента — это рыночная процентная ставка без учета инфляции, отражающая текущую оценку денежных активов.
Реальная ставка процента — это рыночная процентная ставка с учетом инфляции
При пересчете номинальной ставки в реальную и наоборот, целесообразно использовать формулу американского экономиста Фишера, выведенную им еще в 30-е годы:
Rн = Rр + Jинф + Rр * Jинф
Rр = (Rн – Jинф) / (1+ Jинф)
где: Rн — номинальная ставка;
Rр — реальная ставка;
Jинф — годовые темпы прироста инфляции.
Важно отметить, что при использовании номинальных потоков доходов коэффициент капитализации (и ее составные части) должны быть рассчитаны в номинальном выражении, а при реальных потоках доходов — реальном. Для преобразования номинальных потоков доходов в реальные нужно номинальную величину разделить на соответствующий индекс цен, то есть выраженное в процентах отношение уровня цен за тот год, в котором возникнут денежные потоки к уровню цен базового периода.
Например:
Объект недвижимости, сданный на условиях чистой аренды, будет приносить по 1000 долл. ежегодно в течение 2-х лет. Индекс цен в текущем периоде равен 140% и ожидается, что в следующем году он составит 156,7%, а через год 178,5%. Для преобразования номинальных величин в реальные, их необходимо выразить в ценах базисного года. Построим базисный индекс цен для каждого из трех лет. Индексы цен текущего года равны 140/140 = 1, для прогнозного периода: первый год — 156,7/140 = 1,119; второй год — 178,5/140 = 1,275.
Таким образом, реальная величина номинальной 1000 долл., которая будет получена в первом прогнозном году, равна 1000 долл./1,119 = 893,65 долл., во 2-м году (1000 долл./1,275) = 784,31 долл.).
Таким образом, в результате инфляционной корректировки происходит приведение ретроспективной информации, используемой в оценке, к сопоставимому виду, а также учет инфляционного роста цен при составлении прогнозов денежных потоков.
Реактив Фишера для определения количества воды в различных соединениях
Как приобрести реактив Фишера в нашей компании – подробная информация…
Для определения воды как примеси в органических соединениях предложено много способов. Одним из наиболее употребительных является титрование реактивом Фишера.
Метод был разработан в 1935 году немецким химиком Карлом Фишером. Фундаментальный принцип «определения воды по Фишеру» основан на реакции Бунзена между йодом, диоксидом серы и водой. К. Фишер обнаружил, что использование избытка диоксида серы, метанола в качестве растворителя и пиридина как основания (для буферизации раствора) можно использовать для определения воды в неполярных растворителях.
Данный метод — один из наиболее чувствительных и позволяет определять даже небольшое количество воды в органических жидкостях, например, в нефтепродуктах. При этом реактив одновременно выступает и как индикатор.
Наиболее важное преимущество метода титрования по Карлу Фишеру над термическими (потеря веса при прокаливании) — его специфичность для воды. Потери при прокаливании показывают суммарное содержание всех летучих компонентов.
Используемый для титрования реактив Фишера представляет собой раствор иода, двуокиси серы и пиридина, чаще всего в метаноле. Иногда метанол заменяется на метилцеллозольв, диоксан или ледяную уксусную кислоту.
Применяют реактив Фишера и для определения воды, как гигроскопической, так и кристаллизационной в неорганических веществах, хотя мешающие соединения здесь встречаются чаще, чем при анализе органических веществ. Мешают определению сильные окислители и восстановители, которые реагируют с иодом или иодидом.
Реактив Фишера нельзя применять для определения воды в кетонах или альдегидах или в присутствии небольших количеств этих веществ в анализируемых растворителях. Для удаления альдегидов и кетонов прибавляют 2%-ный пиридиновый раствор цианистого водорода, под действием которого альдегиды и кетоны превращаются в циангидрины.
Было предложено использовать реактив Фишера для определения оксикислот, ангидридов кислот, карбонильных производных и гидратационной воды в солях.
В настоящее время используются два варианта метода: кулонометрический и волюметрический (объемный)
Волюметрическое титрование сводится к добавлению йодсодержащего титранта. Для протекания реакции йода с водой также необходимы спирт (или аналог) и органическое основание, которые могут содержаться как в титранте, так и в реактиве для титровальной ячейки.
В зависимости от состава реактивы Фишера бывают однокомпонентные и двухкомпонентные.
Однокомпонентный реактив Фишера (титрант) содержит в себе все реагирующие вещества: йод, имидазол, спирт, двуокись серы. В этом случае в титровальной ячейки необходима лишь среда для протекания реакции: спирт, хлороформ и т.п.
В двухкомпонентных реактивах Фишера — взаимодействующие вещества распределение между титрантом и растворителем. Первый реактив (титрант) — это спиртовой раствор йода, а второй (растворитель) — это спиртовой раствор SO2 с имидазолом или пиридином.
Однокомпонентные реактивы Фишера не устойчивы и за год хранения в запечатанной бутыли титр может измениться на 0,5 мг. Двухкомпонентные реактивы Фишера устойчивы при хранении и дают большую скорость титрования из-за избыточного содержания SO2 в растворителе.
Популярность титрования по методу Карла Фишера обусловлена следующими преимуществами:
— высокая точность и воспроизводимость;
— селективность по воде;
— малые количества необходимых образцов;
— несложная методика подготовки проб;
— малое время анализа;
— практически неограниченный диапазон измерения;
— пригоден для анализа твердых веществ, жидкостей, газов;
— независимость от наличия других летучих веществ.
Важно помнить:
1. Реактив Фишера неустойчив к действию света и влаги. Его стандартизация должна проводиться перед каждым использованием.
2. Рабочий интервал pH для определения воды по Фишеру между 5 и 8, в противном случае высоко кислотные или основные соединения рекомендуется буферизовать.
3. Плохо растворимые в метаноле соединения (например, жиры, углеводороды) следует растворять в высших спиртах или хлороформе, возможны также добавки формамида (для полярных веществ).
4. Следует избегать титрования веществ, которые способны реагировать с компонентами реактива Фишера (например, альдегиды и кетоны, сильные кислоты и основания, окислители и восстановители, соединения, реагирующие с компонентами реактива Фишера с образованием воды).
Критерий Фишера и критерий Стьюдента в эконометрике
С помощью критерия Фишера оценивают качество регрессионной модели в целом и по параметрам.
Для этого выполняется сравнение полученного значения F и табличного F значения. F-критерия Фишера. F фактический определяется из отношения значений факторной и остаточной дисперсий, рассчитанных на одну степень свободы:
где n — число наблюдений;
m — число параметров при факторе х.
F табличный — это максимальное значение критерия под влиянием случайных факторов при текущих степенях свободы и уровне значимости а.
Уровень значимости а — вероятность не принять гипотезу при условии, что она верна. Как правило а принимается равной 0,05 или 0,01.
Если Fтабл > Fфакт то признается статистическая незначимость модели, ненадежность уравнения регрессии.
Таблицы по нахождению критерия Фишера и Стьюдента
Таблицы значений F-критерия Фишера и t-критерия Стьюдента Вы можете посмотреть здесь.
Табличное значение критерия Фишера вычисляют следующим образом:
- Определяют k1, которое равно количеству факторов (Х). Например, в однофакторной модели (модели парной регрессии) k1=1, в двухфакторной k=2.
- Определяют k2, которое определяется по формуле n — m — 1, где n — число наблюдений, m — количество факторов. Например, в однофакторной модели k2 = n — 2.
- На пересечении столбца k1 и строки k2 находят значение критерия Фишера
Для нахождения табличного значения критерия Стьюдента определяют число степеней свободы, которое определяется по формуле n — m — 1 и находят его значение при определенном уровне значимости (0,10, 0,05, 0,01).
Критерии Стьюдента
Для оценки статистической значимости модели по параметрам рассчитывают t-критерии Стьюдента.
Оценка значимости модели с помощью критерия Стьюдента проводится путем сравнения их значений с величиной случайной ошибки:
Случайные ошибки коэффициентов линейной регрессии и коэффициента корреляции определяются по формулам:
Сравнивая фактическое и табличное значения t-статистики и принимается или отвергается гипотеза о значимости модели по параметрам.
Зависимость между критерием Фишера и значением t-статистики Стьюдента определяется так
Как и в случае с оценкой значимости уравнения модели в целом, модель считается ненадежной если tтабл > tфакт
Видео лекциий по расчету критериев Фишера и Стьюдента
Для более подробного изучения расчетов критериев Фишера и Стьюдента советуем посмотреть это видео
Лекция 1. Критерии и Гипотезы
Лекция 2. Критерии и Гипотезы
Лекция 3. Критерии и Гипотезы
Определение доверительных интервалов
Для построения доверительного интервала определяется предельная ошибка А для обоих показателей:
Формулы для нахождения доверительных интервалов выглядят так
Прогнозное значение у определяется с помощью подстановки в
уравнение регрессии прогнозного значения х. Вычисляется средняя стандартная ошибка прогноза
и находится доверительный интервал
Задача регрессионного анализа в предмете эконометрика состоит в анализе дисперсии изучаемого показателя y:
общая сумма квадратов отклонений (TSS)
сумма квадратов отклонений, обусловленная регрессией (RSS)
остаточная сумма квадратов отклонений (ESS)
Долю дисперсии, обусловленную регрессией, в общей дисперсии показателя у характеризует коэффициент детерминации R, который должен превышать 50% (R2 > 0,5). В контрольных по эконометрике в ВУЗах этот показатель рассчитывается всегда.
Методы статистики
Критерии и методы
Рональд Фишер
Точный критерий Фишера – это критерий, который используется для сравнения двух и более относительных показателей, характеризующих частоту определенного признака, имеющего два значения. Исходные данные для расчета точного критерия Фишера обычно группируются в виде четырехпольной таблицы, но могут быть представлены и многопольной таблицей.
1. История разработки критерия
Впервые критерий был предложен Рональдом Фишером в его книге «Проектирование экспериментов». Это произошло в 1935 году. Сам Фишер утверждал, что на эту мысль его натолкнула Муриэль Бристоль. В начале 1920-х годов Рональд, Муриэль и Уильям Роуч находились в Англии на опытной сельскохозяйственной станции. Муриэль утверждала, что может определить, в какой последовательности наливали в ее чашку чай и молоко. На тот момент проверить правильность ее высказывания не представлялось возможным.
Это дало толчок идее Фишера о «нуль гипотезе». Целью стала не попытка доказать, что Муриэль может определить разницу между по-разному приготовленными чашками чая. Решено было опровергнуть гипотезу, что выбор женщина делает наугад. Было определено, что нуль-гипотезу нельзя ни доказать, ни обосновать. Зато ее можно опровергнуть во время экспериментов.
Было приготовлено 8 чашек. В первые четыре налито молоко сначала, в другие четыре – чай. Чашки были помешаны. Бристоль предложили опробовать чай на вкус и разделить чашки по методу приготовления чая. В результате должно было получиться две группы. История говорит, что эксперимент прошел удачно.
Благодаря тесту Фишера вероятность того, что Бристоль действует интуитивно, была уменьшена до 0.01428. То есть, верно определить чашку можно было в одном случае из 70. Но все же нет возможности свести к нулю шансы того, что мадам определяет случайно. Даже если увеличивать число чашек.
Эта история дала толчок развитию «нуль гипотезы». Тогда же был предложен точный критерий Фишера, суть которого в переборе всех возможных комбинаций зависимой и независимой переменных.
2. Для чего используется точный критерий Фишера?
Точный критерий Фишера в основном применяется для сравнения малых выборок. Этому есть две весомые причины. Во-первых, вычисления критерия довольно громоздки и могут занимать много времени или требовать мощных вычислительных ресурсов. Во-вторых, критерий довольно точен (что нашло отражение даже в его названии), что позволяет его использовать в исследованиях с небольшим числом наблюдений.
Особое место отводится точному критерию Фишера в медицине. Это важный метод обработки медицинских данных, нашедший свое применение во многих научных исследованиях. Благодаря ему можно исследовать взаимосвязь определенных фактора и исхода, сравнивать частоту патологических состояний между разными группами пациентов и т.д.
3. В каких случаях можно использовать точный критерий Фишера?
- Сравниваемые переменные должны быть измерены в номинальной шкале и иметь только два значения, например, артериальное давление в норме или повышено, исход благоприятный или неблагоприятный, послеоперационные осложнения есть или нет.
- Критерий подходит для сравнения очень малых выборок: точный критерий Фишера может применяться для анализа четырехпольных таблиц в случае значений ожидаемого явления менее 10, что является ограничением для применения критерия хи-квадрат Пирсона.
- Точный критерий Фишера бывает односторонним и двусторонним. При одностороннем варианте точно известно, куда отклонится один из показателей. Например, во время исследования сравнивают, сколько пациентов выздоровело по сравнению с группой контроля. Предполагают, что терапия не может ухудшить состояние пациентов, а только либо вылечить, либо нет.
Двусторонний тест является предпочтительным, так как оценивает различия частот по двум направлениям. То есть оценивается верятность как большей, так и меньшей частоты явления в экспериментальной группе по сравнению с контрольной группой.
Аналогом точного критерия Фишера является Критерий хи-квадрат Пирсона, при этом точный критерий Фишера обладает более высокой мощностью, особенно при сравнении малых выборок, в связи с чем в этом случае обладает преимуществом.
4. Как рассчитать точный критерий Фишера?
Допустим, изучается зависимость частоты рождения детей с врожденными пороками развития (ВПР) от курения матери во время беременности. Для этого выбраны две группы беременных женщин, одна из которых — экспериментальная, состоящая из 80 женщин, куривших в первом триместре беременности, а вторая — группа сравнения, включающая 90 женщин, ведущих здоровый образ жизни на протяжении всей беременности. Число случаев ВПР плода в экспериментальной группе составило 10, в группе сравнения — 2.
Вначале составляем четырехпольную таблицу сопряженности:
Исход есть (Наличие ВПР) | Исхода нет (Отсутствие ВПР) | Всего | |
Фактор риска есть (Курящие) | A = 10 | B = 70 | (A + B) = 80 |
Фактор риска отсутствует (Некурящие) | C = 2 | D = 88 | (C + D) = 90 |
Всего | (A + C) = 12 | (B + D) = 158 | (A + B + C + D) = 170 |
Точный критерий Фишера рассчитывается по следующей формуле:
где N — общее число исследуемых в двух группах; ! — факториал, представляющий собой произведение числа на последовательность чисел, каждое из которых меньше предыдущего на 1 (например, 4! = 4 · 3 · 2 · 1)
В результате вычислений находим, что P = 0,0137.
5. Как интерпретировать значение точного критерия Фишера?
Достоинством метода является соответствие полученного критерия точному значению уровня значимости p. То есть, полученное в нашем примере значение 0,0137 и есть уровень значимости различий сравниваемых групп по частоте развития ВПР плода. Необходимо лишь сопоставить данное число с критическим уровнем значимости, обычно принимаемым в медицинских исследованиях за 0,05.
- Если значение точного критерия Фишера больше критического, принимается нулевая гипотеза и делается вывод об отсутствии статистически значимых различий частоты исхода в зависимости от наличия фактора риска.
- Если значение точного критерия Фишера меньше критического, принимается альтернативная гипотеза и делается вывод о наличии статистически значимых различий частоты исхода в зависимости от воздействия фактора риска.
В нашем примере P < 0,05, в связи с чем делаем вывод о наличии прямой взаимосвязи курения и вероятности развития ВПР плода. Частота возникновения врожденной патологии у детей курящих женщин статистически значимо выше, чем у некурящих.
Фишер. Статистический вывод
Рональд Фишер — ученый, снабдивший статистику инструментами, благодаря которым она обрела то огромное значение, которое имеет сегодня. Его основной вклад — статистический вывод, инновационный подход, связанный с понятием вероятности, который дал статистике, состоявшей прежде на службе других дисциплин, необходимый импульс для того, чтобы она стала полноправной наукой. Этому британскому математику и биологу мы обязаны статистическим методом, который применяется в планировании научных экспериментов. Он был ярым сторонником евгеники, зародившейся в первой половине XX века, и в этом контексте его исследования касались также генетики и современной эволюционной теории.
По теме см. также Левин. Статистика для менеджеров с использованием Microsoft Excel
Наука. Величайшие теории: выпуск 47: Возможно да, возможно нет. Фишер. Статистический вывод. — М.: Де Агостини, 2015. — 176 с.
Скачать конспект (краткое содержание) в формате Word или pdf, примеры в формате Excel
Глава 1. Статистика до Фишера
Принято считать, что статистика состоит из двух взаимосвязанных, но самостоятельных разделов: описательной статистики, которая занимается эксплоративным (пробным, разведочным) анализом данных, и статистического вывода, отвечающего за предсказания в условиях неопределенности.
Золотая теорема Бернулли, известная сегодня как просто теорема Бернулли, гласит, что относительная частота события стремится к фиксированному числу (вероятности события) по мере того, как увеличивается количество повторений эксперимента (рис. 1).
Рис. 1. Относительная частота выпадения орла после 100 подбрасываний монеты; откройте файл Excel и поэкспериментируйте, нажимая F9; поскольку график основан на формуле, использующей случайные числа, его вид будет постоянно меняться, но стремление к среднему значению будет неизменным
На рубеже XVIII и XIX веков Пьер-Симон де Лаплас (1749–1827) утверждал, что ситуации, связанные со случайностью, бывают двух типов. В ситуациях первого типа случай проявляется в результатах. Известно, что в урне находятся белые и черные шары; вопрос: какой шар мы вытащим? На основе причин (количество белых и черных шаров в урне) вычислим вероятность результата: вытащим ли мы белый или черный шар? В ситуациях второго типа случай проявляется не в результатах, а в причинах. Мы знаем результат опыта (например, мы вытащили черный шар) и хотим вычислить состав содержимого урны, который нам неизвестен (подробнее см. Пьер Симон Лаплас. Опыт философии теории вероятностей).
Представим себе урну, в которой могут быть два разных набора шаров: 2 белых и 3 черных или 3 белых и 2 черных (рис. 2). Вытаскиваем случайный шар, он оказывается черным; какой набор шаров в урне более вероятен? Теорема Байеса обеспечивает численную оценку этой вероятности (подробнее см. Идеи Байеса для менеджеров). Если предположить a priori, что два состава равновероятны (вероятность каждого равна 0,5), после применения формулы Байеса вероятность первого состава увеличивается до 0,6 из-за извлечения черного шара, тогда как вероятность второго состава уменьшается до 0,4. Вероятности a priori (0,5 и 0,5) корректируются a posteriori (0,6 и 0,4). Для Лапласа, как и для Байеса, эта важная теорема означала возможность обучения через опыт.
Рис. 2. Теорема Байеса в действии. Если мы вытащили черный шар, то согласно теореме Байеса, вероятность a posteriori состава на рисунке слева выше, чем вероятность состава на рисунке справа
В 1835 году ученик Лапласа, Симеон Дени Пуассон исследовал вопросы электоральной математики и юриспруденции и сформулировал «закон больших чисел», который обеспечивал лучшую основу применению исчисления вероятностей к социальным проблемам, объясняя статистическую стабильность в социальных изменениях. Большое количество индивидуумов, действуя в рамках системы, определяют регулярность, которая не зависит от их взаимной координации.
Неполнота генетических теорий Чарльза Дарвина (1809–1882) подтолкнула его кузена Фрэнсиса Гальтона (1822–1911) к попытке разрешить проблемы наследования признаков при помощи математического анализа биологических данных. В книге «Наследственность таланта, ее законы и последствия» (1869) он утверждает: «Точно так же, как методами тщательной и умелой селекции в рамках естественных ограничений удается получить стабильную породу собак или лошадей, обладающих особенными способностями к бегу или к чему-либо еще, представляется вполне возможным произвести высокоталантливую расу людей путем рассчитанных браков в течение нескольких последовательных поколений».
Гальтон считал, что союз двух умных людей приведет к рождению еще более умных детей, точно так же как у двух высоких людей рождаются еще более высокие дети. Однако эксперименты с наследованием, которые он проводил в течение всей жизни, привели к открытию другой статистической закономерности, отличной от ожидаемой. В своей книге «Естественное наследование» (1889) он назвал ее «возвращением к посредственности», а позднее — «регрессией к среднему». Сегодня мы знаем, что это не столько биологическая, сколько исключительно статистическая закономерность: более вероятно, что значения нормальной случайной величины будут ближе к ее среднему, ожидаемому значению (подробнее см. Даниэль Канеман. Думай медленно… решай быстро).
Гальтон обнаружил следующее линейное соотношение:
рост ребенка, см = 85 см + 0,5 * рост родителя, см
Это и была одна из прямых регрессии.
Глава 2. Карл Пирсон и биометрическая школа
Карл Пирсон (1857–1936) родился в Лондоне. Его семья принадлежала к среднему классу, что позволило юноше изучать математику в Кембридже, а после окончания в 1879 году продолжить обучение в университетах Гейдельберга и Берлина.
Одним из первых введенных им понятий было «стандартное отклонение». Далее он придумал коэффициент вариации, определив его как связь между стандартным отклонением и средним в абсолютном значении. Пирсон разработал еще два описательных показателя, коэффициент асимметрии и эксцесс. Пирсон сделал использование абстрактной математики в статистике обязательным и анализировал большие наборы данных (более 1000 объектов). Пирсон первым предупреждал об опасности выявления «ложных корреляций»: две переменные могут коррелировать между собой в отсутствие причинно-следственной связи или общей причины. В 1900 году Пирсон вывел критерий хи-квадрат (χ2) для определения качества подгонки наблюдаемого и теоретического, или ожидаемого, распределения.
В 1901 году Уэлдон и Пирсон, при участии Фрэнсиса Гальтона, основали журнал «Биометрика». В 1914 году Пирсон получил для публикации в «Биометрике» статью, подписанную 24-летним студентом по имени Р.А. Фишер. Статья была посвящена выборочным распределениям. Этот вопрос превращался в важную тему для дальнейшего развития статистического вывода, так как позволял количественно оценить надежность предположений, сделанных на основании репрезентативной выборки, имеющей своей целью узнать характеристики генеральной совокупности, набора объектов, который считался слишком большим, чтобы исследовать его полностью.
Характеристики генеральной совокупности, которые нужно было оценить, стали называть параметрами: генеральное среднее μ, стандартное отклонение генеральной совокупности σ или коэффициент корреляции генеральной совокупности ρ. Значения, вычисленные на основе данных выборки для оценки этих параметров, окрестили статистиками: выборочное среднее, стандартное отклонение выборки S.
Уильям Сили Госсет (1876–1937) по образованию был химиком, а познакомился со статистикой после работы в биометрической лаборатории Пирсона. В 1908 году он издал свою знаменитую статью «Возможная ошибка среднего» под псевдонимом Стьюдент. Причиной такой таинственности послужило то, что дублинская пивоварня Гиннесс, где он работал, не разрешала сотрудникам публиковать результаты исследований, проведенных на производстве. Стремясь контролировать качество производимого пива, Стьюдент собирал небольшие выборки (из соображений экономии). Он обнаружил, что один из типов кривых Пирсона соответствовал распределению, удобному для проведения этих экспериментов в малом масштабе. Так, если Стьюдент хотел оценить среднюю кислотность пива, произведенного на заводе за определенный период, он вычислял среднее из уровней кислотности, измеренных в дюжине бочек, которые составляли выборку. Проблема заключалась в том, что Стьюдент не знал возможную ошибку, возникающую при оценке среднего генеральной совокупности по выборочному среднему, то есть число, необходимое для определения точности статистического вывода и допустимого предела кислотности. Чтобы определить это, Стьюденту нужно было знать распределение выборочного среднего. Было замечено, что если выборка достаточно большая — 30 или более объектов, — распределение выборочного среднего приближается к нормальному (благодаря центральной предельной теореме), но если выборка маленькая, то дело обстоит иначе.
Стьюдент вычислил верное распределение, известное сегодня — после того как в 1925 году его отметил Фишер — как распределение t Стьюдента. На самом деле это целое семейство распределений, зависящих от количества степеней свободы; в общем случае оно более плоское, чем нормальное распределение, с более длинными хвостами, что отражает большую неуверенность в выводах. Эта вероятностная модель оказалась неотъемлемой частью статистических методов наших дней ввиду своей надежности, поэтому она используется не только для выводов на основании малых выборок из нормальной генеральной совокупности (для которой среднее и стандартное отклонение неизвестны), но также и для ненормально распределенных данных. Распределение t оказалось практически нечувствительным к гипотезе нормальности (рис. 3).
Рис. 3. Распределение t Стьюдента (серый цвет) обладает более широкими хвостами по отношению к нормальному (черный цвет)
Глава 3. Математические основы статистического вывода
В 1920-е годы Фишер принял эстафету у поколения статистиков, сформировавшегося вокруг Пирсона. Его статья «О математических основах теоретической статистики» и последовавшие за ней книги «Статистические методы для исследователей» и «Планирование экспериментов» ознаменовали становление статистического вывода как математической дисциплины. В них Фишер изложил критерий значимости, дисперсионный анализ и рандомизацию в качестве основных принципов любой работы ученого-натуралиста с фактами.
Статистический вывод — это набор методов, которые позволяют формулировать суждения об общем (генеральная совокупность) на основании частного (выборка), предоставляя меру уверенности в предсказании, вероятность ошибки.
До Фишера статистика, в которой доминировал титан Карл Пирсон, находилась в следующей ситуации. В описательной статистике, хотя и не существовало явного различия между выборкой и генеральной совокупностью, были известны простые графические представления (столбчатая диаграмма, гистограмма, диаграмма рассеяния) и вычислялись основные показатели центральной тенденции (среднее, медиана, мода), дисперсии (стандартное отклонение, хотя оно было и не единственной мерой), позиции (квартили и перцентили) и формы (асимметрия и эксцесс). Переход от разведывательного анализа данных к математической теории вероятностей проходил через подгонку теоретических распределений — нормальной кривой или кривых Пирсона — к наблюдаемым распределениям частот методом наименьших квадратов или методом моментов. Качество подгонки оценивалось с помощью критерия хи-квадрат. Наконец, статистический вывод мог похвастаться только двумя быстрыми методами: предсказаниями, основанными на анализе регрессии и корреляции, и в особенности обратными вероятностными методами, базирующимися на теореме Байеса (байесовский, или субъективный вывод).
Фишер заполнил пробел в этом важнейшем секторе, заложив основы методов оценки и вывода. Если Пирсон учил, как извлекать интересующую информацию из путаницы данных, Фишер показал, как познать целое (генеральную совокупность), наблюдая часть (выборку). Можно сказать, что Фишер довел до абсолюта создание методологического корпуса статистики: выбор теоретической модели на основании эмпирических данных, математическая дедукция свойств этой модели, оценка неизвестных параметров и заключительная проверка модели с помощью эксперимента. Подход, состоящий в сборе информации в ходе эксперимента и подготовке выводов на ее основе, составляет суть статистического вывода, и в отличие от вычисления вероятностей это не индуктивный, а дедуктивный процесс, сопровождающийся определенными ошибками, которые можно оценить количественно.
Сосредоточившись на проблемах теории оценки, Фишер утверждал, что речь шла о выборе наиболее подходящего значения из параметров генеральной совокупности традиционно обозначаемых греческими буквами, например, Θ на основании данных выборки, или, точнее, на базе статистик — обозначаемых латинскими буквами (например, Т), — которые вычисляются по наблюдаемым данным. Теория статистической оценки, разработанная Фишером, определяет, каким критериям должна удовлетворять хорошая оценка.
Сегодня три критерия, приведенные Фишером, претерпели небольшие изменения, хотя их смысл сохранился. Смещение. Оценка Т считается правильной и несмещенной для параметра Θ, если для любого размера выборки среднее значение ее распределения равно Θ, то есть если ожидаемое значение статистики Т является истинным значением Θ. Иначе оценка считается смещенной (рис. 4; подробнее см. СТАНДОТКЛОН.В и СТАНДОТКЛОН.Г: в чем различие?). Эффективность. Эффективность, или точность, оценки определяется как величина, обратная дисперсии ее выборочного распределения: чем больше дисперсия оценки, тем менее точна эта оценка, и наоборот. Это понятие имеет особенное значение для сравнения несмещенных оценок, так как среди них предпочтение должно отдаваться более эффективным, то есть обладающим наименьшей дисперсией.
Рис. 4. Пример смещения и эффективности оценки. Если мы сравним статистические оценки с выстрелами нескольких стрелков, то сможем лучше понять, какими качествами должна обладать хорошая оценка. Выстрелы стрелка А не отклоняются в каком-то одном направлении, но видно, что они очень разбросаны (это соответствует несмещенной, но неэффективной оценке). Выстрелы стрелка В смещены влево и разбросаны (оценка смещенная и неэффективная). Выстрелы стрелка С кучные, но отклоняются (смещенная и эффективная оценка).
Чтобы понять функцию правдоподобия, заново введенную Фишером и одну из самых важных для вывода, следует четко различать два довольно близких понятия. Пусть Θ — неизвестный параметр генеральной совокупности, а X — случайная выборка из этой совокупности. С одной стороны, существует вероятность получения выборки X при условии некоторого значения Θ (предполагается, что оно известно), обозначаемая р(Х|Θ) (X — переменная, Θ фиксирован) и определяющая вероятность появления всякой выборки. В проблеме оценки случай противоположный: наблюдается выборка X, но значение Θ неизвестно. Тем не менее описанная функция остается полезной, так как, заменив X на наблюдаемые значения, из Р(Х|Θ) для каждого Θ у нас будет вероятность получить выборочное значение X.
Изменяя Θ при фиксированном X, можно получить функцию, называемую функцией правдоподобия, L(Θ|Х), где X фиксирован, а Θ — переменная. Следует отметить, что после перемены ролей X и Θ в соответствии со сменой точки зрения, возникающей при выводе, функция правдоподобия не является распределением вероятности, следовательно, не подчиняется правилам исчисления вероятностей (когда подставляются конкретные значения из выборки). Эта функция отражает наши знания о параметре генеральной совокупности. И вместо предположения известного Θ и вычисления вероятности наблюдать различные выборки X допустим, что наблюдается конкретная выборка X, и оценим правдоподобие различных значений Θ.
Во время Второй мировой войны статистики, работавшие на стороне союзников, столкнулись с трудно разрешимой проблемой: как оценить общее количество немецких танков по серийным номерам захваченных у противника машин? Допустим, у захваченных танков следующие серийные номера: 2, 3, 7, 16. Сколько всего выпущено танков? Наилучшей возможной оценкой будет эффективная оценка (несмещенная и с минимальной дисперсией), формула которой для N следующая: m + (m –n)/n, где m — наибольший наблюдаемый серийный номер и n — размер выборки. Эту формулу можно интерпретировать как сумму выборочного максимума и «пустого среднего» выборки. К максимальному значению прибавляется среднее из промежутков между имеющимися наблюдениями, исходя из предположения, что за максимумом находится еще столько же элементов, сколько пропущено между имеющимися значениями. В нашем примере лучшей оценкой для N будет: 16 + (16-4)/4 = 19 танков всего (см. также Малые выборки в конкурентной разведке).
Фишер в качестве оценки Θ предлагал выбирать то значение, которое соответствует максимальной вероятности появления наблюдаемых значений выборки. Другими словами, выбрать такое значение параметра, которое, оказавшись реальным, максимизирует вероятность иметь данные, наблюдаемые в реальности.
Пусть имеется монета, вероятность выпадения сторон которой, р, неизвестна. Монету подбрасывают четыре раза и получают следующую серию: 0Р00 (орел-решка-орел-орел). Из исчисления вероятностей мы знаем, что Р(ОРОО|р) = р3(1 – р). Следовательно, функция правдоподобия: L(p|ОРОО) = р3(1 – р). Функция правдоподобия достигает максимума для значения 0,75. Таким образом, наша оценка, основанная на имеющейся выборке, будет р = 0,75. В сущности, это и есть основа метода оценки параметров с помощью максимального правдоподобия.
В период с 1923 по 1924 год Фишер писал книгу «Статистические методы для исследователей», которая увидела свет в 1925 году и на сегодняшний день выдержала 14 переизданий. Это самый влиятельный и известный труд Фишера. Он больше похож на учебник, чем на научную работу, благодаря убедительному стилю и характерному отсутствию математических доказательств. Возможно, этим и обусловлен его успех.
Фишер пишет, что статистика — не что иное, как математика, примененная к результатам наблюдений. При исследовании доступных выборок статистик делает выводы о полной совокупности, но они должны быть выражены не языком вероятности (как считали сторонники теоремы Байеса, или обратных вероятностных методов), а скорее языком правдоподобия.
Краеугольный камень всего произведения — «критерий значимости». В чем он состоит? В первую очередь в нулевой гипотезе Н0, которая устанавливает, например, что истинное значение неизвестного параметра таково, что Θ = Θ0. Затем выбирается статистика Т и вычисляется ее значение по данным имеющейся выборки X, которое обозначается как Т(Х). Так как распределение статистики Т в выборке известно, определяется вероятность того, что статистика Т примет значение, равное или большее наблюдаемому значению Т(Х), при условии, что нулевая гипотеза верна.
Математически это выглядит так: Р(Т ≥ Т(Х)|Н0). Это число было названо p-значением. Отсюда, если значение р достаточно мало — обычно меньше 0,05,— считается, что критерий оказался значимым и поэтому позволил опровергнуть нулевую гипотезу Н0. В противоположном случае тест оказывается незначимым для заранее заданного уровня значимости α = 0,05, отклонить нулевую гипотезу Н0 оказывается невозможным, ее временно принимают.
Нулевая гипотеза отвергалась только в том случае, когда вероятность наблюдать выборку вроде имеющейся была очень низка. Статистическое рассуждение основывалось на следующей логической дизъюнкции: «Или произошло исключительное событие (очень маловероятное), или нулевая гипотеза неверна». Очень малое р-значение указывало на то, что наблюдаемая выборка сильнее отличается от ожидаемой, чем это можно объяснить чистой случайностью, и поэтому исследователь имеет дело с неправдоподобной нулевой гипотезой, которую следует отвергнуть.
Закрепим эти понятия с помощью простой иллюстрации. Предположим, мы использовали новое удобрение на 20 растениях и наблюдали их рост в течение определенного периода времени, чтобы измерить влияние нового удобрения: увеличение (+) или уменьшение (–) скорости роста по отношению к росту без удобрения. Нашей нулевой гипотезой будет отсутствие всякого положительного эффекта от удобрения, то есть распределение ускорений (+) и замедлений (–) будет полностью случайным, как выпадение «орла» и «решки» при подбрасывании совершенно симметричной монеты. Поэтому, согласно нулевой гипотезе Н0, вероятность «+» будет равна вероятности «–», а именно Θ = 0,5. Представим, что после эксперимента мы видим 16 «+» и только 4 «–». Если мы выберем в качестве статистики Т количество наблюдаемых «+», то выяснится, что вероятность получить 16 или более «+» исходя из предположения, что вероятность положительного эффекта равна 0,5, составляет, как легко вычислить (см. рис. 5), только 0,006. То есть: Р(Т ≥ Т(Х)|Н0) = 0,006. Так как значение р ниже порогового α = 0,05, тест оказался значимым, и мы можем отвергнуть исходную нулевую гипотезу: есть эмпирические данные, противоречащие гипотезе, что удобрение не оказывает эффекта. Наоборот, все указывает на то, что оно стимулирует рост растений.
Рис. 5. Таблица расчета вероятности получить 16 и более «+» согласно формуле биномиальной вероятности:
Фишер предупреждал, что уровень значимости α не должен быть фиксированным, жестким. Впрочем, его предупреждение забыли и уровень 0,05 широко приняли, вплоть до того, что значение р = 0,051 стали считать незначимым, а 0,049 — значимым. Выбор этого граничного значения — вопрос не математический, универсальный, а зависит от прагматического контекста: если проверяется новое лекарство, то уровень значимости 0,05 несет 5%-ный риск того, что неэффективное лекарство будет признано эффективным (в этом случае, как и в некоторых других, уровень 0,01 или 0,001 может оказаться более подходящим).
Фишер описывал критерий значимости как способ отвергнуть нулевую гипотезу, которая никаким образом не может быть доказана и установлена окончательно. Этот подход, связанный с опровержением, соответствовал направлению «фальсификации», инициированному философом Карлом Поппером (1902–1994). Для статистика и философа наука характеризуется постановкой экспериментальных доказательств, которые могут опровергнуть или фальсифицировать теории, описываемые учеными (подробнее см. Карл Поппер. Логика научного исследования).
Методологически подход Фишера был разновидностью фальсификации в приложении к статистике: он состоял в опровержении гипотез, для которых наблюдения были относительно неправдоподобными. Нулевая гипотеза никогда не подтверждалась, но ее можно было опровергнуть. Если критерий оказывался значимым, гипотеза оказывалась неприемлемой в свете имеющихся данных; если нет, это говорило лишь о том, что гипотеза была совместима с данными.
Кроме критерия значимости, книга Фишера описывала дисперсионный анализ — другую инновационную статистическую методику, известную во всем мире по английской аббревиатуре ANOVA (англ. ANalysis Of VAriance; подробнее см. Однофакторный дисперсионный анализ).
В 1935 году Фишер издает книгу «Планирование экспериментов». Биологические исследования требуют проведения контролируемых экспериментов. Пассивного наблюдения недостаточно. Выборочная техника состоит в исследовании репрезентативной выборки из генеральной совокупности и измерении изучаемых показателей. Планирование экспериментов, наоборот, заключается в фиксации одних параметров и наблюдении за другими, с измерением возникающих изменений.
Объекты, получающие «обработку», являются экспериментальными единицами. Каждая обработка должна встречаться как минимум два раза, а лучше — несколько раз. Если мы хотим сравнить обработку А и В, в идеале следует применить их одновременно на множестве участков. Может случиться, что наблюдаемая разница между обработками А и В связана просто с разной плодородностью почвы на разных участках, а не с тем, например, что А эффективнее, чем В. Принцип повторения, сформулированный Фишером, помогал ограничить ошибку эксперимента, то есть случайную вариацию, не контролируемую экспериментатором (как, например, разная плодородность участков, на которых применяются обработки А и В).
Уравнение Фишера — обзор, формула и пример
Что такое уравнение Фишера?
Уравнение Фишера — это экономическая концепция, описывающая взаимосвязь между номинальными и реальными процентными ставками под влиянием инфляции Инфляция Инфляция — это экономическая концепция, которая относится к увеличению уровня цен на товары в течение определенного периода времени. Повышение уровня цен означает, что валюта в данной экономике теряет покупательную способность (то есть за ту же сумму денег можно купить меньше).. Уравнение утверждает, что номинальная процентная ставка равна сумме реальной процентной ставки плюс инфляция.
Уравнение Фишера часто используется в ситуациях, когда инвесторы или кредиторы просят дополнительное вознаграждение для компенсации потерь покупательной способности из-за высокой инфляции.
Концепция широко используется в области финансов и экономики. Он часто используется при расчете окупаемости инвестиций. Возврат инвестиций (ROI). Рентабельность инвестиций (ROI) — это показатель эффективности, используемый для оценки окупаемости инвестиций или сравнения эффективности различных инвестиций.или в прогнозировании поведения номинальных и реальных процентных ставок. Один из примеров — когда инвестор хочет определить фактическую (реальную) процентную ставку, полученную по инвестициям, после учета влияния инфляции.
Один интересный вывод уравнения Фишера связан с денежно-кредитной политикой. Это мощный инструмент. Уравнение показывает, что денежно-кредитная политика движет инфляцию и номинальную процентную ставку в одном направлении.Принимая во внимание, что денежно-кредитная политика обычно не влияет на реальную процентную ставку.
Американский экономист Ирвинг Фишер предложил уравнение.
Формула уравнения Фишера
Уравнение Фишера выражается следующей формулой:
(1 + i) = (1 + r) (1 + π)
Где:
i — номинальная процентная ставка ставка
r — реальная процентная ставка
π — уровень инфляции
Однако можно также использовать приблизительную версию предыдущей формулы:
i ≈ r + π
Пример уравнения Фишера
Предположим, что Сэм владеет инвестиционным портфелем.В прошлом году доходность портфеля составила 3,25%. Однако в прошлом году инфляция составляла около 2%. Сэм хочет определить реальную прибыль, которую он получил от своего портфеля. Чтобы найти реальную норму прибыли, мы используем уравнение Фишера. Уравнение утверждает, что:
(1 + i) = (1 + r) (1 + π)
Мы можем изменить уравнение, чтобы найти реальную процентную ставку:
Следовательно, реальная процентная ставка , или фактическая доходность инвестиций, портфеля равна:
Реальный процент, который инвестиционный портфель Сэма заработал в прошлом году с учетом инфляции, составляет 1.26% .
Дополнительная литература
CFI предлагает аналитика по финансовому моделированию и оценке (FMVA) ® Стать сертифицированным аналитиком финансового моделирования и оценки (FMVA) ® Сертификат CFI по анализу финансового моделирования и оценки (FMVA) ® поможет вам обрести уверенность в себе. необходимость в вашей финансовой карьере. Запишитесь сегодня! программа сертификации для тех, кто хочет вывести свою карьеру на новый уровень. Чтобы продолжать учиться и продвигаться по карьерной лестнице, вам будут полезны следующие ресурсы CFI:
- Эффективная годовая процентная ставка Эффективная годовая процентная ставка Эффективная годовая процентная ставка (EAR) — это процентная ставка, которая корректируется с учетом начисления сложных процентов за определенный период.Проще говоря, эффективная
- плавающая процентная ставка Плавающая процентная ставка Плавающая процентная ставка относится к переменной процентной ставке, которая изменяется в течение срока действия долгового обязательства. Это противоположность фиксированной ставки.
- Премия за рыночный риск Премия за рыночный риск Премия за рыночный риск — это дополнительная прибыль, которую инвестор ожидает от владения рискованным рыночным портфелем вместо безрисковых активов.
- Нормативная экономика Нормативная экономика Нормативная экономика — это школа мысли, которая считает, что экономика как предмет должна передавать оценочные заявления, суждения и мнения по
Формула эффекта Фишера | Финансы
Автор: Джеки Лохри | Рецензент: Райан Кокерхэм, CISI по рынкам капитала и корпоративным финансам | Обновлено 5 февраля 2019 г.
Формула «эффекта Фишера» пытается показать, как ожидание инфляции влияет как на процентные ставки, так и на покупательную способность.Эта формула, разработанная как часть общей экономической теории в 1930 году экономистом-математиком Ирвином Фишером, сегодня настолько широко используется в области экономики и финансов, что многие считают ее стилизованным фактом — термином, используемым для описания столь последовательных наблюдений и выводов. они общеприняты как истинные. Формула эффекта Фишера не только играет важную роль в научных кругах и бизнесе, но и имеет приложения, которые могут принести пользу вам как инвестору.
Совет
Формула, используемая для расчета эффекта Фишера, требует трех важных параметров данных: номинальной процентной ставки, реальной процентной ставки и ожидаемого уровня инфляции.
Компоненты формулы
Формула эффекта Фишера предполагает период в один год и разбивается на три компонента: номинальная процентная ставка, реальная процентная ставка и ожидаемый уровень инфляции. Номинальная процентная ставка — это процент, показывающий цену, которую вы платите за использование денег, без учета инфляции. Реальная процентная ставка — это процент, который корректируется для устранения влияния инфляции и, как следствие, является мерой «реальной» покупательной способности.Ожидаемый уровень инфляции — это процент, который меняется в зависимости от текущих экономических циклов.
Расчет формулы
Расчет эффекта Фишера нетрудно. Технический формат формулы: «Rном = Rreal + E [I]» или номинальная процентная ставка = реальная процентная ставка + ожидаемый уровень инфляции. Более простой способ вычислить формулу и определить покупательную способность — разбить уравнение на два этапа. Используйте уравнение «Rnon — E [I] = Rreal», чтобы получить реальный уровень инфляции.Затем определите реальную покупательную способность, скажем, 100 долларов, умножив 100 на Rreal и вычтя сумму из 100. Если Rreal составляет 2 процента, покупательная способность 100 долларов через год будет только 98 долларов.
Краткосрочные применения
Эффект Фишера имеет краткосрочные практические применения, которые вы можете использовать в повседневной жизни. Поскольку это оценка покупательной способности в соответствии с темпами инфляции, вы можете использовать ее для определения «реальной» нормы прибыли на инвестиции с точки зрения того, что ваша прибыль купит.Например, если вы приобретаете 12-месячный депозитный сертификат на 5000 долларов США с процентной ставкой 0,94 процента, а инфляция остается постоянной, в конце года у вас будет 5 047,22 доллара США.
Однако, если инфляция вырастет на 2 процента в течение года, вы фактически потеряете деньги на вложениях, поскольку ваша прибыль и покупательная способность ваших денег уменьшатся на 53,50 доллара или 1,06 процента, в результате чего вам останется потратить всего 4 993,72 доллара: 2,0 — 0,94 = 1,06 5047,22 * 1,06 = 53,50 5047,22 — 53.50 = 4993,72
Долгосрочное значение
Знание реальной процентной ставки имеет решающее значение для принятия правильных инвестиционных решений. Хотя использование эффекта Фишера для расчета краткосрочных эффектов может оказаться полезным и интересным, результаты становятся еще более значимыми при его применении в долгосрочной перспективе. Сложные расчеты с использованием ваших временных рамок для инвестирования, начальной суммы инвестиций, номинальной процентной ставки и различных темпов инфляции могут быть инструментом для анализа потенциала инвестиций и помочь вам определить, в какой момент экономического цикла имеет смысл продавать.
ИССЛЕДОВАНИЕ СТРУКТУРИРОВАННОГО УРАВНЕНИЯ ФИШЕРА С ПРИМЕНЕНИЕМ В БИОХИМИИ
DOI: 10.1137 / 16M1108546. Epub 2018 19 июня.Принадлежности Расширять
Принадлежность
- 1 Департамент прикладной математики, Университет Колорадо, Боулдер 80309-0526, США.
Элемент в буфере обмена
Джон Т. Нардини и др. SIAM J Appl Math. 2018.
Бесплатная статья PMC Показать детали Показать вариантыПоказать варианты
Формат АннотацияPubMedPMID
DOI: 10.1137 / 16М1108546. Epub 2018 19 июня.Принадлежность
- 1 Департамент прикладной математики, Университет Колорадо, Боулдер 80309-0526, США.
Элемент в буфере обмена
Полнотекстовые ссылки Опции CiteDisplayПоказать варианты
Формат АннотацияPubMedPMID
Абстрактный
Недавние биологические исследования были направлены на понимание того, как биохимические сигнальные пути, такие как семейство митоген-активируемых протеинкиназ (MAPK), влияют на миграцию популяции клеток во время заживления ран.Уравнение Фишера широко использовалось для моделирования экспериментальных анализов заживления ран из-за его простой природы и известных решений с бегущей волной. Однако это уравнение в частных производных с независимыми переменными времени и пространства не может учесть влияние биохимической активности на заживление ран. С этой целью мы выводим структурированное уравнение Фишера с независимыми переменными времени, пространства и уровня активности биохимических путей и доказываем существование самоподобного решения бегущей волны для этого уравнения.Мы показываем, что эти методы также применимы к общему структурированному уравнению реакции-диффузии и уравнению хемотаксиса. Мы также рассматриваем более сложную модель с различными фенотипами, основанную на активации MAPK, и численно исследуем, как различные временные паттерны биохимической активности могут приводить к увеличению или уменьшению скорости миграции населения.
Ключевые слова: Биохимические сигнальные пути; Сценическое строение; Решения бегущей волны; Лечение раны.
Цифры
Фиг.1
Аналитическое решение (9) для г…
Рис.1
Аналитическое решение (9) относительно g ( y ) = (1 — y ) / 2,…
рисунок 1Аналитическое решение (9) для g ( y ) = (1 — y ) / 2, A = 0 и ϕ ( y ) = e — ( y -0,1) 2 /,16 . (а) Графики t = σ ( y ; y̲ ) для y̲ = 0, 0.2, 0,4 и 0,6. (б) Графики y = σ −1 ( t ; y̲ ) для y̲ = 0, 0,2, 0,4 и 0,6. (в) Графики υ ( t ; y̲ ) для y̲ = 0, 0,2, 0,4, 0,6. Кривые на плоскости ( t, y ) обозначают y = σ −1 ( t ; y̲ ). (d) График в логарифмическом масштабе и ( t, y ). Пунктирная линия обозначает минимальный уровень поддержки с течением времени.
Фиг.2
Решения для (5) с f…
Фиг.2
Решения (5) с f ( t ) = βe γt — 1,…
Инжир.2Решения (5) с f ( t ) = βe γt — 1, g ( m ) = αm (1 — m ) при t = 0 , 5 и 10 для α = 1, β = 3, γ = −1/4. В каждом кадре контур в градациях серого изображает численное моделирование и ( t, x, m ) вдоль, красный профиль представляет собой график p ( t, m ) для обозначения профиля активации в момент времени. t , а синий профиль является результатом численного интегрирования u ( t, x, m ) по размеру m , чтобы отобразить w ( t, x ).
Фиг.3
Два изображения изоклины…
Фиг.3
Два изображения изоклины для u = 1 из численного моделирования…
Инжир.3Два изображения изоклины для u = 1 из численного моделирования (7) с g ( m ) = αm (1 — m ) и f ( t ) = β sin ( γt ) для α = 0,5, β = 1, γ = 1,615, D̄ = 100, λ̄ = 1/100 и начальное условие ϕ 1 ( м ) = 10 / 3 I [.05,0.35] ( м ) и ϕ 2 ( x ) = I [ x ≤5] ( x ). Числовая схема обсуждается в разделе 5.3, и используются размеры шага: Δ м = 1 / 80 , Δ x = 1 / 5 , Δ t = 10 −3 . Из (32) моделирование не должно пересекать плоскость м = м крит , которая задана красной плоскостью.На кадре (а) мы видим, что моделирование действительно пересекает плоскость м = м крит из-за числовой диффузии, которая вызывает высокую скорость диффузии вдоль x , видимую на кадре (b).
Фиг.4
Аналитическое решение для…
Рис.4
Аналитическое решение для активационного профиля, п ( т, м ), для…
Рис. 4Аналитическое решение для профиля активации, p ( t, m ), для примера 1 для α = 0,5 и ϕ 1 ( м ) = I (0,05, 0,35) ( м ). Сплошные черные кривые обозначают h ( t ; m̲ ) для m̲ = 0.05, 0,15 и 0,35, а пунктирная линия обозначает м = м крит . Обратите внимание, что шкала журнала используется по p для облегчения наглядности.
Фиг.5
Численное моделирование усредненного…
Рис.5
Численное моделирование усредненного нестационарного уравнения Фишера для примера 1. В (а),…
Рис. 5Численное моделирование усредненного нестационарного уравнения Фишера для примера 1. В (а) мы изображаем моделирование w ( t, x ) с течением времени для α = 0,05. Буквы «P» и «D» обозначают, когда популяция в первую очередь размножается или распространяется, соответственно.В (b) мы изображаем, как профиль для w ( t = 40, x ) изменяется для различных значений α . Описания «Нет активации», «Активация» и «Полная активация» обозначают значения α , для которых популяция полностью находится в неактивной популяции, разделена между активной и неактивной популяциями или полностью входит в активную популяцию при т. = 40 соответственно.
Рис.6
Аналитическое решение для…
Фиг.6
Аналитическое решение для активационного профиля, п ( т, м ), для…
Рис. 6Аналитическое решение для профиля активации, p ( t, m ), для примера 2 для α = 0.5, β = 3, γ = −1/4 и ϕ 1 ( м ) = I (0,05,0,35) ( м ). Сплошные черные кривые обозначают h ( t ; m̲ ) для m̲ = 0,05, 0,15 и 0,35, а пунктирная линия обозначает м = м крит . Обратите внимание, что шкала журнала используется по p для облегчения наглядности.
Рис.7
Численное моделирование усредненного…
Фиг.7
Численное моделирование усредненного нестационарного уравнения Фишера для примера 2. В (а),…
Рис. 7Численное моделирование усредненного нестационарного уравнения Фишера для Примера 2.В (а) мы изображаем моделирование w ( t, x ) с течением времени для α = 1, β = 8, γ = −1. Области, обозначенные буквами «P» или «D», обозначают, когда популяция в первую очередь размножается или диффундирует, соответственно. В (b) мы изображаем, как профиль для w ( t = 30, x ) изменяется для различных значений β . Описания «Нет активации», «Активация» и «Полная активация» обозначают значения β , для которых популяция полностью находится в неактивной популяции, разделена между активной и неактивной популяциями или полностью входит в активную популяцию при t = т макс .
Фиг.8
Аналитическое решение для…
Фиг.8
Аналитическое решение для активационного профиля, п ( т, м ), для…
Инжир.8Аналитическое решение для профиля активации, p ( t, m ), для примера 3 для α = 1/2, β = 4, γ = 1 и ϕ 1 ( м ) = I (0,05,0,35) ( м ). Сплошные черные кривые обозначают h ( t ; m̲ ) для m̲ = 0,05, 0,15 и 0,35, а пунктирная линия обозначает м = м крит .Обратите внимание, что шкала журнала используется по p для облегчения наглядности.
Фиг.9
Численное моделирование усредненного…
Фиг.9
Численное моделирование усредненного нестационарного уравнения Фишера для Примера 3.В (а),…
Рис.9.Численное моделирование усредненного нестационарного уравнения Фишера для примера 3. В (a) мы изображаем моделирование w ( t, x ) с течением времени для α = 0,5, β = 1, и γ = 4 π . Области, обозначенные буквами «P» или «D», обозначают, когда популяция в первую очередь размножается или диффундирует, соответственно. В (b) мы изображаем w ( t = 40, x ) для различных значений γ .Описания «Нет активации», «Активация» и «Полная активация» обозначают значения γ , для которых популяция полностью находится в неактивной популяции, разделена между активной и неактивной популяциями или полностью входит в активную популяцию при т. = т макс .
Все фигурки (9)
Похожие статьи
- Пространственная динамика для модели заживления эпидермальных ран.
Ван Х, Ву С. Ван Х и др. Math Biosci Eng. 2014 Октябрь; 11 (5): 1215-27. DOI: 10.3934 / mbe.2014.11.1215. Math Biosci Eng. 2014 г. PMID: 25347813
- Модель бегущей волны для вторжения клеток-предшественников и дифференцированных клеток.
Trewenack AJ, Landman KA. Trewenack AJ, et al. Bull Math Biol. 2009 Февраль; 71 (2): 291-317. DOI: 10.1007 / s11538-008-9362-х. Epub 2009 9 января. Bull Math Biol. 2009 г. PMID: 145
- Модель бегущей волны для интерпретации анализа миграции раневых клеток мезотелиальных клеток брюшины человека.
Майни П.К., МакЭлвейн Д.Л., Ливсли Д.И. Maini PK, et al. Tissue Eng. 2004 март-апрель; 10 (3-4): 475-82. DOI: 10.1089 / 107632704323061834. Tissue Eng. 2004 г. PMID: 15165464
- Эффекты транспортной памяти и нелинейного затухания в обобщенном уравнении Фишера.
Абрамсон Г., Епископ А.Р., Кенкре В.М. Abramson G, et al. Phys Rev E Stat Nonlin Soft Matter Phys. 2001 декабрь; 64 (6, балл 2): 066615. DOI: 10.1103 / PhysRevE.64.066615. Epub 2001 26 ноября. Phys Rev E Stat Nonlin Soft Matter Phys. 2001 г. PMID: 11736304
- Применение моделирования Монте-Карло в моделировании биохимических процессов.
Тенекеджиев К.И., Николова Н.Д., Колев К.Тенекеджиев К.И., и др. В: Mode CJ, редактор. Применение методов Монте-Карло в биологии, медицине и других областях науки [Интернет]. Риека (HR): InTech; 2011 28 февраля. Глава 4. В: Mode CJ, редактор. Применение методов Монте-Карло в биологии, медицине и других областях науки [Интернет]. Риека (HR): InTech; 2011 28 февраля. Глава 4. PMID: 28045483 Бесплатные книги и документы. Рассмотрение.
Процитировано
2 статей- Влияние численной ошибки на оценку параметров и количественную оценку неопределенности для адвективных моделей PDE.
Нардини JT, Bortz DM. Нардини Дж. Т. и др. Обратная Пробл. 2019 июн; 35 (6): 065003. DOI: 10.1088 / 1361-6420 / ab10bb. Epub 2019 29 мая. Обратная Пробл. 2019. PMID: 34121793 Бесплатная статья PMC.
- Изучение уравнений на основе биологических данных с ограниченными по времени выборками.
Нардини Дж. Т., Лагергрен Дж. Х., Хокинс-Дааруд А., Куртин Л., Моррис Б., Раттер Е. М., Суонсон К. Р., Флорес КБ.Нардини Дж. Т. и др. Bull Math Biol. 2020 9 сентября; 82 (9): 119. DOI: 10.1007 / s11538-020-00794-z. Bull Math Biol. 2020. PMID: 32
- 7 Бесплатная статья PMC.
Как уравнение Фишера является секретом экономики децентрализованного хранения | Тор | ppio
Уравнение Фишера не только помогает понять инфляцию и дефляцию в макроэкономике, но и служит моделью для децентрализованного рынка хранения и распределения.
Эта статья была опубликована на канале PPIO Medium Channel . PPIO — это децентрализованная платформа хранения и доставки для разработчиков, которые ценят скорость, доступность и конфиденциальность. Щелкните ссылку, чтобы просмотреть статьи о блокчейне, технологиях, обновлениях проектов и многом другом.
Основная информация
- Уравнение обмена Фишера — это математическая связь между количеством денег и экономикой.
- Обнаружено экономистом Ирвингом Фишером в рамках академической основы его количественной теории денег.
- Основываясь на уравнении обмена Фишера, мы обнаруживаем, что масштаб рынка децентрализованной полосы пропускания лучше, чем рынок хранилища.
- Инвестиционная стоимость сети блокчейн, в которой отсутствует жизнеспособное использование токена, в конечном итоге вернется к нулю.
Что такое уравнение обмена Фишера?
Концепция уравнения обмена впервые была предложена британским философом Дж. Милл. Впоследствии на его развитие оказал влияние шотландский философ Дэвид Юм.В 1911 году она была наконец известна миру благодаря элегантной алгебраической модели Ирвинга Фишера:
- M (денежная масса) относится к общей сумме номинальной денежной массы в обращении в экономике
- V (скорость обращения денег). относится к скорости обращения валюты и частоте потребления на единицу валюты
- P (Уровень цен) относится к уровню цен (на основе гипотетического метода измерения, такого как упаковка товаров)
- Q (Индекс реальных расходов) представляет индекс реального потребления
Формулу Фишера можно просто объяснить, умножив сумму денег на количество использований валюты.Результат равен экономическому выпуску, умноженному на среднюю цену, которая представляет собой номинальный доход общества, обычно известный как ВВП. С точки зрения традиционной макроэкономики, эта статья объяснит взаимосвязь между общим социальным доходом и количеством и скоростью спроса на деньги . Для нас также поучительно понять некоторые типичные экономические явления в реальном мире, такие как инфляция или дефляция.
В области блокчейнов дизайн экономической модели играет жизненно важную роль, поэтому формула Фишера имеет значение при разработке токеномики.Соучредитель Ethereum Виталик Бутерин в октябре 2017 года в своем блоге, посвященном оценке токенов среды обмена, резюмировал качественный анализ закона функционирования экономической системы, используя простую трансформацию уравнения обмена Фишера.
Об экономической инфляции и дефляции по формуле Фишера
Из формулы Фишера MV = PQ относительно легко вывести экономическую инфляцию и дефляцию в реальном мире.
Что касается инфляции, предполагается, что в течение определенного периода времени из-за относительно постоянных факторов, таких как социальные системы и привычки, скорость V денег в обращении остается в основном неизменной, в то время как в случае относительно стабильных условий занятости, сальдировано количество Q операций с товарами и услугами.Следовательно, можно вывести пропорциональную зависимость между суммой денег M и ценой товара P . В результате относительно легко понять, что рост цен вызван выпуском дополнительной валюты, то есть возникновением инфляции. Например, в Китае 1988, 1993 и 1995 годы были тремя годами высокого предложения валюты, что привело к увеличению индекса потребительских цен на 18,5%, 30% и 15% соответственно. Инфляция — это, по сути, денежное явление.При экономическом росте денежная масса должна увеличиваться, но если увеличившаяся часть превышает нормальный спрос экономического развития, и если нет эффективной политики корректировки цен или денежно-кредитной политики, избыточные деньги станут реальной покупательной способностью, вызывая общую уровень цен будет продолжать расти.
Как использовать уравнение Фишера для объяснения дефляции? Если формулу Фишера преобразовать в P = MV / Q, можно увидеть, что знаменатель — это количество Q операций с товарами и услугами.За исключением крупномасштабных войн, стихийных бедствий или события черного лебедя, Q практически не изменится, тогда как цена P прямо пропорциональна M и V . В некоторых крайних случаях вновь выпущенная валюта M не так быстро, как падение V , поэтому цена P будет снижаться, что является феноменом дефляции. Например, кризис субстандартного кредитования в США в 2008 году был, по сути, кредитным кризисом, который вызвал трудности с окончательной оплатой жилья, а затем сформировал цепной эффект на рынок.
В этой ситуации все ужесточили свои кошельки, что снизило скорость обращения денег V , что в конечном итоге привело к падению цен. Поэтому, если мы хотим сохранить P без изменений или вызвать его рост, самый простой способ борьбы с дефляцией — это напечатать деньги и сохранить большую денежную массу. Пакет мер стимулирования Китая в размере 4 триллионов и 10 триллионов кредитов в конце 2008 года были способом борьбы с дефляцией. Конечно, такое высвобождение денег не безгранично, и денежная инфляция от центрального банка может быть эффективной в краткосрочной перспективе, но токсичной в долгосрочной перспективе.
В идеальном дизайне модели токена экологические заинтересованные стороны сети блокчейн надеются, что цена токена будет продолжать расти, поэтому общий дизайн экономической модели будет учитывать следующие аспекты:
- Уменьшить количество денег в обращении и увеличить ожидания аналогично выкупу ценных бумаг. Это краткосрочное эффективное решение.
- Уменьшите скорость циркуляции. Примером этого является STEEM Power to STEEM, преобразование которого занимает 13 недель.Этот процесс чем-то похож на передачу токенов. Это краткосрочное эффективное решение.
- Повысьте общую емкость рынка, например DApp и сообщества разработчиков крупных публичных сетей. Это долгосрочное эффективное решение.
Что такое рынок пропускной способности децентрализованного хранилища?
В настоящее время основные компании, занимающиеся облачными вычислениями, такие как AWS и Alibaba Cloud, являются централизованными платформами с пропускной способностью хранилища, поскольку стоимость пропускной способности хранилища определяется предприятием в соответствии с рынком.Это дает преимущества высокой эффективности и гарантии коммерческого обслуживания, но его недостатки — высокая стоимость, низкий уровень безопасности и недостаточная защита конфиденциальности.
С другой стороны, рынок пропускной способности децентрализованного хранилища полностью использует неиспользуемые ресурсы пропускной способности хранилища. Сторона спроса работает путем поиска подходящего поставщика на открытом и прозрачном рынке. Механизм децентрализованного рыночного ценообразования может устранить монопольную ценовую власть компании (также известной как посредник) и полностью защитить основные интересы обеих сторон.Сторона предложения может даже снизить цену на стороне спроса. Это в сочетании с технологией распределенного хранения с открытым исходным кодом и сквозным шифрованием может обеспечить высокую безопасность данных, а также эффективно защитить конфиденциальность пользователей.
Конечно, на ранних этапах развития рынка пропускной способности децентрализованного хранилища из-за недостаточного количества ресурсов эффективность использования не будет высокой. Кроме того, сложно предоставить гарантию обслуживания на бизнес-уровне. Однако в долгосрочной перспективе, с увеличением неиспользуемых ресурсов полосы пропускания, развитием технологии блокчейн и эффективной экономической модели для стимулирования согласования спроса и предложения, рынок пропускной способности децентрализованного хранилища в конечном итоге переместится из периферийного в основной поток.
Со стороны спроса рынок децентрализованных хранилищ может стимулировать эффективное, дешевое и безопасное предложение пропускной способности хранилища, в то время как со стороны предложения большое количество недорогих и эффективных систем хранения данных может стимулировать спрос. Эти двое дополняют друг друга. Но когда мы вместе обсуждаем рынок хранилищ и рынок пропускной способности, мы теряем глубокое понимание их уникальных свойств.
На основе формулы Фишера мы проанализируем экономические модели спроса и предложения для рынка пропускной способности и хранилища соответственно, чтобы получить более глубокое и всестороннее представление о рынке пропускной способности децентрализованного хранилища.
Насколько важно уравнение обмена Фишера для разработки модели рыночной экономики децентрализованного хранилища с пропускной способностью
Давайте воспользуемся аналогией, чтобы лучше понять, как оно работает. Если автомобиль используется для представления данных пользователя, то шоссе может использоваться для представления рынка пропускной способности, и, наконец, парковка может использоваться для представления рынка хранения данных. Частота обращения денег, вызванная постоянным потоком транспортных средств на шоссе, будет намного выше, чем скорость обращения денег, приносимая транспортным средством на стоянке.Другими словами, с качественной точки зрения скорость обращения валюты на рынке пропускной способности намного выше, чем на рынке хранилищ. Из формулы Фишера можно сделать вывод, что рынок с той же суммой денег M , но с большей V будет иметь большую стоимость PQ . Поскольку значение PQ относится к номинальному доходу общества или ВВП, мы знаем, что на этом рынке будет более сильный экономический рост. Если мы сначала разработаем рынок пропускной способности, это будет способствовать быстрому запуску и развитию общего экономического объема, а значит, также будет стимулировать развитие рынка хранения данных.В обратном порядке результаты могут быть не такими здоровыми и упорядоченными.
Независимо от рынка пропускной способности или рынка хранения, цена P должна быть относительно стабильной, по крайней мере, с небольшими колебаниями, чтобы коммерческие пользователи могли легко пользоваться услугой или товарами на рынке. Однако, чем больше колебание значения скорости V , тем более своевременная регулировка и управление M потребуются для стабилизации P . Это требует, чтобы модель токена была более надежной и подробной, чтобы избежать нестабильности P из-за нестабильности, вызванной другими факторами.
С помощью формулы Фишера мы проанализировали децентрализованную пропускную способность, а также скорость и стоимость рынка хранилищ и пришли к выводу, что развитие рынка пропускной способности превосходит развитие рынка хранилищ. Мы также осознаем проблему поддержания стабильности цены. Основываясь на нашем текущем понимании, мы обсудим дизайн модели рыночной экономики децентрализованной пропускной способности хранилища.
Обсуждение Проектирование экономической модели пропускной способности децентрализованного хранилища
При разработке модели рыночной экономики децентрализованной пропускной способности хранилища мы сначала поставили несколько основных целей:
- Спрос на первом месте, предложение второе.Спрос стимулирует предложение, а предложение способствует спросу.
- Сторона предложения предлагает относительно стабильные возможности обслуживания и цену для стороны спроса.
- Дизайн токена может полностью побудить сообщество к максимальному развитию экосистемы.
Полезность токена должна быть больше, чем его инвестиционная стоимость, если она была меньше, то инвестиционная стоимость токена в конечном итоге вернется к нулю. Стоимость полезности должна сначала оплачиваться стороной спроса, поэтому мы должны сначала сосредоточиться на стороне спроса, а затем — на стороне предложения.Во-первых, мы ожидаем, что спрос будет стимулировать предложение, способствуя росту. Однако это создает возможность для стороны предложения быть подверженной колебаниям спроса, поэтому нам нужно подумать о том, как мотивировать сторону предложения предоставлять стабильные и обильные товары и услуги, чтобы в конечном итоге стимулировать рост спроса.
Если стоимость использования полосы пропускания децентрализованной платформы хранения напрямую связана с колебаниями цены одного токена, цена использования хранилища или полосы пропускания будет колебаться.Это неизбежно вызовет конфликт с пользователями со стороны спроса, которые стремятся к стабильным эксплуатационным расходам. Благодаря стабильной конструкции служебного токена эффективно реализуется стоимость использования услуг пропускной способности хранилища. Другими словами, цена служебного токена и токена стейкинга разделена, так что стоимость использования на уровне предприятия является управляемой и предсказуемой.
Если взаимная корреляция между токеном стекинга и токеном полезности сильна (например, в соответствии с моделью дизайна токенов VeChain), то токен стекинга может чеканить токен полезности.Однако для этого требуется централизованная модель для контроля стабильности служебного токена, которая не подчиняется духу децентрализации. Если, с другой стороны, взаимная корреляция слабая, то служебный токен может выкупить, сжечь или выплачивать дивиденды на этот токен, таким образом сохраняя децентрализованный подход. Однако в этом случае это может не соответствовать финансовым требованиям и требованиям безопасности. Нахождение баланса между централизацией и децентрализацией — ключ к успешной двухвалютной экономической модели.Мы подробно рассмотрим эту тему в следующих статьях.
Резюме
Уравнение обмена Фишера не только помогает нам понять инфляцию и дефляцию в макроэкономике, но также служит отличным ориентиром для разработки экономической модели блокчейна, особенно модели, используемой для рынка распределения децентрализованного хранилища.
(PDF) Мутимасштабный анализ уравнения Фишера
Мутимасштабный анализ уравнения Фишера
Карло Каттани1 и Алексей Кудрейко2
1diFarma
2D Кафедра математики,
c4, Университет Салерно, 9408, Виа Терно,
(SA), Италия
ccattani @ unisa.it, [email protected]
Аннотация. Вейвлет-подход Галеркина к решению одномерного
-мерного уравнения Фишера с использованием комплексных гармонических вейвлетов был представлен в этой статье до
. Коэффициенты линейной и нелинейной связи
были получены явно. Вейвлет-решение сравнивается с аналогичным аналитическим результатом.
Ключевые слова: уравнение Фишера, гармонические всплески, нелинейные параболические
уравнения, коэффициенты связности, метод Галеркина.
1 Введение
Классическое гармоническое вейвлет-преобразование было разработано Newland в 1993 году.
Подобно обычному дискретному вейвлет-преобразованию, классическое гармоническое вейвлет-преобразование
может также выполнять анализ функции с множественным разрешением. Кроме того,
имеет быстрый алгоритм, основанный на быстром преобразовании Фурье для численной реализации. Явным преимуществом гармонических вейвлетов является то, что они не пересекаются на частоте
.
Есть много статей и других публикаций о применении вейвлетов
для решения PDE, опубликованных с середины 1990-х годов (см.
, например, [1,6] и ссылки в нем). Кроме того, существуют другие исследования
на основе вейвлетов, направленные на решение нелинейных уравнений в частных производных. Один из таких методов был представлен
Бейлкиным и Кейзером [13], которые использовали метод полугрупп для преобразования
данного УЧП в интегральное уравнение перед его дискретизацией.
В этой статье мы используем комплексные гармонические вейвлеты в качестве масштабных функций в
проекции Галеркина уравнения Фишера [1,2]
∂u
∂t = c∂2u
∂x2 − bu2 + au, ( 1)
, где
t∈ [0, T], 0 u — частота нового мутантного гена. Работа А. Кудрейко поддержана Istituto Nazionale di Alta Matem- atica Francesco Severi (Rome-IT) в рамках стипендии U 2007/000394, 07.02.2007. O. Gervasi et al. (Ред.): ICCSA 2008, Часть I, LNCS 5072, стр. 1171–1180, 2008. c Springer-Verlag Berlin Heidelberg 2008 50 стр. Размещено: 29 сентября 2005 г. Дата написания: 15 сентября 2005 г. В этой статье мы изучаем инфляционный риск и временную структуру премий за инфляционный риск в США.S. номинальные процентные ставки через казначейские ценные бумаги защиты от инфляции (TIPS) и аналитическую двухфакторную модель Кокса-Ингерсолла-Росса (CIR) с коррелированными реальной ставкой и инфляцией. Аналитическая формула облегчает оценку параметров модели и повышает точность оценки номинальных ставок и TIPS, и, в частности, позволяет нам оценить временную структуру премий за риск инфляции. Мы используем двухфакторную модель для оценки облигаций с индексом инфляции и изучаем взаимосвязь между реальной ставкой и ожидаемым уровнем инфляции, подразумеваемым номинальными ставками казначейских облигаций с постоянным сроком погашения (CMT) за период с января 1998 года по декабрь 2004 года.Мы используем фильтр Калмана без запаха (UKF) для оценки модели и премии за инфляционный риск. Эмпирические данные показывают, что ожидаемый уровень инфляции, в отличие от уровня, рассчитываемого на основе индексов потребительских цен, очень стабилен, а премии за риск инфляции демонстрируют резкую временную структуру. Ключевые слова: CIR модель временной структуры процентных ставок, TIPS, фильтр Калмана без запаха, премия за инфляционный риск Классификация JEL: C51, G12, E43 Рекомендуемое цитирование: Предлагаемая ссылка Модель Фишера-КПП [1, 2, 3, 4, 5, 6] представляет собой одномерное уравнение реакции-диффузии, объединяющее линейную диффузию с нелинейным логистическим исходным членом, где — плотность населения, зависящая от положения и времени.Размерными параметрами в модели Фишера-КПП являются коэффициент диффузии, скорость распространения и плотность несущей способности. Решения модели Фишера-КПП на полубесконечной области, которые эволюционируют от начальных условий с компактной опорной асимптотой до бегущей волны с минимальной волновой скоростью в длительном пределе времени [1, 2, 3, 4, 5, 6 ]. Модель Фишера-КПП также приводит к решениям бегущей волны с начальными условиями, которые затухают достаточно медленно, как [1, 2, 3, 4, 5, 6], хотя для большинства практических приложений нас интересуют решения бегущей волны, которые движутся с минимальная волновая скорость, поскольку чаще всего актуальны начальные условия с компактной опорой [7, 8, 9, 10]. Модель Фишера-КПП и ее расширения успешно использовались в широком спектре приложений, включая изучение пространственного распространения инвазивных видов в экологии [11, 12, 13, 14]. В клеточной биологии пространственное распространение популяций инвазивных клеток моделировалось с использованием модели Фишера-КПП и ее расширений для ряда приложений, включая экспериментов in vitro, экспериментов по клеточной биологии [15, 16, 17, 18, 19] и . vivo злокачественное распространение [20, 21, 22].Другие области применения включают теорию горения [23, 24] и вторжение лесных пожаров [25]. Некоторые расширения модели Фишера-КПП включают работу с различными геометриями [9, 16], такими как внутреннее и внешнее распространение в геометриях с [26] и без [27] радиальной симметрии. Другие варианты включают: (i) рассмотрение моделей с нелинейным коэффициентом диффузии [28, 29, 30, 31, 32, 33]; (ii) включение различных нелинейных транспортных механизмов [34, 35, 36]; (iii) модели множественных вторгающихся субпопуляций [31, 37]; и (iv) многомерные модели, учитывающие анизотропию [38].Модель Фишера-КПП порождает решения в виде бегущей волны, которые не позволяют раствору исчезнуть. Карикатура такого поведения схематически показана на рис. 1 (a) — (c), где вторгающаяся популяция клеток будет бесконечно распространяться в полубесконечной области. (a) — (c) Эволюция популяции вторгшихся клеток при t = 0, t 1 и t 2 , с t 2 > t 1 > 0.В этом случае популяция вторгается в положительном направлении x на неопределенное время при условии, что область бесконечна. (d) — (f) Эволюция популяции клеток при t = 0, t 1 и t 2 , с t 2 > t 1 > 0. В этом случае популяция стремится к вымиранию. Обратите внимание, что мы намеренно сделали первоначальное распределение клеток в ( a ) и ( d ) одинаковым. Здесь нас интересует модель Фишера-Стефана [39, 40, 41, 42, 43, 44, 45] где — плотность населения, зависящая от местоположения и времени. Параметры в модели Фишера-Стефана такие же, как в модели Фишера-КПП, а также параметр Стефана, который связывает скорость изменения во времени с пространственным градиентом плотности на движущейся границе. Как и в случае с моделью Фишера-КПП, решения модели Фишера-Стефана, (2) — (5), также могут привести к постоянной скорости, бегущим волнам постоянной формы в течение длительного времени, как [39, 40, 41 , 42, 43, 44, 45].Интересно, что для тех же начальных условий, но другого выбора параметра ĩ, модель Фишера-Стефана также может привести к совершенно другому результату, когда популяция стремится к вымиранию, как [39, 40, 41, 42, 43, 44, 45]. Это существенное различие между моделями Фишера-КПП и Фишера-Стефана представляет большой интерес, поскольку модель Фишера-КПП никогда не приводит к вымиранию, независимо от выбора параметров. Один из способов интерпретации этого различия состоит в том, что модель Фишера-Стефана способна уловить и предсказать дополнительные детали, которые представляют практический интерес, поскольку хорошо известно, что многие изначально небольшие перемещенные популяции вымрут [46].Это один из недостатков модели Фишера-КПП, поскольку эта модель подразумевает, что каждая небольшая начальная популяция всегда приводит к успешному вторжению. Модель Фишера-Стефана является адаптацией модели Фишера-КПП, которая включает движущуюся границу, x = L ( t ), вдохновленную классической задачей Стефана [47, 48]. Классическая проблема Стефана представляет собой одномерную модель теплопроводности, которая включает фазовый переход, такой как теплопроводность, связанная с таянием льда в воду [47, 48].Интересной математической и физической особенностью классической задачи Стефана является то, что граница раздела между двумя фазами может перемещаться со временем, что дает начало понятию задачи с движущейся границей [47, 48]. В отличие от классических моделей теплопроводности без какого-либо фазового перехода [49], решение задачи Стефана требует задания двух граничных условий на движущейся границе раздела [47, 48]. Во-первых, на движущейся границе задается температура, при которой происходит фазовый переход.Это аналогично уравнению (3) в модели Фишера-Стефана. Во-вторых, условие Стефана определяет баланс последней тепловой энергии с удельной тепловой энергией на движущейся границе, связывая скорость изменения положения движущейся границы во времени с потоком тепла на границе [47, 48]. Это аналогично уравнению (5) в модели Фишера-Стефана. Задача с подвижной границей (2) — (5) с представляет собой однофазную задачу Стефана, которая имеет начальную область твердого тела, при некоторой начальной температуре, изолированную при.Первоначально интервал занят жидкостью, предположительно уже имеющей температуру плавления. Для этой конкретной постановки задачи граница раздела распространяется в жидкую область, поскольку тепловая энергия, содержащаяся в твердом теле, постоянно используется в качестве скрытой энергии для преобразования жидкости в твердое тело. Процесс продолжается до тех пор, пока простой энергетический баланс не показывает, что К этому результату мы вернемся позже. Более общая двухфазная задача Стефана включает теплопроводность в обеих фазах, опять же разделенных движущейся границей [47, 48].Так же, как модель Фишера-КПП, существует множество расширений классической проблемы Стефана, таких как работа с более высокими измерениями [50, 51, 52, 53], а также изменение подвижного граничного условия [54]. Мы вводим безразмерные переменные и, чтобы изменить масштаб уравнения Фишера-КПП как на x ≥ 0 и t > 0. Полезно отметить, что уравнение (7) не включает свободных параметров и что все решения уравнения (7) с начальными условиями с компактной опорой в конечном итоге приведут к бегущим волнам, которые движутся со скоростью c мин = 2.В этой работе мы всегда указываем начальное условие. где H ( x ) — обычная функция Хевисайда, а α и β — положительные константы, так что мы имеем u ( x , 0) = α для x < β и u ( x , 0) = 0 для x > β . Все численные решения уравнения (7), которые мы рассматриваем, задают ∂u / ∂x = 0 при x = 0 и ∂u / ∂x = 0 при x = x ∞ .Здесь x ∞ выбрано достаточно большим, чтобы мы могли численно аппроксимировать задачу с бесконечной областью на 0 ≤ x <∞ с помощью задачи конечной области 0 ≤ x ≤ x ∞ [31 ]. Полная информация о численном методе, используемом для решения уравнения (7), а также тестовых примерах для подтверждения точности наших численных решений приведены в дополнительном материале. Чтобы сделать модель Фишера-Стефана безразмерной, мы используем те же безразмерные переменные, что и в модели Фишера-КПП с и, так что мы имеем Уместно отметить, что безразмерная модель Фишера-Стефана включает один параметр, κ > 0.Для всех численных решений уравнений (9) — (12) в этой работе мы всегда применяем начальное условие (8) такое, что L (0) = β . Полная информация о численном методе, используемом для решения уравнений (9) — (12), вместе с тестовыми примерами для подтверждения точности наших численных решений, приведены в дополнительном материале. Чтобы проиллюстрировать ключевые особенности моделей Фишера-КПП и Фишера-Стефана, мы представляем численные решения обеих моделей на Рисунке 2. Результаты на Рисунке 2 (а) показывают временную эволюцию решения модели Фишера-КПП от начального состояние с компактной опорой.Построение графика решений при равных промежутках времени предполагает, что решение приближается к бегущей волне, которая движется в положительном направлении x с постоянной скоростью и постоянной формой. Наши численные решения подтверждают, что скорость распространения равна c = 2, как и ожидалось. Результаты на рисунке 2 (b) показывают временную эволюцию решения модели Фишера-Стефана (9) — (12) для тех же начальных условий, что и на рисунке 2 (a), вместе с конкретным выбором κ .Опять же, построение графика решений при равных промежутках времени предполагает, что решение приближается к бегущей волне, которая движется в положительном направлении x с постоянной скоростью и постоянной формой. В этом случае для нашего выбора κ наше численное решение предполагает, что c = 1,2, что меньше минимальной скорости волны для модели Фишера-КПП. Еще одно важное различие между решениями бегущей волны на рис. 2 (a) — (b) состоит в том, что решение бегущей волны модели Фишера-КПП не имеет компактной опоры, поскольку u ( x, t )> 0 для всех x ≥ 0 и t > 0.Эта особенность модели Фишера-КПП ранее подвергалась критике как биологически неправдоподобная [7, 8, 30] (и это наблюдение побудило расширить модель Фишера-КПП, включив в нее различные члены нелинейной диффузии, чтобы результирующие бегущие волны имели компактные размеры). поддержка [7, 8, 30]). Напротив, из-за граничных условий при x = L ( t ) решение бегущей волны Фишера-Стефана имеет компактную поддержку, поскольку у нас есть u ( L ( t ), т ) = 0 для т > 0.Таким образом, мы определили две особенности модели Фишера-Стефана, которые привлекательны по сравнению с моделью Фишера-КПП: (i) модель Фишера-Стефана допускает вымирание популяции, тогда как модель Фишера-КПП подразумевает, что все начальные популяции успешно вторгаются; и (ii) решения бегущей волны модели Фишера-Стефана имеют компактный носитель, тогда как решения бегущей волны модели Фишера-КПП нет. Результаты на рисунке 2 (c) показывают временную эволюцию решения модели Фишера-Стефана для того же начального условия, которое использовалось на рисунке 2 (a) — (b), но на этот раз мы выбрали немного меньшее значение κ .В этом случае нанесение решений на те же интервалы времени, что и на рисунках 2 (a) — (b), указывает на то, что решение не стремится к бегущей волне, а вместо этого, кажется, исчезает. На рисунке 2 (d) показано увеличенное изображение решения на рисунке 2 (c), чтобы мы могли видеть дополнительные детали, поскольку популяция стремится к вымиранию. Из этих увеличенных решений мы видим, что L (0) = 1. По t = 20 мы имеем L (20) ≈ 1,4, и после этого времени решение быстро стремится к нулю.Вместе результаты на рис. 2 (b) — (d) иллюстрируют дихотомию распространения-экстинкции, поскольку для определенных вариантов κ мы наблюдаем распространение как бегущую волну на рис. 2 (b), при сохранении всего остального идентичного, за исключением выбора меньшее значение κ мы наблюдаем вымирание популяции на рис. 2 (c) — (d). Эти первоначальные сравнения на Рисунке 2 показывают, что модели Фишера-КПП и Фишера-Стефана выглядят очень разными, когда они описываются в терминах дифференциальных уравнений в частных производных.Поскольку обе модели Фишера-КПП и Фишера-Стефана поддерживают решения с бегущей волной, мы теперь исследуем эти модели на фазовой плоскости. Численные решения модели Фишера-КПП на рисунке 2 (a) предполагают, что мы ищем решения бегущей волны с координатой бегущей волны z = x — ct , где c > 0 — постоянная скорость распространения в положительном направлении x . В координате бегущей волны уравнение (7) упрощается до нелинейного обыкновенного дифференциального уравнения второго порядка где –∞ < z <∞, при этом U (–∞) = 1 и U (∞) = 0. (a) Численное решение уравнения Фишера-КПП, превращающееся в решение бегущей волны с минимальной скоростью волны, c = 2. (b) Численное решение модели Фишера-Стефана, развивающееся в решение бегущей волны с c = 1.2. (c) Численные решения модели Фишера-Стефана, ведущие к вымиранию. (d) Увеличенная область раствора в (c) , для 0 ≤ x ≤ 2, чтобы четко показать динамику поведения поглощения.Для модели Фишера-Стефана мы устанавливаем κ = 20 в (b) и κ = 0,45 в (c). Численные решения модели Фишера-Стефана получены с Δ ξ = 1 × 10 -4 , тогда как численные решения модели Фишера-КПП получены с Δ x = 1 × 10 -4 . Для моделей Fisher-KPP и Fisher-Stefan мы устанавливаем Δ t = 1 × 10 -3 и ϵ = 1 × 10 -8 . Для всех представленных здесь результатов начальным условием является уравнение (8) с α = 0.5 и β = 1. Наш подход к анализу решений бегущей волны модели Фишера-Стефана очень похож, за исключением того, что мы должны сначала предположить, что наш выбор начального условия и κ в уравнении (9) — Уравнение (12) таково, что образуется решение с бегущей волной, как на рисунке 2 (b). В этом случае запись уравнения (9) в координате бегущей волны приводит к тому же нелинейному обыкновенному дифференциальному уравнению второго порядка, уравнению (13) с U (–∞) = 1.Другое граничное условие трактуется иначе, и чтобы увидеть это различие, мы выражаем условие Стефана (12) через z , давая при z = L ( t ) — ct . Для решения с бегущей волной мы имеем d L ( t ) / d t = c , так что дифференциальное уравнение (13) выполняется при –∞ < z <0. Приведены граничные условия на U (–∞) = 1, и при z = 0.Следовательно, хотя модели Фишера-КПП и Фишера-Стефана, представленные в виде уравнений в частных производных, кажутся очень разными, когда мы ищем решения этих моделей в виде бегущей волны в координатах бегущей волны, мы обнаруживаем, что уравнения, управляющие фазовыми плоскостями для двух моделей одинаковы, с отличиями только в одном граничном условии. Сначала мы исследуем решения бегущей волны модели Фишера-КПП, переписав уравнение (13) как динамическую систему первого порядка. для –∞ < z <∞.Эта динамическая система имеет две точки равновесия: (1, 0) и (0, 0). Решение с бегущей волной соответствует гетероклинической траектории между этими двумя точками равновесия. Линейный анализ показывает, что собственные значения в (1,0) таковы, что локальное поведение в (1, 0) является седловой точкой [5]. Собственные значения в точке (0,0) равны, что означает, что локальное поведение в точке (0, 0) является стабильной спиралью для c <2 и стабильным узлом, если c ≥ 2 [5]. Следовательно, чтобы избежать нефизических отрицательных решений вблизи (0, 0), нам требуется c ≥ 2, что дает хорошо известную минимальную скорость волны для модели Фишера-КПП [5, 55]. Форма гетероклинической орбиты между (0, 0) и (1, 0) задается решением Ни система (16) — (17), ни уравнение (18) не имеют точных решений для произвольного выбора c > 0. Поэтому мы будем рассматривать численные решения этих обыкновенных дифференциальных уравнений, когда мы представляем и визуализируем фазу плоскости в этой работе, с подробностями о численных методах, приведенными в дополнительном материале. Решение уравнения (18) дает V ( U ), а решение системы (16) — (17) дает U ( z ) и V ( z ), и мы будем использовать оба подхода, где это уместно.Все фазовые плоскости, представленные в этой работе, созданы с использованием комбинации точных и численных методов, которые описаны в дополнительном материале. Результаты на рисунке 3 показывают набор фазовых плоскостей для модели Фишера-КПП для диапазона c. Фазовая плоскость на рисунке 3 (а) показывает поток, расположение точек равновесия и гетероклиническую орбиту для c = 10, которая приближается к (0, 0) без спирали. Когда мы строим решение в терминах плотности, U = U ( z ) на рисунке 3 (b), мы видим, что это решение является положительным и монотонно убывающим.Аналогичные результаты представлены на рисунке 3 (c) — (d) для c = 2. Напротив, на рисунке 3 (e) показана фазовая плоскость для c = 0,5, где мы видим, что гетероклиническая орбита приближается к (0, 0) в виде спирали, что указывает на то, что U ( z ) <0 в определенных регионах. Это колебательное поведение часто используется для обоснования условия, согласно которому c ≥ c мин = 2 для модели Фишера-КПП, а возможность бегущих волн с c <2 обычно игнорируется [5]. Поскольку решения бегущей волны моделей Фишера-КПП и Фишера-Стефана управляются одними и теми же фазовыми плоскостями, стоит изучить, как фазовые плоскости на рисунке 3 соотносятся с решениями бегущей волны модели Фишера-Стефана. Как указывалось ранее, решения модели Фишера-Стефана с бегущей волной удовлетворяют другому граничному условию (15). Траектория в фазовой плоскости должна пересекаться и заканчиваться в (0, — c / κ ). Фазовые плоскости и гетероклинические орбиты для модели Фишера-КПП на рисунках 3 (a) — (b) и 3 (c) — (d) показывают, что такое пересечение невозможно для этих вариантов c ≥ 2.В частности, линеаризация около (0, 0) означает, что всякий раз, когда c ≥ 2, стабильный узел в (0, 0) исключает возможность такого пересечения, демонстрируя, что не существует решения бегущей волны для уравнения Фишера-Стефана. модель с c ≥ 2. Напротив, для c <2 траектория пересекает ось V в какой-то момент (0, — c / κ ), как показано оранжевым диском на рисунке 3 (е). Следовательно, это дополнительное граничное условие для модели Фишера-Стефана вместе с линеаризацией около (0, 0) указывает, что решения бегущей волны для модели Фишера-Стефана требуют c <2.В этих условиях траектория на фазовой плоскости, соответствующая решениям бегущей волны, проходит между (1, 0) и (0, - c / κ ), и, поскольку эта траектория никогда не покидает четвертый квадрант, мы избегаем проблемы отрицательного плотности. Это очень интересный результат, потому что стандартный анализ фазовой плоскости решений бегущей волны модели Фишера-КПП обычно отбрасывает любые решения, для которых c <2, исходя из физических соображений [5, 4]. В этой работе мы показываем, что именно эти обычно отбрасываемые решения фактически составляют основу решений бегущей волны модели Фишера-Стефана при условии, что мы рассматриваем только ту часть траектории между (1, 0) и (0, - c / κ ), где U ( z ) ≥ 0.Таким образом, пересматривая решения бегущей волны модели Фишера-КПП, мы обеспечиваем понимание свойств решений бегущей волны модели Фишера-Стефана. (a) Фазовая плоскость и гетероклиническая траектория для c = 10. (b) Соответствующий профиль плотности для гетероклинической траектории в (a ). (c) Фазовая плоскость и гетероклиническая траектория для c = 2. (d) Соответствующий профиль плотности для гетероклинической траектории в (c). (д) Фазовая плоскость и гетероклиническая траектория для c = 0,5. (f) Соответствующий профиль плотности для гетероклинической траектории в (e) . Точки равновесия в (1,0) и (0,0) показаны черными дисками. Синие стрелки показывают поток, связанный с динамической системой, а сплошная оранжевая линия показывает гетероклиническую траекторию, которая проходит между (0, 0) и (1, 0).Оранжевый диск на (c) показывает место, где гетероклиническая траектория пересекается с осью U ( z ) = 0, где V ( z ) <0. Теперь мы предоставляем дополнительные результаты на рисунке 4. сравнение траекторий в фазовой плоскости для более широкого диапазона c . Траектории для c = 10 и c = 2 показывают гетероклиническую орбиту, которая проходит между (1,0) и (0,0) без спирали вокруг начала координат. Эти траектории связаны с решениями бегущей волны модели Фишера-КПП для этих вариантов выбора c .Также показаны дополнительные результаты для c = 1 и c = 0,5, и эти траектории явно закручиваются по спирали вблизи начала координат. Однако обе эти траектории в какой-то момент пересекают ось V , показанную соответствующим цветным диском на рисунке 4, где U ( z ) = 0 и V ( z ) <0, которое удовлетворяет граничному условию Фишера-Стефана. Поэтому траектории на Рисунке 4 с c <2 показаны как комбинация сплошных и пунктирных линий.Те части траекторий, которые показаны сплошным шрифтом, соответствуют решению бегущей волны модели Фишера-Стефана, тогда как штриховые части траекторий не связаны с решением бегущей волны. Наконец, мы также включаем траекторию на рис. 4 для c = 0. В этом случае траектория образует гомоклиническую орбиту с (1, 0). Хотя эта траектория не соответствует никаким бегущим волнам с c > 0, позже мы покажем, что это важно при построении приближенных решений возмущений для c ≪ 1. Интересно напомнить, что все решения модели Фишера-КПП с компактной опорой всегда в конечном итоге будут формировать бегущую волну с минимумом, c = c min = 2. Напротив, вопрос о том, сформируется ли решение бегущей волны для модели Фишера-Стефана для конкретного начального условия, и если да, то с какой скоростью движутся эти решения бегущей волны, зависит от выбора κ .Чтобы исследовать эту взаимосвязь, мы используем комбинацию анализа фазовой плоскости и численных решений модели Фишера-Стефана (9) — (12). Повторно решая уравнения фазовой плоскости (16) для различных значений c <2, мы можем оценить точку, в которой траектория, выходящая из точки (1, 0), сначала пересекает ось V , а затем использовать граничное условие, V = — c / κ при U = 0, чтобы вычислить соответствующие значения κ .Сплошная синяя кривая на рисунке 5 показывает эту взаимосвязь. Мы обнаружили, что по мере того, как мы исследуем возрастающие значения c по направлению к пороговому значению c = 2, κ , кажется, неограниченно растут. Этот числовой результат предполагает, что κ → ∞ как c → 2 — . Гетероклинические траектории между (1, 0) и (0, 0) для c = 10, 2, 1 и 0.5, как указано. Дополнительная траектория с c = 0 образует гомоклиническую траекторию к (1, 0). Точки равновесия в (1, 0) и (0, 0) показаны черными дисками. Для траекторий с 0 < c <2 пересечение с осью U ( z ) = 0 показано диском соответствующего цвета: зеленый для c = 1; синий для c = 0,5; и красный для c = 0. В дополнение к изучению взаимосвязи между κ и c на фазовой плоскости, мы также решаем модель Фишера-Стефана (9) — (12) численно с конкретным выбором начального условия, заданного уравнением (8) с β = 1 и α = 0.5. Мы позволяем таким численным решениям развиваться в течение достаточного времени, чтобы полученное решение, казалось, превращалось в бегущую волну, по которой мы можем оценить скорость, c . Повторение этой процедуры для различных значений κ позволяет нам оценить, как наш выбор κ влияет на c . Дополнительные результаты, показанные в виде красных дисков на Рисунке 5, подтверждают, что решения длинновременной бегущей волны из описания уравнения в частных производных очень хорошо сравниваются с соотношением, подразумеваемым анализом фазовой плоскости. Как мы заявили в разделе 2.2, всякий раз, когда мы работаем на фазовой плоскости, мы делаем неявное предположение, что было сгенерировано решение с бегущей волной. Тем не менее, когда мы сравниваем результаты на рис. 2 (b) — (c), мы знаем, что решения с бегущей волной с длительным временем существования не всегда образуются, поскольку этот результат зависит от выбора κ . Мы исследуем это на рисунке 5, поддерживая постоянное начальное условие при численном решении уравнений (9) — (12) и выбирая последовательность все более малых значений κ .Численные решения предполагают, что для этого начального условия существует пороговое значение, κ крит ≈ 0,48. Если κ > κ крит , мы наблюдаем решения с долгой бегущей волной, а для κ < κ крит популяция в конечном итоге вымирает. Это приблизительное пороговое значение показано на рисунке 5 в виде черного диска, а соотношение между κ и c , полученное на фазовой плоскости, показано сплошной линией для κ > κ crit и как пунктирная линия для κ < κ крит . Теперь мы рассмотрим дихотомию распространения-исчезновения, исследуя зависящие от времени численные решения уравнений (9) — (12) на рисунке 6. Рисунок 6 (a) показывает численные решения для конкретного выбора κ , где мы видим очень интересное поведение. При t = 30 решение, кажется, уменьшается, почти до исчезновения, тогда как к t = 60 и t = 90 решение восстанавливается после начального спада, чтобы в конечном итоге сформировать бегущую волну.Результаты на Рисунке 6 (b) показывают детали при промежуточных значениях т , чтобы четко выделить это начальное снижение, за которым последовало восстановление. На рисунке 6 (c) показаны оценки L ( t ) в зависимости от t , где мы видим, что L ( t ) медленно увеличивается в ранний момент, прежде чем в конечном итоге будет увеличиваться с постоянной скоростью, что соответствует к решению с бегущей волной. Напротив, на рисунке 6 (d) показано решение уравнений (9) — (12) для тех же начальных условий, что и на рисунке 6 (a), с той лишь разницей, что κ уменьшается.Рисунок 6 (d) показывает, что популяция, похоже, почти вымерла при t = 90, а дополнительные результаты, увеличенные на Рисунке 6 (e), показывают, что популяция никогда не восстанавливается, а вместо этого вымирает с увеличением t. Рисунок 6 (f) показывает L ( t ) как функцию от t , где мы видим, что распространяющаяся популяция сначала увеличивает свою область, а затем в конечном итоге останавливается. Ду и его коллеги [39] предоставляют формальное доказательство существования критической длины, L крит = π /2, так что если вообще L ( t )> L крит популяция превратится в бегущую волну, тогда как если L ( t ) никогда не превысит эту критическую длину, то популяция в конечном итоге вымрет.Здесь мы приводим несколько простых физических и математических аргументов, подтверждающих этот результат. Визуальный осмотр зависящих от времени решений уравнений (9) — (12) на рис. 6 (b) и (e) подтверждает, что у нас u ( x, t ) ≪ 1 близко ко времени, когда популяция восстанавливается ( Рисунок 6 (b)) или не восстанавливается после первоначального спада (Рисунок 6 (e)). Это наблюдение предполагает, что мы можем изучить дихотомию распространения-вымирания, используя приближенную линеаризованную модель, где. В этих условиях мы можем аппроксимировать модель Фишера-Стефана с помощью для 0 < x ≤ L , при x = 0 и x = L .В этом приближенном анализе мы рассматриваем длину домена как фиксированную величину, и это позволяет нам записать точное решение линейного уравнения (19) как где γ n = π (2 n — 1) / (2 L ), n = 1,2,3,…, и A n — константы, выбранные таким образом что решение соответствует начальному условию,. Решение линеаризованной модели (20) можно дополнительно упростить, если предположить, что динамику около времени восстановления или сокращения популяции можно аппроксимировать, взяв приближение главного собственного значения, так что Первая строка представляет явление вторжения, а вторая строка представляет вымирание. (a) Эволюция во времени профилей плотности инвазивной клеточной популяции. (b) Увеличенные профили плотности в (a) от x = от 0 до 2 в промежуточные моменты времени. (c) Прогрессия L ( t ) с наложением критической длины π /2. (d) Временная эволюция профилей плотности поглощения. (e) Увеличенные профили плотности в (d) от x = от 0 до 2 в промежуточные моменты времени. (f) Прогрессия L ( t ) с наложением критической длины π /2. Для обоих расчетов Δ ξ = 1 × 10 −4 , Δ t = 1 × 10 −3 , ϵ = 1 × 10 −8 , и начальное условие, заданное уравнением ( 8), где α = 0,5 и β = 1. Результаты в (a) — (c) соответствуют κ = 0.5, а результаты в (г) — (е) соответствуют κ = 0,4. С помощью этого приближенного решения мы формулируем утверждение сохранения, описывающее скорость изменения всего населения в пределах области во времени, где — общая численность населения в домене. Первое слагаемое справа в уравнении (22) — это скорость увеличения общей численности населения из-за исходного члена, а второе слагаемое справа в уравнении (22) — это скорость уменьшения общей численности населения из-за диффузионных потерь. на границе x = L .Установка d n / d t = 0 и подстановка уравнения (21) в (22) дает L = L крит = π /2, подтверждая результаты Du и его коллег [39 ]. Мы интерпретируем это следующим образом. Как только зависящее от времени решение уравнений (9) — (12) развивается так, что L ( t )> π /2, чистое положительное накопление массы в системе означает, что решение с бегущей волной в конечном итоге будет форма, как на Рисунке 6 (a) — (c).В качестве альтернативы, если зависящее от времени решение уравнений (9) — (12) развивается так, что L ( t ) никогда не превышает π /2, чистая потеря массы в системе означает, что популяция всегда будет вымерли, как на Рисунке 6 (d) — (f). Также стоит отметить, что, поскольку результат L крит = π /2 регулируется линеаризованной моделью (19), этот результат будет справедливым для любого обобщения модели Фишера-Стефана, которое может быть линеаризовано. чтобы получить уравнение (19).Например, если мы расширили модель Фишера-Стефана, чтобы рассмотреть обобщенный логистический источник, u (1 — u ) m , где m > 0 — некоторая экспонента [56, 57, 58 , 59], то применимо то же значение L крит = π /2. Теперь, когда мы использовали анализ фазовой плоскости и линеаризацию, чтобы установить условия для формирования решений бегущей волны модели Фишера-Стефана, мы обращаем наше внимание на то, можно ли обеспечить дополнительное математическое понимание деталей формы решений бегущей волны.Хорошо известно, что решения бегущей волны модели Фишера-КПП распространяются со скоростью c ≥ c min = 2, и что можно получить приближенные решения возмущений для описания формы решений бегущей волны в большой c предел [4, 5, 31]. Теперь мы попытаемся провести аналогичный анализ, чтобы описать форму решения бегущей волны модели Фишера-Стефана, для которой c <2, предлагая нам попытаться найти решение возмущения для малых c : Подставляя это разложение в уравнение (18), мы получаем обыкновенные дифференциальные уравнения, определяющие V 0 ( U ) и V 1 ( U ), которые можно точно проинтегрировать.Убедившись, что U (1) = 0, получим Поскольку мы решаем как для V 0 ( U ), так и для V 1 ( U ), мы можем построить как решение для возмущений, заданное как, так и решение для возмущений, заданное как. Чтобы сравнить точность этих решений возмущений для формы кривой V ( U ) на фазовой плоскости, мы генерируем серию численных фазовых плоскостей для c = 0.05,0.1, 0,2 и 0,5 на рисунке 7. Числовые траектории, показанные синим цветом, проходят между точками равновесия (1, 0) и (0, 0) и проходят через точку (0, — c / κ ). В численных решениях выделим (0, — c / κ ) синим диском. На каждой части рисунка 7 мы сравниваем численную траекторию с решением возмущения, выделенным красным цветом. В каждом случае решение возмущения проходит между (1, 0) и сначала пересекает ось V в точке. Мы показываем эту точку красным диском.Сравнение красной и синей траекторий в четвертом квадранте ясно показывает расхождение, которое, как и ожидалось, увеличивается с c . Точно так же на каждом подфигуре рисунка 7 мы также сравниваем численную траекторию с решением возмущения, показанным зеленым цветом. В каждом случае решение возмущения проходит между (1, 0) и сначала пересекает ось V в точке, и мы показываем эту точку зеленым диском. Сравнение зеленой и синей траекторий в четвертом квадранте показывает, что у нас есть отличное совпадение численных решений и решений возмущений.Это сравнение показывает, что решение по возмущениям можно использовать для обеспечения высокоточной аппроксимации формы решений бегущей волны модели Фишера-Стефана для c <0,5. Используя уравнения (23) — (25), мы можем получить дополнительное аналитическое представление о взаимосвязи между c и κ , которую мы ранее численно исследовали на рисунке 5. Поскольку ордината точки пересечения равна V = — κ / c , мы можем построить приблизительные отношения в замкнутой форме между c и κ .Эти зависимости показаны на вставке к рисунку 5 для κ <2 и c <1. Сравнение вычисленных численно соотношений между c и κ с решениями возмущений показывает хорошее совпадение с ожидаемым результатом, что решение возмущения приводит к особенно точному приближению. Поскольку наши решения для возмущений обеспечивают хорошее приближение к кривой V ( U ) в четвертом квадранте фазовой плоскости, показанной на рисунке 7, теперь мы сравним точность решений для возмущений в системе координат бегущей волны.Для возмущающего решения имеем Интегрирование уравнения (26) с U = 0 при z = 0 дает неявное решение Для решения возмущения имеем для которого мы не можем найти точное решение. Тем не менее, уравнение (28) можно интегрировать численно, чтобы получить численное приближение U ( z ). На рисунке 8 показан набор решений бегущей волны, полученных путем решения уравнений (9) — (12) (выделены синим пунктиром), представленных для c = 0.05, 0,1,0,2 и 0,5. В каждом случае решения получаются в течение достаточно длительного времени, чтобы численные решения, зависящие от времени, превратились в бегущую волну. Эти бегущие волны затем сдвигаются так, что u ( x, t ) = 0 при z = 0, где z = x − ct . Для каждого значения c мы накладываем графики решения, заданные уравнением (27) (сплошным красным цветом). Результаты показывают, что решение отлично соответствует форме полных численных решений уравнений (9) — (12) для c = 0.05. Однако для больших c решение дает относительно плохое приближение. Аналогичным образом для каждого значения c мы также наносим на график решение (сплошным зеленым цветом). Здесь мы видим, что решение по возмущениям обеспечивает отличное согласование, неотличимое от полных численных решений уравнений (9) — (12) для c <0,1. На каждой части рисунка 8 мы также предоставляем вставку, показывающую увеличенную область сразу за передним фронтом волнового фронта, чтобы сделать сравнение между численным решением уравнения в частных производных и двумя решениями для возмущений более наглядным. В этой работе мы непосредственно сравниваем особенности решений бегущей волны известной модели Фишера-КПП и решений модели Фишера-Стефана. Ключевой особенностью модели Fisher-KPP является то, что любое положительное начальное состояние с компактной опорой всегда будет развиваться в бегущую волну, которая движется с минимальной волновой скоростью, c min = 2. Следовательно, согласно Fisher-KPP модели, любая начальная популяция приведет к успешному вторжению.Эта особенность является потенциальной слабостью модели Фишера-КПП, поскольку хорошо известно, что небольшие перемещенные популяции не всегда вторгаются и могут вымереть [46]. Напротив, модель Фишера-Стефана является адаптацией модели Фишера-КПП с подвижной границей, x = L ( t ). В модели Фишера-Стефана эволюция движущейся границы регулируется однофазным условием Стефана [39, 40, 41, 42, 43, 44, 45]. Модель Фишера-Стефана может поддерживать решения с бегущей волной со скоростью c <2.Поскольку и модель Фишера-КПП, и модель Фишера-Стефана поддерживают решения с бегущей волной, обе эти модели могут использоваться для моделирования процессов вторжения. Однако, в отличие от модели Фишера-КПП, модель Фишера-Стефана также предсказывает вымирание определенных начальных популяций, что приводит к дихотомии распространения-вымирания [39, 40, 41, 42, 43, 44, 45]. Дихотомия растекания-вымирания, связанная с моделью Фишера-Стефана, была изучена, в основном, с использованием строгих математических подходов, что привело ко многим важным результатам существования [39, 40, 41, 42, 43, 44, 45].Одна из целей этой работы — обеспечить более практическое сравнение моделей Фишера-КПП и Фишера-Стефана с использованием стандартных инструментов прикладной математики и инженерии, чтобы дать представление о сходствах и различиях между этими двумя моделями вторжения. Интересно отметить, что описания дифференциальных уравнений в частных производных моделей Фишера-КПП и Фишера-Стефана сильно различаются, поскольку модель Фишера-КПП связана с фиксированной областью, а модель Фишера-Стефана представляет собой задачу с подвижной границей.Напротив, когда мы анализируем решения бегущей волны обеих моделей, мы обнаруживаем, что уравнения, управляющие траекториями на фазовой плоскости, одинаковы. Интересно, что стандартные аргументы на фазовой плоскости для модели Фишера-КПП приводят нас к выводу, что решения с бегущей волной со скоростью меньше минимальной c min = 2 невозможны на основе физических аргументов и, следовательно, обычно не регистрируются. В этой работе мы показываем, что решения бегущей волны модели Фишера-Стефана требуют, чтобы c <2, и оказывается, что именно эти обычно отбрасываемые решения для модели Фишера-КПП имеют отношение к модели Фишера-КПП. Модель Стефана. В безразмерной модели Фишера-Стефана есть один свободный параметр, κ > 0, который связывает динамику движущейся границы, x = L ( t ), и пространственный градиент функция плотности на движущейся границе. Анализируя траектории на фазовой плоскости, связанные с решениями бегущей волны модели Фишера-Стефана, мы можем прийти к соотношению между κ и c , подтверждая, что c → 2 — как κ → ∞.Однако весь анализ фазовой плоскости модели Фишера-Стефана неявно предполагает, что решение бегущей волны сформировалось, тогда как численные решения полного описания уравнения в частных производных модели Фишера-Стефана показывают, что для фиксированного начального условия существует критическое значение κ cris : для κ > κ crit решение в конечном итоге превращается в бегущую волну, тогда как для κ < κ crit решение в конечном итоге исчезает.Зависящие от времени решения моделей уравнений в частных производных предполагают, что около этого перехода между окончательным исчезновением и возможным образованием бегущей волны мы имеем u ( x, t ) ≪ 1, предполагая, что мы можем получить представление, используя линеаризованную модель . Работая в линеаризованной структуре, мы получаем приближенное решение, из которого мы формируем заявление о сохранении, описывающее чистое накопление общей численности населения в области. В критическом случае, когда нетто-накопление массы в области, мы обнаруживаем, что существует критическая длина, L крит = π /2, и всякий раз, когда L ( t ) превышает π /2 решение всегда будет развиваться в бегущую волну, а если L ( t ) никогда не превысит π /2, плотность всегда в конечном итоге исчезнет.Используя сравнение с задачей Стефана, уравнения (9) — (12) с, мы можем усилить эти результаты так, что если L ( t ) превышает тогда сформируется бегущая волна. Или, если популяция вымирает с L ( t ) → L e как t → ∞, то Хотя хорошо известно, что не существует точных решений в замкнутой форме, описывающих решения бегущей волны уравнения Фишера-КПП для произвольного c , можно получить приближенные решения для возмущений для c ≫ 1 [4, 5].Поскольку решения бегущей волны для модели Фишера-Стефана движутся со скоростью c <2, мы построим решение возмущения для c ≪ 1, что приведет к приближенным решениям в замкнутой форме для формы траектории на фазовой плоскости, которая может быть используется для оценки формы бегущей волны. Мы обнаружили, что решения по возмущениям превосходно совпадают с нашими полными численными решениями для c <0,5, тем самым обеспечивая аналитическое понимание взаимосвязи между скоростью и формой решений бегущей волны модели Фишера-Стефана. Целью данной работы является сравнение моделей вторжения Фишера-Стефана и Фишера-КПП. Хотя мы начинаем нашу работу с указания на то, что модель Фишера-Стефана позволяет нам моделировать вымирание популяций, в то время как модель Фишера-КПП — нет, существуют также ограничения модели Фишера-Стефана, которые требуют признания и обсуждения. Например, зависящие от времени решения модели Фишера-Стефана, которые движутся в положительном направлении x , включая решения с бегущей волной, включают потерю популяции на свободной границе, x = L ( t ), так как ∂u / ∂x <0 при x = L ( t ).Может быть трудно оправдать эту потерю на движущейся границе на основе биологических, экологических или физических оснований и / или откалибровать модель для оценки соответствующего значения κ . Чтобы обратиться к этому моменту, стоит вспомнить, что модель Фишера-Стефана использует однофазное граничное условие Стефана, которое является упрощением более реалистичного двухфазного граничного условия Стефана [47, 48]. В более реалистичных применениях инвазии, таких как инвазия злокачественных клеток в окружающие ткани, на границе раздела будет происходить преобразование израсходованной ткани в злокачественные клетки [34, 60, 61, 62].Один из способов интерпретации этого преобразования из тканей в клетки состоит в том, что происходит потеря одного вида, в данном случае окружающей ткани, которая превращается в другой вид, в данном случае инвазивные клетки. Поэтому, хотя мы понимаем, что практическую интерпретацию потерь на движущейся границе в однофазной модели Фишера-Стефана трудно оправдать, мы ожидаем, что эти потери на движущейся границе будут очень естественными в расширенной структуре, где модель Фишера-Стефана Модель переделана как двухэтапная задача. Доступ к данным В этой статье нет дополнительных данных. Теория и доказательства из советов Рен-Роу Чена, Бо Лю, Сяолинь Чэн :: SSRN
Абстрактные
Возвращение к уравнению Фишера-КПП для интерпретации дихотомии распространения-поглощения
1 Введение
2 Результаты и обсуждение
2.1 Обезразмерение
2.2 Анализ фазовой плоскости
2.3 Связь между
κ и c 2.4 Критическая длина и дихотомия распространения-исчезновения
2.5 Решение возмущений, когда
c ≪ 1 3 Заключение