Home » date » 2010 » Dec » 26 »

paper MR

*The author of this computation has been verified*
R Software Module: /rwasp_multipleregression.wasp (opens new window with default values)
Title produced by software: Multiple Regression
Date of computation: Sun, 26 Dec 2010 22:10:23 +0000
 
Cite this page as follows:
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL http://www.freestatistics.org/blog/date/2010/Dec/26/t129340208496v1818pq5tc140.htm/, Retrieved Sun, 26 Dec 2010 23:21:36 +0100
 
BibTeX entries for LaTeX users:
@Manual{KEY,
    author = {{YOUR NAME}},
    publisher = {Office for Research Development and Education},
    title = {Statistical Computations at FreeStatistics.org, URL http://www.freestatistics.org/blog/date/2010/Dec/26/t129340208496v1818pq5tc140.htm/},
    year = {2010},
}
@Manual{R,
    title = {R: A Language and Environment for Statistical Computing},
    author = {{R Development Core Team}},
    organization = {R Foundation for Statistical Computing},
    address = {Vienna, Austria},
    year = {2010},
    note = {{ISBN} 3-900051-07-0},
    url = {http://www.R-project.org},
}
 
Original text written by user:
 
IsPrivate?
No (this computation is public)
 
User-defined keywords:
 
Dataseries X:
» Textbox « » Textfile « » CSV «
5,2 7,9 8,7 8,9 15,3 15,4 18,1 19,7 13 12,6 6,2 3,5 3,4 0 9,5 8,9 10,4 13,2 18,9 19 16,3 10,6 5,8 3,6 2,6 5 7,3 9,2 15,7 16,8 18,4 18,1 14,6 7,8 7,6 3,8 5,6 2,2 6,8 11,8 14,9 16,7 16,7 15,9 13,6 9,2 2,8 2,5 4,8 2,8 7,8 9 12,9 16,4 21,8 17,8 13,5 10 10,4 5,5 4 6,8 5,7 9,1 13,6 15 20,9 20,4 14 13,7 7,1 0,8 2,1 1,3 3,9 10,7 11,1 16,4 17,1 17,3 12,9 10,9 5,3 0,7 -0,2 6,5 8,6 8,5 13,3 16,2 17,5 21,2 14,8 10,3 7,3 5,1 4,4 6,2 7,7 9,3 15,6 16,3 16,6 17,4 15,3 9,7 3,7 4,6 5,4 3,1 7,9 10,1 15 15,6 19,7 18,1 17,7 10,7 6,2 4,2 4 5,9 7,1 10,5 15,1 16,8 15,3 18,4 16,1 11,3 7,9 5,6 3,4 4,8 6,5 8,5 15,1 15,7 18,7 19,2 12,9 14,4 6,2 3,3 4,6 7,2 7,8 9,9 13,6 17,1 17,8 18,6 14,7 10,5 8,6 4,4 2,3 2,8 8,8 10,7 13,9 19,3 19,5 20,4 15,3 7,9 8,3 4,5 3,2 5 6,6 11,1 12,8 16,3 17,4 18,9 15,8 11,7 6,4 2,9 4,7 2,4 7,2 10,7 13,4 18,5 18,3 16,8 16,6 14, etc...
 
Output produced by software:

Enter (or paste) a matrix (table) containing all data (time) series. Every column represents a different variable and must be delimited by a space or Tab. Every row represents a period in time (or category) and must be delimited by hard returns. The easiest way to enter data is to copy and paste a block of spreadsheet cells. Please, do not use commas or spaces to seperate groups of digits!


Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time14 seconds
R Server'Gwilym Jenkins' @ 72.249.127.135


Multiple Linear Regression - Estimated Regression Equation
Temperatuur[t] = + 3.24210526315791 + 0.110847953216363M1[t] + 0.762134502923958M2[t] + 3.49342105263159M3[t] + 6.43470760233919M4[t] + 10.3809941520468M5[t] + 12.7672807017544M6[t] + 14.7885672514620M7[t] + 14.6898538011696M8[t] + 11.2711403508772M9[t] + 7.3924269005848M10[t] + 3.21371345029240M11[t] + 0.00371345029239764t + e[t]


Multiple Linear Regression - Ordinary Least Squares
VariableParameterS.D.T-STAT
H0: parameter = 0
2-tail p-value1-tail p-value
(Intercept)3.242105263157910.4158357.796600
M10.1108479532163630.5208180.21280.8316470.415824
M20.7621345029239580.520771.46350.1447210.07236
M33.493421052631590.5207276.708700
M46.434707602339190.52068912.358100
M510.38099415204680.52065519.938300
M612.76728070175440.52062624.52300
M714.78856725146200.52060128.406700
M814.68985380116960.5205828.218200
M911.27114035087720.52056421.651800
M107.39242690058480.52055314.201100
M113.213713450292400.5205466.173700
t0.003713450292397640.0015362.41830.0163830.008192


Multiple Linear Regression - Regression Statistics
Multiple R0.958336827443597
R-squared0.91840947483466
Adjusted R-squared0.914096319319311
F-TEST (value)212.932149459136
F-TEST (DF numerator)12
F-TEST (DF denominator)227
p-value0
Multiple Linear Regression - Residual Statistics
Residual Standard Deviation1.64610468633684
Sum Squared Residuals615.092964912282


Multiple Linear Regression - Actuals, Interpolation, and Residuals
Time or IndexActualsInterpolation
Forecast
Residuals
Prediction Error
15.23.356666666666621.84333333333338
27.94.011666666666613.88833333333339
38.76.746666666666581.95333333333342
48.99.69166666666671-0.79166666666671
515.313.64166666666671.65833333333329
615.416.0316666666666-0.631666666666627
718.118.05666666666680.0433333333332414
819.717.96166666666661.73833333333338
91314.5466666666667-1.54666666666666
1012.610.67166666666661.92833333333338
116.26.4966666666667-0.296666666666691
123.53.286666666666680.213333333333323
133.43.40122807017544-0.00122807017544409
1404.05622807017544-4.05622807017544
159.56.791228070175442.70877192982456
168.99.73622807017544-0.836228070175437
1710.413.6862280701754-3.28622807017544
1813.216.0762280701754-2.87622807017544
1918.918.10122807017540.798771929824564
201918.00622807017540.993771929824557
2116.314.59122807017541.70877192982456
2210.610.7162280701754-0.116228070175443
235.86.54122807017544-0.74122807017544
243.63.331228070175440.268771929824558
252.63.44578947368422-0.845789473684218
2654.100789473684210.899210526315787
277.36.835789473684220.464210526315784
289.29.7807894736842-0.58078947368421
2915.713.73078947368421.96921052631579
3016.816.12078947368420.679210526315786
3118.418.14578947368420.254210526315792
3218.118.05078947368420.049210526315787
3314.614.6357894736842-0.0357894736842123
347.810.7607894736842-2.96078947368421
357.66.585789473684211.01421052631579
363.83.375789473684210.424210526315786
375.63.490350877192992.10964912280701
382.24.14535087719298-1.94535087719298
396.86.88035087719299-0.0803508771929872
4011.89.825350877192981.97464912280702
4114.913.77535087719301.12464912280702
4216.716.1653508771930.534649122807013
4316.718.1903508771930-1.49035087719298
4415.918.095350877193-2.19535087719298
4513.614.680350877193-1.08035087719298
469.210.805350877193-1.60535087719299
472.86.63035087719298-3.83035087719298
482.53.42035087719298-0.920350877192985
494.83.534912280701761.26508771929824
502.84.18991228070176-1.38991228070176
517.86.924912280701760.87508771929824
5299.86991228070175-0.869912280701752
5312.913.8199122807018-0.919912280701753
5416.416.20991228070180.190087719298241
5521.818.23491228070173.56508771929825
5617.818.1399122807018-0.339912280701757
5713.514.7249122807018-1.22491228070176
581010.8499122807018-0.849912280701758
5910.46.674912280701753.72508771929825
605.53.464912280701762.03508771929824
6143.579473684210530.420526315789467
626.84.234473684210532.56552631578947
635.76.96947368421053-1.26947368421053
649.19.91447368421052-0.814473684210524
6513.613.8644736842105-0.264473684210526
661516.2544736842105-1.25447368421053
6720.918.27947368421052.62052631578948
6820.418.18447368421052.21552631578947
691414.7694736842105-0.769473684210527
7013.710.89447368421052.80552631578947
717.16.719473684210530.380526315789474
720.83.50947368421053-2.70947368421053
732.13.62403508771930-1.52403508771930
741.34.2790350877193-2.9790350877193
753.97.0140350877193-3.1140350877193
7610.79.95903508771930.740964912280703
7711.113.9090350877193-2.80903508771930
7816.416.29903508771930.100964912280697
7917.118.3240350877193-1.22403508771929
8017.318.2290350877193-0.9290350877193
8112.914.8140350877193-1.91403508771930
8210.910.9390350877193-0.0390350877193015
835.36.7640350877193-1.4640350877193
840.73.5540350877193-2.8540350877193
85-0.23.66859649122807-3.86859649122807
866.54.323596491228072.17640350877193
878.67.058596491228071.54140350877192
888.510.0035964912281-1.50359649122807
8913.313.9535964912281-0.653596491228068
9016.216.3435964912281-0.143596491228074
9117.518.3685964912281-0.868596491228065
9221.218.27359649122812.92640350877193
9314.814.8585964912281-0.0585964912280698
9410.310.9835964912281-0.683596491228073
957.36.808596491228070.49140350877193
965.13.598596491228071.50140350877193
974.43.713157894736850.686842105263153
986.24.368157894736841.83184210526316
997.77.103157894736850.596842105263154
1009.310.0481578947368-0.748157894736839
10115.613.99815789473681.60184210526316
10216.316.3881578947368-0.0881578947368438
10316.618.4131578947368-1.81315789473683
10417.418.3181578947368-0.918157894736846
10515.314.90315789473680.396842105263158
1069.711.0281578947368-1.32815789473685
1073.76.85315789473684-3.15315789473684
1084.63.643157894736840.956842105263156
1095.43.757719298245621.64228070175438
1103.14.41271929824562-1.31271929824562
1117.97.147719298245620.752280701754383
11210.110.09271929824560.00728070175438869
1131514.04271929824560.957280701754388
11415.616.4327192982456-0.832719298245617
11519.718.45771929824561.24228070175439
11618.118.3627192982456-0.262719298245615
11717.714.94771929824562.75228070175439
11810.711.0727192982456-0.372719298245618
1196.26.89771929824561-0.697719298245613
1204.23.687719298245620.512280701754384
12143.802280701754390.197719298245609
1225.94.457280701754391.44271929824561
1237.17.19228070175439-0.0922807017543897
12410.510.13728070175440.362719298245617
12515.114.08728070175441.01271929824562
12616.816.47728070175440.322719298245613
12715.318.5022807017544-3.20228070175438
12818.418.4072807017544-0.00728070175438963
12916.114.99228070175441.10771929824562
13011.311.11728070175440.182719298245612
1317.96.942280701754380.957719298245616
1325.63.732280701754391.86771929824561
1333.43.84684210526316-0.446842105263163
1344.84.501842105263160.298157894736841
1356.57.23684210526316-0.736842105263161
1368.510.1818421052632-1.68184210526315
13715.114.13184210526320.968157894736844
13815.716.5218421052632-0.821842105263161
13918.718.54684210526320.153157894736848
14019.218.45184210526320.74815789473684
14112.915.0368421052632-2.13684210526316
14214.411.16184210526323.23815789473684
1436.26.98684210526316-0.786842105263156
1443.33.77684210526316-0.476842105263162
1454.63.891403508771930.708596491228065
1467.24.546403508771932.65359649122807
1477.87.281403508771930.518596491228067
1489.910.2264035087719-0.326403508771926
14913.614.1764035087719-0.576403508771928
15017.116.56640350877190.53359649122807
15117.818.5914035087719-0.791403508771923
15218.618.49640350877190.103596491228070
15314.715.0814035087719-0.38140350877193
15410.511.2064035087719-0.706403508771932
1558.67.031403508771931.56859649122807
1564.43.821403508771930.57859649122807
1572.33.93596491228071-1.63596491228071
1582.84.5909649122807-1.79096491228070
1598.87.32596491228071.47403508771930
16010.710.27096491228070.429035087719301
16113.914.2209649122807-0.320964912280699
16219.316.61096491228072.68903508771930
16319.518.63596491228070.864035087719304
16420.418.54096491228071.85903508771930
16515.315.12596491228070.174035087719300
1667.911.2509649122807-3.3509649122807
1678.37.07596491228071.2240350877193
1684.53.86596491228070.634035087719298
1693.23.98052631578948-0.780526315789478
17054.635526315789470.364473684210526
1716.67.37052631578948-0.770526315789477
17211.110.31552631578950.78447368421053
17312.814.2655263157895-1.46552631578947
17416.316.6555263157895-0.355526315789475
17517.418.6805263157895-1.28052631578947
17618.918.58552631578950.314473684210523
17715.815.17052631578950.629473684210528
17811.711.29552631578950.404473684210523
1796.47.12052631578947-0.72052631578947
1802.93.91052631578948-1.01052631578948
1814.74.025087719298250.67491228070175
1822.44.68008771929824-2.28008771929824
1837.27.41508771929825-0.215087719298248
18410.710.36008771929820.339912280701758
18513.414.3100877192982-0.910087719298243
18618.516.70008771929821.79991228070175
18718.318.7250877192982-0.425087719298238
18816.818.6300877192982-1.83008771929825
18916.615.21508771929821.38491228070176
19014.111.34008771929822.75991228070175
1916.17.16508771929824-1.06508771929824
1923.53.95508771929825-0.455087719298246
1931.74.06964912280702-2.36964912280702
1942.34.72464912280702-2.42464912280702
1954.57.45964912280702-2.95964912280702
1969.310.4046491228070-1.10464912280701
19714.214.354649122807-0.154649122807015
19817.316.74464912280700.555350877192982
1992318.7696491228074.23035087719299
20016.318.6746491228070-2.37464912280702
20118.415.2596491228073.14035087719298
20214.211.38464912280702.81535087719298
2039.17.209649122807021.89035087719298
2045.93.999649122807021.90035087719298
2057.24.114210526315793.08578947368421
2066.84.769210526315792.03078947368421
20787.50421052631580.495789473684208
20814.310.44921052631583.85078947368421
20914.614.39921052631580.200789473684214
21017.516.78921052631580.71078947368421
21117.218.8142105263158-1.61421052631578
21217.218.7192105263158-1.51921052631579
21314.115.3042105263158-1.20421052631579
21410.511.4292105263158-0.92921052631579
2156.87.25421052631579-0.454210526315787
2164.14.044210526315790.0557894736842105
2176.54.158771929824562.34122807017544
2186.14.813771929824561.28622807017544
2196.37.54877192982456-1.24877192982456
2209.310.4937719298246-1.19377192982456
22116.414.44377192982461.95622807017544
22216.116.8337719298246-0.733771929824561
2231818.8587719298246-0.858771929824554
22417.618.7637719298246-1.16377192982456
2251415.3487719298246-1.34877192982456
22610.511.4737719298246-0.973771929824562
2276.97.29877192982456-0.398771929824558
2282.84.08877192982456-1.28877192982456
2290.74.20333333333334-3.50333333333334
2303.64.85833333333333-1.25833333333333
2316.77.59333333333334-0.893333333333335
23212.510.53833333333331.96166666666667
23314.414.4883333333333-0.0883333333333292
23416.516.8783333333333-0.378333333333334
23518.718.9033333333333-0.203333333333326
23619.418.80833333333330.591666666666665
23715.815.39333333333330.406666666666669
23811.311.5183333333333-0.218333333333334
2399.77.343333333333332.35666666666667
2402.94.13333333333333-1.23333333333333


Goldfeld-Quandt test for Heteroskedasticity
p-valuesAlternative Hypothesis
breakpoint indexgreater2-sidedless
160.9668014483741360.06639710325172830.0331985516258642
170.9583616625359310.08327667492813690.0416383374640685
180.927841214637240.1443175707255220.0721587853627609
190.9386352952263160.1227294095473680.061364704773684
200.9112898969036440.1774202061927120.0887101030963558
210.9699238660637540.06015226787249180.0300761339362459
220.951612487103250.09677502579350080.0483875128967504
230.9295229519903870.1409540960192260.0704770480096132
240.9060011525131940.1879976949736120.0939988474868062
250.8674048353483520.2651903293032970.132595164651648
260.8805114934843770.2389770130312460.119488506515623
270.8401954474114210.3196091051771590.159804552588579
280.8167244497105130.3665511005789740.183275550289487
290.892874772781870.214250454436260.10712522721813
300.9143148495527070.1713703008945860.085685150447293
310.8853955773749520.2292088452500970.114604422625048
320.85475093543720.2904981291255990.145249064562799
330.8154604034140150.3690791931719710.184539596585985
340.8650526457911220.2698947084177550.134947354208878
350.8629013709693650.274197258061270.137098629030635
360.8309602168476420.3380795663047160.169039783152358
370.8461576683058070.3076846633883850.153842331694193
380.836376832592270.3272463348154590.163623167407729
390.8085340288148530.3829319423702930.191465971185146
400.8600909543651360.2798180912697280.139909045634864
410.843480404072020.3130391918559610.156519595927981
420.8309819991459860.3380360017080280.169018000854014
430.817365154855110.3652696902897780.182634845144889
440.8421411527123660.3157176945752680.157858847287634
450.8125311687986670.3749376624026660.187468831201333
460.7846359689117010.4307280621765970.215364031088299
470.8530498914060010.2939002171879990.146950108593999
480.8245735003903230.3508529992193550.175426499609677
490.8100239275579260.3799521448841490.189976072442074
500.7799249832579320.4401500334841360.220075016742068
510.7471343031532650.505731393693470.252865696846735
520.7079657894139160.5840684211721690.292034210586084
530.6683265488412980.6633469023174040.331673451158702
540.6483101325298020.7033797349403970.351689867470198
550.8100061810201570.3799876379596860.189993818979843
560.7759964226555020.4480071546889960.224003577344498
570.7454774790895620.5090450418208750.254522520910438
580.7107664239049280.5784671521901450.289233576095072
590.8857147204682030.2285705590635940.114285279531797
600.896026253666710.2079474926665810.103973746333290
610.8756922559204170.2486154881591660.124307744079583
620.9127835170703780.1744329658592440.0872164829296221
630.914950677396310.1700986452073810.0850493226036907
640.8979170743762330.2041658512475350.102082925623767
650.8767886503382230.2464226993235550.123211349661777
660.8588010964615960.2823978070768080.141198903538404
670.8805592021657640.2388815956684730.119440797834236
680.8950475128925120.2099049742149760.104952487107488
690.875514018758310.2489719624833790.124485981241689
700.9203690130637580.1592619738724840.079630986936242
710.9045082540728540.1909834918542930.0954917459271463
720.9296062519776240.1407874960447510.0703937480223756
730.9307291422215790.1385417155568430.0692708577784214
740.9491516933500660.1016966132998690.0508483066499344
750.9695216349805650.06095673003886920.0304783650194346
760.9654012502750690.06919749944986240.0345987497249312
770.9732337282033880.05353254359322410.0267662717966121
780.9678451054915160.06430978901696720.0321548945084836
790.9648927231958840.07021455360823130.0351072768041157
800.9580052365584090.08398952688318270.0419947634415914
810.9560682099025750.087863580194850.043931790097425
820.9465331494057130.1069337011885740.0534668505942871
830.9406093825246320.1187812349507360.0593906174753678
840.95340162065760.09319675868480.0465983793424
850.9789439682806210.04211206343875730.0210560317193787
860.985668802412680.02866239517463930.0143311975873196
870.9860049984561570.02799000308768630.0139950015438432
880.9844283203173270.03114335936534670.0155716796826733
890.9809996268425860.0380007463148290.0190003731574145
900.977055975882980.04588804823404060.0229440241170203
910.972443622696280.05511275460744080.0275563773037204
920.9842072149571210.03158557008575770.0157927850428789
930.981174963310850.03765007337830080.0188250366891504
940.976925286757010.04614942648597890.0230747132429894
950.9724850179307960.05502996413840770.0275149820692039
960.9737577188686080.05248456226278330.0262422811313916
970.9692312331207370.06153753375852580.0307687668792629
980.9719237768049020.05615244639019680.0280762231950984
990.966250622508750.06749875498249830.0337493774912492
1000.9599617934877470.08007641302450520.0400382065122526
1010.9613086706932080.07738265861358320.0386913293067916
1020.9532315217212620.09353695655747560.0467684782787378
1030.9543184541024120.09136309179517630.0456815458975881
1040.9472711943729480.1054576112541040.052728805627052
1050.939278344042340.1214433119153200.0607216559576602
1060.934399083565420.1312018328691600.0656009164345799
1070.9598869998493820.08022600030123530.0401130001506176
1080.955053437063320.08989312587335940.0449465629366797
1090.955664237844330.08867152431133910.0443357621556696
1100.9517959507990680.09640809840186420.0482040492009321
1110.944031059091850.1119378818163000.0559689409081499
1120.9338113492930080.1323773014139840.0661886507069919
1130.9259768058267680.1480463883464650.0740231941732324
1140.9172245887656580.1655508224686850.0827754112343425
1150.9108837610253330.1782324779493350.0891162389746673
1160.8946569606444280.2106860787111430.105343039355572
1170.923685328569430.1526293428611400.0763146714305701
1180.9107884140326110.1784231719347780.0892115859673888
1190.8996573655020130.2006852689959730.100342634497987
1200.8835299306422390.2329401387155220.116470069357761
1210.8637579754402080.2724840491195840.136242024559792
1220.8585165671199190.2829668657601620.141483432880081
1230.8364773575939310.3270452848121370.163522642406069
1240.8139738569372680.3720522861254640.186026143062732
1250.7969767730909070.4060464538181850.203023226909093
1260.7715492483959770.4569015032080460.228450751604023
1270.843674432692670.3126511346146610.156325567307330
1280.8194113149576850.361177370084630.180588685042315
1290.804106676876260.3917866462474800.195893323123740
1300.7780380387847970.4439239224304070.221961961215203
1310.7582848803577660.4834302392844680.241715119642234
1320.7644889806249380.4710220387501240.235511019375062
1330.735505894696960.528988210606080.26449410530304
1340.7036659223182280.5926681553635430.296334077681772
1350.6753153590655480.6493692818689040.324684640934452
1360.6868520364686760.6262959270626490.313147963531324
1370.6629676556913610.6740646886172770.337032344308639
1380.6447633012288370.7104733975423270.355236698771164
1390.6076618397987670.7846763204024660.392338160201233
1400.5794703045007880.8410593909984230.420529695499212
1410.6144082404502710.7711835190994570.385591759549729
1420.7060956905611650.587808618877670.293904309438835
1430.6877940411752570.6244119176494850.312205958824743
1440.6553756060729580.6892487878540840.344624393927042
1450.6253548283714370.7492903432571260.374645171628563
1460.6909824440943020.6180351118113970.309017555905698
1470.6634687219735010.6730625560529980.336531278026499
1480.635463009711470.7290739805770610.364536990288530
1490.6015083125233260.7969833749533490.398491687476674
1500.5662337201485290.8675325597029420.433766279851471
1510.5374513775251360.9250972449497280.462548622474864
1520.5000425479890320.9999149040219350.499957452010967
1530.4671523627462670.9343047254925340.532847637253733
1540.4363386891767870.8726773783535740.563661310823213
1550.4219293470224670.8438586940449350.578070652977533
1560.3868099720462230.7736199440924450.613190027953777
1570.3840978852245570.7681957704491140.615902114775443
1580.383603998326480.767207996652960.61639600167352
1590.3912115913330880.7824231826661760.608788408666912
1600.3565846559020110.7131693118040220.643415344097989
1610.3201471552527160.6402943105054320.679852844747284
1620.3611467009988040.7222934019976080.638853299001196
1630.3310992897148580.6621985794297160.668900710285142
1640.3689085864407430.7378171728814850.631091413559257
1650.3305129435662810.6610258871325630.669487056433719
1660.4786054007002160.9572108014004320.521394599299784
1670.4511435457299810.9022870914599620.548856454270019
1680.4177348078386860.8354696156773730.582265192161314
1690.3853628391476130.7707256782952260.614637160852387
1700.3509877242233880.7019754484467770.649012275776612
1710.3169205323193340.6338410646386690.683079467680666
1720.282815265632540.565630531265080.71718473436746
1730.2764912556450740.5529825112901490.723508744354926
1740.2460066457986320.4920132915972630.753993354201368
1750.2366641097388080.4733282194776160.763335890261192
1760.2173746016612430.4347492033224850.782625398338757
1770.1875944786717010.3751889573434010.8124055213283
1780.1606911143227110.3213822286454210.839308885677289
1790.1451947633401470.2903895266802940.854805236659853
1800.1288034487943660.2576068975887320.871196551205634
1810.1083782755486600.2167565510973210.89162172445134
1820.1242701962801260.2485403925602520.875729803719874
1830.1040019913315210.2080039826630420.895998008668479
1840.08758635304638030.1751727060927610.91241364695362
1850.08160341573362730.1632068314672550.918396584266373
1860.07600185070970850.1520037014194170.923998149290291
1870.0645838267883770.1291676535767540.935416173211623
1880.05921649065442920.1184329813088580.94078350934557
1890.04949904075296910.09899808150593830.95050095924703
1900.05770814570849290.1154162914169860.942291854291507
1910.05906657198699780.1181331439739960.940933428013002
1920.04747926070689760.09495852141379520.952520739293102
1930.07311939808965070.1462387961793010.92688060191035
1940.1216073420905600.2432146841811190.87839265790944
1950.1770002821397160.3540005642794320.822999717860284
1960.2567743192005040.5135486384010080.743225680799496
1970.257027834215980.514055668431960.74297216578402
1980.2198494644349240.4396989288698470.780150535565076
1990.380735036193470.761470072386940.61926496380653
2000.4540672088993690.9081344177987380.545932791100631
2010.5194092250380960.9611815499238070.480590774961904
2020.5727806331031190.8544387337937630.427219366896881
2030.5259576963480470.9480846073039060.474042303651953
2040.5326986367015930.9346027265968140.467301363298407
2050.6463197914576840.7073604170846330.353680208542316
2060.6412280552106030.7175438895787940.358771944789397
2070.6111217280759720.7777565438480570.388878271924028
2080.7673013639006180.4653972721987630.232698636099382
2090.7116670643418940.5766658713162120.288332935658106
2100.687476331524970.625047336950060.31252366847503
2110.6326661348877540.7346677302244920.367333865112246
2120.5792590578943820.8414818842112360.420740942105618
2130.509520129259880.980959741480240.49047987074012
2140.4316263351547150.863252670309430.568373664845285
2150.3865273299232590.7730546598465180.613472670076741
2160.3343769759433540.6687539518867090.665623024056646
2170.847081896119610.3058362077607790.152918103880390
2180.9345127950117470.1309744099765050.0654872049882527
2190.89565360193990.2086927961201990.104346398060099
2200.917374116142910.1652517677141800.0826258838570902
2210.9854930816450550.02901383670989020.0145069183549451
2220.9724193989268280.05516120214634410.0275806010731720
2230.940781710140440.1184365797191210.0592182898595604
2240.8599514683348330.2800970633303330.140048531665167


Meta Analysis of Goldfeld-Quandt test for Heteroskedasticity
Description# significant tests% significant testsOK/NOK
1% type I error level00OK
5% type I error level100.0478468899521531OK
10% type I error level390.186602870813397NOK
 
Charts produced by software:
http://www.freestatistics.org/blog/date/2010/Dec/26/t129340208496v1818pq5tc140/109qxk1293401408.png (open in new window)
http://www.freestatistics.org/blog/date/2010/Dec/26/t129340208496v1818pq5tc140/109qxk1293401408.ps (open in new window)


http://www.freestatistics.org/blog/date/2010/Dec/26/t129340208496v1818pq5tc140/12pi81293401408.png (open in new window)
http://www.freestatistics.org/blog/date/2010/Dec/26/t129340208496v1818pq5tc140/12pi81293401408.ps (open in new window)


http://www.freestatistics.org/blog/date/2010/Dec/26/t129340208496v1818pq5tc140/22pi81293401408.png (open in new window)
http://www.freestatistics.org/blog/date/2010/Dec/26/t129340208496v1818pq5tc140/22pi81293401408.ps (open in new window)


http://www.freestatistics.org/blog/date/2010/Dec/26/t129340208496v1818pq5tc140/3dgzt1293401408.png (open in new window)
http://www.freestatistics.org/blog/date/2010/Dec/26/t129340208496v1818pq5tc140/3dgzt1293401408.ps (open in new window)


http://www.freestatistics.org/blog/date/2010/Dec/26/t129340208496v1818pq5tc140/4dgzt1293401408.png (open in new window)
http://www.freestatistics.org/blog/date/2010/Dec/26/t129340208496v1818pq5tc140/4dgzt1293401408.ps (open in new window)


http://www.freestatistics.org/blog/date/2010/Dec/26/t129340208496v1818pq5tc140/5dgzt1293401408.png (open in new window)
http://www.freestatistics.org/blog/date/2010/Dec/26/t129340208496v1818pq5tc140/5dgzt1293401408.ps (open in new window)


http://www.freestatistics.org/blog/date/2010/Dec/26/t129340208496v1818pq5tc140/66phw1293401408.png (open in new window)
http://www.freestatistics.org/blog/date/2010/Dec/26/t129340208496v1818pq5tc140/66phw1293401408.ps (open in new window)


http://www.freestatistics.org/blog/date/2010/Dec/26/t129340208496v1818pq5tc140/7ghyz1293401408.png (open in new window)
http://www.freestatistics.org/blog/date/2010/Dec/26/t129340208496v1818pq5tc140/7ghyz1293401408.ps (open in new window)


http://www.freestatistics.org/blog/date/2010/Dec/26/t129340208496v1818pq5tc140/8ghyz1293401408.png (open in new window)
http://www.freestatistics.org/blog/date/2010/Dec/26/t129340208496v1818pq5tc140/8ghyz1293401408.ps (open in new window)


http://www.freestatistics.org/blog/date/2010/Dec/26/t129340208496v1818pq5tc140/9ghyz1293401408.png (open in new window)
http://www.freestatistics.org/blog/date/2010/Dec/26/t129340208496v1818pq5tc140/9ghyz1293401408.ps (open in new window)


 
Parameters (Session):
par1 = 1 ; par2 = Include Monthly Dummies ; par3 = Linear Trend ;
 
Parameters (R input):
par1 = 1 ; par2 = Include Monthly Dummies ; par3 = Linear Trend ;
 
R code (references can be found in the software module):
library(lattice)
library(lmtest)
n25 <- 25 #minimum number of obs. for Goldfeld-Quandt test
par1 <- as.numeric(par1)
x <- t(y)
k <- length(x[1,])
n <- length(x[,1])
x1 <- cbind(x[,par1], x[,1:k!=par1])
mycolnames <- c(colnames(x)[par1], colnames(x)[1:k!=par1])
colnames(x1) <- mycolnames #colnames(x)[par1]
x <- x1
if (par3 == 'First Differences'){
x2 <- array(0, dim=c(n-1,k), dimnames=list(1:(n-1), paste('(1-B)',colnames(x),sep='')))
for (i in 1:n-1) {
for (j in 1:k) {
x2[i,j] <- x[i+1,j] - x[i,j]
}
}
x <- x2
}
if (par2 == 'Include Monthly Dummies'){
x2 <- array(0, dim=c(n,11), dimnames=list(1:n, paste('M', seq(1:11), sep ='')))
for (i in 1:11){
x2[seq(i,n,12),i] <- 1
}
x <- cbind(x, x2)
}
if (par2 == 'Include Quarterly Dummies'){
x2 <- array(0, dim=c(n,3), dimnames=list(1:n, paste('Q', seq(1:3), sep ='')))
for (i in 1:3){
x2[seq(i,n,4),i] <- 1
}
x <- cbind(x, x2)
}
k <- length(x[1,])
if (par3 == 'Linear Trend'){
x <- cbind(x, c(1:n))
colnames(x)[k+1] <- 't'
}
x
k <- length(x[1,])
df <- as.data.frame(x)
(mylm <- lm(df))
(mysum <- summary(mylm))
if (n > n25) {
kp3 <- k + 3
nmkm3 <- n - k - 3
gqarr <- array(NA, dim=c(nmkm3-kp3+1,3))
numgqtests <- 0
numsignificant1 <- 0
numsignificant5 <- 0
numsignificant10 <- 0
for (mypoint in kp3:nmkm3) {
j <- 0
numgqtests <- numgqtests + 1
for (myalt in c('greater', 'two.sided', 'less')) {
j <- j + 1
gqarr[mypoint-kp3+1,j] <- gqtest(mylm, point=mypoint, alternative=myalt)$p.value
}
if (gqarr[mypoint-kp3+1,2] < 0.01) numsignificant1 <- numsignificant1 + 1
if (gqarr[mypoint-kp3+1,2] < 0.05) numsignificant5 <- numsignificant5 + 1
if (gqarr[mypoint-kp3+1,2] < 0.10) numsignificant10 <- numsignificant10 + 1
}
gqarr
}
bitmap(file='test0.png')
plot(x[,1], type='l', main='Actuals and Interpolation', ylab='value of Actuals and Interpolation (dots)', xlab='time or index')
points(x[,1]-mysum$resid)
grid()
dev.off()
bitmap(file='test1.png')
plot(mysum$resid, type='b', pch=19, main='Residuals', ylab='value of Residuals', xlab='time or index')
grid()
dev.off()
bitmap(file='test2.png')
hist(mysum$resid, main='Residual Histogram', xlab='values of Residuals')
grid()
dev.off()
bitmap(file='test3.png')
densityplot(~mysum$resid,col='black',main='Residual Density Plot', xlab='values of Residuals')
dev.off()
bitmap(file='test4.png')
qqnorm(mysum$resid, main='Residual Normal Q-Q Plot')
qqline(mysum$resid)
grid()
dev.off()
(myerror <- as.ts(mysum$resid))
bitmap(file='test5.png')
dum <- cbind(lag(myerror,k=1),myerror)
dum
dum1 <- dum[2:length(myerror),]
dum1
z <- as.data.frame(dum1)
z
plot(z,main=paste('Residual Lag plot, lowess, and regression line'), ylab='values of Residuals', xlab='lagged values of Residuals')
lines(lowess(z))
abline(lm(z))
grid()
dev.off()
bitmap(file='test6.png')
acf(mysum$resid, lag.max=length(mysum$resid)/2, main='Residual Autocorrelation Function')
grid()
dev.off()
bitmap(file='test7.png')
pacf(mysum$resid, lag.max=length(mysum$resid)/2, main='Residual Partial Autocorrelation Function')
grid()
dev.off()
bitmap(file='test8.png')
opar <- par(mfrow = c(2,2), oma = c(0, 0, 1.1, 0))
plot(mylm, las = 1, sub='Residual Diagnostics')
par(opar)
dev.off()
if (n > n25) {
bitmap(file='test9.png')
plot(kp3:nmkm3,gqarr[,2], main='Goldfeld-Quandt test',ylab='2-sided p-value',xlab='breakpoint')
grid()
dev.off()
}
load(file='createtable')
a<-table.start()
a<-table.row.start(a)
a<-table.element(a, 'Multiple Linear Regression - Estimated Regression Equation', 1, TRUE)
a<-table.row.end(a)
myeq <- colnames(x)[1]
myeq <- paste(myeq, '[t] = ', sep='')
for (i in 1:k){
if (mysum$coefficients[i,1] > 0) myeq <- paste(myeq, '+', '')
myeq <- paste(myeq, mysum$coefficients[i,1], sep=' ')
if (rownames(mysum$coefficients)[i] != '(Intercept)') {
myeq <- paste(myeq, rownames(mysum$coefficients)[i], sep='')
if (rownames(mysum$coefficients)[i] != 't') myeq <- paste(myeq, '[t]', sep='')
}
}
myeq <- paste(myeq, ' + e[t]')
a<-table.row.start(a)
a<-table.element(a, myeq)
a<-table.row.end(a)
a<-table.end(a)
table.save(a,file='mytable1.tab')
a<-table.start()
a<-table.row.start(a)
a<-table.element(a,hyperlink('http://www.xycoon.com/ols1.htm','Multiple Linear Regression - Ordinary Least Squares',''), 6, TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Variable',header=TRUE)
a<-table.element(a,'Parameter',header=TRUE)
a<-table.element(a,'S.D.',header=TRUE)
a<-table.element(a,'T-STAT<br />H0: parameter = 0',header=TRUE)
a<-table.element(a,'2-tail p-value',header=TRUE)
a<-table.element(a,'1-tail p-value',header=TRUE)
a<-table.row.end(a)
for (i in 1:k){
a<-table.row.start(a)
a<-table.element(a,rownames(mysum$coefficients)[i],header=TRUE)
a<-table.element(a,mysum$coefficients[i,1])
a<-table.element(a, round(mysum$coefficients[i,2],6))
a<-table.element(a, round(mysum$coefficients[i,3],4))
a<-table.element(a, round(mysum$coefficients[i,4],6))
a<-table.element(a, round(mysum$coefficients[i,4]/2,6))
a<-table.row.end(a)
}
a<-table.end(a)
table.save(a,file='mytable2.tab')
a<-table.start()
a<-table.row.start(a)
a<-table.element(a, 'Multiple Linear Regression - Regression Statistics', 2, TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a, 'Multiple R',1,TRUE)
a<-table.element(a, sqrt(mysum$r.squared))
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a, 'R-squared',1,TRUE)
a<-table.element(a, mysum$r.squared)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a, 'Adjusted R-squared',1,TRUE)
a<-table.element(a, mysum$adj.r.squared)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a, 'F-TEST (value)',1,TRUE)
a<-table.element(a, mysum$fstatistic[1])
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a, 'F-TEST (DF numerator)',1,TRUE)
a<-table.element(a, mysum$fstatistic[2])
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a, 'F-TEST (DF denominator)',1,TRUE)
a<-table.element(a, mysum$fstatistic[3])
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a, 'p-value',1,TRUE)
a<-table.element(a, 1-pf(mysum$fstatistic[1],mysum$fstatistic[2],mysum$fstatistic[3]))
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a, 'Multiple Linear Regression - Residual Statistics', 2, TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a, 'Residual Standard Deviation',1,TRUE)
a<-table.element(a, mysum$sigma)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a, 'Sum Squared Residuals',1,TRUE)
a<-table.element(a, sum(myerror*myerror))
a<-table.row.end(a)
a<-table.end(a)
table.save(a,file='mytable3.tab')
a<-table.start()
a<-table.row.start(a)
a<-table.element(a, 'Multiple Linear Regression - Actuals, Interpolation, and Residuals', 4, TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a, 'Time or Index', 1, TRUE)
a<-table.element(a, 'Actuals', 1, TRUE)
a<-table.element(a, 'Interpolation<br />Forecast', 1, TRUE)
a<-table.element(a, 'Residuals<br />Prediction Error', 1, TRUE)
a<-table.row.end(a)
for (i in 1:n) {
a<-table.row.start(a)
a<-table.element(a,i, 1, TRUE)
a<-table.element(a,x[i])
a<-table.element(a,x[i]-mysum$resid[i])
a<-table.element(a,mysum$resid[i])
a<-table.row.end(a)
}
a<-table.end(a)
table.save(a,file='mytable4.tab')
if (n > n25) {
a<-table.start()
a<-table.row.start(a)
a<-table.element(a,'Goldfeld-Quandt test for Heteroskedasticity',4,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'p-values',header=TRUE)
a<-table.element(a,'Alternative Hypothesis',3,header=TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'breakpoint index',header=TRUE)
a<-table.element(a,'greater',header=TRUE)
a<-table.element(a,'2-sided',header=TRUE)
a<-table.element(a,'less',header=TRUE)
a<-table.row.end(a)
for (mypoint in kp3:nmkm3) {
a<-table.row.start(a)
a<-table.element(a,mypoint,header=TRUE)
a<-table.element(a,gqarr[mypoint-kp3+1,1])
a<-table.element(a,gqarr[mypoint-kp3+1,2])
a<-table.element(a,gqarr[mypoint-kp3+1,3])
a<-table.row.end(a)
}
a<-table.end(a)
table.save(a,file='mytable5.tab')
a<-table.start()
a<-table.row.start(a)
a<-table.element(a,'Meta Analysis of Goldfeld-Quandt test for Heteroskedasticity',4,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Description',header=TRUE)
a<-table.element(a,'# significant tests',header=TRUE)
a<-table.element(a,'% significant tests',header=TRUE)
a<-table.element(a,'OK/NOK',header=TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'1% type I error level',header=TRUE)
a<-table.element(a,numsignificant1)
a<-table.element(a,numsignificant1/numgqtests)
if (numsignificant1/numgqtests < 0.01) dum <- 'OK' else dum <- 'NOK'
a<-table.element(a,dum)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'5% type I error level',header=TRUE)
a<-table.element(a,numsignificant5)
a<-table.element(a,numsignificant5/numgqtests)
if (numsignificant5/numgqtests < 0.05) dum <- 'OK' else dum <- 'NOK'
a<-table.element(a,dum)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'10% type I error level',header=TRUE)
a<-table.element(a,numsignificant10)
a<-table.element(a,numsignificant10/numgqtests)
if (numsignificant10/numgqtests < 0.1) dum <- 'OK' else dum <- 'NOK'
a<-table.element(a,dum)
a<-table.row.end(a)
a<-table.end(a)
table.save(a,file='mytable6.tab')
}
 





