能否通过贝叶斯的credible interval比较不同水平间的差异

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

各位坛友好！


这里有一个分层模型的问题，想请教一下大家。

小弟的试验中，响应变量是RB，固定因素是treatment（四个水平为ck tr1 tr2 tr3），随机因素是year和block。我想分别比较ck和tr1 tr2 tr3的差异，即不关注tr1 tr2 tr3 三者之间的差异。

我用R做的，不过我想对于结果的解释应该跟软件无关。

我想请教一下：
（1）能否直接根据summary的结果，得出“tr1和ck的差异不显著，tr2 tr3和ck的差异显著“的结论？
（2）我能使用附图这样的方式来表示这种差异么？
（3）方法中我能用类似下面的话描述这种统计方法么？“ In the model, the RB under "ck" was set as the intercept, so for the other three conditions ("tr1","tr2" and "tr3"), 95% credibility interval that not overlapping zero indicate a significant difference from that in "ck" "

1. library(MCMCglmm)
   
2. fit1<-MCMCglmm(RB~Treatment,random=~Year+Block,data=tempRB)
   
3. 
   
4. summary(fit1)
   
5. 
   
6. Iterations = 3001:12991
   
7. Thinning interval  = 10
   
8. Sample size  = 1000
   
9. 
   
10. DIC: -14.23551
    
11. 
    
12. G-structure:  ~Year
    
13. 
    
14. post.mean l-95% CI u-95% CI eff.samp
    
15. Year    0.3482  6.4e-56   0.1674      436
    
16. 
    
17.            ~Block
    
18. 
    
19. post.mean  l-95% CI  u-95% CI eff.samp
    
20. Block 0.0001119 1.156e-66 3.083e-08    798.7
    
21. 
    
22. R-structure:  ~units
    
23. 
    
24. post.mean l-95% CI u-95% CI eff.samp
    
25. units   0.02743  0.01134    0.048      114
    
26. 
    
27. Location effects: RB ~ Treatment
    
28. 
    
29.               **post.mean  l-95% CI  u-95% CI eff.samp  pMCMC   
    
30. (Intercept)    0.660959  0.391469  0.878285     1000  0.024 *  
    
31. Treatment_tr1 -0.004188 -0.207064  0.179150     1000  0.966   
    
32. Treatment_tr2 -0.387922 -0.586966 -0.186081     1085 <0.001 ***
    
33. Treatment_tr3 -0.467356 -0.646885 -0.261730     1000 <0.001 *****
    
34. ---
    
35. Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

复制代码

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：Interval inter Inte 不同水平 cred