Home » date » 2010 » Dec » 24 »

MR

*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, 24 Dec 2010 13:32:30 +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/24/t1293197449u4ehrs4m9jxe9hr.htm/, Retrieved Fri, 24 Dec 2010 14:31:01 +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/24/t1293197449u4ehrs4m9jxe9hr.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 «
4 1 27 5 26 49 35 4 1 36 4 25 45 34 5 1 25 4 17 54 13 2 1 27 3 37 36 35 3 2 25 3 35 36 28 5 2 44 3 15 53 32 4 1 50 4 27 46 35 4 1 41 4 36 42 36 4 1 48 5 25 41 27 4 2 43 4 30 45 29 5 2 47 2 27 47 27 4 2 41 3 33 42 28 3 1 44 2 29 45 29 4 2 47 5 30 40 28 3 2 40 3 25 45 30 3 2 46 3 23 40 25 4 1 28 3 26 42 15 3 1 56 3 24 45 33 4 2 49 4 35 47 31 2 2 25 4 39 31 37 4 2 41 4 23 46 37 3 2 26 3 32 34 34 4 1 50 5 29 43 32 4 1 47 4 26 45 21 3 1 52 2 21 42 25 3 2 37 5 35 51 32 2 2 41 3 23 44 28 4 1 45 4 21 47 22 5 2 26 4 28 47 25 4 1 0 3 30 41 26 2 1 52 4 21 44 34 5 1 46 2 29 51 34 4 1 58 3 28 46 36 3 1 54 5 19 47 36 4 1 29 3 26 46 26 2 2 50 3 33 38 26 3 1 43 2 34 50 34 3 2 30 3 33 48 33 3 2 47 2 40 36 31 5 1 45 3 24 51 33 0 2 48 1 35 35 22 4 2 48 3 35 49 29 4 2 26 4 32 38 24 4 1 46 5 20 47 37 2 2 0 3 35 36 32 4 2 50 3 35 47 23 3 1 25 4 21 46 29 4 1 47 2 33 43 35 1 2 47 2 40 53 20 2 1 41 3 22 55 28 2 2 45 2 35 39 26 4 2 41 4 20 55 36 3 2 45 5 28 41 26 4 2 40 3 46 33 33 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 time14 seconds
R Server'George Udny Yule' @ 72.249.76.132


Multiple Linear Regression - Estimated Regression Equation
Teamwork33[t] = -0.627811222048486 + 0.000559943079122198geslacht[t] -0.00134553113394833leeftijd[t] + 0.096198517648116opleiding[t] + 0.0166751551093951Neuroticisme[t] + 0.0635975737699658Extraversie[t] + 0.0110399029138406Openheid[t] + e[t]


Multiple Linear Regression - Ordinary Least Squares
VariableParameterS.D.T-STAT
H0: parameter = 0
2-tail p-value1-tail p-value
(Intercept)-0.6278112220484860.816352-0.7690.4428320.221416
geslacht0.0005599430791221980.1491910.00380.9970090.498505
leeftijd-0.001345531133948330.006477-0.20770.8356650.417832
opleiding0.0961985176481160.0787211.2220.2232320.111616
Neuroticisme0.01667515510939510.0108541.53630.1261430.063072
Extraversie0.06359757376996580.0124895.09231e-060
Openheid0.01103990291384060.0146230.7550.4512170.225608


Multiple Linear Regression - Regression Statistics
Multiple R0.386354784279872
R-squared0.149270019335946
Adjusted R-squared0.122119062506242
F-TEST (value)5.49778117479238
F-TEST (DF numerator)6
F-TEST (DF denominator)188
p-value2.87937910941061e-05
Multiple Linear Regression - Residual Statistics
Residual Standard Deviation1.00036654065481
Sum Squared Residuals188.137844544393


