問題タブ [broom]

R で plot_ly() の周りにカスタム関数を書きたいと思います。そうすれば、同じフォーマットとスタイルで一連の散布図を作成できますが、コードは重複しません。このページをガイドとして使用しました。このコードはエラーを再現します:

問題の行は次のとおりです。

as.formula()を使ってみましたが、エラーメッセージは似ています。

エラーメッセージは次のとおりです。

model.frame.default(formula = y ~ x, data = mtcars, drop.unused.levels = TRUE) のエラー: オブジェクトは行列ではありません

関数ではない場合、コードは機能します。

ここで説明されているように、ネストされたデータ フレーム内で回帰を実行しようとしています。私の目的では、多くのレベルの固定効果があるためfelm、パッケージから使用しています。lfe

上記のリンクの例を のfelm代わりに使用してやり直すlmと、 を使用するまではほとんどの場合機能しますbroom::augment。

コードの最後の行で数式の代わりに関数を使用する必要があったことを除いて、すべてがこの時点まで機能します。purrr::mapおそらく別のfelm癖です。

broomモデル出力を抽出するために使用しようとすると、glanceとでは機能しますが、 では機能しtidyませんaugment。

使用しようとするaugmentと、次のエラー メッセージが表示されます。

グループ化変数によってタグ付けされた 3 つの従属変数と 6 つの独立変数のセットを含むデータ フレームがあるとします。この形式の例は、以下のサンプル コードで生成されます。

x1 から x6 までのセットの y1、y2、y3 のそれぞれを回帰させたい場合は、次のようなものを使用できます。

(lm() ヘルプの次の行を利用することにより: 「応答が行列の場合、行列の各列に最小二乗法によって線形モデルが個別に適合されます。」)

ただし、最初にグループ化変数でグループ化し、次に lm 関数を適用する必要がある場合、その方法がよくわかりません。以下を試してみましたが、両方のグループで同じ係数セットが生成されます。