Home » date » 2009 » Dec » 28 »

Multiple Regression 2

*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: Mon, 28 Dec 2009 10:27:53 -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/28/t1262021332ghmyfyued1fxe2e.htm/, Retrieved Mon, 28 Dec 2009 18:29:05 +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/28/t1262021332ghmyfyued1fxe2e.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 zonder lineaire trend
 
Dataseries X:
» Textbox « » Textfile « » CSV «
4.6 11.7 4.5 11.4 4.4 11.2 4.4 11.1 4.3 10.8 4.1 10.4 3.9 10.1 3.7 9.8 3.6 9.7 3.9 10.3 4.2 10.9 4.2 10.8 4.1 10.6 4.1 10.4 4.1 10.3 4.1 10.2 4.1 10 4 9.7 3.9 9.4 3.8 9.2 3.8 9.1 4 9.6 4.4 10.2 4.6 10.2 4.6 10 4.6 9.9 4.7 9.9 4.8 9.9 4.8 9.7 4.7 9.5 4.7 9.4 4.7 9.3 4.6 9.3 5 9.9 5.4 10.5 5.5 10.6 5.6 10.6 5.6 10.5 5.8 10.6 6 10.8 6.1 10.8 6.1 10.7 6 10.6 6 10.6 6.1 10.8 6.5 11.4 7.1 12.2 7.4 12.4 7.4 12.4 7.5 12.3 7.6 12.4 7.8 12.5 7.8 12.5 7.7 12.4 7.6 12.3 7.5 12.2 7.3 12.1 7.6 12.6 8 13.2 8 13.4 7.9 13.2 7.8 12.9 7.7 12.8 7.8 12.7 7.7 12.6 7.5 12.4 7.3 12.1 7.1 12 7 11.9 7.3 12.5 7.8 13.2 7.9 13.4 7.9 13.3 7.8 13 7.8 12.9 7.9 13 7.8 12.9 7.6 12.6 7.4 12.4 7.2 12.1 6.9 11.9 7.1 12.3 7.5 13 7.6 13 7.4 12.6 7.3 12.2 7.2 12.1 7.3 12 7.2 11.8 7.1 11.6 7 11.4 6.9 11.2 6.8 11.2 7.2 11.8 7.6 12.5 7.7 12.6 7.6 12.4 7.5 12.1 7.5 12 7.6 12 7.6 11.9 7.6 11.8 7.5 11.5 7.3 11.3 7.2 11.2 7.4 11.6 8 12.2 8.2 12.2 8 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] = + 4.99679000275911 + 0.180172709660390X[t] -0.164189637420380M1[t] -0.177840415167166M2[t] -0.107840415167166M3[t] + 0.00815094934981393M4[t] + 0.0097716745090601M5[t] -0.0651851480758428M6[t] -0.147439380015839M7[t] -0.274603110987156M8[t] -0.42690052034225M9[t] -0.215228325490940M10[t] -0.059008635483022M11[t] + e[t]


Multiple Linear Regression - Ordinary Least Squares
VariableParameterS.D.T-STAT
H0: parameter = 0
2-tail p-value1-tail p-value
(Intercept)4.996790002759110.5545469.010600
X0.1801727096603900.0463313.88890.0001326.6e-05
M1-0.1641896374203800.378872-0.43340.6651620.332581
M2-0.1778404151671660.379199-0.4690.6395270.319764
M3-0.1078404151671660.379199-0.28440.7763710.388185
M40.008150949349813930.3791080.02150.9828650.491433
M50.00977167450906010.379350.02580.9794720.489736
M6-0.06518514807584280.380113-0.17150.8639920.431996
M7-0.1474393800158390.381304-0.38670.6993620.349681
M8-0.2746031109871560.382526-0.71790.4735770.236788
M9-0.426900520342250.382622-1.11570.2657210.13286
M10-0.2152283254909400.37935-0.56740.5710290.285515
M11-0.0590086354830220.378869-0.15570.8763690.438184


