p値とベータから標準誤差を計算するR

Question

私は2つのモデルの出力を持っています - 最初のモデルはglmを使って実行し、2番目のモデルはlmekinを使って実行しました。両方のモデルの入力は同一であったが、lmekinモデルは親密行列に適合した。

私はglmモデルから係数を取り、lmekinで実行されたモデルからのp値を使用してこれらの値のSEを計算します。

これは可能ですか？スケーリングされていない共分散行列の対角要素の平方根を使って、SEを手動で計算することができます。したがって、次のコードは私が達成したいことで成功するでしょうか？

beta <- glm_model$coefficients 
nvar <- length(beta) 
nfrail <- nrow(lmekin_model$var) - nvar 
se <- sqrt(diag(lmekin_model$var)[nfrail + 1:nvar]) 

Answer 1

これは私にとってうまくいくようです。異なる方法と追加の情​​報を使用する別のモデルのあるモデルの推定された（そして計算された）SEを使用することは、うまくいかないようです。

しかし、最初にいくつかの人に不快感を与えたことの多くは、正確で有用であると証明されています。

あなたのモデルに適したデータをシミュレートすることをお勧めしますが、「真実」が何であるかをシミュレートしています。それから、あなたの分析を行い、あなたのコードを試して、それが何を表示しているかを見てくださいその後、プロセス全体を数回繰り返して、分析間のばらつきを確認します。あなたのメソッドがシミュレートされたデータに対して機能しない場合は、実際のデータに対して信頼できないことがわかります。しかし、シミュレートされたデータを処理すると、実際のデータが現実的であるという確信が得られます。