2016-05-25 26 views
0

私はメタアナリシスの研究を行っており、Rのメタファパッケージを本当に好んでいます。 しかし、ネストされたdata.frameの多くの異なるモデルを管理するなど、Rでできることをすべて活用するにはHadley WickhamとDavid Robinsonのおかげです)私はmetaforパッケージrma() - broomのモデルからの出力を使用できるようにする必要があります。rma()出力を箒パッケージで使用する方法はありますか?

誰でもこれを行う方法はありますか。私はそれがまだ箒の箱からサポートされていない(まだ)ことを知っているが、それまでの回避策を非常に感謝するだろう。

多くの感謝!ここで

+1

### Example from metafor package: data(dat.bcg) dat <- escalc(measure="RR", ai=tpos, bi=tneg, ci=cpos, di=cneg, data=dat.bcg) mod <-rma(yi, vi, data=dat, method="REML") ### And call tidy() tidy(mod) 

がフィット統計の素敵なdata.frameを提供します:普通のようなtidy機能を使用して

adv-r.had.co.nz/OO-essentials.html)のHadley WickhamのAdvanced Rの本。彼は、特定の関数の追加メソッドを作成し、ディスパッチする方法について説明します。これにより、 'rma'オブジェクトを扱う' tidy.rma'メソッドを作成することができます。 –

答えて

1

rmaクラスのオブジェクトを扱うことができるtidyのための新しい方法を作成するための非常に基本的な例です。

tidy.rma <- function(x) { 
    return(data.frame(summary(x)$fit.stats)) 
} 

これは、関数tidy方法rmaを作成し、我々はフィットを含むdata.frameを返すように定義しますrmaオブジェクトの統計情報

代わりに、tidy.rma関数をdata.frameタイプのオブジェクトを、適切と思われるデータで返すようにすることができます。あなたが気にしているのは、あなたが好きなよ、お母さんの冗談を返してもらうことができ、それは入力のrmaモデルオブジェクトとはまったく関係がありません。 tidy.rmaメソッドは、それが何をするのかを決めるのはあなた自身です。 //:私は、[このセクション](HTTPのチェックアウトをお勧めし

   ML  REML 
ll -12.68777 -12.20237 
dev 37.16141 24.40474 
AIC 29.37554 28.40474 
BIC 30.50544 29.37456 
AICc 30.57554 29.73808 
関連する問題