R で縦断的項目応答理論 (IRT) モデルを適合させようとしています。複数の測定機会で実施されたテストがあります。段階的反応モデル (GRM) から、因子スコア (つまり、能力レベル) の個人の成長曲線を調べたいと思います。R の ltm パッケージを使用して、IRT の断面 GRM モデルに適合させましたが、モデルを拡張して同じ項目の繰り返し測定を処理する方法 (または ltm で可能かどうか) が明確ではありません。経時的な能力レベルの平均/分散の変化を確認するために、経時的な GRM 因子スコアに成長曲線を当てはめるにはどうすればよいですか? これが ltm パッケージで不可能な場合、これを許可するパッケージ/機能は何ですか? 特定のコード例は特に高く評価されます。
これは、私がやろうとしていることに似た経験的な例です（例では、多項データにGRMを使用しているバイナリアイテムにRaschモデルを使用しています）：