Copyright

Creative Commons License

This work is licensed under a Creative Commons Attribution-Noncommercial-Share Alike 3.0 License.

Software written by Ed van Stee & Patrick Wessa


Disclaimer

Information provided on this web site is provided "AS IS" without warranty of any kind, either express or implied, including, without limitation, warranties of merchantability, fitness for a particular purpose, and noninfringement. We use reasonable efforts to include accurate and timely information and periodically update the information, and software without notice. However, we make no warranties or representations as to the accuracy or completeness of such information (or software), and we assume no liability or responsibility for errors or omissions in the content of this web site, or any software bugs in online applications. Your use of this web site is AT YOUR OWN RISK. Under no circumstances and under no legal theory shall we be liable to you or any other person for any direct, indirect, special, incidental, exemplary, or consequential damages arising from your access to, or use of, this web site.


Privacy Policy

We may request personal information to be submitted to our servers in order to be able to:

  • personalize online software applications according to your needs
  • enforce strict security rules with respect to the data that you upload (e.g. statistical data)
  • manage user sessions of online applications
  • alert you about important changes or upgrades in resources or applications

We NEVER allow other companies to directly offer registered users information about their products and services. Banner references and hyperlinks of third parties NEVER contain any personal data of the visitor.

We do NOT sell, nor transmit by any means, personal information, nor statistical data series uploaded by you to third parties.

