残差分析代码

残差分析图

figure(2002)

Y2=y2';

X2=[ones(length(y2),1),t2',h2',p2'];

[b2,bint2,r2,rint2,stats2]=regress(Y2,X2)

rcoplot(r2,rint2)

figure(2003)

Y3=y3';

X3=[ones(length(y3),1),t3',h3',p3'];

[b3,bint3,r3,rint3,stats3]=regress(Y3,X3)

rcoplot(r3,rint3)

figure(2004)

Y4=y4';

X4=[ones(length(y4),1),t4',h4',p4'];

[b4,bint4,r4,rint4,stats4]=regress(Y4,X4)

rcoplot(r4,rint4)

figure(2005)

Y5=y5';

X5=[ones(length(y5),1),t5',h5',p5'];

[b5,bint5,r5,rint5,stats5]=regress(Y5,X5)

rcoplot(r5,rint5)

figure(2006)

Y6=y6';

X6=[ones(length(y6),1),t6',h6',p6'];

[b6,bint6,r6,rint6,stats6]=regress(Y6,X6)

rcoplot(r6,rint6)

figure(2007)

Y7=y7';

X7=[ones(length(y7),1),t7',h7',p7'];

[b7,bint7,r7,rint7,stats7]=regress(Y7,X7)

rcoplot(r7,rint7)

figure(2008)

Y8=y8';

X8=[ones(length(y8),1),t8',h8',p8'];

[b8,bint8,r8,rint8,stats8]=regress(Y8,X8)

rcoplot(r8,rint8)

figure(2009)

Y9=y9';

X9=[ones(length(y9),1),t9',h9',p9'];

[b9,bint9,r9,rint9,stats9]=regress(Y9,X9)

rcoplot(r9,rint9)

b2 =

1.0e+003 *

2.1285

-0.0022

-0.0027

-0.0018

bint2 =

1.0e+003 *

-4.8718

9.1287

-0.0080

0.0036

-0.0055

0.0002

-0.0085

0.0048

r2 =

-23.9155

2.1988

0.8720

5.3871

0.2913

6.7232

-5.1657

-2.4146

-0.4199

-8.6214

8.8987

16.1659

rint2 =

-35.8274

-12.0036

-21.3997

25.7972

-24.6686

26.4127

-20.0774

30.8516

-25.7866

26.3692

-15.1348

28.5813

-28.1630

17.8317

-27.5002

22.6711

-22.8906

22.0508

-31.7191

14.4763

-12.4282

30.2255

10.0514

22.2805

stats2 =

0.6529

5.0168

b3 =

1.0e+003 *

-5.1606

0.0019

0.0028

0.0049

bint3 =

1.0e+004 *

0.0303

137.4243

-1.1668

0.1347

-0.0003

0.0007

-0.0001

0.0006

-0.0001

0.0011

r3 =

16.9876

-13.4058

-17.2214

2.4798

-5.1463

-4.0615

18.2898

8.8141

-21.6383

-3.8351

13.4352

5.3018

rint3 =

-10.1025

44.0776

-43.2222

16.4106

-46.8617

12.4189

-32.7948

37.7543

-39.4735

29.1810

-35.9056

27.7827

-8.2274

44.8071

-22.8084

40.4366

-45.1638

1.8872

-21.6098

13.9396

-11.8217

38.6921

-27.0912

37.6949

stats3 =

0.6191

4.3345

b4 =

-611.2928

-1.2365

-1.4924

0.7731

bint4 =

1.0e+003 *

-6.2039

4.9814

-0.0059

0.0035

-0.0043

0.0013

-0.0045

0.0061

r4 =

14.3686

0.0431

239.7180

-22.1351

-14.8314

7.1045

1.2160

3.7356

1.7434

-3.3339

-8.5361

-17.9624

34.0498

4.5810

rint4 =

-19.4004

48.1376

-50.2215

5.9514

-53.7754

24.1126

-33.9306

48.1395

-35.5116

37.9435

-32.1544

39.6256

-35.2939

38.7807

-41.7155

35.0478

-43.3929

26.3208

-48.3442

12.4193

13.4904

54.6092

-17.9228

27.0848

stats4 =

0.6552

5.0676

b5 =

1.0e+003 *

-2.1121

0.0006

-0.0017

0.0022

bint5 =

1.0e+003 *

-9.1799

4.9556

-0.0055

0.0068

-0.0054

0.0020

-0.0045

0.0089

r5 =

12.5140

-7.9311

-25.4398

-12.2676

-3.8587

19.6319

4.7886

9.6521

-31.8747

0.0296

321.5417