1

ツイン モデリングを行うために数年前に記述されたOpenMxスクリプトの例があります。

OpenMx バージョン 1.0 用に作成されました (スクリプトはここにリンクされています) 。

実行すると、適合関数と目的関数の更新に関するいくつかの警告が表示されます。OpenMx 2.0 fit 関数呼び出しを使用するには、どのように更新すればよいですか?

4

1 に答える 1

2

OpenMx 1.0 から 2.0 以降では、いくつかの変更があります。ほぼすべてのスクリプトは問題なく動作しますが、OpenMx 2.x に更新すると、一部の 2012 年以前のスクリプトや機能が改善されます。

参考例はこちら

ユーザーは次のような問題を抱えていました。

1. ヘルパー関数へのパスがない

これは、R コードなどのより一般的な堅牢性の問題です。ディスクベースのファイル パスではなく、Web URL を含めることをお勧めします。

source("http://www.vipbg.vcu.edu/~vipbg/Tc24/GenEpiHelperFunctions.R")

より良い解決策は、 umxのような CRAN ベースのヘルパー パッケージです。これらは、最新の状態に保ちやすく、アクセスしやすいものです。

2. 古いスタイルの目標 (期待値と適合関数の代わり)

このような呼び出しは非推奨です。

objMZ<- mxFIMLObjective(covariance="expCovMZ", means="expMean", dimnames=selVars)

+ への呼び出しに置き換えmxFIMLObjectiveて、スクリプトのスタック全体でこれらを簡単に更新できます。mxExpectationNormalmxFitFunctionML

さらに、次のような古いスタイルの複数グループ目標では:

minus2ll <- mxAlgebra( expression = MZ.objective + DZ.objective, name="m2LL")
obj  <- mxAlgebraObjective("m2LL")

と交換する必要がmxAlgebraObjectiveありますmxFitFunctionAlgebra

ただし、OpenMx 2 には、これを 1 行で処理し、識別チェック、参照モデルの生成などを可能にするきちんとしたマルチグループ機能があります。

したがって、全体を（たとえば）次のように置き換えてください。

mxFitFunctionMultigroup(c("MZ", "DZ"))}
于 2016-02-23T13:50:14.360 に答える