Home » date » 2009 » Dec » 29 »

Multiple Regression 4

*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: Tue, 29 Dec 2009 04:08:12 -0700
 
Cite this page as follows:
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL http://www.freestatistics.org/blog/date/2009/Dec/29/t126208500330p1w5t28nm4rkh.htm/, Retrieved Tue, 29 Dec 2009 12:10:16 +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/2009/Dec/29/t126208500330p1w5t28nm4rkh.htm/},
    year = {2009},
}
@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 = {2009},
    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:
Met Seasonal Dummies Met Lineair trend Met 4lags
 
Dataseries X:
» Textbox « » Textfile « » CSV «
4.3 10.8 4.4 4.4 4.5 4.6 4.1 10.4 4.3 4.4 4.4 4.5 3.9 10.1 4.1 4.3 4.4 4.4 3.7 9.8 3.9 4.1 4.3 4.4 3.6 9.7 3.7 3.9 4.1 4.3 3.9 10.3 3.6 3.7 3.9 4.1 4.2 10.9 3.9 3.6 3.7 3.9 4.2 10.8 4.2 3.9 3.6 3.7 4.1 10.6 4.2 4.2 3.9 3.6 4.1 10.4 4.1 4.2 4.2 3.9 4.1 10.3 4.1 4.1 4.2 4.2 4.1 10.2 4.1 4.1 4.1 4.2 4.1 10 4.1 4.1 4.1 4.1 4 9.7 4.1 4.1 4.1 4.1 3.9 9.4 4 4.1 4.1 4.1 3.8 9.2 3.9 4 4.1 4.1 3.8 9.1 3.8 3.9 4 4.1 4 9.6 3.8 3.8 3.9 4 4.4 10.2 4 3.8 3.8 3.9 4.6 10.2 4.4 4 3.8 3.8 4.6 10 4.6 4.4 4 3.8 4.6 9.9 4.6 4.6 4.4 4 4.7 9.9 4.6 4.6 4.6 4.4 4.8 9.9 4.7 4.6 4.6 4.6 4.8 9.7 4.8 4.7 4.6 4.6 4.7 9.5 4.8 4.8 4.7 4.6 4.7 9.4 4.7 4.8 4.8 4.7 4.7 9.3 4.7 4.7 4.8 4.8 4.6 9.3 4.7 4.7 4.7 4.8 5 9.9 4.6 4.7 4.7 4.7 5.4 10.5 5 4.6 4.7 4.7 5.5 10.6 5.4 5 4.6 4.7 5.6 10.6 5.5 5.4 5 4.6 5.6 10.5 5.6 5.5 5.4 5 5.8 10.6 5.6 5.6 5.5 5.4 6 10.8 5.8 5.6 5.6 5.5 6.1 10.8 6 5.8 5.6 5.6 6.1 10.7 6.1 6 5.8 5.6 6 10.6 6.1 6.1 6 5.8 6 10.6 6 6.1 6.1 6 6.1 10.8 6 6 6.1 6.1 6.5 11.4 6.1 6 6 6.1 7.1 12 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 time9 seconds
R Server'Gwilym Jenkins' @ 72.249.127.135


Multiple Linear Regression - Estimated Regression Equation
Y[t] = + 0.0335157175898164 + 0.0317667969328308X[t] + 1.49494432782640Y1[t] -0.669453640778338Y2[t] -0.126862020467570Y3[t] + 0.246168085780859Y4[t] -0.169060236440471M1[t] -0.143500017446115M2[t] -0.101294922082687M3[t] -0.179767159035877M4[t] -0.185797235933362M5[t] + 0.26915662315277M6[t] -0.114909765968540M7[t] -0.171035776651497M8[t] -0.00890299086717574M9[t] -0.068787914595903M10[t] + 0.0463262564646009M11[t] + 0.000843338901157726t + e[t]


