Home » date » 2009 » Dec » 29 »

Multiple Regression 6

*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:48:04 -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/t12620873577p07r06hcxgqmpr.htm/, Retrieved Tue, 29 Dec 2009 12:49:30 +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/t12620873577p07r06hcxgqmpr.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:
 
Dataseries X:
» Textbox « » Textfile « » CSV «
10.8 4.3 11.1 11.2 11.4 11.7 10.4 4.1 10.8 11.1 11.2 11.4 10.1 3.9 10.4 10.8 11.1 11.2 9.8 3.7 10.1 10.4 10.8 11.1 9.7 3.6 9.8 10.1 10.4 10.8 10.3 3.9 9.7 9.8 10.1 10.4 10.9 4.2 10.3 9.7 9.8 10.1 10.8 4.2 10.9 10.3 9.7 9.8 10.6 4.1 10.8 10.9 10.3 9.7 10.4 4.1 10.6 10.8 10.9 10.3 10.3 4.1 10.4 10.6 10.8 10.9 10.2 4.1 10.3 10.4 10.6 10.8 10 4.1 10.2 10.3 10.4 10.6 9.7 4 10 10.2 10.3 10.4 9.4 3.9 9.7 10 10.2 10.3 9.2 3.8 9.4 9.7 10 10.2 9.1 3.8 9.2 9.4 9.7 10 9.6 4 9.1 9.2 9.4 9.7 10.2 4.4 9.6 9.1 9.2 9.4 10.2 4.6 10.2 9.6 9.1 9.2 10 4.6 10.2 10.2 9.6 9.1 9.9 4.6 10 10.2 10.2 9.6 9.9 4.7 9.9 10 10.2 10.2 9.9 4.8 9.9 9.9 10 10.2 9.7 4.8 9.9 9.9 9.9 10 9.5 4.7 9.7 9.9 9.9 9.9 9.4 4.7 9.5 9.7 9.9 9.9 9.3 4.7 9.4 9.5 9.7 9.9 9.3 4.6 9.3 9.4 9.5 9.7 9.9 5 9.3 9.3 9.4 9.5 10.5 5.4 9.9 9.3 9.3 9.4 10.6 5.5 10.5 9.9 9.3 9.3 10.6 5.6 10.6 10.5 9.9 9.3 10.5 5.6 10.6 10.6 10.5 9.9 10.6 5.8 10.5 10.6 10.6 10.5 10.8 6 10.6 10.5 10.6 10.6 10.8 6.1 10.8 10.6 10.5 10.6 10.7 6.1 10.8 10.8 10. 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.406327472245542 + 0.132002787047602X[t] + 1.33719141843655Y1[t] -0.523857828846532Y2[t] -0.238580032736566Y3[t] + 0.322342059837373Y4[t] -0.170474709866218M1[t] -0.123311871032308M2[t] -0.0849884362519037M3[t] -0.0917642589554215M4[t] + 0.0375192409404774M5[t] + 0.662186968591818M6[t] + 0.136244777286711M7[t] -0.0875822905375944M8[t] + 0.171234680238685M9[t] + 0.0642899482917782M10[t] + 0.160873466453772M11[t] -0.00279915567983406t + e[t]


Multiple Linear Regression - Ordinary Least Squares
VariableParameterS.D.T-STAT
H0: parameter = 0
2-tail p-value1-tail p-value
(Intercept)0.4063274722455420.1393592.91570.003920.00196
X0.1320027870476020.0221465.960600
Y11.337191418436550.06447620.739300
Y2-0.5238578288465320.110877-4.72474e-062e-06
Y3-0.2385800327365660.109924-2.17040.0310570.015528
Y40.3223420598373730.0593455.431700
M1-0.1704747098662180.063764-2.67350.0080740.004037
M2-0.1233118710323080.065029-1.89630.0592490.029625
M3-0.08498843625190370.064575-1.31610.1895150.094758
M4-0.09176425895542150.063158-1.45290.147680.07384
M50.03751924094047740.0627310.59810.5503960.275198
M60.6621869685918180.06167910.73600
M70.1362447772867110.0717411.89910.058870.029435
M8-0.08758229053759440.081154-1.07920.2816860.140843
M90.1712346802386850.0813542.10480.0364540.018227
M100.06428994829177820.0717960.89550.3715340.185767
M110.1608734664537720.0641662.50710.0129010.00645
t-0.002799155679834060.000478-5.852700