Multiple Linear Regression - Actuals, Interpolation, and Residuals
Time or IndexActualsInterpolation
Forecast
Residuals
Prediction Error
143.753643718211550.246356281788452
243.363230067254870.636769932745128
353.585169871592181.41483012840782
422.91790493010918-0.917904930109179
532.810526304840530.189473695159474
653.576776476852391.42322352314761
743.452380418282190.547619581717807
843.371216202306260.628783797693738
943.111612595818860.888387404181138
1043.382547553374130.61745244662587
1153.239858269926171.76014173007383
1243.137232939098360.862767060901642
1333.17156988875543-0.171569888755432
1443.144336174722780.855663825277217
1533.21804975649472-0.218049756494724
1632.803438876053210.196561123946788
1742.893920077114871.10607992288513
1833.21240586890456-0.212405868904555
1943.607125095485030.392874904514971
2022.75479670030096-0.754796700300959
2143.420429326956950.579570673043048
2232.698199578321500.301800421678495
2343.358016816097680.64198318390232
2443.221585642010910.77841435798909
2532.792452065843420.207547934156575
2633.98490018473423-0.984900184734229
2723.09767653554434-1.09767653554434
2843.27913597918560.720864020814397
2953.455106808317031.54489319168297
3043.056136927585280.943863072414716
3123.21140337490415-1.21140337490415
3253.605663783676531.39433621632347
3343.373132709585730.626867290414273
3433.48443304720316-0.484433047203162
3543.268403773113030.731596226886968
3622.84865305798528-0.848653057985278
3733.62947857885539-0.629478578855389
3833.58881873876079-0.588818738760788
3932.801221586534040.198778413465955
4053.608792153997781.39120784600222
4102.45734506221047-2.45734506221047
4243.617387450683110.382612549316889
4342.938389361911081.06161063808892
4443.522912354297980.477087645702016
4522.8883241948446-0.888324194844597
4643.421261823392240.57873817660776
4733.31972834849149-0.319728348491487
4843.173278185734280.82672181426572
4913.76094140857122-2.76094140857122
5023.78001474882545-1.78001474882545
5122.85613007999566-0.85613007999566
5243.931742122644620.0682578773553817
5333.15519469471417-0.155194694714174
5442.838176837293951.16182316270605
5533.60174658959194-0.601746589591941
5633.16666188904272-0.166661889042723
5754.041436534792420.958563465207578
5833.32070047122651-0.320700471226513
5922.58354337094263-0.583543370942629
6013.53321023251919-2.53321023251919
6123.23052946146232-1.23052946146232
6253.149606880030051.85039311996995
6343.695257085424490.304742914575507
6443.485694311530610.514305688469387
6533.64495012833058-0.644950128330581
6643.635796226424990.364203773575012
6743.278641412238360.721358587761636
6822.95253250674872-0.952532506748719
6933.07083049868552-0.0708304986855244
7043.001057263096050.998942736903952
7132.945204156830540.0547958431694561
7222.40128891354019-0.401288913540193
7342.525060368455611.47493963154439
7443.410996123175890.589003876824111
7532.878277269114630.121722730885366
7654.136734036002110.863265963997893
7712.61092628989287-1.61092628989287
7833.73832572235812-0.738325722358117
7934.06327627005132-1.06327627005132
8053.767228883744381.23277111625562
8123.38412736487840-1.38412736487840
8232.555808096911260.444191903088743
8333.74968470747056-0.749684707470564
8443.635138138276830.364861861723174
8522.99149131654039-0.991491316540394
8643.784266879454170.215733120545827
8732.658768559755980.341231440244022
8833.09770696706006-0.0977069670600576
8933.22805600383106-0.228056003831056
9023.76185677584723-1.76185677584723
9133.63001282029384-0.630012820293843
9223.09090756665300-1.09090756665300
9342.870541524502211.12945847549779
9443.464548247693350.535451752306653
9523.42974985231713-1.42974985231713
9612.52580681682443-1.52580681682443
9753.25946080422371.7405391957763
9843.1212040487590.878795951241001
9943.586576017612640.413423982387363
10043.73957841154130.260421588458703
10133.00860882731488-0.00860882731488219
10233.25408229833061-0.254082298330611
10312.49944736682413-1.49944736682413
10453.179838944325001.82016105567500
10533.35390698895573-0.353906988955732
10633.37153522425452-0.371535224254524
10722.79459903959432-0.794599039594317
10843.492532926379180.507467073620822
10942.971333598429301.02866640157070
11033.73719630685801-0.737196306858013
11142.839477451068761.16052254893124
11243.258580289943180.741419710056822
11322.53987277880064-0.539872778800642
11432.956085371013530.0439146289864661
11533.7157707752796-0.715770775279602
11633.62629693309924-0.626296933099238
11743.184362064335010.815637935664988
11852.384889047456772.61511095254323
11933.28065199518145-0.280651995181448
12033.26004845610290-0.260048456102898
12122.49907239445033-0.499072394450325
12233.28377844945346-0.283778449453464
12311.30155966924746-0.301559669247457
12443.112484597185270.887515402814727
12543.440253949817640.55974605018236
12642.613583988168061.38641601183194
12733.63077817620244-0.630778176202442
12853.871946071594451.12805392840555
12923.59680481356183-1.59680481356183
13023.29102684628615-1.29102684628615
13133.01470838295106-0.0147083829510622
13233.37380048411051-0.373800484110505
13323.51245628936445-1.51245628936445
13413.3093877721919-2.30938777219190
13533.08790880754384-0.0879088075438377
13653.619141913239111.38085808676089
13743.217230014135470.782769985864535
13843.306707547734220.693292452265785
13943.234957002084180.76504299791582
14033.16176574378456-0.161765743784560
14152.995229790167842.00477020983216
14233.26072750645867-0.260727506458674
14332.976188445187620.0238115548123835
14433.17958283646325-0.179582836463246
14533.38795725170416-0.387957251704158
14643.329503981825210.670496018174786
14723.06873231644301-1.06873231644301
14821.676437597233560.32356240276644
14942.665705779723291.33429422027671
15032.776385254809160.223614745190840
15133.55472342364738-0.554723423647376
15222.92830572717015-0.928305727170146
15333.19357053542211-0.193570535422113
15433.40556090593604-0.405560905936042
15543.513280258163890.486719741836113
15612.70783055026909-1.70783055026909
15713.19984171439962-2.19984171439962
15853.268221569019021.73177843098098
15943.451456875661270.54854312433873
16033.23016669561742-0.230166695617418
16133.25380016812417-0.253800168124167
16243.425493495248230.574506504751773
16333.48901795556219-0.489017955562189
16422.60178464202460-0.601784642024597
16512.77973896192742-1.77973896192742
16613.03803355103566-2.03803355103566
16754.006805766556620.99319423344338
16843.72275529775250.277244702247502
16933.20311190091086-0.203111900910859
17043.604464404482510.395535595517493
17153.60350547384431.3964945261557
17243.596091864191560.403908135808436
17342.405007793462121.59499220653788
17423.34211939143489-1.34211939143489
17533.30298915510907-0.302989155109067
17643.312014248639650.687985751360353
17733.04548299925180-0.0454829992518048
17843.375684448568830.624315551431166
17933.3783165802092-0.378316580209203
18043.051427269282320.948572730717685
18113.17262874550607-2.17262874550607
18223.32924153260512-1.32924153260512
18333.74283318327879-0.742833183278789
18433.35989410893563-0.359894108935633
18553.82943606927161.17056393072840
18644.00527605502908-0.0052760550290807
18733.8861961539305-0.886196153930499
18833.38824036747852-0.388240367478520
18932.921250206376010.0787497936239947
19033.1698491899915-0.169849189991498
19143.084228248116170.915771751883827
19233.15805426554439-0.158054265544392
19322.62202281304329-0.622022813043293
19443.044716957177480.955283042822517
19523.40359518010176-1.40359518010176