Multiple Linear Regression - Ordinary Least Squares
VariableParameterS.D.T-STAT
H0: parameter = 0
2-tail p-value1-tail p-value
(Intercept)0.03351571758981640.1323170.25330.8002760.400138
X0.03176679693283080.0191351.66010.0983290.049165
Y11.494944327826400.06694722.330300
Y2-0.6694536407783380.121043-5.530700
Y3-0.1268620204675700.121326-1.04560.2968890.148445
Y40.2461680857808590.066213.7180.0002550.000128
M1-0.1690602364404710.057317-2.94960.0035290.001765
M2-0.1435000174461150.058663-2.44620.015230.007615
M3-0.1012949220826870.058463-1.73260.0845760.042288
M4-0.1797671590358770.058673-3.06390.002460.00123
M5-0.1857972359333620.059449-3.12530.0020180.001009
M60.269156623152770.0588424.57428e-064e-06
M7-0.1149097659685400.060531-1.89840.0589680.029484
M8-0.1710357766514970.065159-2.62490.009280.00464
M9-0.008902990867175740.061737-0.14420.8854690.442734
M10-0.0687879145959030.058749-1.17090.2429240.121462
M110.04632625646460090.0579020.80010.4245330.212266
t0.0008433389011577260.0005191.62590.1054080.052704


Multiple Linear Regression - Regression Statistics
Multiple R0.989982053773046
R-squared0.980064466792699
Adjusted R-squared0.97850986099213
F-TEST (value)630.426354020904
F-TEST (DF numerator)17
F-TEST (DF denominator)218
p-value0
Multiple Linear Regression - Residual Statistics
Residual Standard Deviation0.175250034029212
Sum Squared Residuals6.69534122513832


