Порядковый номер MCMCglmm в R - PullRequest
Новая
       18

Порядковый номер MCMCglmm в R

0 голосов
/ 02 февраля 2020

Я пытаюсь объяснить изменения жизнеспособности дерева с 1 до 3 (1 = зеленый, 2 = урон, 3 = dry), используя климатические переменные c в модели MCMCglmm. К сожалению, я борюсь с двумя вопросами:

1. Как мне интерпретировать сводку модели MCMCglmm?

Я вижу, какие переменные значимы, но что это значит? Означает ли это, что с увеличением количества осадков жизнеспособность достигает 1 (зеленый)? Я построил заднюю часть модели и, похоже, что осадки и диаметр не влияют на жизнеспособность дерева, но виды реагируют по-разному.

2. Как проверить правильность модели?

Структура данных: вид, диаметр (DBH), жизнеспособность 2018 (вит_2018) и количество осадков (стр. 18). Были сделаны участки, 1 участок содержит 10 замеренных деревьев, количество осадков составляет одно значение на участок. 
###MCMCglmm model, family= ordinal

    prior1<-list(R=list(V=diag(1),nu=0.002))
    m1 <-MCMCglmm(VIT_2018~ pp18 +DBH+Species, 
                                 family = "ordinal", data = comsp,prior=prior1,pr = TRUE, 
                                 nitt = 60000, burnin =30000, thin = 50)
    summary(m1)

    #Output:
    Iterations = 30001:59951
     Thinning interval  = 50
     Sample size  = 600 

     DIC: -146008.7 

     R-structure:  ~units

          post.mean l-95% CI u-95% CI eff.samp
    units      2697     1016     5304    6.075

     Location effects: VIT_2018 ~ pp18 + DBH + Species 

                            post.mean  l-95% CI  u-95% CI eff.samp   pMCMC   
    (Intercept)              60.64229  31.71644  92.78733    11.49 < 0.002 **
    pp18                     -0.07020  -0.11518  -0.03556    15.47 < 0.002 **
    DBH                       0.17357  -0.01797   0.37281    58.26 0.06000 . 
    SpeciesBetula pendula   -12.88510 -29.42191   1.87807    58.40 0.08667 . 
    SpeciesCarpinus betulus  15.57570   0.96439  30.89528    45.04 0.02000 * 
    SpeciesCorylus avellana  -3.81337 -18.31965  13.12771   600.00 0.59000   
    SpeciesCrataegus spec.  -14.90077 -36.86880   4.24729   286.00 0.10333   
    SpeciesFagus sylvatica  -15.03559 -29.04547  -2.07809    60.56 0.00667 **
    SpeciesFrangula alnus    20.10817  -0.10598  38.66354    73.88 0.01333 * 
    SpeciesQuercus spec.     -9.09458 -24.52595   6.42502   293.33 0.26000   
    SpeciesSambucus nigra    24.29894   2.58339  46.29901    49.66 0.02333 * 
    SpeciesSorbus aucuparia  39.56930  22.97282  63.15175    10.13 < 0.002 **
    ---
    Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

     Cutpoints: 
                             post.mean l-95% CI u-95% CI eff.samp
    cutpoint.traitVIT_2018.1     80.46    55.38    117.9    4.204
#### Estimating Credible Intervals

    HPDinterval(mcmc(randomprior1$Sol[,"(Intercept)"]))
    #        lower    upper
    #var1 31.71644 92.78733
    #attr(,"Probability")
    #[1] 0.95

Данные и R скрипт можно найти здесь

https://drive.google.com/drive/folders/1LmgEAssR5FfFw1CkYjygsaawg84dDAwk?usp=sharing Здесь прилагаются части данных, когда я запускаю модель со всеми данными, плотность графики выглядят лучше.

PS Я также попробовал функцию clm из порядкового пакета в R, но псевдо r ^ 2 по Макфаддену было 0,06 = не подходит. Поэтому я проверил данные об осадках и местоположение участков с помощью теста Морана I, и была пространственная автокорреляция. Я использую MCMCglmm, потому что я прочитал, что он может иметь дело с пространственной автокорреляцией. До сих пор я следовал нескольким учебникам или вопросам, которые я нашел в Интернете

Любые предложения будут действительно полезны.

Спасибо за ваше время и помощь!

...