We carefully protect your data from loss, misuse, alteration, and destruction. However, at any time, and under any circumstance you are solely responsible for managing your passwords, and keeping them secret.

We store a unique ANONYMOUS USER ID in the form of a small 'Cookie' on your computer. This allows us to track your progress when using this website which is necessary to create state-dependent features. The cookie is used for NO OTHER PURPOSE. At any time you may opt to disallow cookies from this website - this will not affect other features of this website.

We examine cookies that are used by third-parties (banner and online ads) very closely: abuse from third-parties automatically results in termination of the advertising contract without refund. We have very good reason to believe that the cookies that are produced by third parties (banner ads) do NOT cause any privacy or security risk.

FreeStatistics.org is safe. There is no need to download any software to use the applications and services contained in this website. Hence, your system's security is not compromised by their use, and your personal data - other than data you submit in the account application form, and the user-agent information that is transmitted by your browser - is never transmitted to our servers.

As a general rule, we do not log on-line behavior of individuals (other than normal logging of webserver 'hits'). However, in cases of abuse, hacking, unauthorized access, Denial of Service attacks, illegal copying, hotlinking, non-compliance with international webstandards (such as robots.txt), or any other harmful behavior, our system engineers are empowered to log, track, identify, publish, and ban misbehaving individuals - even if this leads to ban entire blocks of IP addresses, or disclosing user's identity.


FreeStatistics.org is powered by