如何計算GLMM fitme（）的p值 - 模型spaMM-Package？

Question

我使用了spaMM-package的fitme（）函數的GLMM。 該模型運行正常，但是沒有p值發射（見下文）：

Model1<-fitme(arten~ SGroesse*KF*DeckKF +(1|Hof) + Matern(1|Koordinaten_1 + Koordinaten_11), 
     data=kopfschlag1,family=negbin(),HLmethod="ML") 

summary(Model1) 
formula: arten ~ SGroesse * KF * DeckKF + (1 | Hof) + Matern(1 | Koordinaten_1 + 
    Koordinaten_11) 
Estimation of lambda, NB_shape, nu and rho by ML approximation (p_v). 
Estimation of fixed effects by ML approximation (p_v). 
Estimation of lambda, NB_shape, nu and rho by 'outer' ML, maximizing p_v. 
Family: Neg.binomial(shape=442.5) (link = log) 
------------ Fixed effects (beta) ------------ 
           Estimate Cond. SE t-value 
(Intercept)     1.219337 1.479982 0.82389 
SGroesse      0.680895 0.353169 1.92796 
KFMais      2.792657 2.509501 1.11283 
[the list goes on...]   ... 
--------------- Random effects --------------- 
Family: gaussian (link = identity) 
Correlation parameters: 
     nu  rho 
    16.66667 3494.08409 
Outer estimate of lambda (Hof): 1e-06 
Outer estimate of lambda (Koordinat.): 0.04556 
# of obs: 200; # of groups: Hof, 15; Koordinat., 200 
------------- Likelihood values ------------- 
         logLik 
p_v(h) (marginal L): -577.0237 
    p_beta,v(h) (ReL): -577.0237 

任何人可以幫助我，如何讓p值這個模型（lmerTest沒有工作，加上p值，似乎不能用spaMM包運行）？

另外，您是否有一個想法，如何計算fitme（）模型的事後測試（例如tukey-test）？

問候， A.

Answer 1

p值的混合效應模型往往是一個巨大爭論的話題/問題。如果你必須有一個p值，你有一個估計值，你有一個標準誤差，你有一個t值...選擇一個你可以出售的df並計算你的p值... https://stats.stackexchange.com/questions/45153/manually-calculating-p-value-from-t-value-in-t-test