平行線モデルで相対効力を計算するにはどうすればいいですか?こんにちは、
近似曲線に平行線モデル LogY=Sample+Log_Dose を近似しました。
並列性テストとモデルパラメータは問題ありません。RやPLAソフトウェアで得られたものと全く同じです。しかし、相対的な効力は全く異なります。計算に関する技術的な詳細を知りたい方はいらっしゃいますか？ありがとうございます！
高度な統計モデリング基本的なデータ分析とモデリングRe: 平行線モデルで相対効力を計算するにはどうすればいいですか?小さなテスト データ テーブル、使用されたさまざまなコード、および結果の概要があれば、調査が容易になります。
Re: 平行線モデルで相対効力を計算するにはどうすればいいですか?予測および特殊モデリング > 曲線の適合 > モデル適合オプション > テストの並列性にはいくつかの情報が含まれています
Re: 平行線モデルで相対効力を計算するにはどうすればいいですか?PLAの例を使用しました。RのRPは2.170981で、PLAソフトウェアで得た値と同じです。しかし、JMPの結果は5.9です。
推定パラメータは同じなので、両方とも同じモデルを使用していると思いますが、JMP の RP が大きすぎるだけです。
データ <- read_excel("PLA Example_parallel line.xlsx", sheet = "Sheet1")
データ$Log_Dose = log(データ$Dose)
データ$LogY = log(データ$Y)
fit <- lm(LogY ~ Sample+Log_Dose , data = data) #平行線モデル: 共通傾き
要約（フィット）
係数 <- coef(fit)
db <- 係数["サンプルT"] #
ベータ <- 係数["Log_Dose"]
LogRP <- db / ベータ
RP<-exp(logRP)
JMPでは、特殊モデリング→曲線のフィット、多項式→線形フィット、並列性のテストでモデルをフィットしました。
Re: 平行線モデルで相対効力を計算するにはどうすればいいですか?こんにちは@アクトン 、
違いが見られる理由は、JMP の平行線テストの推定値が S の「ベースライン」との差を取得せず、代わりに各項を表示しているためです。スクリーンショットと同じ値を取得できます (サンプル グループにダミー コーディングを使用しています)。
- 切片 = 5.71347 (Sの場合)
- サンプルT = 0.70423 → Tの切片 = 5.71347 + 0.70423 = 6.4177 (JMP表と一致)
- サンプルU = 0.51262 → Uの切片 = 6.2261
- サンプルV = 0.61602 → Vの切片 = 6.3295
標準誤差の推定値は、JMP が共分散行列から推定値を取得するため異なりますが、スクリーンショットの推定値は生の推定値のように見えます。
ありがとう、
ベン
Re: 平行線モデルで相対効力を計算するにはどうすればいいですか?こんにちは、ベン。
JMPの切片と傾きの推定値には問題ないのですが、相対効力が正しくないようです。平行性検定をクリックして相対効力の部分を表示しましたが、平行線モデルの相対効力なのか分かりません。この分析のJMPヘルプファイルが見つかりません。
ありがとう、
アクトン
