Home » date » 2009 » Dec » 29 »

Multiple Regression 5

*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:15:55 -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/t1262085442u332kmk510s44k4.htm/, Retrieved Tue, 29 Dec 2009 12:17:35 +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/t1262085442u332kmk510s44k4.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 2lags
 
Dataseries X:
» Textbox « » Textfile « » CSV «
4.4 11.2 4.5 4.6 4.4 11.1 4.4 4.5 4.3 10.8 4.4 4.4 4.1 10.4 4.3 4.4 3.9 10.1 4.1 4.3 3.7 9.8 3.9 4.1 3.6 9.7 3.7 3.9 3.9 10.3 3.6 3.7 4.2 10.9 3.9 3.6 4.2 10.8 4.2 3.9 4.1 10.6 4.2 4.2 4.1 10.4 4.1 4.2 4.1 10.3 4.1 4.1 4.1 10.2 4.1 4.1 4.1 10 4.1 4.1 4 9.7 4.1 4.1 3.9 9.4 4 4.1 3.8 9.2 3.9 4 3.8 9.1 3.8 3.9 4 9.6 3.8 3.8 4.4 10.2 4 3.8 4.6 10.2 4.4 4 4.6 10 4.6 4.4 4.6 9.9 4.6 4.6 4.7 9.9 4.6 4.6 4.8 9.9 4.7 4.6 4.8 9.7 4.8 4.7 4.7 9.5 4.8 4.8 4.7 9.4 4.7 4.8 4.7 9.3 4.7 4.7 4.6 9.3 4.7 4.7 5 9.9 4.6 4.7 5.4 10.5 5 4.6 5.5 10.6 5.4 5 5.6 10.6 5.5 5.4 5.6 10.5 5.6 5.5 5.8 10.6 5.6 5.6 6 10.8 5.8 5.6 6.1 10.8 6 5.8 6.1 10.7 6.1 6 6 10.6 6.1 6.1 6 10.6 6 6.1 6.1 10.8 6 6 6.5 11.4 6.1 6 7.1 12.2 6.5 6.1 7.4 12.4 7.1 6.5 7.4 12.4 7.4 7.1 7.5 12.3 7.4 7.4 7.6 12.4 7.5 7.4 7.8 12.5 7.6 7.5 7.8 12.5 7.8 7.6 7.7 12.4 7.8 7.8 7.6 12.3 7.7 7.8 7.5 12.2 7.6 7.7 7.3 12.1 7.5 7.6 7.6 12.6 7.3 7.5 8 13.2 7.6 7.3 8 13.4 8 7.6 7.9 13.2 8 8 7.8 12.9 7.9 8 7.7 12.8 7.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 time8 seconds
R Server'Gwilym Jenkins' @ 72.249.127.135


Multiple Linear Regression - Estimated Regression Equation
Y[t] = -0.230404073016308 + 0.063966790884237X[t] + 1.41748644722477Y1[t] -0.518663070635735Y2[t] + 0.151072539337104M1[t] + 0.138136528285187M2[t] -0.0168673683779487M3[t] -0.0294935784935763M4[t] + 0.0153822127343551M5[t] -0.0159516121276719M6[t] -0.00533356626561581M7[t] + 0.424661480534343M8[t] + 0.0636908019214387M9[t] -0.0461693671096149M10[t] + 0.00929006307398224M11[t] + 0.00180052061064903t + e[t]


Multiple Linear Regression - Ordinary Least Squares
VariableParameterS.D.T-STAT
H0: parameter = 0
2-tail p-value1-tail p-value
(Intercept)-0.2304040730163080.137337-1.67770.0948210.04741
X0.0639667908842370.0191563.33920.0009850.000493
Y11.417486447224770.05839124.275600
Y2-0.5186630706357350.057404-9.035300
M10.1510725393371040.0598032.52620.0122280.006114
M20.1381365282851870.0598992.30620.0220240.011012
M3-0.01686736837794870.060112-0.28060.7792810.38964
M4-0.02949357849357630.05995-0.4920.6232260.311613
M50.01538221273435510.0606540.25360.8000360.400018
M6-0.01595161212767190.060902-0.26190.7936230.396811
M7-0.005333566265615810.060819-0.08770.9301970.465099
M80.4246614805343430.0616236.891300
M90.06369080192143870.0654740.97280.3317260.165863
M10-0.04616936710961490.06233-0.74070.4596470.229823
M110.009290063073982240.0607540.15290.8786060.439303
t0.001800520610649030.0005173.4850.0005930.000296