Multiple Linear Regression - Regression Statistics
Multiple R0.994382149367365
R-squared0.98879585898046
Adjusted R-squared0.987922141561505
F-TEST (value)1131.71128047690
F-TEST (DF numerator)17
F-TEST (DF denominator)218
p-value0
Multiple Linear Regression - Residual Statistics
Residual Standard Deviation0.189540022317192
Sum Squared Residuals7.83174157308037


Multiple Linear Regression - Actuals, Interpolation, and Residuals
Time or IndexActualsInterpolation
Forecast
Residuals
Prediction Error
110.810.8278723794691-0.0278723794690997
210.410.4480772511634-0.0480772511634411
310.110.038871345440.0611286545600049
49.89.85062131949202-0.0506213194920155
59.79.71863470326973-0.0186347032697343
610.310.3461795040518-0.0461795040518444
710.910.68661101899750.213388981002476
810.810.8751403345698-0.0751403345698372
910.610.49454180618430.105458193815727
1010.410.22000263401540.179997365984642
1110.310.3683835177556-0.0683835177555948
1210.210.19124512011120.00875487988878097
13109.9198854901860.0801145098139981
149.79.695385985138850.00461401486115476
159.49.41294792306292-0.0129479230629162
169.29.161654389661380.0383456103386212
179.19.18496339669759-0.084963396697587
189.69.77915634187403-0.179156341874026
1910.29.975210990407150.224789009592847
2010.210.2747608522574-0.074760852257376
211010.0649397476939-0.0649397476938813
229.99.705780586656580.194219413343422
239.99.97722288767157-0.0772228876715738
249.99.9268523336747-0.0268523336746943
259.79.71296805943483-0.0129680594348255
269.59.44445897421310.0555410257869075
279.49.317316535395660.0826834646043384
289.39.32650998748528-0.0265099874852751
299.39.34170828861742-0.0417082886174169
309.910.0281533495988-0.1281533495988
3110.510.34615176578470.153848234215254
3210.610.58849176875560.0115082312443625
3310.610.53396628745060.0660337125493615
3410.510.42209383319970.077906166800272
3510.610.57810684387650.0218931561234788
3610.810.65917390986450.140826090135521
3710.810.73801082709950.0619891729004987
3810.710.62151073522690.0784892647731225
3910.610.49463379321550.105366206784546
4010.610.46819386784060.131806132159422
4110.810.68412227691970.115877723080288
4211.411.6178540446875-0.217854044687489
4312.211.833623449240.366376550760007
4412.412.35432049274410.0456795072558651
4512.412.38001072077620.0199892792237989
4612.312.18123675579810.118763244201916
4712.412.36465989646390.0353401035360656
4812.512.47796116843560.0220388315643688
4912.512.40987866512220.0901213348777636
5012.412.33256407742950.0674359225704956
5112.312.22954513869170.070454861308259
5212.212.15767072862840.0423292713716335
5312.112.2002791597496-0.100279159749564
5412.612.7720390061663-0.172039006166271
5513.213.00870406339320.191295936606783
5613.413.31408757381770.0859124261823443
5713.213.3585044742367-0.158504474236713
5812.912.88137346872530.0186265312746604
5912.812.8112609220962-0.011260922096193
6012.712.7964112039924-0.0964112039924424
6112.612.53570929863610.0642907013638791
6212.412.39949445074410.000505549255879314
6312.112.1851894689224-0.0851894689224362
641211.84445187065780.155548129342176
6511.911.996655943543-0.09665594354301
6612.512.5838975905233-0.0838975905232926
6713.212.90301365633120.296986343668823
6813.413.30293080441940.0970691955806218
6913.313.28430419738490.0156958026151152
701312.94926853642730.0507314635727493
7112.912.87220469160190.0277953083980552
721312.83349697022450.166503029775463
7312.912.87246755453930.0275324454607355
7412.612.6314811408780-0.0314811408779579
7512.412.23574101066530.164258989334698
7612.112.1455767490965-0.0455767490964766
7711.911.9754142012738-0.0754142012738348
7812.312.4644157842176-0.164415784217614
791312.63522928304910.364770716950865
801313.0993075882128-0.099307588212824
8112.612.8023239406451-0.202323940645077
8212.212.10643400795830.0935659920416899
8312.112.08732409778590.0126759022141359
841212.1081077571466-0.108107757146604
8511.811.8067954430965-0.00679544309646715
8611.611.51782752608180.0821724739181641
8711.411.36910860584960.0308913941504401
8811.211.19914843140700.000851568592975974
8911.211.13301337358020.0669866264198363
9011.811.8957022207199-0.095702220719855
9112.512.20532443419570.294675565804276
9212.612.54914937302650.05085062697347
9312.412.4158375514274-0.0158375514273565
9412.111.99946853151070.100531468489277
951211.99864847284370.00135152715627272
961211.95156454875620.0484354512437615
9711.911.83778206394830.0622179360516651
9811.811.67558199058120.124418009418799
9911.511.5843384260343-0.0843384260342731
10011.311.22344925086870.0765507491312545
10111.211.2180761786366-0.0180761786366228
10211.611.8767375357805-0.276737535780532
10312.211.96547359987950.234526400120474
10412.212.3174092446144-0.117409244614403
10511.712.1050455859150-0.405045585915045
10611.211.2728939572488-0.0728939572487603
1071111.1534167608383-0.153416760838336
10810.911.1167250645137-0.216725064513731
10910.810.8726226093429-0.0726226093429273
11010.510.7089976314619-0.208997631461869
11110.210.3287393018128-0.128739301812791
1121010.0403874864328-0.0403874864328324
1139.910.0695301420433-0.169530142043262
11410.310.6373225298102-0.337322529810178
11510.710.66005720038540.0399427996146291
11610.410.6653528892044-0.265352889204375
11710.110.1566033707434-0.0566033707433625
1189.79.80996365967042-0.109963659670419
1199.49.72653963718784-0.326539637187837
1208.99.3593243916364-0.459324391636401
1218.48.64694100325993-0.246941003259930
1228.18.17447396268598-0.0744739626859784
1238.38.040556014276880.259443985723120
1248.18.37409481857012-0.274094818570116
12587.999171457117570.00082854288243362
1268.78.513675222040.186324777959997
1279.29.138339184179120.0616608158208758
12899.15979750230211-0.159797502302106
1298.98.674007611683880.225992388316117
1308.58.62846529479591-0.128465294795912
1318.18.43544563054935-0.335445630549346
1327.57.90582916388592-0.405829163885916
1337.17.20298138592744-0.10298138592744
1346.97.00647866687925-0.106478666879253
1357.17.03791982565240.0620801743475994
13677.2893811952131-0.289381195213105
1376.77.0565531783143-0.356553178314293
13877.24386624653422-0.243866246534223
1397.37.38816664638486-0.0881666463848558
1407.77.4976814184140.202318581586007
1418.48.042343496687480.357656503312524
1428.48.67101879985319-0.271018799853187
1438.88.399373286999360.400626713000642
1449.18.692907197145440.407092802854556
14598.923686788773140.0763132112268579
1468.68.581741968334980.0182580316650219
1477.98.20533855576432-0.305338555764322
1487.77.576633058534110.123366941465890
1497.87.82597657025205-0.0259765702520479
1509.28.764005884931440.435994115068557
1519.49.87702330553416-0.477023305534163
1529.29.135708826205340.0642911737946601
1538.78.83654146034055-0.136541460340554
1548.48.60613714260414-0.206137142604136
1558.68.68607769119815-0.0860776911981478
15698.988622027101860.0113779728981428
1579.19.12945558565388-0.0294555856538811
1588.78.94037637590971-0.240376375909713
1598.28.37087498232862-0.170874982328620
1607.97.99412596822214-0.0941259682221396
1617.98.08264746299935-0.182647462999352
1629.18.891627412172430.208372587827567
1639.49.86471846850887-0.464718468508868
1649.49.32711793667340.0728820633265968
1659.19.17928319994629-0.0792831999462903
16698.983618348772460.0163816512275383
1679.39.223944093425660.076055906574341
1689.99.59858896823340.301411031766599
1699.89.97122743300824-0.171227433008240
1709.39.39774766768223-0.0977476676822338
1718.38.67821466782514-0.378214667825135
17287.77103958849030.228960411509695
1738.58.120480425111240.379519574888757
17410.49.777513844820470.622486155179535
17511.111.3295403432443-0.229540343244333
17610.910.84082588222260.0591741177773727
1771010.1177727863394-0.117772786339362
1789.29.32837152125493-0.128371521254929
1799.29.110430522087960.0895694779120363
1809.59.54249833793654-0.0424983379365393
1819.69.67113807025707-0.0711380702570696
1829.59.434189898730940.0658101012690597
1839.19.17243440716795-0.0724344071679533
1848.98.75321325897220.146786741027809
18598.825093545477910.174906454522086
18610.19.801451746992310.298548253007688
18710.310.6100143600153-0.310014360015279
18810.29.99945667193090.200543328069104
1899.69.79998022809265-0.199980228092652
1909.29.24716753156243-0.0471675315624352
1919.39.235116996628550.0648830033714542
1929.49.55202101894493-0.152021018944933
1939.49.375307568254840.0246924317451564
1949.29.21449064131566-0.0144906413156608
19598.9909528394390.0090471605609987
19698.81134451300710.188655486992895
19799.01111475731123-0.0111147573112289
1989.89.563429809043530.236570190956468
1991010.0003723487261-0.000372348726073684
2009.89.635298424536780.164701575463224
2019.39.35464292139688-0.0546429213968742
20298.904432810348510.0955671896514872
20398.984373358942520.0156266410574788
2049.19.045879968568440.0541200314315553
2059.18.903527946063570.196472053936430
2069.18.772402670972260.32759732902774
2079.28.797269225503930.402730774496065
2088.88.92724703753845-0.127247037538453
2098.38.4136679166762-0.113667916676207
2108.48.59222674380868-0.192226743808678
2118.18.57359939346426-0.473599393464261
2127.77.85738259656832-0.157382596568317
2137.97.5374518810470.362548118952998
2147.97.99529734574613-0.0952973457461306
21588.02264037371668-0.0226403737166777
2167.97.855634899058740.044365100941256
2177.67.560724520751860.03927547924814
2187.17.17965744316693-0.0796574431669274
2196.86.72023473484630.0797652651537024
2206.56.5711702130782-0.0711702130782014
2216.96.515842714948610.384157285051386
2228.27.885351248603340.314648751396655
2238.78.87348728462182-0.173487284621817
2248.38.38950584697058-0.0895058469705781
2257.97.561898732008470.338101267991527
2267.57.386975233851750.113024766148253
2277.87.464830318330250.335169681669745
2288.37.997155950768740.302844049231257
2298.48.381017307175250.0189826928247552
2308.28.123060941403310.0769390585966898
2317.77.78977319810532-0.0897731981053238
2327.27.31408626680376-0.114086266803758
2337.37.127054307460630.172945692539374
2348.18.26539393362366-0.165393933623663
2358.58.72533920365766-0.225339203657662
2368.48.41627397255381-0.0162739725538122