Multiple Linear Regression - Actuals, Interpolation, and Residuals
Time or IndexActualsInterpolation
Forecast
Residuals
Prediction Error
14.34.40203335242442-0.102033352424419
24.14.25430515223282-0.154305152232818
33.94.03116323735205-0.131163237352046
43.73.7915923648573-0.0915923648573
53.63.61888640527351-0.0188864052735082
63.94.05427876373087-0.154278763730867
74.24.181683241033510.0183167589664943
84.24.33432368056343-0.134323680563429
94.14.22743493891048-0.127434938910483
104.14.048337381507690.0516626184923065
114.14.30191400158816-0.201914001588164
124.14.26594060637819-0.165940606378193
134.14.066753540874230.0332464591257717
1444.08362705968989-0.0836270596898938
153.93.96765102209199-0.0676510220919917
163.83.80111969594859-0.00111969594858663
173.83.722893411600930.0771065883990729
1844.24958876560114-0.249588765601138
194.44.172484052574630.227515947425368
204.64.556671575189640.0433284248103599
214.64.71912934564898-0.119129345648983
224.64.521509161941610.0784908380583934
234.74.7105615021221-0.0105615021220973
244.84.86380663449747-0.0638066344974667
254.84.771785446276390.0282145537236084
264.74.71220407866075-0.0122040786607490
274.74.614512006980740.0854879930192584
284.74.625268601891350.0747313981086547
294.64.63276806594178-0.0327680659417752
3054.933514100728040.0664858992719626
315.45.234274223875980.165725776124024
325.55.52505070865344-0.0250507086534416
335.65.494378193045110.105621806954890
345.65.562431423354380.0375685766456217
355.85.700401281197080.099598718802924
3665.972191195116810.0278088048831929
376.15.993689243565190.106310756434808
386.16.007147422300880.0928525776991196
3966.003935025857-0.00393502585700686
4065.813359110131750.186640889868250
416.15.906087904177910.193912095822092
426.56.54312581515429-0.0431258151542925
437.16.691731760955040.408268239044956
447.47.259301386897560.140698613102444
457.47.44296062585501-0.0429606258550102
467.57.202256291132460.297743708867543
477.67.580527158898290.0194728411017143
487.87.694620415467190.105379584532811
497.87.745760817368560.0542391826314357
507.77.647027573946460.0529724260535435
517.67.536649300219690.0633506997803081
527.57.422528270925740.0774717290742593
537.37.34430198657808-0.0443019865780823
547.67.572008475013010.0279915249869868
5587.778288922924820.221711077075184
5688.12725684494207-0.127256844942070
577.97.9288059306332-0.0288059306332038
587.87.733845491490380.0661545085096245
597.77.86254448736629-0.162544487366291
607.87.744022023451520.0559779765484836
617.77.77713763654806-0.0771376365480636
627.57.56881743166521-0.0688174316652107
637.37.33298931473766-0.0329893147376583
647.17.12438861020757-0.0243886102075737
6576.951682650623780.0483173493762211
667.37.38707500908114-0.0870750090811387
677.87.517656166077060.282343833922938
687.97.97881551025211-0.078815510252111
697.97.890707152919420.00929284708057888
707.87.765609580434640.0343904195653585
717.87.83929381876405-0.0392938187640516
727.97.888549753549810.0114502464501885
737.87.87933681114661-0.0793368111466127
747.67.65515372452372-0.0551537245237173
757.47.44711909586753-0.0471190958675341
767.27.23216503195088-0.0321650319508847
776.97.0562823926738-0.156282392673806
787.17.18633252617932-0.0863325261793207
797.57.301309978548270.198690021451729
807.67.69893929872546-0.0989392987254639
817.47.63169885128134-0.231698851281342
827.37.192505127006670.107494872993328
837.27.37546328491366-0.175463284913664
847.37.294243831623730.00575616837626721
857.27.29956595644891-0.0995659564489112
867.17.091245751566060.00875424843394182
8777.00808874711443-0.00808874711442774
886.96.878860431595870.0211395684041336
896.86.77919401836340.0208059816365949
907.27.159571619274260.0404283807257425
917.67.451577815584150.148422184415849
927.77.71773749178353-0.0177374917835284
937.67.68071161678863-0.0807116167886335
947.57.443422622146060.0565773778539448
957.57.55943541597521-0.0594354159752149
967.67.61820087311445-0.0182008731144489
977.67.584371122133160.0156288778668375
987.67.516035827679470.0839641723205256
997.57.53686802081745-0.0368680208174539
1007.37.3280081391743-0.0280081391743017
1017.27.087601219997240.112398780002755
1027.47.55318763417745-0.153187634177453
10387.555714487275540.444285512724460
1048.28.2269602689245-0.0269602689244969
10588.22138046357024-0.221380463570235
1067.77.686692291430940.0133077085690649
1077.77.604032319188780.0959676808112202
1087.87.83081483541524-0.0308148354152411
1097.87.797740679949380.00225932005061674
1107.77.673818408952960.0261815910470431
1117.57.5451561693083-0.0451561693082975
1127.37.253747218960340.0462527810396612
1137.17.092971865907870.0070281340921267
1147.17.39713324077411-0.297133240774112
1157.27.13664642442010.0633535755798995
1166.87.19746693327844-0.397466933278436
1176.66.6367563065195-0.0367563065195005
1186.46.5211143916181-0.121114391618097
1196.46.53780534185541-0.137805341855413
1206.56.53723492376239-0.0372349237623927
1216.36.47876784747665-0.178767847476646
1225.96.08047351949303-0.180473519493025
1235.55.65310210812253-0.15310210812253
1245.25.28891278853631-0.0889127885363065
1254.95.10135871984097-0.201358719840969
1265.45.284023043441510.115976956558492
1275.85.8045830196624-0.00458301966239773
1285.75.97040607964144-0.270406079641435
1295.65.575648199571610.0243518004283855
1305.55.493690061969510.00630993803049341
1315.45.62554522081233-0.225545220812331
1325.45.46652254985306-0.0665225498530555
1335.45.340613691087110.0593863089128854
1345.55.348733283064730.151266716935266
1355.85.523012700920440.276987299079560
1365.75.82374505744521-0.123745057445209
1375.45.44601155330613-0.0460115533061346
1385.65.5163590585410.0836409414589957
1395.85.729027632980490.0709723670195055
1406.25.864991614943630.335008385056375
1416.86.41506867062920.384931329370794
1426.77.0090734392488-0.309073439248796
1436.76.585059859703090.114940140296909
1446.46.63840236732913-0.238402367329132
1456.36.178912545263890.121087454736111
1466.36.219334235259050.0806657647409535
1476.46.345149881888760.0548501181112442
1486.36.3494978335453-0.0494978335452978
14966.10643116980369-0.106431169803692
1506.36.21247774718010.0875222528198968
1516.36.52223045755278-0.222230457552780
1526.66.27320013171310.326799868286901
1537.56.756867124405550.743132875594448
1547.87.90675972904265-0.106759729042649
1557.97.836987013898020.0630129861019814
1567.87.712543762970290.0874562370297083
1577.67.514556419326290.0854435806737116
1587.57.357373980648730.142626019351275
1597.67.406238322444770.193761677555232
1607.57.53627477768879-0.0362747776887867
1617.37.278100827722570.0218991722774283
1627.67.502671669916970.0973283300830303
1637.57.7486556959051-0.248655695905095
1647.67.343798094622580.256201905377416
1657.97.626391753792240.273608246207757
1667.98.0322480513225-0.132248051322492
1678.17.919596497505660.180403502494341
1688.28.1787207261050.0212792738949928
16988.09678127923364-0.0967812792336396
1707.57.71599480492611-0.215994804926113
1716.87.15024242160975-0.350242421609751
1726.56.401338488060160.09866151193984
1736.66.446366791804780.153633208195221
1747.67.278570800417460.321429199582542
17588.21132441789252-0.211324417892520
1768.17.991675849295360.108324150704639
1777.77.90552962132311-0.205529621323112
1787.57.351574781334720.148425218665283
1797.67.522105914308020.0778940856919793
1807.87.84489981352785-0.0448998135278474
1817.87.83880826695043-0.0388082669504336
1827.87.666224597794070.133775402205931
1837.57.6958107177701-0.195810717770094
1847.57.212578779139750.287421220860252
1857.17.4114048130702-0.311404813070205
1867.57.342226362693320.157773637306680
1877.57.75726543356737-0.257265433567370
1887.67.481769433967980.118230566032019
1897.77.625967870777030.0740321292229708
1907.77.73523587019348-0.0352358701934781
1917.97.774738493723830.125261506276169
1928.18.043351727950280.0566482720497196
1938.28.064849776398660.135150223601336
1948.28.075131275441070.12486872455893
1958.28.068742199303910.131257800696086
1967.98.0276607163613-0.127660716361295
1977.37.59860748859514-0.298607488595136
1986.97.38368761966635-0.483687619666351
1996.66.84857098830948-0.248570988309482
2006.76.608499901650820.0915000983491846
2016.97.00896710960453-0.108967109604534
20277.11203035901249-0.112030359012489
2037.17.1570549458201-0.0570549458201079
2047.27.196542181139320.00345781886067522
2057.17.14742176741423-0.0474217674142335
2066.96.9693161349806-0.0693161349806013
20776.795428353982350.204571646017645
2086.87.02578090872034-0.22578090872034
2096.46.63953213812991-0.239532138129911
2106.76.572499193632670.127500806367334
2116.66.94600007166352-0.346000071663518
2126.46.5291913471233-0.129191347123301
2136.36.32995148925529-0.0299514892552874
2146.26.34184282758176-0.141842827581759
2156.56.379183544047330.120816455952673
2166.86.80940519410694-0.00940519410693876
2176.86.86737485707086-0.0673748570708648
2186.46.6143835095481-0.214383509548105
2196.16.085715993196270.0142840068037293
2205.85.89170563976206-0.0917056397620612
2216.15.702323222611480.397676777388523
2227.26.78832801902080.411671980979206
2237.37.8287992160486-0.528799216048607
2246.97.06199622154561-0.161996221545613
2256.16.4816447354695-0.381644735469499
2265.85.73982111823120.0601788817687991
2276.26.027749898312580.172250101687423
2287.16.799986584641320.300013415358676
2297.77.55373894304330.146261056956702
2307.97.74365222762630.156347772373705
2317.77.652425360414270.047574639585727
2327.47.271467535097110.128532464902892
2337.57.077193353976810.422806646023186
23487.98334053577620.0166594642238047
2358.18.28217599314864-0.182175993148639
23687.951947626286010.0480523737139876


