How to use chi2gof in system identification white residual test?

1 次查看（过去 30 天）
Tarek Hajj Shehadi
Tarek Hajj Shehadi 2022-5-22
回答： Rajiv Singh 2022-5-31
I have a BJ model that I am trying to validate using flexibility tests. The chi2gof sounds like a good option.
Goal: I want to apply white null hypotheses test to check whether the residual is white or not.
I will show you how I applied it:
First create the BJ model using training data
BJ_Mod=bj(train_dat,[8 7 1 1 0]);
Now we compute the correlation of the residual on the validating data
res_BJ = resid(val_dat,BJ_Mod); %val_dat (contains 411 points)
res_out_BJ = res_BJ.OutputData;
[res_corr_BJ, lags_BJ] = xcorr(res_out_BJ, 'unbiased');
corr_0_BJ = max(res_corr_BJ); %Compute residual at origin
Now define I wish to now apply chi2 test on
[chi_gof,p,stat] = chi2gof(res_corr_BJ/corr_0_BJ,'Alpha',0.05);
I received the following as output:
chi_gof = 1
p=9.0250e-05
chi2stat: 18.6259
df: 2
edges: [-0.2853 -0.1567 -0.0282 0.1003 0.2288 1.0000]
O: [68 222 372 154 5]
E: [63.4314 255.3488 341.3373 142.0305 18.8521]
Question: I am afraid that the results might be incorrect based on the way I implemented the algorithm. What I received from the output (as shown above) that my model fails this hypotheses. The main concern is whether this algorithm makes sense. I hope someone can check this code and confirm whether this is the way to do white null test using chi2 criteria on a system identification model.
Furthermore, do we apply this test on just the validating points or on all the data points of the set?

请先登录，再回答此问题。

采纳的回答

Rajiv Singh
Rajiv Singh 2022-5-31
The approach looks fine to me in principle, hard to say more without looking at the data. One thing you can try is call RESID without output arguments. This will plot the residual with confidence bounds based on chi2 stats on the prediction errors. If the residual curve is within the bounds, you basiclaly have uncorrelated, white residuals.

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Linear Model Identification 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by