Здравствуйте, статистические гуру и мастера программирования R,
Я заинтересован в моделировании захвата животных в зависимости от условий окружающей среды и дня года. Как часть другого исследования, у меня есть подсчеты по ~ 160 дней за три года. В каждый из этих дней у меня есть температура, количество осадков, скорость ветра, относительная влажность и т. Д. Поскольку данные неоднократно собирались с одних и тех же 5 графиков, я использую график как случайный эффект.
Насколько я понимаю, nlme может легко учитывать временную автокорреляцию в невязках, но не обрабатывает функции негауссовой связи, такие как lme4 (которая не может обрабатывать автокорреляцию?). В настоящее время я думаю, что это может работать с использованием пакета nlme в R в журнале (количество). Так что мое решение сейчас было бы запустить что-то вроде:
m1 <- lme(lcount ~ AirT + I(AirT^2) + RainAmt24 + I(RainAmt24^2) + RHpct + windspeed +
sin(2*pi/360*DOY) + cos(2*pi/360*DOY), random = ~1|plot, correlation =
corARMA(p = 1, q = 1, form = ~DOY|plot), data = Data)
где DOY = день года. В окончательной модели может быть больше взаимодействий, но это моя общая идея. Я также мог бы потенциально попытаться смоделировать структуру дисперсии с чем-то вроде
weights = v1Pow
Я не уверен, есть ли лучший способ сделать с регрессией смешанной модели Пуассона или что-нибудь? Я только что нашел математическое обсуждение в главе 4 «Модели регрессии для анализа временных рядов» Кедема и Фокианоса. На данный момент это было немного за пределами меня, особенно в приложении (кодирование в R). Я также видел решение MCMC в Zuur et al. Книга «Модели смешанных эффектов» (гл. 23) на языке BUGS (с использованием winBUGS или JAG). Это мой лучший вариант? Есть ли в R простой пакет MCMC, который бы справился с этим? Я не очень знаком с методами GAMM или GEE, но был бы готов изучить эти возможности, если бы люди думали, что они дадут лучшее понимание.Моя главная цель - создать модель для прогнозирования отлова животных с учетом условий окружающей среды. Во-вторых, я хотел бы объяснить, на что реагируют животные с точки зрения их активности.
Мы будем благодарны за любые мысли о том, как лучше (философски) действовать, как закодировать это в R или в BUGS. Я довольно плохо знаком с R и BUGS (winBUGS), но учусь. Это также первый раз, когда я пытался решить проблему временной автокорреляции.
Спасибо Дэн