Goldfeld-Quandt test for Heteroskedasticity
p-valuesAlternative Hypothesis
breakpoint indexgreater2-sidedless
210.006442658301847510.01288531660369500.993557341698152
220.0009109792072158070.001821958414431610.999089020792784
230.0002090779504048750.0004181559008097490.999790922049595
240.0005761631396115660.001152326279223130.999423836860388
250.0006563734466978320.001312746893395660.999343626553302
260.0003485169475872570.0006970338951745140.999651483052413
278.96399266819688e-050.0001792798533639380.999910360073318
282.19963601556938e-054.39927203113876e-050.999978003639844
291.90620351699062e-053.81240703398124e-050.99998093796483
301.05721725709989e-052.11443451419979e-050.99998942782743
315.3936416795443e-061.07872833590886e-050.99999460635832
322.33416672105400e-064.66833344210801e-060.999997665833279
338.23652789176766e-071.64730557835353e-060.99999917634721
343.74829896684443e-077.49659793368886e-070.999999625170103
352.47940732551697e-074.95881465103394e-070.999999752059267
361.19818998548701e-072.39637997097402e-070.999999880181001
374.02971109759671e-088.05942219519341e-080.999999959702889
381.08805663802705e-082.17611327605411e-080.999999989119434
391.48298202169896e-082.96596404339791e-080.99999998517018
405.5706228643657e-091.11412457287314e-080.999999994429377
411.60029280362787e-093.20058560725574e-090.999999998399707
424.36967195772226e-108.73934391544452e-100.999999999563033
431.38418523647804e-092.76837047295608e-090.999999998615815
441.04401022361943e-092.08802044723886e-090.99999999895599
454.8726843748224e-109.7453687496448e-100.999999999512732
461.79902142067903e-103.59804284135806e-100.999999999820098
472.39484405261704e-104.78968810523407e-100.999999999760516
487.37724730767207e-111.47544946153441e-100.999999999926227
491.42318908806446e-102.84637817612891e-100.999999999857681
501.71891674028303e-103.43783348056606e-100.999999999828108
511.76019130221396e-103.52038260442793e-100.99999999982398
521.1643601186109e-102.3287202372218e-100.999999999883564
536.4358788732335e-101.2871757746467e-090.999999999356412
542.28364562113425e-104.56729124226851e-100.999999999771635
551.03645872998783e-102.07291745997565e-100.999999999896354
567.24603842663935e-111.44920768532787e-100.99999999992754
573.69552030193347e-117.39104060386694e-110.999999999963045
581.99240930861477e-113.98481861722953e-110.999999999980076
591.23897708372943e-112.47795416745887e-110.99999999998761
607.26203329402113e-121.45240665880423e-110.999999999992738
612.99773834583325e-125.9954766916665e-120.999999999997002
621.11504846877296e-122.23009693754592e-120.999999999998885
636.26302348018417e-131.25260469603683e-120.999999999999374
642.46068597871230e-134.92137195742461e-130.999999999999754
659.48347131057598e-141.89669426211520e-130.999999999999905
663.59455669605720e-147.18911339211439e-140.999999999999964
674.80927123469128e-149.61854246938256e-140.999999999999952
682.49645869872536e-144.99291739745072e-140.999999999999975
698.88982862405407e-151.77796572481081e-140.99999999999999
702.32135593830841e-144.64271187661682e-140.999999999999977
719.4347775885208e-151.88695551770416e-140.99999999999999
723.20870221717152e-156.41740443434305e-150.999999999999997
731.39625200389364e-152.79250400778729e-150.999999999999999
746.22368160450782e-161.24473632090156e-151
752.5135518799057e-165.0271037598114e-161
769.03702151149203e-171.80740430229841e-161
772.58167891738803e-165.16335783477606e-161
789.14945655799429e-171.82989131159886e-161
798.87482458878098e-171.77496491775620e-161
806.44346270996427e-171.28869254199285e-161
811.31439497495713e-162.62878994991426e-161
825.57884789996516e-171.11576957999303e-161
834.24751601058713e-178.49503202117425e-171
843.05687217111534e-176.11374434223067e-171
851.17320359285788e-172.34640718571577e-171
866.70664631433826e-181.34132926286765e-171
872.43235739040793e-184.86471478081586e-181
888.2604445982458e-191.65208891964916e-181
892.9258413488026e-195.8516826976052e-191
901.59541182734444e-193.19082365468888e-191
911.27524982189468e-192.55049964378935e-191
924.68490512581827e-209.36981025163655e-201
932.74270424833523e-205.48540849667047e-201
941.27432309322228e-202.54864618644456e-201
954.30334837639926e-218.60669675279852e-211
961.43506818057316e-212.87013636114632e-211
975.84711841273181e-221.16942368254636e-211
985.94798521485604e-221.18959704297121e-211
992.02163173018307e-224.04326346036615e-221
1003.23741406802124e-226.47482813604247e-221
1011.28303797341589e-222.56607594683179e-221
1021.31993860796538e-222.63987721593075e-221
1034.21214497247513e-208.42428994495026e-201
1041.77997077883085e-203.55994155766171e-201
1051.47436612273315e-192.94873224546630e-191
1066.94955036668124e-181.38991007333625e-171
1074.58921393179378e-189.17842786358756e-181
1082.22221057514245e-184.4444211502849e-181
1091.24342180935936e-182.48684361871872e-181
1104.82770576921595e-199.6554115384319e-191
1112.97608174649378e-195.95216349298756e-191
1121.55425948393775e-193.1085189678755e-191
1137.62313672321588e-201.52462734464318e-191
1141.15855883277740e-182.31711766555480e-181
1154.54925163677502e-189.09850327355004e-181
1161.75159115727003e-163.50318231454005e-161
1171.56318856505945e-153.12637713011891e-150.999999999999998
1188.91055516190888e-161.78211103238178e-151
1191.21267840915774e-152.42535681831548e-150.999999999999999
1206.07908193155314e-161.21581638631063e-151
1218.14771676626216e-161.62954335325243e-151
1221.59211672814437e-153.18423345628873e-150.999999999999998
1231.15378198966875e-152.3075639793375e-150.999999999999999
1245.81364074988943e-161.16272814997789e-151
1254.66453539165664e-169.32907078331327e-161
1262.52313652180084e-145.04627304360168e-140.999999999999975
1272.03202265689459e-144.06404531378919e-140.99999999999998
1281.34327339214077e-132.68654678428155e-130.999999999999866
1296.94813703636337e-141.38962740727267e-130.99999999999993
1303.4906241227469e-146.9812482454938e-140.999999999999965
1318.87321085081556e-141.77464217016311e-130.999999999999911
1325.82096617596449e-141.16419323519290e-130.999999999999942
1336.35837069352818e-141.27167413870564e-130.999999999999936
1342.88162774234658e-135.76325548469316e-130.999999999999712
1354.47545258928639e-128.95090517857277e-120.999999999995524
1361.06248457091986e-112.12496914183971e-110.999999999989375
1371.46200551952317e-112.92401103904635e-110.99999999998538
1388.06099276667398e-121.61219855333480e-110.999999999991939
1398.23986583425835e-121.64797316685167e-110.99999999999176
1401.02154600892063e-092.04309201784125e-090.999999998978454
1414.71048866497935e-089.4209773299587e-080.999999952895113
1423.25013103991964e-076.50026207983927e-070.999999674986896
1432.01571406589029e-074.03142813178057e-070.999999798428593
1441.83802903346237e-063.67605806692475e-060.999998161970967
1451.50435325145885e-063.00870650291769e-060.999998495646748
1461.18283031988839e-062.36566063977679e-060.99999881716968
1479.95875401689736e-071.99175080337947e-060.999999004124598
1487.80689499212362e-071.56137899842472e-060.9999992193105
1498.62449491141593e-071.72489898228319e-060.99999913755051
1506.44915039090882e-071.28983007818176e-060.999999355084961
1512.06534475054701e-064.13068950109403e-060.99999793465525
1525.99661998861596e-061.19932399772319e-050.999994003380011
1530.005646040407268160.01129208081453630.994353959592732
1540.004882507023718440.009765014047436870.995117492976282
1550.003644885585284060.007289771170568120.996355114414716
1560.002978101268243930.005956202536487850.997021898731756
1570.002338598303410310.004677196606820610.99766140169659
1580.001940303422291840.003880606844583670.998059696577708
1590.001923020123581750.00384604024716350.998076979876418
1600.001462796394553760.002925592789107510.998537203605446
1610.001027272334565250.002054544669130490.998972727665435
1620.0007434962731484050.001486992546296810.999256503726852
1630.001420930023363600.002841860046727190.998579069976636
1640.002037416016124660.004074832032249320.997962583983875
1650.004953020382558950.00990604076511790.99504697961744
1660.00380393400397160.00760786800794320.996196065996028
1670.004041413688822850.00808282737764570.995958586311177
1680.002990516352311460.005981032704622920.997009483647689
1690.002275652917716620.004551305835433240.997724347082283
1700.002697668345894270.005395336691788550.997302331654106
1710.007460479788899650.01492095957779930.9925395202111
1720.006032505920013230.01206501184002650.993967494079987
1730.005655990798601040.01131198159720210.994344009201399
1740.01037641564831080.02075283129662160.989623584351689
1750.01177661458445500.02355322916891010.988223385415545
1760.01223473708781940.02446947417563880.98776526291218
1770.01273691002996040.02547382005992080.98726308997004
1780.01813920733792520.03627841467585040.981860792662075
1790.01412328896132030.02824657792264060.98587671103868
1800.01109315769906580.02218631539813170.988906842300934
1810.008085706282988060.01617141256597610.991914293717012
1820.007272636238833420.01454527247766680.992727363761167
1830.008606942579959950.01721388515991990.99139305742004
1840.02437763605377700.04875527210755390.975622363946223
1850.04553788129769530.09107576259539060.954462118702305
1860.06926275763997160.1385255152799430.930737242360028
1870.06721886815312070.1344377363062410.93278113184688
1880.07647859600967860.1529571920193570.923521403990321
1890.0934944088920030.1869888177840060.906505591107997
1900.07968581273924460.1593716254784890.920314187260755
1910.08397389469566540.1679477893913310.916026105304335
1920.07053948504243620.1410789700848720.929460514957564
1930.08662665963719280.1732533192743860.913373340362807
1940.1131043772919890.2262087545839780.886895622708011
1950.1340993502590470.2681987005180940.865900649740953
1960.1592254709159560.3184509418319130.840774529084044
1970.1535103016323500.3070206032647010.84648969836765
1980.4168683417733260.8337366835466510.583131658226674
1990.3865783480561130.7731566961122270.613421651943887
2000.3428170744943500.6856341489886990.65718292550565
2010.2927654569949030.5855309139898070.707234543005097
2020.2587145429758380.5174290859516750.741285457024162
2030.2201723006801830.4403446013603650.779827699319817
2040.1908556847526550.3817113695053090.809144315247345
2050.1490049017413370.2980098034826740.850995098258663
2060.1078471629967170.2156943259934350.892152837003283
2070.1215065804792200.2430131609584410.87849341952078
2080.08637888445771690.1727577689154340.913621115542283
2090.5812556498516660.8374887002966670.418744350148334
2100.8779997806187010.2440004387625980.122000219381299
2110.9600598725048270.0798802549903450.0399401274951725
2120.936323688287010.1273526234259810.0636763117129903
2130.9261348762495390.1477302475009220.0738651237504612
2140.8598370347936860.2803259304126290.140162965206314
2150.7423960347179440.5152079305641120.257603965282056


