Home » date » 2008 » Dec » 05 »

*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: Fri, 05 Dec 2008 02:59:14 -0700
 
Cite this page as follows:
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL http://www.freestatistics.org/blog/date/2008/Dec/05/t1228471629zz53asfw3z5xpue.htm/, Retrieved Fri, 05 Dec 2008 10:07:09 +0000
 
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/2008/Dec/05/t1228471629zz53asfw3z5xpue.htm/},
    year = {2008},
}
@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 = {2008},
    note = {{ISBN} 3-900051-07-0},
    url = {http://www.R-project.org},
}
 
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
 
Feedback Forum:
2008-11-27 13:41:43 [a2386b643d711541400692649981f2dc] [reply
test

Post a new message
 
Original text written by user:
 
IsPrivate?
No (this computation is public)
 
User-defined keywords:
 
Dataseries X:
» Textbox « » Textfile « » CSV «
0,24 0,23 0,23 0,24 0,23 0,23 0,25 0,21 0,26 0,25 0,24 0,24 0,27 0,25 0,26 0,29 0,24 0,26 0,24 0,26 0,25 0,26 0,24 0,21 0,20 0,22 0,20 0,21 0,20 0,19 0,20 0,20 0,21 0,24 0,22 0,19 0,23 0,23 0,23 0,22 0,23 0,25 0,25 0,22 0,25 0,25 0,24 0,19 0,24 0,26 0,24 0,24 0,25 0,23 0,27 0,24 0,26 0,27 0,29 0,28 0,32 0,29 0,27 0,26 0,28 0,31 0,29 0,31 0,31 0,32 0,32 0,26 0,31 0,31 0,31 0,31 0,29 0,27 0,30 0,27 0,27 0,30 0,28 0,24 0,28 0,28 0,33 0,28 0,29 0,25 0,31 0,29 0,37 0,31 0,29 0,28 0,30 0,32 0,31 0,28 0,29 0,29 0,28 0,26 0,28 0,30 0,33 0,31 0,37 0,36 0,37 0,37 0,36 0,33 0,33 0,40 0,32 0,39 0,39 0,37 0,37 0,30 0,33 0,33 0,34 0,35 0,34 0,37 0,37 0,37 0,36 0,32 0,33 0,35 0,36 0,35 0,37 0,35 0,32 0,33 0,28 0,32 0,35 0,30 0,32 0,32 0,32 0,32 0,36 0,31 0,26 0,33 0,31 0,34 0,33 0,38 0,32 0,30 0,32 0,33 0,34 0,29 0,33 0,36 0,32 0,32 0,32 0,31 0,30 0,34 0,34 0,30 0,28 0,25 0,27 0,33 0,28 0,33 0,32 0,27 0,27 0,28 0,27 0,27 0,25 0,25 0,22 0,27
 
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'Sir Ronald Aylmer Fisher' @ 193.190.124.24


Multiple Linear Regression - Estimated Regression Equation
X[t] = + 0.232134597594820 + 0.00612345339962998Q1[t] + 0.00578442992599445Q2[t] + 0.00629647028214616Q3[t] + 0.000551789431082331t + e[t]


Multiple Linear Regression - Ordinary Least Squares
VariableParameterS.D.T-STAT
H0: parameter = 0
2-tail p-value1-tail p-value
(Intercept)0.2321345975948200.00744331.187200
Q10.006123453399629980.0078860.77650.4384340.219217
Q20.005784429925994450.0078850.73360.4641170.232059
Q30.006296470282146160.0078840.79860.4255490.212775
t0.0005517894310823315.1e-0510.740800


Multiple Linear Regression - Regression Statistics
Multiple R0.622847105294973
R-squared0.387938516574328
Adjusted R-squared0.374560123493985
F-TEST (value)28.9973926049715
F-TEST (DF numerator)4
F-TEST (DF denominator)183
p-value0
Multiple Linear Regression - Residual Statistics
Residual Standard Deviation0.0382194765515835
Sum Squared Residuals0.267313294981499


Multiple Linear Regression - Actuals, Interpolation, and Residuals
Time or IndexActualsInterpolation
Forecast
Residuals
Prediction Error
10.240.2388098404255310.0011901595744685
20.230.239022606382979-0.00902260638297874
30.230.240086436170213-0.0100864361702128
40.240.2343417553191490.00565824468085105
50.230.241016998149861-0.0110169981498612
60.230.241229764107308-0.0112297641073080
70.250.2422935938945420.0077064061054579
80.210.236548913043478-0.0265489130434783
90.260.2432241558741910.0167758441258094
100.250.2434369218316370.00656307816836263
110.240.244500751618871-0.00450075161887143
120.240.2387560707678080.0012439292321924
130.270.245431313598520.0245686864014801
140.250.2456440795559670.0043559204440333
150.260.2467079093432010.0132920906567993
160.290.2409632284921370.0490367715078631
170.240.247638471322849-0.00763847132284924
180.260.2478512372802960.0121487627197040
190.240.24891506706753-0.00891506706753008
200.260.2431703862164660.0168296137835338
210.250.2498456290471790.000154370952821448
220.260.2500583950046250.00994160499537466
230.240.251122224791859-0.0111222247918594
240.210.245377543940796-0.0353775439407956
250.20.252052786771508-0.0520527867715079
260.220.252265552728955-0.0322655527289547
270.20.253329382516189-0.0533293825161887
280.210.247584701665125-0.0375847016651249
290.20.254259944495837-0.0542599444958372
300.190.254472710453284-0.064472710453284
310.20.255536540240518-0.055536540240518
320.20.249791859389454-0.0497918593894542
330.210.256467102220167-0.0464671022201665
340.240.256679868177613-0.0166798681776133
350.220.257743697964847-0.0377436979648474
360.190.251999017113784-0.0619990171137835
370.230.258674259944496-0.0286742599444958
380.230.258887025901943-0.0288870259019426
390.230.259950855689177-0.0299508556891767
400.220.254206174838113-0.0342061748381129
410.230.260881417668825-0.0308814176688252
420.250.261094183626272-0.0110941836262720
430.250.262158013413506-0.0121580134135060
440.220.256413332562442-0.0364133325624422
450.250.263088575393155-0.0130885753931545
460.250.263301341350601-0.0133013413506013
470.240.264365171137835-0.0243651711378354
480.190.258620490286772-0.0686204902867715
490.240.265295733117484-0.0252957331174838
500.260.265508499074931-0.00550849907493061
510.240.266572328862165-0.0265723288621647
520.240.260827648011101-0.0208276480111008
530.250.267502890841813-0.0175028908418131
540.230.26771565679926-0.0377156567992599
550.270.2687794865864940.00122051341350603
560.240.26303480573543-0.0230348057354302
570.260.269710048566143-0.00971004856614247
580.270.2699228145235897.71854764107524e-05
590.290.2709866443108230.0190133556891767
600.280.2652419634597590.0147580365402405
610.320.2719172062904720.0480827937095282
620.290.2721299722479190.0178700277520814
630.270.273193802035153-0.00319380203515262
640.260.267449121184089-0.0074491211840888
650.280.2741243640148010.0058756359851989
660.310.2743371299722480.0356628700277521
670.290.2754009597594820.014599040240518
680.310.2696562789084180.0403437210915819
690.310.2763315217391300.0336684782608696
700.320.2765442876965770.0434557123034228
710.320.2776081174838110.0423918825161887
720.260.271863436632747-0.0118634366327475
730.310.278538679463460.0314613205365402
740.310.2787514454209070.0312485545790934
750.310.2798152752081410.0301847247918594
760.310.2740705943570770.0359294056429232
770.290.2807458371877890.00925416281221089
780.270.280958603145236-0.0109586031452359
790.30.282022432932470.0179775670675300
800.270.276277752081406-0.00627775208140609
810.270.282952994912118-0.0129529949121184
820.30.2831657608695650.0168342391304348
830.280.284229590656799-0.00422959065679924
840.240.278484909805735-0.0384849098057354
850.280.285160152636448-0.00516015263644772
860.280.285372918593895-0.00537291859389452
870.330.2864367483811290.0435632516188714
880.280.280692067530065-0.000692067530064733
890.290.2873673103607770.00263268963922291
900.250.287580076318224-0.0375800763182239
910.310.2886439061054580.0213560938945421
920.290.2828992252543940.0071007747456059
930.370.2895744680851060.0804255319148936
940.310.2897872340425530.0202127659574468
950.290.290851063829787-0.00085106382978726
960.280.285106382978723-0.00510638297872338
970.30.2917816258094360.00821837419056427
980.320.2919943917668830.0280056082331175
990.310.2930582215541170.0169417784458834
1000.280.287313540703053-0.0073135407030527
1010.290.293988783533765-0.00398878353376506
1020.290.294201549491212-0.00420154949121186
1030.280.295265379278446-0.0152653792784459
1040.260.289520698427382-0.0295206984273821
1050.280.296195941258094-0.0161959412580943
1060.30.2964087072155410.00359129278445882
1070.330.2974725370027750.0325274629972248
1080.310.2917278561517110.0182721438482886
1090.370.2984030989824240.0715969010175763
1100.360.2986158649398700.0613841350601295
1110.370.2996796947271050.0703203052728955
1120.370.2939350138760410.0760649861239593
1130.360.3006102567067530.059389743293247
1140.330.30082302266420.0291769773358002
1150.330.3018868524514340.0281131475485662
1160.40.296142171600370.10385782839963
1170.320.3028174144310820.0171825855689177
1180.390.3030301803885290.0869698196114709
1190.390.3040940101757630.0859059898242368
1200.370.2983493293246990.0716506706753006
1210.370.3050245721554120.0649754278445883
1220.30.305237338112858-0.00523733811285847
1230.330.3063011679000930.0236988320999075
1240.330.3005564870490290.0294435129509713
1250.340.3072317298797410.0327682701202590
1260.350.3074444958371880.0425555041628122
1270.340.3085083256244220.0314916743755782
1280.370.3027636447733580.067236355226642
1290.370.3094388876040700.0605611123959297
1300.370.3096516535615170.0603483464384829
1310.360.3107154833487510.0492845166512488
1320.320.3049708024976870.0150291975023127
1330.330.3116460453284000.0183539546716004
1340.350.3118588112858460.0381411887141535
1350.360.3129226410730800.0470773589269195
1360.350.3071779602220170.0428220397779833
1370.370.3138532030527290.056146796947271
1380.350.3140659690101760.0359340309898242
1390.320.315129798797410.0048702012025902
1400.330.3093851179463460.0206148820536540
1410.280.316060360777058-0.0360603607770583
1420.320.3162731267345050.00372687326549492
1430.350.3173369565217390.0326630434782609
1440.30.311592275670675-0.0115922756706753
1450.320.3182675185013880.00173248149861239
1460.320.3184802844588340.00151971554116559
1470.320.3195441142460680.000455885753931549
1480.320.3137994333950050.00620056660499538
1490.360.3204746762257170.039525323774283
1500.310.320687442183164-0.0106874421831637
1510.260.321751271970398-0.0617512719703978
1520.330.3160065911193340.0139934088806661
1530.310.322681833950046-0.0126818339500463
1540.340.3228945999074930.0171054000925070
1550.330.3239584296947270.00604157030527291
1560.380.3182137488436630.0617862511563367
1570.320.324888991674376-0.00488899167437558
1580.30.325101757631822-0.0251017576318224
1590.320.326165587419056-0.00616558741905642
1600.330.3204209065679930.00957909343200742
1610.340.3270961493987050.0129038506012951
1620.290.327308915356152-0.0373089153561517
1630.330.3283727451433860.00162725485661426
1640.360.3226280642923220.0373719357076781
1650.320.329303307123034-0.00930330712303423
1660.320.329516073080481-0.00951607308048103
1670.320.330579902867715-0.0105799028677151
1680.310.324835222016651-0.0148352220166512
1690.30.331510464847364-0.0315104648473636
1700.340.3317232308048100.00827676919518966
1710.340.3327870605920440.00721293940795562
1720.30.327042379740981-0.0270423797409806
1730.280.333717622571693-0.0537176225716929
1740.250.33393038852914-0.0839303885291397
1750.270.334994218316374-0.0649942183163737
1760.330.329249537465310.000750462534690122
1770.280.335924780296022-0.0559247802960222
1780.330.336137546253469-0.00613754625346899
1790.320.337201376040703-0.0172013760407031
1800.270.331456695189639-0.0614566951896392
1810.270.338131938020352-0.0681319380203515
1820.280.338344703977798-0.0583447039777983
1830.270.339408533765032-0.0694085337650324
1840.270.333663852913969-0.0636638529139685
1850.250.340339095744681-0.0903390957446809
1860.250.340551861702128-0.0905518617021277
1870.220.341615691489362-0.121615691489362
1880.270.335871010638298-0.0658710106382979


Goldfeld-Quandt test for Heteroskedasticity
p-valuesAlternative Hypothesis
breakpoint indexgreater2-sidedless
80.06566491058986850.1313298211797370.934335089410132
90.0476397665378910.0952795330757820.952360233462109
100.01964565931975350.03929131863950700.980354340680247
110.007418751337653830.01483750267530770.992581248662346
120.002508521254988880.005017042509977750.997491478745011
130.001038079543531980.002076159087063960.998961920456468
140.0003121525129825360.0006243050259650730.999687847487017
159.16155367664095e-050.0001832310735328190.999908384463234
160.0006237864238401140.001247572847680230.99937621357616
170.001148281985932870.002296563971865750.998851718014067
180.0004601670126732790.0009203340253465570.999539832987327
190.0003657339904231990.0007314679808463990.999634266009577
200.0001461616854307480.0002923233708614970.99985383831457
218.22895054797856e-050.0001645790109595710.99991771049452
223.14970949477591e-056.29941898955181e-050.999968502905052
232.08217765189514e-054.16435530379029e-050.999979178223481
240.0001522413042993320.0003044826085986650.9998477586957
250.0009054493830664920.001810898766132980.999094550616934
260.0007984889963929350.001596977992785870.999201511003607
270.001268820053670300.002537640107340610.99873117994633
280.001041741408676890.002083482817353780.998958258591323
290.001099265010514700.002198530021029390.998900734989485
300.001603936601994830.003207873203989660.998396063398005
310.001283096188352400.002566192376704810.998716903811648
320.0009753390353982150.001950678070796430.999024660964602
330.0006254148248920610.001250829649784120.999374585175108
340.0005136422231954160.001027284446390830.999486357776805
350.0003280564331328060.0006561128662656130.999671943566867
360.0003091371118543930.0006182742237087870.999690862888146
370.0002348069316080590.0004696138632161170.999765193068392
380.0001622373670532040.0003244747341064090.999837762632947
390.0001244548351044310.0002489096702088620.999875545164896
408.50402062743735e-050.0001700804125487470.999914959793726
416.36295203248859e-050.0001272590406497720.999936370479675
426.94634244616076e-050.0001389268489232150.999930536575538
438.41114372859404e-050.0001682228745718810.999915888562714
445.92415227028426e-050.0001184830454056850.999940758477297
456.41945934508712e-050.0001283891869017420.99993580540655
465.58600710650061e-050.0001117201421300120.999944139928935
474.43578005508643e-058.87156011017285e-050.99995564219945
486.84677792884433e-050.0001369355585768870.999931532220712
495.74493168326042e-050.0001148986336652080.999942550683167
506.39210211300437e-050.0001278420422600870.99993607897887
515.39418707050595e-050.0001078837414101190.999946058129295
525.46408656245064e-050.0001092817312490130.999945359134376
535.11149244281496e-050.0001022298488562990.999948885075572
544.45560433127285e-058.9112086625457e-050.999955443956687
557.39610014853121e-050.0001479220029706240.999926038998515
567.33795190630281e-050.0001467590381260560.999926620480937
577.77686510458294e-050.0001555373020916590.999922231348954
588.82753346579807e-050.0001765506693159610.999911724665342
590.0002055105751990330.0004110211503980660.9997944894248
600.0004062268799779280.0008124537599558560.999593773120022
610.002022214894553020.004044429789106050.997977785105447
620.002241764788889260.004483529577778530.99775823521111
630.001957844335514750.003915688671029510.998042155664485
640.001812074697219550.003624149394439100.99818792530278
650.00160594295491870.00321188590983740.998394057045081
660.002293235403680290.004586470807360590.99770676459632
670.002186946338268370.004373892676536740.997813053661732
680.003869519571951630.007739039143903260.996130480428048
690.004276396086032220.008552792172064440.995723603913968
700.005124999215916950.01024999843183390.994875000784083
710.006204586835674260.01240917367134850.993795413164326
720.005694980481226140.01138996096245230.994305019518774
730.005141283583698020.01028256716739600.994858716416302
740.004327787587320720.008655575174641440.99567221241268
750.003747749690249080.007495499380498170.99625225030975
760.003763398986450680.007526797972901360.99623660101355
770.002882163482361630.005764326964723260.997117836517638
780.00276642919736170.00553285839472340.997233570802638
790.002115771384271920.004231542768543840.997884228615728
800.001927401492895590.003854802985791190.998072598507104
810.001980081997249620.003960163994499250.99801991800275
820.001482052195320960.002964104390641920.99851794780468
830.001338906147801600.002677812295603190.998661093852198
840.002980546011488670.005961092022977350.997019453988511
850.002878191087681930.005756382175363860.997121808912318
860.002812885230437030.005625770460874050.997187114769563
870.002754226633736950.00550845326747390.997245773366263
880.002780113239905010.005560226479810030.997219886760095
890.002555773055712080.005111546111424170.997444226944288
900.006837905971899230.01367581194379850.9931620940281
910.005758407126412510.01151681425282500.994241592873587
920.005963133377664860.01192626675532970.994036866622335
930.01266612835131470.02533225670262940.987333871648685
940.01078436447752590.02156872895505180.989215635522474
950.01138990039928720.02277980079857440.988610099600713
960.01472601617919910.02945203235839830.9852739838208
970.01415730654443640.02831461308887280.985842693455564
980.01218334170821310.02436668341642620.987816658291787
990.01103838111104420.02207676222208840.988961618888956
1000.01712578312567550.0342515662513510.982874216874325
1010.02195804807766330.04391609615532660.978041951922337
1020.02946030596105260.05892061192210520.970539694038947
1030.05354499750174970.1070899950034990.94645500249825
1040.1724620999977940.3449241999955880.827537900002206
1050.2996061070337560.5992122140675130.700393892966244
1060.3788779146165720.7577558292331440.621122085383428
1070.3912604981891440.7825209963782870.608739501810856
1080.5026405071191230.9947189857617540.497359492880877
1090.5420351157073890.9159297685852230.457964884292611
1100.5500206572921580.8999586854156840.449979342707842
1110.5699276040638620.8601447918722760.430072395936138
1120.6225941281590130.7548117436819730.377405871840987
1130.6082396838471660.7835206323056670.391760316152834
1140.6064132974308360.7871734051383280.393586702569164
1150.606895678043850.78620864391230.39310432195615
1160.7282672885765480.5434654228469030.271732711423452
1170.7459572441345110.5080855117309770.254042755865489
1180.7836810720479530.4326378559040940.216318927952047
1190.814467723872040.3710645522559190.185532276127959
1200.8097015938070480.3805968123859050.190298406192952
1210.795161827450350.4096763450993000.204838172549650
1220.8613136894503760.2773726210992480.138686310549624
1230.852785068334730.2944298633305410.147214931665270
1240.8548055982512730.2903888034974540.145194401748727
1250.8335711973043660.3328576053912680.166428802695634
1260.8055005528633420.3889988942733160.194499447136658
1270.7805862665123590.4388274669752820.219413733487641
1280.7627535734004550.474492853199090.237246426599545
1290.74392675823640.51214648352720.2560732417636
1300.7234775532009490.5530448935981020.276522446799051
1310.6885105064146450.6229789871707090.311489493585355
1320.7023926210607540.5952147578784930.297607378939246
1330.6769472275401780.6461055449196450.323052772459822
1340.633361708830250.73327658233950.36663829116975
1350.5981147188561750.803770562287650.401885281143825
1360.5517434419163580.8965131161672840.448256558083642
1370.5504522883525890.8990954232948220.449547711647411
1380.5092164843015050.981567031396990.490783515698495
1390.4947317639332350.989463527866470.505268236066765
1400.4635952425947740.9271904851895480.536404757405226
1410.6423299532241360.7153400935517280.357670046775864
1420.6271890566895750.745621886620850.372810943310425
1430.58656584207170.82686831585660.4134341579283
1440.686111924889280.627776150221440.31388807511072
1450.6721254244346310.6557491511307370.327874575565369
1460.6560758611679630.6878482776640740.343924138832037
1470.63745393403210.72509213193580.3625460659679
1480.6560107929006560.6879784141986880.343989207099344
1490.6323246673123670.7353506653752650.367675332687633
1500.6437829011586290.7124341976827410.356217098841371
1510.9308801661403370.1382396677193260.0691198338596628
1520.9330361417428150.1339277165143710.0669638582571855
1530.9400483811133790.1199032377732430.0599516188866215
1540.9208845014148020.1582309971703960.0791154985851978
1550.9067686349997820.1864627300004350.0932313650002177
1560.9041607970341080.1916784059317840.0958392029658918
1570.886595445529530.2268091089409420.113404554470471
1580.913271422363640.1734571552727200.0867285776363601
1590.9057813247042760.1884373505914480.0942186752957242
1600.8911363478319360.2177273043361280.108863652168064
1610.8635068494080560.2729863011838890.136493150591944
1620.9294887220549160.1410225558901670.0705112779450837
1630.9070472570682770.1859054858634450.0929527429317225
1640.8873156545433380.2253686909133240.112684345456662
1650.8538742746931560.2922514506136880.146125725306844
1660.8177697226753250.3644605546493500.182230277324675
1670.7710120968146310.4579758063707380.228987903185369
1680.7453106956404380.5093786087191230.254689304359562
1690.6984394098518680.6031211802962650.301560590148132
1700.6612528701591490.6774942596817010.338747129840851
1710.6738328798224750.652334240355050.326167120177525
1720.6199024283959780.7601951432080440.380097571604022
1730.5705354726475070.8589290547049860.429464527352493
1740.8670578213738010.2658843572523980.132942178626199
1750.9403238510132360.1193522979735290.0596761489867643
1760.8959605477625680.2080789044748630.104039452237432
1770.8637453654588760.2725092690822480.136254634541124
1780.8199048948485780.3601902103028440.180095105151422
1790.9042047401915820.1915905196168350.0957952598084175
1800.9067935021349570.1864129957300860.0932064978650429


Meta Analysis of Goldfeld-Quandt test for Heteroskedasticity
Description# significant tests% significant testsOK/NOK
1% type I error level740.427745664739884NOK
5% type I error level920.531791907514451NOK
10% type I error level940.543352601156069NOK
 
Charts produced by software:
http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2008/Dec/05/t1228471629zz53asfw3z5xpue/108yk31228471140.png (open in new window)
http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2008/Dec/05/t1228471629zz53asfw3z5xpue/108yk31228471140.ps (open in new window)


http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2008/Dec/05/t1228471629zz53asfw3z5xpue/1i0y81228471139.png (open in new window)
http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2008/Dec/05/t1228471629zz53asfw3z5xpue/1i0y81228471139.ps (open in new window)


http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2008/Dec/05/t1228471629zz53asfw3z5xpue/24rzm1228471139.png (open in new window)
http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2008/Dec/05/t1228471629zz53asfw3z5xpue/24rzm1228471139.ps (open in new window)


http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2008/Dec/05/t1228471629zz53asfw3z5xpue/3k7op1228471139.png (open in new window)
http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2008/Dec/05/t1228471629zz53asfw3z5xpue/3k7op1228471139.ps (open in new window)


http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2008/Dec/05/t1228471629zz53asfw3z5xpue/4stkb1228471139.png (open in new window)
http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2008/Dec/05/t1228471629zz53asfw3z5xpue/4stkb1228471139.ps (open in new window)


http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2008/Dec/05/t1228471629zz53asfw3z5xpue/5nkru1228471139.png (open in new window)
http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2008/Dec/05/t1228471629zz53asfw3z5xpue/5nkru1228471139.ps (open in new window)


http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2008/Dec/05/t1228471629zz53asfw3z5xpue/68hes1228471139.png (open in new window)
http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2008/Dec/05/t1228471629zz53asfw3z5xpue/68hes1228471139.ps (open in new window)


http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2008/Dec/05/t1228471629zz53asfw3z5xpue/79eij1228471140.png (open in new window)
http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2008/Dec/05/t1228471629zz53asfw3z5xpue/79eij1228471140.ps (open in new window)


http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2008/Dec/05/t1228471629zz53asfw3z5xpue/8gwx41228471140.png (open in new window)
http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2008/Dec/05/t1228471629zz53asfw3z5xpue/8gwx41228471140.ps (open in new window)


http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2008/Dec/05/t1228471629zz53asfw3z5xpue/96pmu1228471140.png (open in new window)
http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2008/Dec/05/t1228471629zz53asfw3z5xpue/96pmu1228471140.ps (open in new window)


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