Multiple Linear Regression - Regression Statistics
Multiple R0.988826313218369
R-squared0.977777477713031
Adjusted R-squared0.976275955936885
F-TEST (value)651.191007180978
F-TEST (DF numerator)15
F-TEST (DF denominator)222
p-value0
Multiple Linear Regression - Residual Statistics
Residual Standard Deviation0.186199809715552
Sum Squared Residuals7.6968219486599


Multiple Linear Regression - Actuals, Interpolation, and Residuals
Time or IndexActualsInterpolation
Forecast
Residuals
Prediction Error
14.44.63173593242199-0.231735932421992
24.44.52432142523339-0.124321425233390
34.34.4037943189792-0.103794318979201
44.14.22563326839806-0.125633268398057
53.94.02148856058998-0.121488560589985
63.73.79300054375553-0.093000543755531
73.63.619257755822-0.0192577558220039
83.94.05141736716782-0.151417367167825
94.24.20773952492712-0.00773952492711593
104.24.362930210395-0.162930210395000
114.14.25179788182168-0.151797881821678
124.14.089766336459020.0102336635409818
134.14.28810902438192-0.188109024381922
144.14.27057685485223-0.17057685485223
154.14.1045801206229-0.00458012062289699
1644.07456439385265-0.0745643938526464
173.93.96030202370348-0.0603020237034797
183.83.82809302361635-0.0280930236163505
193.83.744232573341730.055767426658272
2044.25987784325803-0.259877843258028
214.44.222585049231270.177414950768731
224.64.577787365573630.0222126344263732
234.64.69828601938168-0.0982860193816847
244.64.580667183702780.0193328162972189
254.74.73354024365053-0.0335402436505339
264.84.86415339793174-0.0641533979317437
274.84.788039001361310.0119609986386873
284.74.71255364661591-0.0125536466159133
294.74.611084634643590.0889153653564066
304.74.627020958367360.072979041632635
314.64.63943952484007-0.0394395248400707
3254.967866522058740.0321334779412574
335.45.265937324540510.134062675459488
345.55.52380370584415-0.0238037058441463
355.65.515347073106580.0846529268934248
365.65.591343189213720.00865681078627844
375.85.698746621186320.101253378813675
3865.983901778366860.0160982216331417
396.16.010463077632180.0895369223678201
406.16.031256739634110.0687432603658929
4166.01967006532069-0.0196700653206902
4265.848388116346840.151611883653164
436.15.925466348059960.174533651940038
446.56.53739063472359-0.0373906347235882
457.16.744522181255060.355477818744941
467.47.292282531092070.107717468907931
477.47.46359057367231-0.0635905736723066
487.57.294105430929830.205894569070171
497.67.595123814688480.00487618531151684
507.87.680267340994540.119732659005459
517.87.758694947323440.0413050526765637
527.77.637739964602890.0622600353971133
537.67.536270952630570.0637290473694326
547.57.410458631631860.0895413683681394
557.37.32659818135724-0.0265981813572393
567.67.558746161828580.0412538381714152
5787.766934626651450.233065373348550
5888.08306399410708-0.0830639941070812
597.97.92006535847019-0.0200653584701852
607.87.75163713401910.0483628659808955
617.77.80823117721953-0.108231177219529
627.87.700816670030930.0991833299690657
637.77.73483156667607-0.0348315666760743
647.57.5175975672082-0.0175975672081985
657.37.31345285940013-0.0134528594001268
667.17.097758200742520.00224179925748228
6776.924015412808990.0759845871910085
687.37.35617502415481-0.0561750241548124
697.87.518893861002530.281106138997473
707.97.97677187318064-0.0767718731806352
717.97.91005225429107-0.0100522542910680
727.87.8315063674989-0.0315063674988899
737.87.83623410363574-0.0362341036357417
747.97.883361599346470.0166384006535300
757.87.86551018892804-0.0655101889280386
767.67.64187951037174-0.0418795103717374
777.47.44413148165209-0.0441314816520891
787.27.21564346481763-0.0156434648176338
796.97.03550399779568-0.135503997795684
807.17.1713729615197-0.0713729615197032
817.57.296075767772090.203924232227913
827.67.65127808411444-0.0512780841144451
837.47.61723473502318-0.217234735023178
847.37.248794879697620.0512051203023759
857.27.35725522996162-0.157255229961622
867.37.249840722773030.0501592772269729
877.27.27745894032974-0.0774589403297435
887.17.060224940861870.0397750591381321
8977.0042255568647-0.004225556864696
906.96.872016556777570.0279834432224327
916.86.794552785591370.00544721440863035
927.27.174846089873620.0251539101263849
937.67.479313571443810.120686428556191
947.77.73717995274744-0.0371799527474412
957.67.71592996183302-0.115929961833024
967.57.495635430318370.00436456968163192
977.57.55222947351879-0.0522294735187948
987.67.59296029014110.00703970985890022
997.67.575108879722670.0248911202773332
1007.67.506020204065690.0939797959343088
1017.57.533506478639-0.0335064786390002
1027.37.3494311714883-0.0494311714882984
1037.27.12382207649120.076177923508802
1047.47.54318832966017-0.143188329660171
10587.557761842696990.442238157303014
1068.28.19646144848430.00353855151570254
10788.19403745089994-0.194037450899937
1087.77.76733460942239-0.0673346094223852
1097.77.585900991153010.114099008846994
1107.87.723967742814040.076032257185965
1117.87.70611633239560.0938836676043984
1127.77.624234298561780.0757657014382217
1137.57.50997192841261-0.00997192841261137
1147.37.2360142836030.0639857163969952
1157.17.062271495669480.0377285043305208
1167.17.33988910411597-0.239889104115974
1177.27.110038276594560.0899617234054397
1186.87.12453723563136-0.324537235631363
1196.66.543746263206850.0562537367931451
1206.46.43463794319917-0.034637943199166
1216.46.388556290563840.0114437094361579
1226.56.44917001880760.050829981192398
1236.36.40573189203547-0.105731892035474
1245.96.0403525687567-0.140352568756696
1255.55.63656027400936-0.136560274009363
1265.25.23470426094552-0.0347042609455224
1274.95.02294544241667-0.122945442416667
1285.45.229870750469530.17012924953047
1295.85.76702613271250.0329738672875002
1305.75.95383616968729-0.253836169687288
1315.65.65548556841634-0.0554855684163403
1325.55.53252697194041-0.0325269719404072
1335.45.56993097787556-0.169930977875563
1345.45.43053307524485-0.0305330752448483
1355.45.303609289902240.096390710097759
1365.55.279990242220420.220009757779585
1375.85.481208556958320.31879144304168
1385.75.81865820072238-0.118658200722376
1395.45.51453916401661-0.114539164016613
1405.65.592145141588630.00785485841136534
1415.85.691261231487320.108738768512676
1426.25.788552974738420.411447025261579
1436.86.35385164391440.446148356085604
1446.76.98938874153163-0.289388741531631
1456.76.71490203072916-0.0149020307291611
1466.46.77482288461674-0.374822884616737
1476.36.18997689530840.110023104691604
1486.36.167414765917960.132585234082035
1496.46.221180631201150.178819368798847
1506.36.32060261349541-0.0206026134954062
15166.14580290727048-0.145802907270483
1526.36.293772354815170.00622764518483442
1536.36.52824041034791-0.228240410347909
1546.66.251788482559940.348211517440063
1557.56.70231097207950.797689027920504
1567.87.795770273662470.00422972633753365
1577.97.91988586238234-0.0198858623823358
1587.87.92048681182652-0.120486811826520
1597.67.58006516307640.0199348369235944
1607.57.312021774836350.187978225163649
1617.67.288698660637480.311301339362516
1627.57.433590270906880.0664097290931154
1637.37.252393885593540.0476061144064601
1647.67.529318619683850.0706813803161488
1657.57.71831704724144-0.218317047241445
1667.67.312909832907840.287090167092157
1677.97.544594698222870.355405301777132
1687.97.90408810377497-0.00408810377496998
1698.17.920552279797270.179447720202726
1708.28.2312941533315-0.0312941533315014
17188.10971012878592-0.109710128785921
1727.57.7315374473303-0.231537447330297
1736.87.1092363587994-0.309236358799402
1746.56.327604039543280.172395960456721
1756.66.309824216735690.290175783264313
1767.67.160504252739540.439495747260458
17788.21173098851745-0.211730988517451
1788.18.13920949017437-0.0392094901743731
1797.78.07318274564099-0.373182745640988
1807.57.395658884516780.104341115483216
1817.67.472499883273880.127500116726124
1827.87.72603568894750.0739643110524969
1837.87.81085997436482-0.0108599743648217
1847.87.689904991644270.110095008355728
1857.57.71099458712916-0.210994587129158
1867.57.24342199053350.256578009466499
1877.17.41783615728535-0.317836157285351
1887.57.353000615778710.146999384221290
1897.57.7810836230975-0.281083623097506
1907.67.459162067334380.140837932665617
1917.77.619790588320560.0802094116794367
1927.77.676596667162440.0234033328375604
1937.97.784000099135040.115999900864957
1948.18.062758577227150.0372414227728468
1958.28.089319876492470.110680123507526
1968.28.103716859405980.0962831405940225
1978.28.085733506004140.114266493995863
1987.98.05620020175276-0.156200201752758
1997.37.64337283405803-0.343372834058034
2006.97.43144888703189-0.531448887031887
2016.66.82927535069801-0.229275350698013
2026.76.490641638187620.209358361812378
2036.96.813265759452950.0867342405470516
20477.01821716210573-0.0182171621057252
2057.17.2091062526488-0.109106252648809
2067.27.29424977895487-0.0942497789548663
2077.17.23092874056128-0.130928740561285
2086.97.02648809927025-0.126488099270254
20976.847930107815880.152069892184121
2106.87.03829134606043-0.23829134606043
2116.46.68336292058249-0.283362920582488
2126.76.658293202318760.0417067976812409
2136.66.91264416947296-0.312644169472958
2146.46.48165023878566-0.0816502387856596
2156.36.32007256537537-0.0200725653753743
2166.26.27456699231671-0.0745669923167094
2176.56.343954393693980.156045606306016
2186.86.8035344653953-0.00353446539529628
2196.86.90078806505425-0.100788065054250
2206.46.70238005891643-0.302380058916432
2216.16.16287175459983-0.0628717545998337
2225.85.89636770717005-0.0963677071700458
2236.15.664725977019740.435274022980264
2247.26.7605232279380.439476772061997
2257.37.8369726361344-0.536972636134396
2266.97.27454553838346-0.374545538383463
2276.16.68735788687053-0.587357886870534
2285.85.727757698527980.072242301472019
2296.25.889505318082160.310494681917837
2307.16.632946723163640.467053276836357
2317.77.554412600447580.145587399552421
2327.97.91448865752846-0.0144886575284553
2337.77.90148102098843-0.201481020988431
2347.47.45273441772283-0.0527344177228324
2357.57.150036343243680.349963656756322
23687.930352909274870.0696470907251278
2378.18.25364638417512-0.153646384175125
23888.0216071660709-0.0216071660709046


