http サーバーからの断続的なエラー応答をテストして定量化しようとしています。必要な要求を行うことができる curl 呼び出しを作成しましたが、上向き矢印を入力して戻るのにうんざりしています。

エラーが発生するまで繰り返しリクエストを作成し、エラーが発生する前の成功応答の数を出力するようにcurlに指示するにはどうすればよいですか?

現在、次の構造を持つ r のデータセットに対して lmer() を実行しています。

固定効果に関する次の要約を受け取ります。

次に、深さの中央値 (14) のこの方程式を使用して、対数変数とスケーリングされた変数から実際の値の推定値を予測します...これは、下のグラフの正しい実際の値の線です。

lmer 実質値の推定値

しかし、これらの式を使用して +/- 1 つの標準誤差行を追加すると、次のようになります。

マイナス SE

プラス SE

私はこのグラフを持っています: lmer 実際の見積もり +/- 1 SE

すべての線の交点は x 値 0.068341709 にあります。

SE ラインが交差するのはなぜですか? ここで何が起こっているのか、どうすれば修正できますか?

私は、1990 年から 2010 年の間に米国の国勢調査区の空間パネル モデル (splm を使用) を、距離または隣接性に基づく隣接行列で実行しています。しかし、国勢調査区が含まれる都市の標準誤差もクラスター化したいと考えています。

パネル モデル (plm) 内でこれを行う方法があります (事後推定関数「vcovHC」と「coeftest」(こちらを参照) を使用)。 splm オブジェクト。同じコードを実行すると、エラー"Error in terms.default(object) : no terms component nor attribute."が表示されます。他のモデルでのクラスタリングについては多くの質問がありますが、splm パッケージでは質問がないので、助けていただければ幸いです。

これは、XSLME パッケージを使用して Stata で実行でき、呼び出しの最後に vce(cluster "group_var") を含めますが、3 つのパネルで 37k 以上のケースがあると、Stata がクラッシュします。

同じ問題の架空の例を以下に示します。

線形回帰から変換された変数の標準誤差を計算したいと思います。つまり、2 つの変数を除算し、この変数から標準誤差を取得します。

deltamethodパッケージの関数を使用していmsmますが、正確な標準誤差を取得できません。

例えば：

データのシミュレーション:

線形回帰:

coef と vcov を取得します (ここでは、実際のデータの一部の係数が NA であり、ループで多くの回帰を計算するため、欠落を取り除く必要があります)

変数 x1 を x3 で割った Deltamethod (つまり、msm パッケージに従って x2 と x4 を使用する必要があります):

これにより、変換された変数 (x1/x3) の標準誤差は 3.21 になりますが、この回帰によるすべての標準誤差は約 0.03 です。

ここでなぜ/何が間違っているのですか？それを計算するための他の提案も大歓迎です。

誰かが次の問題で私を助けてくれることを望んでいました:

バイナリ変数 (性別) に対して記録された 3 つの異なる連続変数 (体温、長さ、質量) の平均と標準誤差を示す結合棒グラフを作成しようとしています。

各変数の平均値をプロットできましたが、試したコードを使用してこれら 3 つの変数の標準誤差をうまく計算できないようです。私は多くのことを試しましたが、これで正しい軌道に乗っていたと思います：

その後、エラーが発生し続けました：

集約（test4 ...）の後に3つの変数をターゲットにするように再コーディングしようとしましたが、機能しませんでした...次に、結果のデータフレームでサブセット化して性別を除外しましたが、機能しませんでした。次に、それを行列またはベクトルとして定義しようとしましたが、それでもうまくいきませんでした。

最終的なグラフで、y 軸 = 平均値、x 軸 = 変数 (3 つのサブグループ (Tb、Mass、Length)、比較のために男性と女性の値を示す 2 つのバーを並べて表示したいと思います。

誰かが提供できる助けや方向性は大歓迎です!!

よろしくお願いします！:)

Fama-Macbeth 回帰を推定しています。このサイトからコードを取得しました

fpmg を推定した後、二重クラスタリングでロバストな SE を推定します。

次のエラーが表示されます。

このエラーを解決する方法を提案してください。

更新: 「multiwayvcov」パッケージも試しましたが、同じエラーが表示されます。これらのパッケージ (Sandwich、multiwayvcov など) ではオブジェクト クラスが許可されていないようです。Rは本質的に私のすべての労力を無駄にするようで、私は行き止まりにぶつかりました。上記をPythonで行う方法を見つけました（コードを意味します）が、それについての知識はありません。

Rで問題を解決する方法はありませんか?

下の表 (平均値 ± se) には、出版物に掲載される真のデータが含まれています。ご覧のとおり、データには非常に広い範囲があります。R を使用して、すべての数値を有効数字 3 桁に丸めましたが、テーブルを整然と保つために、末尾のゼロと小数点の左側のゼロを有効なものとして含めたいと考えています。したがって、基本的には、小数点の両側にあるゼロを含むすべての数値を有効と見なして、数値に関係なく3桁のみが出力されるようにします。これは可能ですか？、、、およびを試しましsignifたが、成功しませんでした。これらの結果は、roundsprintfoptions()formatCx$summary = paste(signif(x$Value, digits=3), "\u00b1", signif(x$se, digits=3))