Multiple Linear Regression - Regression Statistics
Multiple R0.288332841742365
R-squared0.0831358276272278
Adjusted R-squared0.0346672370172134
F-TEST (value)1.71525160069438
F-TEST (DF numerator)12
F-TEST (DF denominator)227
p-value0.0646419472974461
Multiple Linear Regression - Residual Statistics
Residual Standard Deviation1.19806691188221
Sum Squared Residuals325.827701853762


Multiple Linear Regression - Actuals, Interpolation, and Residuals
Time or IndexActualsInterpolation
Forecast
Residuals
Prediction Error
14.66.94062106836534-2.34062106836534
24.56.87291847772039-2.37291847772039
34.46.90688393578832-2.50688393578832
44.47.00485802933926-2.60485802933926
54.36.95242694160039-2.65242694160039
64.16.80540103515133-2.70540103515133
73.96.66909499031322-2.76909499031322
83.76.48787944644377-2.78787944644377
93.66.31756476612265-2.71756476612265
103.96.6373405867702-2.73734058677019
114.26.90166390257435-2.70166390257435
124.26.94265526709133-2.74265526709133
134.16.74243108773887-2.64243108773887
144.16.69274576806-2.59274576806000
154.16.74472849709396-2.64472849709397
164.16.8427025906449-2.74270259064491
174.16.80828877387207-2.70828877387207
1846.67928013838905-2.67928013838905
193.96.54297409355094-2.64297409355094
203.86.37977582064755-2.57977582064755
213.86.20946114032641-2.40946114032641
2246.51121969000792-2.51121969000792
234.46.77554300581207-2.37554300581207
244.66.83455164129509-2.23455164129509
254.66.63432746194263-2.03432746194263
264.66.60265941322981-2.00265941322981
274.76.67265941322981-1.97265941322981
284.86.78865077774679-1.98865077774679
294.86.75423696097396-1.95423696097396
304.76.64324559645698-1.94324559645698
314.76.54297409355094-1.84297409355094
324.76.39779309161358-1.69779309161358
334.66.24549568225849-1.64549568225849
3456.56527150290604-1.56527150290604
355.46.82959481871019-1.42959481871019
365.56.90662072515925-1.40662072515925
375.66.74243108773887-1.14243108773887
385.66.71076303902604-1.11076303902604
395.86.79878030999208-0.99878030999208
4066.95080621644114-0.95080621644114
416.16.95242694160039-0.852426941600385
426.16.85945284804944-0.759452848049446
4366.75918134514341-0.759181345143409
4466.63201761417209-0.632017614172091
456.16.51575474674908-0.415754746749075
466.56.83553056739662-0.33553056739662
477.17.13588842513285-0.0358884251328513
487.47.230931602547950.169068397452050
497.47.066741965127570.333258034872431
507.57.035073916414740.464926083585255
517.67.123091187380780.476908812619216
527.87.25709982286380.542900177136197
537.87.258720548023050.541279451976951
547.77.165746454472110.534253545527891
557.67.065474951566070.534525048433927
567.56.920293949628720.579706050371284
577.36.749979269307580.550020730692418
587.67.051737818989090.548262181010911
5987.316061134793240.683938865206758
6087.411104312208340.58889568779166
617.97.210880132855880.689119867144119
627.87.143177542210980.65682245778902
637.77.195160271244940.504839728755061
647.87.293134364795880.50686563520412
657.77.276737818989090.423262181010913
667.57.165746454472110.334253545527891
677.37.0294404096340.270559590366005
687.16.884259407696640.215740592303362
6976.71394472737550.286055272624496
707.37.033720548023050.26627945197695
717.87.316061134793240.483938865206758
727.97.411104312208340.48889568779166
737.97.228897403821920.67110259617808
747.87.161194813177020.638805186822982
757.87.213177542210980.586822457789022
767.97.3471861776940.552813822306002
777.87.33078963188720.469210368112795
787.67.201780996404190.398219003595813
797.47.083492222532110.316507777467889
807.26.902276678662680.297723321337323
816.96.71394472737550.186055272624496
827.16.997686006090970.102313993909028
837.57.280026592861160.219973407138837
847.67.339035228344180.260964771655815
857.47.102776507059650.297223492940353
867.37.01705664544870.282943354551294
877.27.069039374482670.130960625517334
887.37.167013468033610.132986531966392
897.27.132599651260770.0674003487392249
907.17.02160828674380.078391713256203
9176.903319512871720.0966804871282786
926.96.740121239968330.159878760031675
936.86.587823830613230.212176169386769
947.26.907599651260780.292400348739224
957.67.189940238030970.410059761969031
967.77.266966144480030.433033855519972
977.67.066741965127570.53325803487243
987.56.999039374482670.500960625517333
997.57.051022103516630.448977896483373
1007.67.167013468033610.432986531966392
1017.67.150616922226810.449383077773185
1027.67.057642828675880.542357171324125
1037.56.921336783837760.57866321616224
1047.36.758138510934360.541861489065635
1057.26.587823830613230.612176169386769
1067.46.87156510932870.528434890671302
10787.135888425132850.864111574867149
1088.27.194897060615871.00510293938413
10986.94062106836531.05937893163470
1107.76.836883935788320.863116064211685
1117.76.870849393856240.829150606143763
1127.86.968823487407180.831176512592821
1137.86.952426941600390.847573058399615
1147.76.823418306117370.876581693882633
1157.56.687112261279250.812887738720747
1167.36.523913988375860.776086011624143
1177.16.353599308054720.746400691945276
1187.16.637340586770190.462659413229809
1197.26.865629360642270.334370639357735
1206.86.87058618322717-0.0705861832271697
1216.66.65234473290867-0.0523447329086721
1226.46.56662487129773-0.16662487129773
1236.46.58257305839961-0.182573058399613
1246.56.6084780680864-0.108478068086398
1256.36.52001243841545-0.220012438415449
1265.96.39100380293243-0.49100380293243
1275.56.34478411292451-0.844784112924512
1285.26.18158584002112-0.981585840021116
1294.96.01127115969998-1.11127115969998
1305.46.34906425131357-0.949064251313566
1315.86.59537029615168-0.79537029615168
1325.76.61834438970262-0.918344389702623
1335.66.4361374813162-0.836137481316204
1345.56.35041761970526-0.850417619705262
1355.46.3483485358411-0.948348535841105
1365.46.35623627456185-0.956236274561852
1375.46.28578791585694-0.885787915856941
1385.56.17479655133996-0.674796551339963
1395.86.12857686133204-0.328576861332044
1405.75.98339585939469-0.283395859394687
1415.45.77704663714148-0.377046637141475
1425.66.0427706448909-0.442770644890904
1435.86.25304214779694-0.453042147796939
1446.26.38411986714412-0.184119867144116
1456.86.346051126486010.453948873513991
1466.76.332400348739220.367599651260777
1476.76.474469432603380.225530567396621
1486.46.64451261001848-0.244512610018476
1496.36.62811606421168-0.328116064211683
1506.36.48109015776263-0.181090157762626
1516.46.272715029060360.127284970939644
1526.36.109516756156960.190483243843040
15365.97523661776790.0247633822320956
1546.36.43915060614376-0.139150606143762
1556.36.63140483808376-0.331404838083759
1566.66.6543789316347-0.0543789316347017
1577.56.400102939384131.09989706061587
1587.86.332400348739221.46759965126078
1597.96.43843489067131.4615651093287
1607.86.626495339052441.17350466094756
1617.66.646133335177720.953866664822278
1627.56.499107428728661.00089257127134
1637.66.326766841958471.27323315804153
1647.56.145551298089041.35444870191096
1657.35.993253888733941.30674611126606
1667.66.421133335177721.17886666482228
1677.56.631404838083760.868595161916241
1687.66.690413473566780.90958652643322
1697.96.472172023248281.42782797675172
1707.96.440503974535461.45949602546454
1718.16.564555787433571.53544421256643
1728.26.788650777746791.41134922225321
17386.772254231941.22774576806001
1747.56.60721105452490.892788945475102
1756.86.344784112924510.455215887075488
1766.56.163568569055080.336431430944923
1776.66.101357514530180.498642485469822
1787.66.655357857736230.94464214226377
17986.937698444506421.06230155549358
1808.16.960672538057361.13932746194264
1817.76.634327461942631.06567253805737
1827.56.476538516467541.02346148353247
1837.66.546538516467531.05346148353247
1847.86.716581693882631.08341830611737
1857.86.736219690007921.06378030999208
1867.86.643245596456981.15675440354302
1877.56.488922280652821.01107771934718
1887.56.325724007749431.17427599225057
1897.16.191443869360370.908556130639627
1907.56.601306044838110.898693955161887
1917.56.793560276778110.70643972322189
1927.66.834551641295090.765448358704909
1937.76.562258378078481.13774162192152
1947.76.476538516467541.22346148353247
1957.96.564555787433571.33544421256643
1968.16.698564422916591.40143557708341
1978.26.700185148075841.49981485192416
1988.26.589193783558861.61080621644114
1998.26.470905009686781.72909499031321
2007.96.343741278715471.55625872128453
2017.36.191443869360371.10855613063963
2026.96.547254231940.352745768060004
2036.66.73950846387999-0.139508463879993
2046.76.76248255743094-0.0624825574309356
2056.96.508206565180360.391793434819640
20676.440503974535460.559496025464543
2077.16.510503974535460.589496025464543
2087.26.644512610018480.555487389981524
2097.16.646133335177720.453866664822278
2106.96.571176512592820.328823487407179
21176.506939551618860.493060448381137
2126.86.307706736783390.49229326321661
2136.46.06532297259810.334677027401901
2146.76.295012438415450.404987561584551
2156.66.397180315525250.202819684474748
2166.46.384119867144120.0158801328558843
2176.36.255964771655810.0440352283441861
2186.26.24231399390903-0.0423139939090281
2196.56.330331264875070.169668735124933
2206.86.4283053584260.371694641573992
2216.86.375874270687140.424125729312863
2226.46.210831093272040.189168906727960
2236.16.074525048433930.025474951566073
2245.85.89330950456449-0.093309504564492
2256.15.813081179073550.286918820926446
2267.26.258977896483370.941022103516629
2277.36.505283941321490.794716058678514
2286.96.492223492940350.40777650705965
2296.16.25596477165581-0.155964771655814
2305.86.17024491004487-0.370244910044872
2316.26.29429672294299-0.0942967229429881
2327.16.500374442290160.599625557709835
2337.76.520012438415451.17998756158455
2347.96.409021073898471.49097892610153
2357.76.236680487128281.46331951287172
2367.46.019430401326771.38056959867324
2377.55.885150262937711.61484973706229
23886.240960625517331.75903937448267
2398.16.469249399389411.63075060061059
24086.510240763906391.48975923609361