Goldfeld-Quandt test for Heteroskedasticity
p-valuesAlternative Hypothesis
breakpoint indexgreater2-sidedless
190.002688269616597960.005376539233195930.997311730383402
200.000595943647725890.001191887295451780.999404056352274
210.0002032529275160150.0004065058550320290.999796747072484
220.001580776938148610.003161553876297220.998419223061851
230.0007726307873023630.001545261574604730.999227369212698
240.0003349457359412840.0006698914718825680.999665054264059
258.71794458079156e-050.0001743588916158310.999912820554192
263.19607863468172e-056.39215726936344e-050.999968039213653
271.72979681605025e-053.45959363210049e-050.99998270203184
289.08787875022665e-061.81757575004533e-050.99999091212125
292.37285125293653e-064.74570250587306e-060.999997627148747
305.87692787045742e-071.17538557409148e-060.999999412307213
312.21842628124059e-064.43685256248119e-060.999997781573719
327.10599998690887e-071.42119999738177e-060.999999289400001
332.01153697158545e-074.0230739431709e-070.999999798846303
341.07558271057395e-072.1511654211479e-070.999999892441729
353.79418822330369e-087.58837644660739e-080.999999962058118
361.40454389109792e-082.80908778219585e-080.99999998595456
375.15508940663013e-091.03101788132603e-080.99999999484491
381.79121162491179e-093.58242324982358e-090.999999998208788
394.66933835102045e-109.3386767020409e-100.999999999533066
401.20852544034007e-102.41705088068013e-100.999999999879147
412.00848560160978e-104.01697120321956e-100.999999999799151
426.76003748340286e-111.35200749668057e-100.9999999999324
432.15569125048425e-114.31138250096850e-110.999999999978443
445.72234345148926e-121.14446869029785e-110.999999999994278
457.73611835587272e-121.54722367117454e-110.999999999992264
465.40944105343834e-121.08188821068767e-110.99999999999459
474.42865746314888e-128.85731492629776e-120.999999999995571
481.36256031523922e-122.72512063047843e-120.999999999998637
498.74689250001744e-131.74937850000349e-120.999999999999125
502.48661089926106e-134.97322179852212e-130.999999999999751
513.46095648554570e-136.92191297109139e-130.999999999999654
525.41059622235629e-131.08211924447126e-120.999999999999459
535.49361195189866e-131.09872239037973e-120.99999999999945
546.38756718643511e-131.27751343728702e-120.99999999999936
551.63678382950217e-113.27356765900434e-110.999999999983632
566.27711131690048e-121.25542226338010e-110.999999999993723
573.8028754904658e-127.6057509809316e-120.999999999996197
581.29978723278069e-112.59957446556139e-110.999999999987002
596.75095115175425e-121.35019023035085e-110.99999999999325
605.01851822544748e-121.00370364508950e-110.999999999994981
616.86060300809421e-121.37212060161884e-110.99999999999314
622.9652183624627e-125.9304367249254e-120.999999999997035
632.05853389586361e-124.11706779172721e-120.999999999997941
648.34507228163878e-131.66901445632776e-120.999999999999166
653.2567669834301e-136.5135339668602e-130.999999999999674
661.24777290531126e-132.49554581062252e-130.999999999999875
675.33283942666763e-141.06656788533353e-130.999999999999947
682.14194891764195e-144.2838978352839e-140.999999999999979
693.09726022483851e-146.19452044967702e-140.99999999999997
701.15015555938443e-142.30031111876887e-140.999999999999988
714.8992399401395e-159.798479880279e-150.999999999999995
722.8777040492615e-155.755408098523e-150.999999999999997
739.96238666748623e-161.99247733349725e-150.999999999999999
743.57724707600502e-167.15449415201004e-161
751.61405157901376e-163.22810315802751e-161
765.72187595470803e-171.14437519094161e-161
771.96004694346778e-173.92009388693555e-171
787.27070077079389e-181.45414015415878e-171
791.92765870139418e-173.85531740278837e-171
806.8483105510783e-181.36966211021566e-171
814.15645319239058e-188.31290638478115e-181
821.53781021376390e-183.07562042752779e-181
833.79548298232347e-187.59096596464693e-181
841.30012264216112e-182.60024528432225e-181
858.71584117080547e-191.74316823416109e-181
864.11489879170549e-198.22979758341097e-191
871.80526493086474e-193.61052986172949e-191
889.19905386594297e-201.83981077318859e-191
893.32910044286033e-206.65820088572065e-201
901.08696490480745e-202.17392980961489e-201
913.49791330017655e-216.9958266003531e-211
922.89982690434283e-215.79965380868566e-211
931.50434973971522e-213.00869947943045e-211
944.9792872822386e-229.9585745644772e-221
952.12494090154077e-224.24988180308154e-221
967.62536516844007e-231.52507303368801e-221
972.63624774044348e-235.27249548088696e-231
988.14995894318067e-241.62999178863613e-231
992.87042320623896e-245.74084641247793e-241
1002.78095750860788e-245.56191501721575e-241
1019.63452615927495e-251.92690523185499e-241
1027.9857965429004e-251.59715930858008e-241
1032.56604254012549e-255.13208508025098e-251
1044.00585118685726e-258.01170237371453e-251
1051.95517893083641e-233.91035786167281e-231
1066.37323758314536e-241.27464751662907e-231
1071.84274584796949e-223.68549169593898e-221
1082.41079032549471e-214.82158065098941e-211
1091.21159880241212e-212.42319760482425e-211
1104.71412773754827e-229.42825547509654e-221
1111.74926161340500e-223.49852322680999e-221
1126.94820094116785e-231.38964018823357e-221
1135.09606625671806e-231.01921325134361e-221
1142.40317625929207e-234.80635251858413e-231
1151.2298785011732e-232.4597570023464e-231
1162.77228818867019e-225.54457637734039e-221
1178.77634569398608e-221.75526913879722e-211
1183.37391465258046e-206.74782930516092e-201
1191.25220018231002e-192.50440036462004e-191
1204.92349916916862e-209.84699833833723e-201
1212.94379187533154e-205.88758375066308e-201
1221.17789060203440e-202.35578120406880e-201
1231.14134633225625e-202.28269266451250e-201
1241.30551229138812e-202.61102458277623e-201
1258.03585153742279e-211.60717030748456e-201
1263.41632182182146e-216.83264364364291e-211
1272.06158148570339e-214.12316297140677e-211
1281.75784521581407e-193.51569043162814e-191
1291.37530544728563e-192.75061089457126e-191
1305.95801155368229e-191.19160231073646e-181
1316.21245942564001e-191.24249188512800e-181
1322.99731956497135e-195.99463912994269e-191
1335.29699144985095e-191.05939828997019e-181
1342.76067775231877e-195.52135550463754e-191
1352.76489396599679e-195.52978793199358e-191
1361.01634926912176e-182.03269853824352e-181
1371.49318948398553e-172.98637896797105e-171
1385.048034868513e-171.0096069737026e-161
1398.66524093694104e-171.73304818738821e-161
1406.56970040543e-171.313940081086e-161
1415.30390196706113e-171.06078039341223e-161
1421.40081376242263e-142.80162752484525e-140.999999999999986
1431.27506727698361e-122.55013455396722e-120.999999999998725
1443.90991888260276e-117.81983776520552e-110.9999999999609
1452.85349340292087e-115.70698680584175e-110.999999999971465
1461.16884642595101e-092.33769285190202e-090.999999998831154
1478.36137468808692e-101.67227493761738e-090.999999999163863
1486.1429711864961e-101.22859423729922e-090.999999999385703
1494.93184943864322e-109.86369887728645e-100.999999999506815
1504.20691008288364e-108.41382016576728e-100.999999999579309
1511.05317243554020e-092.10634487108041e-090.999999998946828
1521.46484465425378e-092.92968930850756e-090.999999998535155
1539.54553504484628e-091.90910700896926e-080.999999990454465
1543.24187854226299e-086.48375708452598e-080.999999967581215
1557.31987528339226e-050.0001463975056678450.999926801247166
1565.81635580422054e-050.0001163271160844110.999941836441958
1574.97195554458709e-059.94391108917418e-050.999950280444554
1585.12931469380408e-050.0001025862938760820.999948706853062
1593.37706308361398e-056.75412616722796e-050.999966229369164
1602.86921217795681e-055.73842435591362e-050.99997130787822
1614.34288315880779e-058.68576631761557e-050.999956571168412
1622.81133576469533e-055.62267152939066e-050.999971886642353
1631.80781447881686e-053.61562895763371e-050.999981921855212
1641.20803365582477e-052.41606731164955e-050.999987919663442
1652.53011855637146e-055.06023711274293e-050.999974698814436
1663.46271711097364e-056.92543422194728e-050.99996537282889
1670.0001698306029073890.0003396612058147780.999830169397093
1680.0001135257169307270.0002270514338614540.99988647428307
1690.0001035145900805890.0002070291801611790.99989648540992
1706.83757515630365e-050.0001367515031260730.999931624248437
1715.07967942329003e-050.0001015935884658010.999949203205767
1726.18036302818895e-050.0001236072605637790.999938196369718
1730.0001315674808595190.0002631349617190380.99986843251914
1740.0001175020468295990.0002350040936591990.99988249795317
1750.0001644051089036980.0003288102178073970.999835594891096
1760.000461215426710930.000922430853421860.99953878457329
1770.0005320117079571880.001064023415914380.999467988292043
1780.0003648419091997310.0007296838183994620.9996351580908
1790.0006774522620806650.001354904524161330.999322547737919
1800.0005855435353070740.001171087070614150.999414456464693
1810.0004446347561056080.0008892695122112170.999555365243894
1820.0003045624656898100.0006091249313796210.99969543753431
1830.0001999192446264240.0003998384892528480.999800080755374
1840.0001847802006600980.0003695604013201960.99981521979934
1850.0001767115261608990.0003534230523217970.99982328847384
1860.000451833585176260.000903667170352520.999548166414824
1870.000866761354295120.001733522708590240.999133238645705
1880.0007185054848514460.001437010969702890.999281494515148
1890.0007535011227721880.001507002245544380.999246498877228
1900.0006656738000781160.001331347600156230.999334326199922
1910.000814302332465930.001628604664931860.999185697667534
1920.0006244983892535590.001248996778507120.999375501610746
1930.0005347689264550320.001069537852910060.999465231073545
1940.0003938479188212150.000787695837642430.999606152081179
1950.0006994014838366480.001398802967673300.999300598516163
1960.00291530931512530.00583061863025060.997084690684875
1970.01270761736807790.02541523473615580.987292382631922
1980.03040480407281380.06080960814562770.969595195927186
1990.02727416643132340.05454833286264680.972725833568677
2000.0881300284032520.1762600568065040.911869971596748
2010.07375119969733550.1475023993946710.926248800302665
2020.07575534002835010.15151068005670.92424465997165
2030.14642617952390.29285235904780.8535738204761
2040.1655564246033310.3311128492066630.834443575396669
2050.1353514443899930.2707028887799860.864648555610007
2060.1022900296798320.2045800593596630.897709970320168
2070.07554215172444880.1510843034488980.924457848275551
2080.05536302664664880.1107260532932980.944636973353351
2090.08470857391707620.1694171478341520.915291426082924
2100.0729493297930170.1458986595860340.927050670206983
2110.1339149717628670.2678299435257340.866085028237133
2120.1732030854106650.3464061708213290.826796914589335
2130.1398944813800310.2797889627600630.860105518619969
2140.1385112113172340.2770224226344680.861488788682766
2150.6468339986846550.7063320026306910.353166001315346
2160.6978209043522150.6043581912955710.302179095647785
2170.7213914541471820.5572170917056370.278608545852818
2180.7342242706027570.5315514587944870.265775729397243
2190.7491410097662370.5017179804675260.250858990233763


