Home » date » 2010 » Dec » 08 »

*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: Wed, 08 Dec 2010 19:13:34 +0000
 
Cite this page as follows:
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL http://www.freestatistics.org/blog/date/2010/Dec/08/t1291835543aetyp8l44oscs7s.htm/, Retrieved Wed, 08 Dec 2010 20:12:34 +0100
 
BibTeX entries for LaTeX users:
@Manual{KEY,
    author = {{YOUR NAME}},
    publisher = {Office for Research Development and Education},
    title = {Statistical Computations at FreeStatistics.org, URL http://www.freestatistics.org/blog/date/2010/Dec/08/t1291835543aetyp8l44oscs7s.htm/},
    year = {2010},
}
@Manual{R,
    title = {R: A Language and Environment for Statistical Computing},
    author = {{R Development Core Team}},
    organization = {R Foundation for Statistical Computing},
    address = {Vienna, Austria},
    year = {2010},
    note = {{ISBN} 3-900051-07-0},
    url = {http://www.R-project.org},
}
 
Original text written by user:
 
IsPrivate?
No (this computation is public)
 
User-defined keywords:
 
Dataseries X:
» Textbox « » Textfile « » CSV «
1214,82 1192,72 1206,41 1176,86 1145,11 1199,07 1214,82 1192,72 1206,41 1176,86 1157,47 1199,07 1214,82 1192,72 1206,41 1100,1 1157,47 1199,07 1214,82 1192,72 1095,63 1100,1 1157,47 1199,07 1214,82 1105,63 1095,63 1100,1 1157,47 1199,07 1137,79 1105,63 1095,63 1100,1 1157,47 1124,72 1137,79 1105,63 1095,63 1100,1 1152,6 1124,72 1137,79 1105,63 1095,63 1211,85 1152,6 1124,72 1137,79 1105,63 1239,62 1211,85 1152,6 1124,72 1137,79 1244,13 1239,62 1211,85 1152,6 1124,72 1198,42 1244,13 1239,62 1211,85 1152,6 1227,99 1198,42 1244,13 1239,62 1211,85 1304,92 1227,99 1198,42 1244,13 1239,62 1340,26 1304,92 1227,99 1198,42 1244,13 1307,32 1340,26 1304,92 1227,99 1198,42 1356,51 1307,32 1340,26 1304,92 1227,99 1383,29 1356,51 1307,32 1340,26 1304,92 1437,87 1383,29 1356,51 1307,32 1340,26 1494,56 1437,87 1383,29 1356,51 1307,32 1521,42 1494,56 1437,87 1383,29 1356,51 1498,76 1521,42 1494,56 1437,87 1383,29 1488,75 1498,76 1521,42 1494,56 1437,87 1524,62 1488,75 1498,76 1521,42 1494,56 1439,27 1 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 time16 seconds
R Server'Gwilym Jenkins' @ 72.249.127.135


Multiple Linear Regression - Estimated Regression Equation
BEL_20[t] = + 46.3927350630606 + 1.28080975239333Y1[t] -0.353226463695596Y2[t] + 0.198103762183654Y3[t] -0.142505461121703Y4[t] + e[t]


Multiple Linear Regression - Ordinary Least Squares
VariableParameterS.D.T-STAT
H0: parameter = 0
2-tail p-value1-tail p-value
(Intercept)46.392735063060624.2295481.91470.0569070.028454
Y11.280809752393330.06878918.619500
Y2-0.3532264636955960.111746-3.1610.0018080.000904
Y30.1981037621836540.1118611.7710.0780350.039018
Y4-0.1425054611217030.068535-2.07930.0388210.019411


Multiple Linear Regression - Regression Statistics
Multiple R0.991510551722486
R-squared0.98309317417703
Adjusted R-squared0.98276647222876
F-TEST (value)3009.14389823190
F-TEST (DF numerator)4
F-TEST (DF denominator)207
p-value0
Multiple Linear Regression - Residual Statistics
Residual Standard Deviation120.02917403868
Sum Squared Residuals2982249.54242439


Multiple Linear Regression - Actuals, Interpolation, and Residuals
Time or IndexActualsInterpolation
Forecast
Residuals
Prediction Error
11214.821217.86016984901-3.04016984901188
21199.071252.33115344679-53.2611534467927
31157.471217.42901811850-59.9590181184959
41100.11176.03964212915-75.9396421291544
51095.631110.98430257890-15.3543025789030
61105.631119.52702971375-13.8970297137485
71137.791128.477063876599.31293612341217
81124.721173.42565536419-48.7056553641922
91152.61147.943745861014.65625413898837
101211.851193.2153540188518.6346459811517
111239.621252.08318623891-12.4631862389056
121244.131274.10828435544-29.9782843554448
131198.421277.84023309522-79.4202330952208
141227.991214.7592608664313.2307391335663
151304.921265.7148582123339.2051417876715
161340.261344.10462333340-3.84462333339493
171307.321374.56658100652-67.2465810065165
181356.511330.9199204751025.5900795249024
191383.291401.59627374094-18.3062737409363
201437.871406.9594682384730.9105317615273
211494.561481.8455137774912.7144862225049
221521.421533.47089337087-12.0508933708687
231498.761554.84524218439-56.0852421843949
241488.751519.78698459047-31.0369845904672
251524.621512.2126230976212.4073769023840
261439.271553.37433788075-114.104337880747
271423.111432.63314735077-9.52314735077478
281466.851450.6156020438716.2343979561264
291425.831490.32653327407-64.4965332740682
301363.451431.29907601870-67.8490760186979
311389.181376.8594600148412.3205399851644
321395.891417.48955655501-21.5995565550100
331368.431410.48313441288-42.0531344128777
341349.031386.92864950652-37.8986495065166
351299.881369.44314973276-69.563149732758
361365.411306.9478028446358.4621971553689
371451.041408.3103335856442.7296664143565
381433.751487.86694855155-54.116948551546
391464.651455.460848796449.189151203562
401475.571508.77039799117-33.2003979911698
411471.161496.21418607510-25.0541860751040
421429.121495.29390775776-66.173907757763
431452.461440.7662688064311.6937311935705
441538.091483.0802117343755.0097882656258
451631.591576.8118120905154.7781879094948
461665.51676.93541324795-11.4354132479503
471690.61700.97854528928-10.3785452892757
481711.741727.46891981875-15.7289198187506
491734.11739.07269170635-4.97269170635476
501748.091760.38443457152-12.2944345715173
511703.451771.01584573767-67.565845737674
521745.741710.3158948380535.4241051619517
531751.011779.83441812840-28.8244181284016
541795.651760.8093350288634.8406649711434
551852.131830.8624307992421.2675692007614
561877.11882.45198715091-5.3519871509135
571989.311902.5759241624186.7340758375856
582097.762042.3029783836555.4570216163496
592154.872138.4691970370016.4008029630034
602152.182191.97969379881-39.7996937988105
612250.272173.8553674395776.414632560432
622346.92296.2991638388350.6008361611696
632525.562376.74444038376148.815559616238
642409.362591.25697528246-181.896975282460
652394.362384.483847908889.87615209112049
662401.332427.93953214795-26.6095321479495
672354.322393.68549022782-39.3654902278204
682450.412344.60021346544105.809786534561
692504.672487.7967637704916.8732362295108
702661.392513.04584911457148.34415088543
712880.42720.34325782509160.056742174913
723064.422942.55151068328121.868489316723
733141.123124.2004687936816.9195312063165
743327.73178.49109204184149.208907958163
753564.953395.6170391547169.332960845297
763403.133622.55486291757-219.424862917566
773149.93357.52328135369-207.623281353690
783006.843110.75438275233-103.914382752333
793230.662951.10270552890279.557294471103
803361.133261.2005402268199.9294597731898
813484.743356.99457521908127.745424780925
823411.133533.95642731407-122.826427314065
833288.183389.96472380682-101.784723806823
843280.373264.3850832736715.9849167263318
853173.953265.61363483526-91.6636348352577
863165.263118.2015290997147.058470900290
873092.713160.63550868016-67.9255086801571
883053.053050.813064393312.23693560668846
893181.963039.08663903370142.873360966296
902999.933205.07073027562-205.140730275618
913249.572928.87248360864320.697516391363
923210.523344.1009659538-133.580965953799
933030.293151.47468390238-121.184683902381
942803.473009.82272791536-206.352727915356
952767.632739.6604502016627.9695497983353
962882.62743.73565236977138.864347630234
972863.362884.39984998074-21.0398499807447
982897.062844.3696736685852.6903263314226
993012.612922.2124247505990.3975752494082
1003142.953038.11089056352104.839109436476
1013032.933173.65421766801-140.724217668014
1023045.783004.7884471121440.9915528878638
1033110.523069.4631662965941.0568337034147
1043013.243107.47429188999-94.2342918899924
1052987.12978.233322094198.866677905814
1062995.552990.108867943295.44113205671213
1072833.182981.66771257377-148.48771257377
1082848.962779.4023683738869.5576316261244
1092794.832862.36599672107-67.5359967210699
1102845.262754.0915722146691.1684277853362
1112915.022864.0676455972950.9523544027085
1122892.632922.63163053658-30.0016305365788
1132604.422887.01741541053-282.597415410527
1142641.652532.41714524096109.232854759045
1152659.812667.52836722112-7.71836722112442
1162638.532623.7324630567614.7975369432362
1172720.252638.5091409611181.7408590388943
1182745.882748.98565907782-3.10565907782441
1192735.72746.14359918522-10.4435991852235
1202811.72743.2733172996668.4266827003411
1212799.432837.64255702387-38.2125570238746
1222555.282789.41269885356-234.132698853564
1232304.982497.54367803646-192.563678036455
1242214.952249.93608991644-34.9860899164423
1252065.812176.41888024230-110.608880242303
1261940.492002.40722895517-61.9172289551715
12720421912.41018079017129.589819209831
1281995.372069.97609075866-74.6060907586645
1291946.811970.82281466966-24.0128146696587
1301765.91963.0659403826-197.165940382599
1311635.251724.80401736509-89.5540173650926
1321833.421618.39353372254215.026466277459
1331910.431889.4417534115820.9882465884183
1341959.671917.9764305750741.6935694249317
1351969.62011.71809386120-42.1180938612034
1362061.411994.0593271253767.350672874627
1372093.482106.92321539705-13.4432153970478
1382120.882110.5192639772610.3607360227397
1392174.562151.0583056792623.5016943207373
1402196.722203.40352935012-6.6835293501233
1412350.442213.68296983764136.757030162361
1422440.252409.4711068593330.7788931406693
1432408.642466.94294493947-58.3029449394651
1442472.812422.0278692662250.7821307337752
1452407.62511.26867899281-103.668678992808
1462454.622386.0200574779368.5999425220725
1472448.052486.49454577843-38.4445457784335
1482497.842439.4079956100758.432004389933
1492645.642524.10783106583121.532168934167
1502756.762687.8222183426768.937781657327
1512849.272788.7387738931060.5312261068951
1522921.442890.1603485826531.2796514173491
1532981.852950.8703911564630.9796088435419
1543080.583005.2431266133975.3368733866063
1553106.223111.47303110376-5.25303110376167
1563119.313111.121773538828.18822646117862
1573061.263129.78087620253-68.520876202525
1583097.313041.8159519521655.4940480478393
1593161.693107.4332779672954.2567220327073
1603257.163163.7926759293193.367324070695
1613277.013278.74494590241-1.73494590240833
1623295.323278.0630873343417.2569126656557
1633363.993304.2416331849759.7483668150332
1643494.173376.05462563761118.115374362394
1653667.033519.33292442451147.69707557549
1663813.063705.74518753534107.314812464659
1673917.963847.7264069087670.2335930912414
1683895.513946.19454484359-50.6845448435939
1693801.063884.68250824288-83.6225082428776
1703570.123771.61097340476-201.490973404757
1713701.613489.7867563504211.823243649599
1723862.273724.2628974824138.007102517603
1733970.13851.30160255483118.798397445174
1744138.523991.62082937904146.899170620961
1754199.754182.3357056463617.4142943536390
1764290.894199.7358870622491.1541129377581
1774443.914312.83910327751131.070896722493
1784502.644464.7646752839037.8753247160960
1794356.984495.2654860682-138.285486068202
1804591.274305.28363728446285.986362715544
1814696.964646.6439691667950.3160308332056
1824621.44662.03018398665-40.6301839866545
1834562.844595.09077005682-32.2507700568223
1844202.524534.32632469250-331.806324692495
1854296.494063.4797739676233.010226032403
1864435.234310.27878212768124.951217872325
1874105.184391.75000859452-286.570008594522
1884116.683989.97548852775126.704511472252
1893844.494135.38087280676-290.890872806757
1903720.983697.5398075855823.4401924144227
1913674.43684.80382692911-10.4038269291122
1923857.623613.21003336201244.409966637994
1933801.063878.65405066987-77.5940506698651
1943504.373749.86647465682-245.496474656819
1953032.63432.7759936922-400.175993692203
1963047.032896.01253694362151.017463056377
1972962.343030.42097412711-68.080974127106
1982197.822865.67267170060-667.852671700603
1992014.451986.4711876929327.9788123070674
2001862.832002.82403799781-139.994037997805
2011905.411734.01329922554171.39670077446
2021810.991914.72836317312-103.738363173122
2031670.071774.84865751159-104.778657511587
2041864.441657.75052611551206.689473884490
2052052.021931.70535119224120.314648807758
2062029.62088.84160026986-59.2416002698595
2072070.832052.4549233980918.3750766019117
2082293.412122.64356403750170.766435962495
2092443.272361.9900108816881.2799891183245
2102513.172486.6738046391526.4961953608454
2112466.922561.48632370682-94.566323706816
2122502.662475.5273071106827.1326928893235


Goldfeld-Quandt test for Heteroskedasticity
p-valuesAlternative Hypothesis
breakpoint indexgreater2-sidedless
80.03751679466611780.07503358933223560.962483205333882
90.00986393914222380.01972787828444760.990136060857776
100.004464147827701990.008928295655403980.995535852172298
110.001252693687230390.002505387374460770.99874730631277
120.0002921482610636970.0005842965221273940.999707851738936
139.08039189769603e-050.0001816078379539210.999909196081023
140.0001572867337614810.0003145734675229620.999842713266239
150.0002165679430440180.0004331358860880360.999783432056956
166.49223378343536e-050.0001298446756687070.999935077662166
172.13037220063385e-054.2607444012677e-050.999978696277994
182.62684977893709e-055.25369955787417e-050.99997373150221
198.47801178924865e-061.69560235784973e-050.99999152198821
203.98854876607724e-067.97709753215448e-060.999996011451234
211.30506795987980e-062.61013591975959e-060.99999869493204
223.97036031661677e-077.94072063323353e-070.999999602963968
232.24174799027775e-074.48349598055551e-070.999999775825201
246.8355205633633e-081.36710411267266e-070.999999931644794
251.98157043972371e-083.96314087944742e-080.999999980184296
261.71634169432778e-073.43268338865556e-070.99999982836583
276.04502764317469e-081.20900552863494e-070.999999939549724
282.26640699767705e-084.5328139953541e-080.99999997733593
291.03699011289261e-082.07398022578522e-080.999999989630099
304.33165318919887e-098.66330637839775e-090.999999995668347
311.50255618648793e-093.00511237297585e-090.999999998497444
324.48900445541823e-108.97800891083646e-100.9999999995511
331.32381626840329e-102.64763253680658e-100.999999999867618
344.18907649725226e-118.37815299450453e-110.99999999995811
352.79854373570801e-115.59708747141601e-110.999999999972015
365.46777619156715e-111.09355523831343e-100.999999999945322
374.0887107197132e-118.1774214394264e-110.999999999959113
381.29790467910232e-112.59580935820463e-110.99999999998702
396.64159066245203e-121.32831813249041e-110.999999999993358
402.32318519977718e-124.64637039955437e-120.999999999997677
417.05228678593372e-131.41045735718674e-120.999999999999295
423.5759603035836e-137.1519206071672e-130.999999999999642
431.41650994671543e-132.83301989343086e-130.999999999999858
441.64403121101281e-133.28806242202562e-130.999999999999836
452.97917849129757e-135.95835698259513e-130.999999999999702
461.01896267503596e-132.03792535007192e-130.999999999999898
473.42477261994628e-146.84954523989256e-140.999999999999966
481.03384562816490e-142.06769125632981e-140.99999999999999
493.16878379698653e-156.33756759397305e-150.999999999999997
509.31002741257307e-161.86200548251461e-151
515.61462735371734e-161.12292547074347e-151
523.71778001328558e-167.43556002657117e-161
531.31425842875449e-162.62851685750897e-161
548.57190773701958e-171.71438154740392e-161
553.14962450759021e-176.29924901518042e-171
569.19431195216647e-181.83886239043329e-171
572.90590591768971e-175.81181183537942e-171
581.46272437321875e-172.92544874643751e-171
594.42479229763283e-188.84958459526567e-181
602.32467134561679e-184.64934269123358e-181
611.79475879104012e-183.58951758208023e-181
625.57395446226966e-191.11479089245393e-181
634.26036532225548e-188.52073064451095e-181
646.72335125856328e-151.34467025171266e-140.999999999999993
652.60651395534987e-155.21302791069974e-150.999999999999997
661.46716262646871e-152.93432525293742e-150.999999999999999
671.56422858184180e-153.12845716368361e-150.999999999999998
681.53129955955555e-153.06259911911111e-150.999999999999998
695.50079251875685e-161.10015850375137e-151
701.54883520844503e-153.09767041689005e-150.999999999999998
719.68870594048656e-151.93774118809731e-140.99999999999999
727.57259078370395e-151.51451815674079e-140.999999999999992
733.5090410149616e-157.0180820299232e-150.999999999999996
744.55834250365189e-159.11668500730377e-150.999999999999995
755.50079879539159e-151.10015975907832e-140.999999999999994
763.49601962743435e-116.99203925486871e-110.99999999996504
771.13896898923423e-092.27793797846847e-090.99999999886103
781.16113236449152e-092.32226472898303e-090.999999998838868
791.04424061715381e-082.08848123430761e-080.999999989557594
809.33912139915766e-091.86782427983153e-080.999999990660879
817.19031034100283e-091.43806206820057e-080.99999999280969
821.08991437459891e-082.17982874919781e-080.999999989100856
831.07860256531573e-082.15720513063146e-080.999999989213974
845.69260939072117e-091.13852187814423e-080.99999999430739
851.57910650780063e-083.15821301560125e-080.999999984208935
868.6662648607969e-091.73325297215938e-080.999999991333735
879.05722094051436e-091.81144418810287e-080.99999999094278
885.12707599004084e-091.02541519800817e-080.999999994872924
896.47657999629994e-091.29531599925999e-080.99999999352342
901.20353317501257e-072.40706635002513e-070.999999879646682
914.40809678957179e-068.81619357914358e-060.99999559190321
924.3263143665611e-068.6526287331222e-060.999995673685633
931.08932956399464e-052.17865912798928e-050.99998910670436
942.38291413386166e-054.76582826772333e-050.999976170858661
951.53160934361940e-053.06321868723880e-050.999984683906564
961.42091530447942e-052.84183060895884e-050.999985790846955
971.14250284753858e-052.28500569507716e-050.999988574971525
987.59378005602739e-061.51875601120548e-050.999992406219944
996.98923596072866e-061.39784719214573e-050.99999301076404
1006.03189305720379e-061.20637861144076e-050.999993968106943
1018.27920643471923e-061.65584128694385e-050.999991720793565
1025.64649065355321e-061.12929813071064e-050.999994353509347
1033.85267143079439e-067.70534286158879e-060.99999614732857
1044.17842749710755e-068.3568549942151e-060.999995821572503
1052.58908846951288e-065.17817693902576e-060.99999741091153
1061.58756849551228e-063.17513699102456e-060.999998412431504
1072.95046538246123e-065.90093076492246e-060.999997049534618
1082.04832556325021e-064.09665112650041e-060.999997951674437
1091.54464486242535e-063.0892897248507e-060.999998455355138
1101.12737166694604e-062.25474333389208e-060.999998872628333
1117.46629866775443e-071.49325973355089e-060.999999253370133
1124.80069463555475e-079.6013892711095e-070.999999519930536
1135.92907018491662e-061.18581403698332e-050.999994070929815
1145.37963821260629e-061.07592764252126e-050.999994620361787
1153.4009078686071e-066.8018157372142e-060.999996599092131
1162.23953406007660e-064.47906812015321e-060.99999776046594
1171.74693940607209e-063.49387881214417e-060.999998253060594
1181.07030712574746e-062.14061425149492e-060.999998929692874
1196.54062549104463e-071.30812509820893e-060.99999934593745
1204.79213706200488e-079.58427412400976e-070.999999520786294
1213.02591770214801e-076.05183540429602e-070.99999969740823
1221.56384810846496e-063.12769621692991e-060.999998436151892
1233.79740778164393e-067.59481556328785e-060.999996202592218
1242.55553876415477e-065.11107752830954e-060.999997444461236
1253.03264864774938e-066.06529729549877e-060.999996967351352
1262.37592961377598e-064.75185922755195e-060.999997624070386
1272.401213298195e-064.80242659639e-060.999997598786702
1281.91874882951051e-063.83749765902102e-060.99999808125117
1291.22762646047791e-062.45525292095582e-060.99999877237354
1303.04361793970907e-066.08723587941814e-060.99999695638206
1312.70160001853254e-065.40320003706509e-060.999997298399981
1326.5008020924498e-061.30016041848996e-050.999993499197908
1334.07971715488944e-068.15943430977888e-060.999995920282845
1342.60903172502592e-065.21806345005185e-060.999997390968275
1351.79881688528872e-063.59763377057744e-060.999998201183115
1361.28230569698621e-062.56461139397242e-060.999998717694303
1377.9974496644959e-071.59948993289918e-060.999999200255034
1384.8501519174436e-079.7003038348872e-070.999999514984808
1392.98853048920567e-075.97706097841133e-070.99999970114695
1401.79393556748714e-073.58787113497428e-070.999999820606443
1411.95633161331963e-073.91266322663926e-070.999999804366839
1421.17976850372866e-072.35953700745732e-070.99999988202315
1438.04312437755413e-081.60862487551083e-070.999999919568756
1445.18776691636772e-081.03755338327354e-070.999999948122331
1454.76867368929986e-089.53734737859973e-080.999999952313263
1463.06023077633849e-086.12046155267699e-080.999999969397692
1471.89460618546246e-083.78921237092492e-080.999999981053938
1481.13968097411689e-082.27936194823377e-080.99999998860319
1491.06210378232880e-082.12420756465760e-080.999999989378962
1506.75595164572333e-091.35119032914467e-080.999999993244048
1514.15443331781901e-098.30886663563802e-090.999999995845567
1522.35303015216498e-094.70606030432996e-090.99999999764697
1531.29987706904900e-092.59975413809801e-090.999999998700123
1548.46028661106091e-101.69205732221218e-090.999999999153971
1554.41804086198182e-108.83608172396364e-100.999999999558196
1562.26492371375355e-104.52984742750709e-100.999999999773508
1571.45789896744512e-102.91579793489024e-100.99999999985421
1588.11716676081666e-111.62343335216333e-100.999999999918828
1594.44475174825296e-118.88950349650593e-110.999999999955552
1603.02001824962491e-116.04003649924981e-110.9999999999698
1611.46293582230891e-112.92587164461782e-110.99999999998537
1627.02699785702397e-121.40539957140479e-110.999999999992973
1633.83177971638247e-127.66355943276494e-120.999999999996168
1643.23494858846526e-126.46989717693053e-120.999999999996765
1653.88825463940257e-127.77650927880513e-120.999999999996112
1662.97152347702674e-125.94304695405349e-120.999999999997028
1671.69374519475089e-123.38749038950178e-120.999999999998306
1688.68411068277513e-131.73682213655503e-120.999999999999132
1695.48539075779318e-131.09707815155864e-120.999999999999451
1701.76138516961273e-123.52277033922546e-120.999999999998239
1715.87038534034956e-121.17407706806991e-110.99999999999413
1726.33226655242789e-121.26645331048558e-110.999999999993668
1736.22116588278321e-121.24423317655664e-110.99999999999378
1747.84817810366307e-121.56963562073261e-110.999999999992152
1753.61446932337359e-127.22893864674717e-120.999999999996386
1762.54992678811299e-125.09985357622598e-120.99999999999745
1772.97804311245442e-125.95608622490884e-120.999999999997022
1781.71654277992977e-123.43308555985955e-120.999999999998283
1791.5379860837443e-123.0759721674886e-120.999999999998462
1804.84260235669389e-119.68520471338779e-110.999999999951574
1814.15787457658919e-118.31574915317838e-110.999999999958421
1823.03360554045220e-116.06721108090441e-110.999999999969664
1831.57910393511849e-113.15820787023697e-110.999999999984209
1842.54144233601949e-105.08288467203898e-100.999999999745856
1852.40090952412455e-094.8018190482491e-090.99999999759909
1867.07746041020643e-091.41549208204129e-080.99999999292254
1872.12677079981612e-084.25354159963223e-080.999999978732292
1886.03377260179185e-081.20675452035837e-070.999999939662274
1891.58449565856126e-073.16899131712252e-070.999999841550434
1901.98487668564259e-073.96975337128519e-070.99999980151233
1911.16928195387702e-072.33856390775405e-070.999999883071805
1922.34197186914023e-054.68394373828046e-050.999976580281309
1935.87286118646515e-050.0001174572237293030.999941271388135
1948.18004422836658e-050.0001636008845673320.999918199557716
1950.0002464304749311760.0004928609498623510.999753569525069
1960.002410161097669780.004820322195339560.99758983890233
1970.02095012587934200.04190025175868390.979049874120658
1980.3336624719442370.6673249438884730.666337528055763
1990.2507679993284180.5015359986568360.749232000671582
2000.3889458334088210.7778916668176410.611054166591179
2010.7100546624810910.5798906750378190.289945337518909
2020.7738082846451850.4523834307096310.226191715354815
2030.8146589649905270.3706820700189460.185341035009473
2040.6869437448521290.6261125102957420.313056255147871


Meta Analysis of Goldfeld-Quandt test for Heteroskedasticity
Description# significant tests% significant testsOK/NOK
1% type I error level1870.949238578680203NOK
5% type I error level1890.959390862944162NOK
10% type I error level1900.964467005076142NOK
 
Charts produced by software:
http://www.freestatistics.org/blog/date/2010/Dec/08/t1291835543aetyp8l44oscs7s/10nmu71291835596.png (open in new window)
http://www.freestatistics.org/blog/date/2010/Dec/08/t1291835543aetyp8l44oscs7s/10nmu71291835596.ps (open in new window)


http://www.freestatistics.org/blog/date/2010/Dec/08/t1291835543aetyp8l44oscs7s/1ylxd1291835596.png (open in new window)
http://www.freestatistics.org/blog/date/2010/Dec/08/t1291835543aetyp8l44oscs7s/1ylxd1291835596.ps (open in new window)


http://www.freestatistics.org/blog/date/2010/Dec/08/t1291835543aetyp8l44oscs7s/28cey1291835596.png (open in new window)
http://www.freestatistics.org/blog/date/2010/Dec/08/t1291835543aetyp8l44oscs7s/28cey1291835596.ps (open in new window)


http://www.freestatistics.org/blog/date/2010/Dec/08/t1291835543aetyp8l44oscs7s/38cey1291835596.png (open in new window)
http://www.freestatistics.org/blog/date/2010/Dec/08/t1291835543aetyp8l44oscs7s/38cey1291835596.ps (open in new window)


http://www.freestatistics.org/blog/date/2010/Dec/08/t1291835543aetyp8l44oscs7s/48cey1291835596.png (open in new window)
http://www.freestatistics.org/blog/date/2010/Dec/08/t1291835543aetyp8l44oscs7s/48cey1291835596.ps (open in new window)


http://www.freestatistics.org/blog/date/2010/Dec/08/t1291835543aetyp8l44oscs7s/51ldj1291835596.png (open in new window)
http://www.freestatistics.org/blog/date/2010/Dec/08/t1291835543aetyp8l44oscs7s/51ldj1291835596.ps (open in new window)


http://www.freestatistics.org/blog/date/2010/Dec/08/t1291835543aetyp8l44oscs7s/61ldj1291835596.png (open in new window)
http://www.freestatistics.org/blog/date/2010/Dec/08/t1291835543aetyp8l44oscs7s/61ldj1291835596.ps (open in new window)


http://www.freestatistics.org/blog/date/2010/Dec/08/t1291835543aetyp8l44oscs7s/7cvum1291835596.png (open in new window)
http://www.freestatistics.org/blog/date/2010/Dec/08/t1291835543aetyp8l44oscs7s/7cvum1291835596.ps (open in new window)


http://www.freestatistics.org/blog/date/2010/Dec/08/t1291835543aetyp8l44oscs7s/8cvum1291835596.png (open in new window)
http://www.freestatistics.org/blog/date/2010/Dec/08/t1291835543aetyp8l44oscs7s/8cvum1291835596.ps (open in new window)


http://www.freestatistics.org/blog/date/2010/Dec/08/t1291835543aetyp8l44oscs7s/9nmu71291835596.png (open in new window)
http://www.freestatistics.org/blog/date/2010/Dec/08/t1291835543aetyp8l44oscs7s/9nmu71291835596.ps (open in new window)


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