Goldfeld-Quandt test for Heteroskedasticity
p-valuesAlternative Hypothesis
breakpoint indexgreater2-sidedless
100.1120368803851570.2240737607703130.887963119614843
110.06258467452553380.1251693490510680.937415325474466
120.02357307482251950.04714614964503910.97642692517748
130.03539536651229810.07079073302459610.964604633487702
140.01804666314540410.03609332629080810.981953336854596
150.04048790297909770.08097580595819550.959512097020902
160.02382427095137990.04764854190275990.97617572904862
170.01781647217812090.03563294435624190.98218352782188
180.01441924873420750.0288384974684150.985580751265792
190.01437958421164540.02875916842329080.985620415788355
200.007791346500305480.01558269300061100.992208653499695
210.004066791329649460.008133582659298910.99593320867035
220.003328251106190690.006656502212381380.99667174889381
230.001746498994956720.003492997989913440.998253501005043
240.0009016232147158820.001803246429431760.999098376785284
250.0004410418111623880.0008820836223247760.999558958188838
260.005300863571799260.01060172714359850.9946991364282
270.02847953439102520.05695906878205040.971520465608975
280.01962122876044480.03924245752088970.980378771239555
290.02315139865896690.04630279731793390.976848601341033
300.01867333321079590.03734666642159170.981326666789204
310.03553604616212980.07107209232425970.96446395383787
320.04109228750928390.08218457501856790.958907712490716
330.03206362279182760.06412724558365520.967936377208172
340.02484610541261930.04969221082523860.97515389458738
350.01755825679952590.03511651359905190.982441743200474
360.02099106111291690.04198212222583380.979008938887083
370.02611128652795910.05222257305591820.97388871347204
380.02739184696048990.05478369392097990.97260815303951
390.02068034555863230.04136069111726460.979319654441368
400.02312991725481000.04625983450962010.97687008274519
410.1148810941265680.2297621882531360.885118905873432
420.09091813301053460.1818362660210690.909081866989465
430.08720468547351510.1744093709470300.912795314526485
440.06882863432279820.1376572686455960.931171365677202
450.07007064978208710.1401412995641740.929929350217913
460.0554451755073270.1108903510146540.944554824492673
470.05576054872243620.1115210974448720.944239451277564
480.05716086964532580.1143217392906520.942839130354674
490.3272900293584610.6545800587169220.672709970641539
500.5195029421793040.9609941156413910.480497057820696
510.4919076325429590.9838152650859190.508092367457041
520.4490730080952290.8981460161904570.550926991904771
530.411307122316270.822614244632540.58869287768373
540.4765892950554450.953178590110890.523410704944555
550.4813703910159270.9627407820318530.518629608984073
560.4481291432859830.8962582865719660.551870856714017
570.4304252550326240.8608505100652480.569574744967376
580.3885596668840370.7771193337680740.611440333115963
590.3646395628689320.7292791257378630.635360437131068
600.5990425697210520.8019148605578960.400957430278948
610.6094187844046680.7811624311906630.390581215595332
620.7074245358471670.5851509283056660.292575464152833
630.6735781726568780.6528436546862440.326421827343122
640.6415185241371190.7169629517257620.358481475862881
650.6359442623749690.7281114752500620.364055737625031
660.600422294579220.799155410841560.39957770542078
670.5809426643565780.8381146712868430.419057335643422
680.5730011827204820.8539976345590370.426998817279519
690.5295296532718580.9409406934562850.470470346728142
700.5192756093411270.9614487813177470.480724390658873
710.476674288602970.953348577205940.52332571139703
720.4484646449023240.8969292898046480.551535355097676
730.4752111811081690.9504223622163370.524788818891831
740.445883742667910.891767485335820.55411625733209
750.4042263844965270.8084527689930550.595773615503473
760.3864233926065080.7728467852130150.613576607393492
770.4814525125977940.9629050251955870.518547487402206
780.4667317976622570.9334635953245130.533268202337743
790.4705896190448840.941179238089770.529410380955116
800.4898475585621150.979695117124230.510152441437885
810.5553140860613540.8893718278772910.444685913938646
820.5248475797215120.9503048405569760.475152420278488
830.5082141547210440.9835716905579130.491785845278956
840.4704143707539960.9408287415079910.529585629246004
850.4631527241406580.9263054482813150.536847275859342
860.4263120766132830.8526241532265660.573687923386717
870.394430030037020.788860060074040.60556996996298
880.3568867198799910.7137734397599830.643113280120009
890.3238352886720470.6476705773440950.676164711327953
900.4085694876038820.8171389752077640.591430512396118
910.3909574968717910.7819149937435810.609042503128209
920.400852264806630.801704529613260.59914773519337
930.4176116110957620.8352232221915240.582388388904238
940.3896383361982020.7792766723964040.610361663801798
950.4292144955040570.8584289910081140.570785504495943
960.4871488386397070.9742976772794130.512851161360293
970.5686993461116280.8626013077767440.431300653888372
980.5585999818348190.8828000363303620.441400018165181
990.5270027652845210.9459944694309580.472997234715479
1000.4869918857358660.9739837714717320.513008114264134
1010.4456978281881580.8913956563763170.554302171811842
1020.4083067427763450.816613485552690.591693257223655
1030.4629221712973350.9258443425946710.537077828702665
1040.5586200297876780.8827599404246450.441379970212323
1050.521333202907770.957333594184460.47866679709223
1060.4848845411002650.969769082200530.515115458899735
1070.4721058916368940.9442117832737880.527894108363106
1080.4502008972696630.9004017945393260.549799102730337
1090.452398143076710.904796286153420.54760185692329
1100.4289610893805620.8579221787611240.571038910619438
1110.4383398874997670.8766797749995340.561660112500233
1120.4190813757479140.8381627514958280.580918624252086
1130.3919878836497520.7839757672995040.608012116350248
1140.3518288267513460.7036576535026920.648171173248654
1150.3302741275453930.6605482550907860.669725872454607
1160.3069588694670000.6139177389339990.693041130533
1170.2968217223558370.5936434447116740.703178277644163
1180.5487507130481370.9024985739037270.451249286951863
1190.5077411711185430.9845176577629130.492258828881457
1200.4683097571590540.9366195143181080.531690242840946
1210.4393834849837830.8787669699675670.560616515016217
1220.4005383712993930.8010767425987860.599461628700607
1230.3642807652981030.7285615305962060.635719234701897
1240.3542897974476860.7085795948953730.645710202552314
1250.3287555801183410.6575111602366820.671244419881659
1260.3535685247914520.7071370495829050.646431475208548
1270.3302002773553540.6604005547107070.669799722644646
1280.3615615914568490.7231231829136980.638438408543151
1290.4250157845753080.8500315691506160.574984215424692
1300.448301547149730.896603094299460.55169845285027
1310.4053521795136820.8107043590273630.594647820486318
1320.3671114747692030.7342229495384060.632888525230797
1330.4234586948430620.8469173896861230.576541305156938
1340.6170688151651950.765862369669610.382931184834805
1350.5725992635666370.8548014728667270.427400736433363
1360.5952210068131780.8095579863736450.404778993186822
1370.571593911481490.856812177037020.42840608851851
1380.5679333472113530.8641333055772930.432066652788647
1390.5563831256684660.8872337486630670.443616874331534
1400.510737584279540.978524831440920.48926241572046
1410.6612328596468830.6775342807062330.338767140353116
1420.6176836775539940.7646326448920130.382316322446006
1430.5701949472708330.8596101054583340.429805052729167
1440.5365576415983390.9268847168033220.463442358401661
1450.4937514808783020.9875029617566030.506248519121698
1460.4918792110824340.9837584221648680.508120788917566
1470.4673350846954710.9346701693909430.532664915304529
1480.4230376398664920.8460752797329850.576962360133508
1490.5042734532909830.9914530934180350.495726546709017
1500.48625573457510.97251146915020.5137442654249
1510.4598862073042720.9197724146085450.540113792695728
1520.438209009442280.876418018884560.56179099055772
1530.4096109544850280.8192219089700570.590389045514972
1540.3638524968191190.7277049936382380.636147503180881
1550.3213589279367770.6427178558735540.678641072063223
1560.3584797531334140.7169595062668280.641520246866586
1570.5367545933089710.9264908133820570.463245406691029
1580.6149227923817280.7701544152365440.385077207618272
1590.5764416500527020.8471166998945960.423558349947298
1600.5246195738611890.9507608522776220.475380426138811
1610.4688411948400630.9376823896801270.531158805159937
1620.4183584581700660.8367169163401310.581641541829934
1630.3672906028538710.7345812057077430.632709397146129
1640.3205045178527800.6410090357055590.67949548214722
1650.3866325042832450.7732650085664890.613367495716755
1660.6566674909384640.6866650181230720.343332509061536
1670.683087661971240.633824676057520.31691233802876
1680.6364724196344280.7270551607311440.363527580365572
1690.5886558525771490.8226882948457020.411344147422851
1700.526137048520950.94772590295810.47386295147905
1710.5700592237625550.859881552474890.429940776237445
1720.5738527343073810.8522945313852380.426147265692619
1730.6340853658303260.7318292683393480.365914634169674
1740.7195957746138040.5608084507723910.280404225386196
1750.6633465422297880.6733069155404240.336653457770212
1760.6404110067982320.7191779864035370.359588993201768
1770.6069555417075670.7860889165848670.393044458292433
1780.6658383751346290.6683232497307410.334161624865371
1790.7481810899251570.5036378201496860.251818910074843
1800.877633842410860.2447323151782790.122366157589139
1810.8409866704933590.3180266590132830.159013329506641
1820.7850985134533570.4298029730932850.214901486546643
1830.67113849009310.65772301981380.3288615099069
1840.6090558821818240.7818882356363520.390944117818176
1850.5054879596200580.9890240807598850.494512040379942


