こんにちは、
ウィンドウでGraphBuilderとBivariateを使用してグラフを生成するスクリプトがあります。 機能を以下に示します。 紛らわしい場合は申し訳ありませんが、実際の列とデータテーブルの名前を非表示にする必要がありました。 基本的に、関数で呼び出されるデータテーブルには、行に異なるロット番号、列に異なるテストが含まれます。この関数は、時間に対する各テストパラメーターのグラフを生成し、線形回帰統計を含むグラフごとに2変量分析も生成します(線の方程式、Rの二乗、p値)。
setup = Function({dt}, {Default Local},
dt = Open ( "Product_example.JMP" );
dt << Select Where(:"Example_column" == "Example_value" );
cfg1 = EvalList({dt << Data View, { "List of columns would be here" }});
EvalList({cfg1});
);
build_graphs = Function({cfgs}, {Default Local},
win = New Window( "Product", vlb = V List Box() );
ncfg = NItems(cfgs);
for(i = 1, i <= ncfg, i++,
{this_dt, this_cols} = cfgs[i];
ncols = NItems(this_cols);
for(j = 1, j <= ncols, j++,
vlb << Append(H List Box (
this_dt << Graph Builder(
Show Control Panel( 0 ),
Size( 600, 400),
Variables(
X( :"Time" ),
Y( Column(this_dt, this_cols[j]) ),
Overlay( :"Lot Number" )
),
Elements(
Points( X, Y, Legend(12) ),
Line Of Fit( X, Y, Legend(10), Confidence of Fit(0), F Test(1))
),
SendToReport(
Dispatch(
{},
"Time",
ScaleBox,
{Format( "Fixed Dec", 10, 0 ), Min(0), Max(50), Inc(6),
Minor Ticks(0)}
)
)
),
this_dt << Bivariate(
invisible,
Y( Column(this_dt, this_cols[j])),
X( :"Time" ),
By( :"Lot Number" ),
Fit Line( {Line Color( {213, 72, 87} )} )
)));
)
);
win
);
dt = Open ( "Product_example.JMP" );
cfgs = setup(dt);
Show(cfgs);
win = build_graphs(cfgs);
Show(win);
win << Save Journal("Product_example.jrn", embed data(1));
各ロット番号のすべての線形回帰統計を要約する要約テーブルに似たものを作成できるかどうかを確認しようとしています。 たとえば、グラフビルダーを使用すると、F検定(p値)などの要素をグラフの上に表示できますが、この値をグラフの上ではなくテーブルに保存することはできますか? グラフビルダーから統計を抽出する方法がわからなかったため、現在、これらの統計を取得するために2変量解析を生成しています。
私が達成しようとしていることの例を示すために、別のソフトウェア「GraphpadPrism」からの画像を添付しました。 GraphPadでデータの分析を実行すると、一次方程式、Rの二乗、P値などを含む結果テーブルが生成されます。
これが必要かどうかはわかりませんが、フィットモデルを使用し、変数ごとにロット番号を使用すると、ロット番号ごとに回帰モデルが得られます。いずれかのモデルのパラメーター推定テーブルを右クリックして、[結合データテーブルに作成]を選択します。次に、グラフビルダーを使用して、さまざまな回帰モデル全体でパラメーター値、信頼区間、p値などをグラフ化できます。
うーん、これがうまくいくかどうかはわかりません。 Fit Modelを使用してみましたが、X変数のオプションがなく、Y byX分析を実行しようとしています。