2015-06-02 17 views
0

我已經繪製從mcmcglmm包單一曲線中的R

> str(mcmc.list(results3[,2])) 
List of 2 
$ :Class 'mcmc' atomic [1:4000] 1.637 0.869 0.44 1.088 0.652 ... 
    .. ..- attr(*, "mcpar")= num [1:3] 1e+05 5e+05 1e+02 
$ :Class 'mcmc' atomic [1:4000] 1.068 0.472 0.851 1.549 0.927 ... 
    .. ..- attr(*, "mcpar")= num [1:3] 1e+05 5e+05 1e+02 
- attr(*, "class")= chr "mcmc.list" 
> 

和我可以繪製的曲線圖中的R由mcmcglmm包獲得以下對象看起來像這樣:

plot(results3[,2]) 

enter image description here

我希望能夠繪製右側的圖(分佈圖),但無法弄清楚(上週我管理的是什麼樣的真氣 - 沒有問題 - 但已經放棄瞭如何,並在過去一小時內嘗試過很多變體)。你知道這個解決方案嗎?

回答

0

跨解決方案迷迷糊糊的試驗和錯誤:

plot(density(c(results3[,2][[1]],results3[,2][[2]]))) 

(更優雅的答案將被接納爲答案)