Вероятно, это более общий вопрос о том, как кодировать уравнения в R, но я применю концепцию к тому, что было в моем списке дел.Мы изучаем сезонность дождевых осадков месячного временного ряда осадков, который показывает, как один месяц в данном году отличается от осадков, если бы каждый месяц получал 1/12 годового количества осадков.Уравнение можно найти здесь: http://leddris.aegean.gr/ses-parameters/293-rainfall-seasonality.html, но переписано здесь: сезонность осадков
Где Si - годовой индекс сезонности, Ri - общее годовое количество осадков для конкретногогод изучается, и Xin - это фактическое месячное количество осадков за месяц n.
Наши данные имеют следующую структуру и продолжаются ежемесячно в течение всего водного года 2014
head(monthlyRain)
MonthofYr WaterYear monthlyRain
1 1983 289.2
2 1983 150.8
3 1983 134.3
4 1983 16.1
5 1983 490.3
6 1983 328.9
Как упоминалось выше, это болееобщий вопрос о том, как решить слегка сложные уравнения в R, используя пример, который я имею под рукой.Дайте мне знать, если это кросс-пост или есть другие ссылки для решения без лишних здесь.Спасибо!