Meta Analysis of Goldfeld-Quandt test for Heteroskedasticity
Description# significant tests% significant testsOK/NOK
1% type I error level50.0284090909090909NOK
5% type I error level210.119318181818182NOK
10% type I error level290.164772727272727NOK
 
Charts produced by software:
http://www.freestatistics.org/blog/date/2010/Dec/24/t1293197449u4ehrs4m9jxe9hr/10sczj1293197533.png (open in new window)
http://www.freestatistics.org/blog/date/2010/Dec/24/t1293197449u4ehrs4m9jxe9hr/10sczj1293197533.ps (open in new window)


http://www.freestatistics.org/blog/date/2010/Dec/24/t1293197449u4ehrs4m9jxe9hr/14b2p1293197533.png (open in new window)
http://www.freestatistics.org/blog/date/2010/Dec/24/t1293197449u4ehrs4m9jxe9hr/14b2p1293197533.ps (open in new window)


http://www.freestatistics.org/blog/date/2010/Dec/24/t1293197449u4ehrs4m9jxe9hr/24b2p1293197533.png (open in new window)
http://www.freestatistics.org/blog/date/2010/Dec/24/t1293197449u4ehrs4m9jxe9hr/24b2p1293197533.ps (open in new window)


http://www.freestatistics.org/blog/date/2010/Dec/24/t1293197449u4ehrs4m9jxe9hr/34b2p1293197533.png (open in new window)
http://www.freestatistics.org/blog/date/2010/Dec/24/t1293197449u4ehrs4m9jxe9hr/34b2p1293197533.ps (open in new window)