Meta Analysis of Goldfeld-Quandt test for Heteroskedasticity
Description# significant tests% significant testsOK/NOK
1% type I error level1780.885572139303483NOK
5% type I error level1790.890547263681592NOK
10% type I error level1810.90049751243781NOK
 
Charts produced by software:
http://www.freestatistics.org/blog/date/2009/Dec/29/t1262085442u332kmk510s44k4/103na41262085345.png (open in new window)
http://www.freestatistics.org/blog/date/2009/Dec/29/t1262085442u332kmk510s44k4/103na41262085345.ps (open in new window)


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


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


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


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


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


http://www.freestatistics.org/blog/date/2009/Dec/29/t1262085442u332kmk510s44k4/6ss8v1262085345.png (open in new window)
http://www.freestatistics.org/blog/date/2009/Dec/29/t1262085442u332kmk510s44k4/6ss8v1262085345.ps (open in new window)


http://www.freestatistics.org/blog/date/2009/Dec/29/t1262085442u332kmk510s44k4/7qsbr1262085345.png (open in new window)
http://www.freestatistics.org/blog/date/2009/Dec/29/t1262085442u332kmk510s44k4/7qsbr1262085345.ps (open in new window)


http://www.freestatistics.org/blog/date/2009/Dec/29/t1262085442u332kmk510s44k4/8yrk11262085345.png (open in new window)
http://www.freestatistics.org/blog/date/2009/Dec/29/t1262085442u332kmk510s44k4/8yrk11262085345.ps (open in new window)


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