Как проверить правильность коэффициента смертности

Этот вопрос связан с вопросом Я спрашивал некоторое время назад.

Я пытаюсь создать популяцию в NetLogo. Поэтому я пытаюсь начать со смертности. Однако я не уверен, как проверить, работает ли мой код. Возраст самых старых умирающих черепах, кажется, в порядке, однако, когда я печатаю средний возраст смерти, он слишком высок, и я не знаю, как это решить.

Модель доступна здесь.

Проблема может заключаться в расчете вероятности смерти, который я сделал следующим образом:

  if (gender = "female") [
    let die-val mortality-female age ; mortality prob included for the age of the turtle
    if (die-val >= random-float 1)[
       set prob-die? true
    ]
  ]

person Hannah H.    schedule 08.07.2018    source источник


Ответы (1)


Похоже, ваш код работает нормально, но ваша реализация может быть нереалистичной. Все ваши черепахи в конце концов умрут, и поскольку вероятность вашей смертности увеличивается с возрастом, ваши старшие черепахи вносят гораздо больший вклад в список возрастов смерти, чем более молодые черепахи. Например, после бега до тех пор, пока всем черепахам не исполнится 90, в живых осталось 52 черепахи, поэтому все они (более половины) внесут свою старость в список death-ages. Итак, ваш список death-ages сильно перекошен влево, что-то вроде:

введите здесь описание изображения

Где ваш средний death-age более 100 повторений составляет около 90,5 лет. Я думаю, это связано с вашим процессом старения — попробуйте увеличить возраст черепахи только на 1 вместо 5 — посмотрите, даст ли это результат больше, чем вы ожидаете. Это немедленно изменяет распределение, чтобы оно выглядело более нормальным, где ваше среднее значение death-age для 100 повторений составляет около 81,2 года:

введите здесь описание изображения

Изменить:

Вы можете быстро настроить гистограмму с помощью виджета графика с такими настройками:

введите здесь описание изображения

Обратите внимание, что код

if length death-ages > 0 [
   set-plot-x-range ( min death-ages - 5 ) ( max death-ages + 5)
]

включен, потому что автоматическое масштабирование не работает для x-шкалы гистограмм. Вы можете исключить этот код в поле «Команды обновления графика» и вручную установить x min и max, если хотите, но вы пропустите черепах, которые умирают после любого значения x max, которое вы выберете.

Вам также необходимо изменить настройки пера с линии на полосу — щелкните значок карандаша справа от поля пера, чтобы открыть меню ниже, и выберите «Панель» в раскрывающемся меню «Режим». Также обратите внимание, что я установил свой «Интервал» на 5, чтобы получить графики, как показано выше.

введите здесь описание изображения

person Luke C    schedule 09.07.2018
comment
Спасибо. Завтра проверю, работает ли. - person Hannah H.; 10.07.2018
comment
Можно поинтересоваться, как вы рисовали сюжет? Вы использовали гистограмму? - person Hannah H.; 12.07.2018
comment
@Hannah - Да, точно - я отредактировал ответ, включив в него настройку гистограммы. - person Luke C; 12.07.2018