現在、2つの異なるリンテスト値を相互に変換する方法に取り組んでいます。
バックグラウンド
土壌中の植物の利用可能なリンを測定するための多くの(抽出)方法が存在します。国によって適用される方法が異なるため、国ごとの出生率を比較するには、P検定値yに基づいてP検定値xを計算する必要があり、逆も同様です。したがって、応答と共変量は交換可能です。
抽出剤1のP量= [mg / 100g土壌]のP_CAL
抽出剤2のP量= P_DL [mg /土壌100g]
このような「変換方程式」を確立するために、136の土壌サンプルのP含有量をCALおよびDL抽出物で分析しました。土壌pH、全有機炭素、全窒素、粘土、炭酸塩などの追加パラメーターも測定されました。目的は、単純な回帰モデルを導出することです。2番目のステップでは、複数のモデルも使用します。
データの概要を説明するために、単純な線形(OLS)回帰直線を使用した2つの散布図を示します。
質問:
私の理解では、応答変更(y)変数と説明変数(x)の両方に(測定)エラーがあり、互換性がある場合は、デミング回帰が適しています。デミング回帰は、分散比が既知であることを前提としています。P抽出測定の精度に関する詳細はありませんが、分散比を決定する別の方法はありますか?ここでどの差異が意味されますか?計算されていないと思いますvar(DL_P)/var(CAL_P)
か?
Q1:デミング回帰の分散比を決定するにはどうすればよいですか?
デミング回帰の特殊なケースは、直交回帰です。分散比= 1と仮定します。
Q2:仮定δ= 1が「ほぼ」正しいか、または(誤った)仮定が高い予測エラーを伴うかどうかを診断する方法はありますか?
δ= 1と仮定すると、直交回帰は次の(丸められた)出力を提供します
library(MethComp)
deming <- Deming(y=P_CAL, x=P_DL, vr=1)
切片:0.75; 勾配:0.71; シグマP_DL:3.17; シグマP_CAL:3.17
上記のプロットでデミング回帰線をプロットすると、デミング回帰はa)CAL-P = f(DL-P)回帰に非常に近いが、b)DL-P = f(CAL-P)とは非常に異なることがわかります方程式。
Q3:正解です。直交回帰では、CAL-P = f(DL-P)とDL-P = f(CAL-P)は同じ式で表されますか?そうでない場合、どのようにして両方の正しい方程式を導出できますか?ここで何が恋しいですか?
両方の抽出ソリューションの特性により、DL-P値はCAL-P値よりも約25%高くなる傾向があるため、CAL-P = f(DL-P)はDL-P = f(CAL -P)。ただし、勾配が1つしかない場合、これはデミング回帰では表現されません。最後の質問です。
Q4:デミング回帰は私の目的にとって有効なアプローチですか?