Meta Analysis of Goldfeld-Quandt test for Heteroskedasticity
Description# significant tests% significant testsOK/NOK
1% type I error level1480.758974358974359NOK
5% type I error level1640.841025641025641NOK
10% type I error level1660.851282051282051NOK
 
Charts produced by software:
http://www.freestatistics.org/blog/date/2009/Dec/29/t126208500330p1w5t28nm4rkh/104ook1262084880.png (open in new window)
http://www.freestatistics.org/blog/date/2009/Dec/29/t126208500330p1w5t28nm4rkh/104ook1262084880.ps (open in new window)


http://www.freestatistics.org/blog/date/2009/Dec/29/t126208500330p1w5t28nm4rkh/1wl241262084880.png (open in new window)
http://www.freestatistics.org/blog/date/2009/Dec/29/t126208500330p1w5t28nm4rkh/1wl241262084880.ps (open in new window)


http://www.freestatistics.org/blog/date/2009/Dec/29/t126208500330p1w5t28nm4rkh/2mc5h1262084880.png (open in new window)
http://www.freestatistics.org/blog/date/2009/Dec/29/t126208500330p1w5t28nm4rkh/2mc5h1262084880.ps (open in new window)


http://www.freestatistics.org/blog/date/2009/Dec/29/t126208500330p1w5t28nm4rkh/3cert1262084880.png (open in new window)
http://www.freestatistics.org/blog/date/2009/Dec/29/t126208500330p1w5t28nm4rkh/3cert1262084880.ps (open in new window)