Goldfeld-Quandt test for Heteroskedasticity
p-valuesAlternative Hypothesis
breakpoint indexgreater2-sidedless
160.0001041311632676820.0002082623265353640.999895868836732
171.58960369558944e-053.17920739117888e-050.999984103963044
184.9494050823811e-069.8988101647622e-060.999995050594918
193.77801047775776e-067.55602095551551e-060.999996221989522
202.97076593538797e-065.94153187077594e-060.999997029234065
213.43403494175652e-066.86806988351304e-060.999996565965058
221.47070546263599e-062.94141092527197e-060.999998529294537
231.03530260862348e-062.07060521724695e-060.999998964697391
241.91286255831437e-063.82572511662874e-060.999998087137442
251.77367472289486e-063.54734944578971e-060.999998226325277
261.38134666218810e-062.76269332437621e-060.999998618653338
271.86922526969480e-063.73845053938961e-060.99999813077473
283.33602206710436e-066.67204413420872e-060.999996663977933
295.87754872802942e-061.17550974560588e-050.999994122451272
301.14701683381148e-052.29403366762295e-050.999988529831662
314.01184540750343e-058.02369081500687e-050.999959881545925
320.0002042089690067120.0004084179380134240.999795791030993
330.0006215944759590280.001243188951918060.99937840552404
340.002373027535442480.004746055070884960.997626972464557
350.00734685303501930.01469370607003860.99265314696498
360.01890706663533760.03781413327067520.981092933364662
370.04477050410461410.08954100820922820.955229495895386
380.09390658134359610.1878131626871920.906093418656404
390.2113448881606740.4226897763213480.788655111839326
400.4227285380525920.8454570761051840.577271461947408
410.6585799514493520.6828400971012960.341420048550648
420.8308173594997730.3383652810004540.169182640500227
430.9140480354031750.1719039291936500.0859519645968248
440.9528757575848230.09424848483035330.0471242424151767
450.9697657880828470.06046842383430510.0302342119171526
460.978406891944080.04318621611184170.0215931080559208
470.9824932106763320.0350135786473370.0175067893236685
480.9837817037122590.03243659257548280.0162182962877414
490.9859196390369930.02816072192601420.0140803609630071
500.986852482856680.02629503428663870.0131475171433194
510.9858157491682630.02836850166347330.0141842508317366
520.9840278956692960.03194420866140790.0159721043307040
530.980047380572130.03990523885574120.0199526194278706
540.9740822078978780.0518355842042430.0259177921021215
550.966266180689030.06746763862193960.0337338193109698
560.9564894294592650.08702114108147020.0435105705407351
570.9446371117752230.1107257764495550.0553628882247775
580.9304251703444450.1391496593111100.0695748296555548
590.91392752713740.1721449457251990.0860724728625993
600.8962119845383040.2075760309233920.103788015461696
610.8744789777607850.251042044478430.125521022239215
620.850412900583940.2991741988321190.149587099416059
630.8232114116831020.3535771766337950.176788588316898
640.7950646221783250.409870755643350.204935377821675
650.7630276720531380.4739446558937240.236972327946862
660.7294292341366270.5411415317267460.270570765863373
670.6947423898013190.6105152203973630.305257610198681
680.6606920250845290.6786159498309420.339307974915471
690.6226889698568960.7546220602862080.377311030143104
700.5848137028338120.8303725943323750.415186297166188
710.5432895484528870.9134209030942250.456710451547112
720.505563650747810.988872698504380.49443634925219
730.4623703510842850.924740702168570.537629648915715
740.4198600721184490.8397201442368980.580139927881551
750.3789693616058570.7579387232117140.621030638394143
760.3396431850124260.6792863700248510.660356814987574
770.3042988599754190.6085977199508370.695701140024581
780.2716229118470150.5432458236940290.728377088152985
790.2439520994014670.4879041988029340.756047900598533
800.216624440866280.433248881732560.78337555913372
810.1930289169276720.3860578338553440.806971083072328
820.1717227569862070.3434455139724140.828277243013793
830.1503673166387730.3007346332775470.849632683361227
840.1297245723010140.2594491446020270.870275427698986
850.1151042872600220.2302085745200450.884895712739978
860.1053994676845520.2107989353691050.894600532315448
870.09744906137725130.1948981227545030.902550938622749
880.09332298481369120.1866459696273820.906677015186309
890.09104483326266980.1820896665253400.90895516673733
900.08913586436616790.1782717287323360.910864135633832
910.08849149824444560.1769829964888910.911508501755554
920.08864862376661540.1772972475332310.911351376233385
930.08571609549678130.1714321909935630.914283904503219
940.08224623896541540.1644924779308310.917753761034585
950.07542966282275110.1508593256455020.92457033717725
960.0685090508591790.1370181017183580.93149094914082
970.06708290777463830.1341658155492770.932917092225362
980.0668186359499210.1336372718998420.933181364050079
990.06883893219550640.1376778643910130.931161067804494
1000.07134635189313960.1426927037862790.92865364810686
1010.07397509312856330.1479501862571270.926024906871437
1020.07551605390076830.1510321078015370.924483946099232
1030.08165951674242560.1633190334848510.918340483257574
1040.08757191415917970.1751438283183590.91242808584082
1050.09385103888199840.1877020777639970.906148961118002
1060.1028704251381010.2057408502762020.8971295748619
1070.1197358639173100.2394717278346190.88026413608269
1080.1496850304590810.2993700609181610.85031496954092
1090.2295814308198420.4591628616396840.770418569180158
1100.3314983140597730.6629966281195460.668501685940227
1110.4618538693660830.9237077387321650.538146130633917
1120.6000642095865410.7998715808269170.399935790413459
1130.7139491130892390.5721017738215220.286050886910761
1140.8114084127115570.3771831745768850.188591587288443
1150.881322073855830.2373558522883410.118677926144170
1160.9252541706343630.1494916587312730.0747458293656366
1170.9515032801769550.09699343964608940.0484967198230447
1180.96644079677060.0671184064588020.033559203229401
1190.9759586869839820.04808262603203610.0240413130160181
1200.9833739449807710.03325211003845710.0166260550192286
1210.9898819288597870.02023614228042620.0101180711402131
1220.9941990139172490.01160197216550210.00580098608275104
1230.9969139923466230.00617201530675490.00308600765337745
1240.9985483191164040.002903361767192580.00145168088359629
1250.9993343537745820.001331292450836050.000665646225418024
1260.9996608920300460.000678215939907320.00033910796995366
1270.9998497273506840.0003005452986325960.000150272649316298
1280.9999435644596870.0001128710806264635.64355403132313e-05
1290.9999802186103023.95627793954677e-051.97813896977339e-05
1300.9999910631389041.78737221929681e-058.93686109648403e-06
1310.9999949233876651.01532246700737e-055.07661233503685e-06
1320.9999975614605494.87707890207662e-062.43853945103831e-06
1330.9999991356602271.72867954623256e-068.6433977311628e-07
1340.9999997246027445.50794512774769e-072.75397256387385e-07
1350.9999999144194481.71161102965506e-078.55805514827529e-08
1360.9999999620335867.59328284885785e-083.79664142442892e-08
1370.9999999804599723.9080056646137e-081.95400283230685e-08
1380.9999999887496672.25006664516526e-081.12503332258263e-08
1390.999999993141611.37167790695222e-086.8583895347611e-09
1400.9999999955414758.91705030010686e-094.45852515005343e-09
1410.9999999966938196.6123623684032e-093.3061811842016e-09
1420.9999999971857325.62853562950225e-092.81426781475112e-09
1430.999999997102395.79521788105109e-092.89760894052554e-09
1440.999999996839746.32051950177211e-093.16025975088606e-09
1450.9999999967877136.42457351750567e-093.21228675875283e-09
1460.9999999965428536.91429436415779e-093.45714718207889e-09
1470.999999996600116.79977832742143e-093.39988916371072e-09
1480.9999999982119513.57609736929854e-091.78804868464927e-09
1490.9999999992499261.50014890342658e-097.5007445171329e-10
1500.9999999995682648.63470972962869e-104.31735486481434e-10
1510.999999999587418.25180711605138e-104.12590355802569e-10
1520.9999999995712978.57405379651369e-104.28702689825684e-10
1530.9999999995934868.1302802778123e-104.06514013890615e-10
1540.9999999997563124.87375919085769e-102.43687959542885e-10
1550.999999999846533.06940197896375e-101.53470098948187e-10
1560.999999999832913.3417948874349e-101.67089744371745e-10
1570.9999999998579472.8410529696251e-101.42052648481255e-10
1580.9999999999427051.14589793944765e-105.72948969723823e-11
1590.9999999999684576.30867613941891e-113.15433806970946e-11
1600.999999999962737.45405950520101e-113.72702975260051e-11
1610.9999999999447811.10437339687670e-105.52186698438352e-11
1620.9999999999197221.60555350867906e-108.02776754339528e-11
1630.9999999999165161.6696870769279e-108.3484353846395e-11
1640.9999999999265761.46847962961838e-107.3423981480919e-11
1650.9999999999305721.38856908173012e-106.94284540865061e-11
1660.9999999999044321.91135284405105e-109.55676422025526e-11
1670.9999999998379333.24133285554516e-101.62066642777258e-10
1680.9999999997322485.35504869899392e-102.67752434949696e-10
1690.9999999998204193.59161509637496e-101.79580754818748e-10
1700.9999999998745822.50835984507921e-101.25417992253960e-10
1710.9999999998969292.06142708782503e-101.03071354391252e-10
1720.9999999998487683.02463174131204e-101.51231587065602e-10
1730.9999999997314815.37037296886987e-102.68518648443493e-10
1740.9999999995065679.86865015645416e-104.93432507822708e-10
1750.99999999919271.61460217310105e-098.07301086550524e-10
1760.9999999987513452.49730889198201e-091.24865444599101e-09
1770.9999999978790974.24180626966584e-092.12090313483292e-09
1780.9999999961349457.73010916406653e-093.86505458203327e-09
1790.999999992457341.50853212165345e-087.54266060826724e-09
1800.9999999857270722.85458557249366e-081.42729278624683e-08
1810.9999999761166474.77667054598958e-082.38833527299479e-08
1820.9999999632948477.34103053968371e-083.67051526984186e-08
1830.9999999396425141.20714972820842e-076.03574864104208e-08
1840.9999998893608872.21278226832623e-071.10639113416311e-07
1850.9999997924576234.15084753354108e-072.07542376677054e-07
1860.9999996176726457.64654710343927e-073.82327355171963e-07
1870.9999992918283221.41634335590861e-067.08171677954304e-07
1880.9999987263957632.54720847308598e-061.27360423654299e-06
1890.9999976508863724.69822725580204e-062.34911362790102e-06
1900.9999958534175068.29316498842452e-064.14658249421226e-06
1910.9999924535421171.50929157665068e-057.54645788325342e-06
1920.999985877348552.82453028984513e-051.41226514492257e-05
1930.99998232328073.53534386002181e-051.76767193001090e-05
1940.9999823549066333.52901867346475e-051.76450933673237e-05
1950.9999821612263543.56775472920966e-051.78387736460483e-05
1960.9999779942958694.40114082625543e-052.20057041312772e-05
1970.9999712952382785.74095234444861e-052.87047617222430e-05
1980.9999680493407126.39013185753663e-053.19506592876831e-05
1990.9999711954064615.76091870773296e-052.88045935386648e-05
2000.9999703643304025.92713391961715e-052.96356695980858e-05
2010.9999467498906530.0001065002186946175.32501093473083e-05
2020.9999244980634670.0001510038730665177.55019365332585e-05
2030.9999403115414030.0001193769171931475.96884585965734e-05
2040.9999389762081730.0001220475836533456.10237918266727e-05
2050.9998728341237730.0002543317524529920.000127165876226496
2060.9997839262305590.0004321475388830360.000216073769441518
2070.9996008045276640.0007983909446716530.000399195472335826
2080.9991966356033060.001606728793388280.000803364396694142
2090.998632060928330.002735878143340820.00136793907167041
2100.9981504592566640.003699081486672420.00184954074333621
2110.9976855042398430.004628991520314030.00231449576015701
2120.9983172245929920.003365550814016810.00168277540700840
2130.9998222688897040.0003554622205926550.000177731110296327
2140.9999536766482829.26467034366133e-054.63233517183067e-05
2150.9998809654382340.0002380691235315220.000119034561765761
2160.9996481407454070.0007037185091868590.000351859254593429
2170.999037635007790.001924729984422170.000962364992211085
2180.9974780184314370.005043963137125280.00252198156856264
2190.993408460851910.01318307829617810.00659153914808903
2200.9849140327289980.03017193454200380.0150859672710019
2210.968042596512570.06391480697485820.0319574034874291
2220.9368773005467560.1262453989064880.0631226994532441
2230.8677590678015470.2644818643969070.132240932198453
2240.7388165721652090.5223668556695820.261183427834791