Goldfeld-Quandt test for Heteroskedasticity
p-valuesAlternative Hypothesis
breakpoint indexgreater2-sidedless
210.004680342361808400.009360684723616810.995319657638192
220.000644688282222940.001289376564445880.999355311717777
237.80028630655313e-050.0001560057261310630.999921997136934
243.12667336796914e-056.25334673593828e-050.99996873326632
255.28354791017982e-061.05670958203596e-050.99999471645209
263.23463165899134e-056.46926331798268e-050.99996765368341
277.35783727773474e-050.0001471567455546950.999926421627223
282.22316962273037e-054.44633924546074e-050.999977768303773
295.34870308814221e-061.06974061762844e-050.999994651296912
301.95018691783463e-063.90037383566926e-060.999998049813082
315.09035872959235e-071.01807174591847e-060.999999490964127
329.74467451521214e-071.94893490304243e-060.999999025532548
334.77537701670907e-079.55075403341814e-070.999999522462298
341.33602445881224e-072.67204891762448e-070.999999866397554
355.879860624564e-081.1759721249128e-070.999999941201394
362.34444215390459e-084.68888430780917e-080.999999976555578
379.80478249458747e-091.96095649891749e-080.999999990195217
382.50067321666301e-095.00134643332602e-090.999999997499327
399.516947918251e-101.9033895836502e-090.999999999048305
402.76578087878964e-105.53156175757928e-100.999999999723422
418.95961758708269e-111.79192351741654e-100.999999999910404
423.37850306033759e-106.75700612067517e-100.99999999966215
431.62051236684262e-103.24102473368524e-100.999999999837949
447.68050663952966e-111.53610132790593e-100.999999999923195
452.02520476066788e-114.05040952133576e-110.999999999979748
464.03074415966445e-118.0614883193289e-110.999999999959693
471.20299304160804e-112.40598608321607e-110.99999999998797
484.73926986160861e-129.47853972321722e-120.99999999999526
492.63896246742893e-125.27792493485786e-120.99999999999736
507.43467017208196e-131.48693403441639e-120.999999999999257
512.03010102522181e-134.06020205044361e-130.999999999999797
526.32402253814971e-141.26480450762994e-130.999999999999937
536.27601636470225e-141.25520327294045e-130.999999999999937
541.82564010080091e-143.65128020160181e-140.999999999999982
555.88838389805992e-151.17767677961198e-140.999999999999994
562.24262539714358e-134.48525079428717e-130.999999999999776
571.13121135550702e-132.26242271101404e-130.999999999999887
584.00785131564075e-148.0157026312815e-140.99999999999996
591.16118696069211e-142.32237392138423e-140.999999999999988
603.60823411182666e-157.21646822365331e-150.999999999999996
618.34607142398277e-151.66921428479655e-140.999999999999992
624.6134985171474e-159.2269970342948e-150.999999999999995
631.39553470614415e-152.79106941228829e-150.999999999999999
644.42703071681576e-158.85406143363153e-150.999999999999996
651.39468496765448e-152.78936993530896e-150.999999999999999
662.31695508402027e-154.63391016804054e-150.999999999999998
672.62992088491929e-155.25984176983859e-150.999999999999997
685.72248805620978e-151.14449761124196e-140.999999999999994
691.94694304780188e-153.89388609560375e-150.999999999999998
701.17254480296331e-152.34508960592662e-150.999999999999999
713.81562942914398e-167.63125885828796e-161
726.20103929615346e-161.24020785923069e-151
732.31483605454497e-164.62967210908994e-161
747.93467664357183e-171.58693532871437e-161
754.3148360285712e-178.6296720571424e-171
764.00835313012862e-178.01670626025725e-171
771.42922128516097e-172.85844257032194e-171
786.51434631473837e-181.30286926294767e-171
793.95289723796886e-177.90579447593773e-171
801.60559614024445e-173.21119228048891e-171
813.16519133726147e-176.33038267452293e-171
821.95821816940309e-173.91643633880618e-171
836.87077562999184e-181.37415512599837e-171
842.78574180786551e-185.57148361573103e-181
859.58405649259454e-191.91681129851891e-181
865.1577359593891e-191.03154719187782e-181
871.83941362200003e-193.67882724400006e-191
886.07113389630245e-201.21422677926049e-191
893.78123423572719e-207.56246847145438e-201
901.46292169254294e-202.92584338508588e-201
912.74696444671242e-205.49392889342483e-201
921.02167098360439e-202.04334196720879e-201
933.32433310226576e-216.64866620453152e-211
941.82740064493995e-213.65480128987991e-211
956.39731326913074e-221.27946265382615e-211
962.45309643334812e-224.90619286669624e-221
971.09627903086292e-222.19255806172583e-221
981.23297055789773e-222.46594111579546e-221
991.75057998714631e-223.50115997429262e-221
1007.92095415288039e-231.58419083057608e-221
1013.83043134266089e-237.66086268532179e-231
1026.35053586570662e-231.27010717314132e-221
1032.90529155003505e-225.81058310007009e-221
1043.45296172845869e-226.90592345691738e-221
1053.15202202478562e-196.30404404957123e-191
1065.25006508869064e-191.05001301773813e-181
1072.83288246512103e-195.66576493024206e-191
1081.14439400248882e-192.28878800497763e-191
1095.96939808867683e-201.19387961773537e-191
1105.61985725528464e-201.12397145105693e-191
1112.43067273881754e-204.86134547763509e-201
1121.02563571675040e-202.05127143350081e-201
1134.09523729579408e-218.19047459158816e-211
1146.94171562343563e-211.38834312468713e-201
1153.91074920343275e-207.8214984068655e-201
1161.17343671730961e-192.34687343461922e-191
1171.50723420564293e-193.01446841128586e-191
1186.06494033762475e-201.21298806752495e-191
1195.63454425498692e-201.12690885099738e-191
1201.17473502404606e-172.34947004809211e-171
1216.75285814805208e-181.35057162961042e-171
1227.46345832432192e-181.49269166486438e-171
1231.41313030296940e-132.82626060593879e-130.999999999999859
1242.4130502456517e-134.8261004913034e-130.999999999999759
1253.53057690551168e-137.06115381102337e-130.999999999999647
1268.30402208141236e-131.66080441628247e-120.99999999999917
1273.10963475797579e-126.21926951595157e-120.99999999999689
1281.98196772997160e-123.96393545994319e-120.999999999998018
1295.8548169918591e-121.17096339837182e-110.999999999994145
1304.46372376682009e-128.92744753364018e-120.999999999995536
1311.13780421022587e-112.27560842045174e-110.999999999988622
1325.08419553604622e-101.01683910720924e-090.99999999949158
1332.99998194910352e-105.99996389820704e-100.999999999700002
1342.17428583288045e-104.34857166576091e-100.999999999782571
1351.03843478318746e-092.07686956637491e-090.999999998961565
1362.69841022232531e-095.39682044465063e-090.99999999730159
1372.62756562656667e-085.25513125313335e-080.999999973724344
1381.07226764141233e-072.14453528282466e-070.999999892773236
1392.08393799210078e-074.16787598420157e-070.9999997916062
1402.74734474546320e-065.49468949092641e-060.999997252655255
1413.08767140877996e-056.17534281755992e-050.999969123285912
1425.5325392515099e-050.0001106507850301980.999944674607485
1430.0006452756022389860.001290551204477970.99935472439776
1440.002195089764416820.004390179528833640.997804910235583
1450.001715540004899420.003431080009798840.9982844599951
1460.001382721269420450.002765442538840890.99861727873058
1470.004114213435776550.00822842687155310.995885786564223
1480.004393072342684540.008786144685369090.995606927657315
1490.003674521268302760.007349042536605530.996325478731697
1500.04502766988750160.09005533977500310.954972330112498
1510.1232690764974230.2465381529948460.876730923502577
1520.1063504840465300.2127009680930600.89364951595347
1530.1037802151487000.2075604302974000.8962197848513
1540.0944373248734860.1888746497469720.905562675126514
1550.07987104099625370.1597420819925070.920128959003746
1560.06900660003879750.1380132000775950.930993399961202
1570.05591450690940440.1118290138188090.944085493090596
1580.06995619287983280.1399123857596660.930043807120167
1590.06725788354588070.1345157670917610.932742116454119
1600.05565831544788730.1113166308957750.944341684552113
1610.07325131237020970.1465026247404190.92674868762979
1620.08374885734922130.1674977146984430.916251142650779
1630.1463790483665160.2927580967330330.853620951633484
1640.1247174502756050.2494349005512100.875282549724395
1650.1130514719416190.2261029438832390.88694852805838
1660.09418012972464460.1883602594492890.905819870275355
1670.08208463923816560.1641692784763310.917915360761834
1680.1128448825430540.2256897650861070.887155117456946
1690.1115144526139410.2230289052278820.88848554738606
1700.1071218763419370.2142437526838740.892878123658063
1710.2158881815600010.4317763631200030.784111818439999
1720.2249708191214770.4499416382429530.775029180878523
1730.3076803022070680.6153606044141360.692319697792932
1740.5971873461992850.8056253076014310.402812653800715
1750.5681365150637450.863726969872510.431863484936255
1760.5250493170817210.9499013658365580.474950682918279
1770.4970105331100730.9940210662201460.502989466889927
1780.4496152792116250.899230558423250.550384720788375
1790.4209167221488950.841833444297790.579083277851105
1800.3760603300757650.752120660151530.623939669924235
1810.3564212625701050.7128425251402090.643578737429895
1820.3175582210959590.6351164421919180.682441778904041
1830.3082531798911120.6165063597822250.691746820108888
1840.3031151896579910.6062303793159820.696884810342009
1850.2807919847780260.5615839695560510.719208015221974
1860.3137839793615110.6275679587230220.686216020638489
1870.2948032891620830.5896065783241650.705196710837917
1880.3237399249500540.6474798499001080.676260075049946
1890.3969696932941670.7939393865883350.603030306705833
1900.3517574453247440.7035148906494880.648242554675256
1910.3050874439069930.6101748878139850.694912556093007
1920.3003467069718580.6006934139437160.699653293028142
1930.2517893174007420.5035786348014840.748210682599258
1940.2154076141393220.4308152282786430.784592385860678
1950.182680445917470.365360891834940.81731955408253
1960.2453345226661380.4906690453322750.754665477333862
1970.2302654402608140.4605308805216280.769734559739186
1980.2698356018066280.5396712036132570.730164398193372
1990.2639598838994080.5279197677988160.736040116100592
2000.3089639623977480.6179279247954960.691036037602252
2010.4065256546796150.813051309359230.593474345320385
2020.5460892926008050.907821414798390.453910707399195
2030.5362595045356980.9274809909286030.463740495464302
2040.4781029109394070.9562058218788140.521897089060593
2050.4487330655951610.8974661311903220.551266934404839
2060.5508837378045910.8982325243908180.449116262195409
2070.8613797764230970.2772404471538070.138620223576903
2080.8756932270722740.2486135458554530.124306772927726
2090.8160207240131940.3679585519736120.183979275986806
2100.8009851812882440.3980296374235130.199014818711756
2110.7425182997683060.5149634004633870.257481700231694
2120.7780483015004640.4439033969990720.221951698499536
2130.8386533371692670.3226933256614650.161346662830733
2140.8256496207162440.3487007585675120.174350379283756
2150.6805638665940930.6388722668118140.319436133405907


Meta Analysis of Goldfeld-Quandt test for Heteroskedasticity
Description# significant tests% significant testsOK/NOK
1% type I error level1290.661538461538462NOK
5% type I error level1290.661538461538462NOK
10% type I error level1300.666666666666667NOK
 
Charts produced by software:
http://www.freestatistics.org/blog/date/2009/Dec/29/t12620873577p07r06hcxgqmpr/1044vj1262087273.png (open in new window)
http://www.freestatistics.org/blog/date/2009/Dec/29/t12620873577p07r06hcxgqmpr/1044vj1262087273.ps (open in new window)


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


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


http://www.freestatistics.org/blog/date/2009/Dec/29/t12620873577p07r06hcxgqmpr/30tck1262087273.png (open in new window)
http://www.freestatistics.org/blog/date/2009/Dec/29/t12620873577p07r06hcxgqmpr/30tck1262087273.ps (open in new window)


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


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


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


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


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


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