http://www.freestatistics.org/blog/date/2010/Dec/24/t1293197449u4ehrs4m9jxe9hr/4ek1a1293197533.png (open in new window)
http://www.freestatistics.org/blog/date/2010/Dec/24/t1293197449u4ehrs4m9jxe9hr/4ek1a1293197533.ps (open in new window)


http://www.freestatistics.org/blog/date/2010/Dec/24/t1293197449u4ehrs4m9jxe9hr/5ek1a1293197533.png (open in new window)
http://www.freestatistics.org/blog/date/2010/Dec/24/t1293197449u4ehrs4m9jxe9hr/5ek1a1293197533.ps (open in new window)


http://www.freestatistics.org/blog/date/2010/Dec/24/t1293197449u4ehrs4m9jxe9hr/6ek1a1293197533.png (open in new window)
http://www.freestatistics.org/blog/date/2010/Dec/24/t1293197449u4ehrs4m9jxe9hr/6ek1a1293197533.ps (open in new window)


http://www.freestatistics.org/blog/date/2010/Dec/24/t1293197449u4ehrs4m9jxe9hr/7pujv1293197533.png (open in new window)
http://www.freestatistics.org/blog/date/2010/Dec/24/t1293197449u4ehrs4m9jxe9hr/7pujv1293197533.ps (open in new window)


http://www.freestatistics.org/blog/date/2010/Dec/24/t1293197449u4ehrs4m9jxe9hr/8iliy1293197533.png (open in new window)
http://www.freestatistics.org/blog/date/2010/Dec/24/t1293197449u4ehrs4m9jxe9hr/8iliy1293197533.ps (open in new window)


http://www.freestatistics.org/blog/date/2010/Dec/24/t1293197449u4ehrs4m9jxe9hr/9iliy1293197533.png (open in new window)
http://www.freestatistics.org/blog/date/2010/Dec/24/t1293197449u4ehrs4m9jxe9hr/9iliy1293197533.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