Meta Analysis of Goldfeld-Quandt test for Heteroskedasticity
Description# significant tests% significant testsOK/NOK
1% type I error level1150.550239234449761NOK
5% type I error level1310.626794258373206NOK
10% type I error level1400.669856459330144NOK
 
Charts produced by software:
http://www.freestatistics.org/blog/date/2009/Dec/28/t1262021332ghmyfyued1fxe2e/107dpi1262021262.png (open in new window)
http://www.freestatistics.org/blog/date/2009/Dec/28/t1262021332ghmyfyued1fxe2e/107dpi1262021262.ps (open in new window)


http://www.freestatistics.org/blog/date/2009/Dec/28/t1262021332ghmyfyued1fxe2e/1apew1262021262.png (open in new window)
http://www.freestatistics.org/blog/date/2009/Dec/28/t1262021332ghmyfyued1fxe2e/1apew1262021262.ps (open in new window)


http://www.freestatistics.org/blog/date/2009/Dec/28/t1262021332ghmyfyued1fxe2e/2pelp1262021262.png (open in new window)
http://www.freestatistics.org/blog/date/2009/Dec/28/t1262021332ghmyfyued1fxe2e/2pelp1262021262.ps (open in new window)


http://www.freestatistics.org/blog/date/2009/Dec/28/t1262021332ghmyfyued1fxe2e/3ijxu1262021262.png (open in new window)
http://www.freestatistics.org/blog/date/2009/Dec/28/t1262021332ghmyfyued1fxe2e/3ijxu1262021262.ps (open in new window)


