問題タブ [longitudinal]

私は、処理のための次の変数 (栄養、肥料) を含むデータセットを持っています。これは、水中の藻類の成長を時間 (t0、t1...t10) にわたって記録します。「窒素」と記された肥料と連続して、ｔ５日後に窒素を添加した。「なし」とマークされたシリーズでは、窒素は添加されませんでした。

グループ間の勾配を比較し、肥料介入後に勾配が変化するかどうかを確認したいと考えています。コントロールとして肥料を使用したくない (つまり、(良い、なし) または (悪い、なし) コントロールとして)

このデータを、「栄養」、「肥料」、「時間」、「複製」、「成長」などの列ヘッダーを持つ長い形式に変換します。t5 より前の期間と t5 より後の期間を区別するために、加算と呼ばれる新しい列を作成します。t5 前の時間 -> 0、t5 の後 -> 1

各列が次の構造を持つ縦方向の分析を実行します。

lmer(成長 ~ 栄養 + 肥料 + 時間 + 追加 + (1|複製))

固定効果モデルのランクが不足しているため、x 列数を削除するというエラー メッセージが表示されます。とにかくこの問題の周りにありますか? モデルの記述方法を改善するための提案はありますか?

基本的にデータフレーム内の別の変数の開始値である新しい変数を作成しようとしています。サンプルデータ:

私が取得したいのは、 for の値に等しく、また に依存する new_var (またはその他のもの) という名前のoutcome新しいtime == 0列id = idですAttitude。したがって、私が拡張したいのdataframeは次のとおりです。

その後、適切なコーディングを行います。SASでは、関数を使用してこれを実行できることを知っていlagます。「回避策」ではないソリューションを本当に感謝しています。SASのようなものですが、適切なrソリューションです。最終的にはrも強くなりたいです。

関連: SAS としての R の保持と遅延機能 しかし、私は、インデックスまたは「通常の」r の方法に基づくいくつかのソリューションを好みます。そしてここでも、他の条件に依存しません。

したがって、ここで重要なのは、コーディングがさまざまなidsレベルattitude/変数 (A1、A2、...) に対して機能し、outcome valueattime == 0が基本的に にコピーされることnew_varです。

メッセージを明確に伝えることができれば幸いです。そうでない場合は、コード例の小さな部分とそれをどのように拡張したいかが十分に明確であると思います。提案を楽しみにしています。

EDIT @jogo回答の別のコード例。