http://www.freestatistics.org/blog/date/2009/Dec/29/t126208500330p1w5t28nm4rkh/4fo7r1262084880.png (open in new window)
http://www.freestatistics.org/blog/date/2009/Dec/29/t126208500330p1w5t28nm4rkh/4fo7r1262084880.ps (open in new window)


http://www.freestatistics.org/blog/date/2009/Dec/29/t126208500330p1w5t28nm4rkh/5y4ef1262084880.png (open in new window)
http://www.freestatistics.org/blog/date/2009/Dec/29/t126208500330p1w5t28nm4rkh/5y4ef1262084880.ps (open in new window)


http://www.freestatistics.org/blog/date/2009/Dec/29/t126208500330p1w5t28nm4rkh/67ih71262084880.png (open in new window)
http://www.freestatistics.org/blog/date/2009/Dec/29/t126208500330p1w5t28nm4rkh/67ih71262084880.ps (open in new window)


http://www.freestatistics.org/blog/date/2009/Dec/29/t126208500330p1w5t28nm4rkh/70olg1262084880.png (open in new window)
http://www.freestatistics.org/blog/date/2009/Dec/29/t126208500330p1w5t28nm4rkh/70olg1262084880.ps (open in new window)


http://www.freestatistics.org/blog/date/2009/Dec/29/t126208500330p1w5t28nm4rkh/85fce1262084880.png (open in new window)
http://www.freestatistics.org/blog/date/2009/Dec/29/t126208500330p1w5t28nm4rkh/85fce1262084880.ps (open in new window)


http://www.freestatistics.org/blog/date/2009/Dec/29/t126208500330p1w5t28nm4rkh/92r4j1262084880.png (open in new window)
http://www.freestatistics.org/blog/date/2009/Dec/29/t126208500330p1w5t28nm4rkh/92r4j1262084880.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