http://www.freestatistics.org/blog/date/2009/Dec/28/t1262021332ghmyfyued1fxe2e/45o681262021262.png (open in new window)
http://www.freestatistics.org/blog/date/2009/Dec/28/t1262021332ghmyfyued1fxe2e/45o681262021262.ps (open in new window)


http://www.freestatistics.org/blog/date/2009/Dec/28/t1262021332ghmyfyued1fxe2e/5irm91262021262.png (open in new window)
http://www.freestatistics.org/blog/date/2009/Dec/28/t1262021332ghmyfyued1fxe2e/5irm91262021262.ps (open in new window)


http://www.freestatistics.org/blog/date/2009/Dec/28/t1262021332ghmyfyued1fxe2e/6ashq1262021262.png (open in new window)
http://www.freestatistics.org/blog/date/2009/Dec/28/t1262021332ghmyfyued1fxe2e/6ashq1262021262.ps (open in new window)


http://www.freestatistics.org/blog/date/2009/Dec/28/t1262021332ghmyfyued1fxe2e/7205c1262021262.png (open in new window)
http://www.freestatistics.org/blog/date/2009/Dec/28/t1262021332ghmyfyued1fxe2e/7205c1262021262.ps (open in new window)


http://www.freestatistics.org/blog/date/2009/Dec/28/t1262021332ghmyfyued1fxe2e/8issd1262021262.png (open in new window)
http://www.freestatistics.org/blog/date/2009/Dec/28/t1262021332ghmyfyued1fxe2e/8issd1262021262.ps (open in new window)


http://www.freestatistics.org/blog/date/2009/Dec/28/t1262021332ghmyfyued1fxe2e/9tei01262021262.png (open in new window)
http://www.freestatistics.org/blog/date/2009/Dec/28/t1262021332ghmyfyued1fxe2e/9tei01262021262.ps (open in new window)


 
Parameters (Session):
par1 = 1 ; par2 = Include Monthly Dummies ; par3 = No Linear Trend ;
 
Parameters (R input):
par1 = 1 ; par2 = Include Monthly Dummies ; par3 = No 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