これにはいくつか問題があります。線形フィット パラメータを引き出そうとするとエラーが発生し続けます。
標準の YvsX データセットがありますが、4 つの異なる「グループ化」列があります。 抽出しようとすると (さまざまな方法を試しました)、次のようなエラーが発生し続けます。
これで何が起こっているのですか？ 4つの列を1つの列に連結してから、その1つの連結で分割しようとしましたが、それでもエラーが発生します。
//make report
biv = Bivariate(
Y( :Y ),
X( :X ),
By( :A, :B, :C, :D</img> ),
Fit Line( {Line Color( "Red" )} )
);
//list the reports
rbiv = biv << report;
//count the number of reports
nBy = N Items( rbiv );
//array of as many rows as reports. 1col, all set to zero.
vc = J( nBy, 1, 0 );
//Create output table and fill with estimates
dtOut = New Table( "Regression parameters",
New Column( "Group", Character ),
New Column( "Intercept", numeric ),
New Column( "Slope", numeric ),
addrows( nBy ));
//grab the values
For( j = 1, j <= nBy, j++,
Estimates = rbiv[j]["Linear Fit"]["Parameter Estimates"][1][3] << get as matrix;
:Group[j] = groups[j];
:Intercept[j] = Estimates[1];
:Slope[j] = Estimates[2]);
