#EJERCICIO 1 #Literal 1/1 Calcule las matrices A, P, M
options(scipen= 99999999)
load("C:/Users/HP/Desktop/REBE UNIVERSIDAD/ECONOMETRIA/GUIA ECONOMETRIA/EJERICIO N1/modelo_ventas.RData")
matriz_X_1<-model.matrix(modelo_ventas)
matriz_XX_1<-t(matriz_X_1)%*%matriz_X_1
matriz_A_1<-solve(matriz_XX_1)%*%t(matriz_X_1)
print(matriz_A_1)## 1 2 3 4
## (Intercept) -0.01128647020 0.014103779728 0.03506391877 0.00042833810
## tv -0.00006704103 0.000030949142 -0.00061201930 -0.00021206420
## periodico 0.00139818182 -0.001907246902 -0.00254688165 0.00022932432
## radio -0.00058002134 0.000648666541 -0.00010932839 -0.00018997103
## tv:radio 0.00001481990 0.000001945347 0.00003260788 0.00001485919
## 5 6 7 8
## (Intercept) -0.02802264005 0.04966825571 0.006721406058 0.013685771163
## tv 0.00156830834 -0.00096215598 0.000111663780 -0.000326611039
## periodico 0.00023142788 -0.00327078177 -0.001147087178 0.000018738210
## radio 0.00061595623 -0.00033083611 0.000772490837 -0.000100744979
## tv:radio -0.00003937064 0.00004642785 -0.000008511094 0.000004091601
## 9 10 11 12
## (Intercept) 0.06237349027 0.019005753345 0.0190774881507 0.00360558894
## tv -0.00140215257 -0.000194663316 0.0000131809233 -0.00039756395
## periodico -0.00123590258 0.000225171010 -0.0005248899956 0.00094786986
## radio -0.00131344307 -0.000735913630 -0.0003341352189 -0.00006782492
## tv:radio 0.00004044251 0.000007159682 0.0000007922246 0.00000221642
## 13 14 15 16
## (Intercept) 0.014899113818 0.03634139729 -0.011844700457 -0.00371574112
## tv 0.000213286058 -0.00085844430 0.000196044767 -0.00048543648
## periodico -0.002038352833 -0.00023578263 0.000811932802 0.00097903989
## radio 0.000336815220 -0.00089476990 0.000067852348 -0.00029853987
## tv:radio 0.000005251734 0.00002304953 -0.000002557309 0.00002032708
## 17 18 19 20
## (Intercept) 0.01533281090 -0.01698619508 0.013579431193 0.003694513315
## tv 0.00003888414 -0.00013321239 -0.000102411708 -0.000053437811
## periodico -0.00123652539 0.00195099121 -0.000580428810 0.000138506238
## radio -0.00083378557 -0.00044690226 0.000145541883 0.000151937838
## tv:radio 0.00003057413 0.00001102821 -0.000000845325 -0.000003078161
## 21 22 23 24
## (Intercept) -0.015735702231 0.010614027039 0.0039565024 -0.001153283948
## tv 0.000494790312 -0.000052925056 0.0010010281 0.000090362439
## periodico 0.000793792331 0.000598233661 -0.0020168690 0.000758430507
## radio 0.000066848849 -0.000626475545 0.0008051669 -0.000180595683
## tv:radio -0.000008570347 0.000002653757 -0.0000242403 -0.000003687696
## 25 26 27 28
## (Intercept) 0.020932126859 0.017662003516 0.001906479731 0.000459166814
## tv -0.000212519874 -0.000294160295 -0.000060899474 -0.000027522912
## periodico -0.000575150168 0.000571758575 -0.000044632121 0.000892483305
## radio -0.000222315654 -0.000842174697 0.000474802108 -0.000270548547
## tv:radio 0.000004605371 0.000009579634 -0.000006627237 -0.000001017345
## 29 30 31 32
## (Intercept) -0.009229296933 -0.00074255704 -0.020737133237 -0.00169324433
## tv 0.000025600123 0.00064847607 0.000302455839 0.00053023898
## periodico 0.001140034360 -0.00063161328 0.001724818959 -0.00031085451
## radio 0.000072298919 0.00034965909 -0.000097870752 0.00025494172
## tv:radio -0.000005118001 -0.00001618543 -0.000006563958 -0.00001343427
## 33 34 35 36
## (Intercept) 0.008864197705 0.01027581678 0.04236431872 0.03033028408
## tv 0.000349355310 -0.00066871791 -0.00106304463 -0.00092566223
## periodico -0.000041683492 0.00119509517 0.00006961331 0.00082003595
## radio -0.000335141773 -0.00044960780 -0.00136411967 -0.00130298655
## tv:radio -0.000007465636 0.00001074002 0.00003137340 0.00002639173
## 37 38 39 40
## (Intercept) -0.03691371138 0.01189641422 0.008312017572 -0.016421791419
## tv 0.00035119573 -0.00034519692 0.000282477316 0.000055370346
## periodico 0.00199878941 -0.00124947589 -0.001302070784 0.001259943073
## radio 0.00110703560 0.00045954540 0.000562172245 0.000254407845
## tv:radio -0.00002740139 0.00001270793 -0.000008200132 -0.000004675563
## 41 42 43 44
## (Intercept) -0.007098668942 -0.006645942619 -0.007827522911 0.00562418953
## tv 0.000222581337 0.000155364334 -0.000299207325 0.00010847218
## periodico 0.000747334684 0.000268389768 0.001648736739 0.00051002693
## radio 0.000029797138 0.000282618270 0.000056961391 -0.00039368776
## tv:radio -0.000007123256 -0.000004170304 -0.000003001041 -0.00000234203
## 45 46 47 48
## (Intercept) 0.00884005605 -0.002874908479 0.00100977476 -0.02451539730
## tv 0.00046057305 0.000228756456 0.00056234506 0.00014158543
## periodico -0.00171033165 0.000265781998 -0.00025382772 0.00152356607
## radio 0.00063601703 0.000136973068 0.00008691167 0.00066324823
## tv:radio -0.00001057989 -0.000007309017 -0.00001413098 -0.00001407076
## 49 50 51 52
## (Intercept) -0.01878764724 0.00272615988 -0.00178738864 0.03598768964
## tv 0.00095991001 0.00058696509 0.00060244753 -0.00096600682
## periodico 0.00053544810 -0.00060557497 0.00034275541 -0.00005860236
## radio 0.00026402722 0.00021734621 -0.00018745312 -0.00090579432
## tv:radio -0.00002327068 -0.00001482293 -0.00001462919 0.00002476176
## 53 54 55 56
## (Intercept) -0.015681810129 0.00074396845 -0.011168389305 0.00002541883
## tv -0.000067159493 -0.00051349075 -0.000062236867 -0.00074494491
## periodico 0.001345269333 0.00072784252 0.001418155003 0.00125889173
## radio 0.000103698401 -0.00043431222 0.000114865241 -0.00068987402
## tv:radio 0.000001767024 0.00002405979 -0.000005442162 0.00003168547
## 57 58 59 60
## (Intercept) 0.01831083030 0.010653820682 -0.015790159420 -0.005853894390
## tv 0.00037974671 -0.000185452442 -0.000211903555 -0.000115582009
## periodico -0.00266494448 0.000025266826 0.001251591776 0.000891572787
## radio 0.00085114752 -0.000065272194 0.000267070291 0.000307153563
## tv:radio -0.00000890655 0.000001426743 0.000004245123 -0.000006380628
## 61 62 63 64
## (Intercept) 0.024987561800 -0.01302879544 -0.002445181178 0.004949872745
## tv -0.000173057706 -0.00028160660 0.000128664831 -0.000096027177
## periodico -0.000434643049 0.00173643412 0.000887450274 -0.000331339214
## radio -0.000610378781 -0.00043427213 -0.000232545627 0.000581957346
## tv:radio 0.000006783625 0.00001501912 -0.000004268503 -0.000007163908
## 65 66 67 68
## (Intercept) -0.00821712853 0.04273553577 0.0229473149842 0.01820829672
## tv 0.00006580937 -0.00110033999 -0.0003618915503 -0.00052137470
## periodico 0.00001079112 -0.00041586509 -0.0012710528427 0.00037421605
## radio 0.00072776383 -0.00093213639 0.0004476366163 -0.00050459832
## tv:radio -0.00000712428 0.00002837397 0.0000006298279 0.00001117128
## 69 70 71 72
## (Intercept) -0.00587330896 -0.019410327467 -0.010929740347 0.001720050399
## tv -0.00015525538 0.000059280863 0.000208663285 0.000322859874
## periodico 0.00114264022 0.001140588393 0.000751392556 0.000022233140
## radio 0.00012267595 0.000570942606 0.000152655637 0.000009188374
## tv:radio -0.00000395234 -0.000008096259 -0.000005496981 -0.000008621258
## 73 74 75 76
## (Intercept) 0.014155913033 0.004183081585 0.000096498081 0.04133114660
## tv 0.000108055467 0.000387042857 -0.000182478045 -0.00071691915
## periodico -0.001983762315 0.000107969277 0.000778771733 -0.00259939839
## radio 0.001013932711 -0.000213930377 0.000032463947 -0.00064778051
## tv:radio -0.000009976884 -0.000009149524 -0.000001450933 0.00004497655
## 77 78 79 80
## (Intercept) 0.029052219752 0.004019870279 0.02568587940 0.014215520541
## tv -0.000213088043 -0.000054586137 -0.00002185503 -0.000049182026
## periodico -0.000740382123 -0.000223275356 -0.00276030449 0.000039133548
## radio -0.000591894264 0.000460017434 0.00111509664 -0.000423521968
## tv:radio 0.000008015092 -0.000005969221 -0.00000879480 0.000001929415
## 81 82 83 84
## (Intercept) 0.008345104731 -0.0067821677 0.004921867503 0.005786401831
## tv 0.000056112718 0.0007153471 0.000298223994 -0.000016296534
## periodico -0.000802770542 0.0005249677 -0.000613853773 -0.001299317483
## radio 0.000474863427 -0.0001189448 0.000286351652 0.000828294453
## tv:radio -0.000005880627 -0.0000177861 -0.000008593964 -0.000001622635
## 85 86 87 88
## (Intercept) -0.015196826163 -0.02826803967 0.009478256801 0.00697422257
## tv -0.000019107862 0.00135194466 -0.000035948293 -0.00019524819
## periodico 0.001025333548 0.00045642150 -0.000779061930 -0.00042841743
## radio 0.000362986505 0.00046608800 0.000520009708 -0.00015610155
## tv:radio -0.000002475481 -0.00003065541 -0.000005520579 0.00001652708
## 89 90 91 92
## (Intercept) -0.01345947077 0.00933095016 0.0324700557 0.01006201929
## tv 0.00101277100 -0.00042350545 -0.0008446069 0.00055051470
## periodico -0.00055310023 -0.00060621191 0.0003259917 -0.00069641935
## radio 0.00040642574 0.00009972553 -0.0011152082 -0.00006211493
## tv:radio -0.00001719827 0.00001793239 0.0000238975 -0.00001293978
## 93 94 95 96
## (Intercept) -0.012561667584 -0.01359268266 0.0226058265 -0.008511490317
## tv 0.000232457597 0.00002954089 -0.0004968032 0.000309545708
## periodico 0.000875090545 0.00148757677 -0.0001259276 0.000288419189
## radio -0.000128556632 -0.00059940482 -0.0004127950 0.000113139890
## tv:radio 0.000001429484 0.00001306496 0.0000108430 -0.000003180198
## 97 98 99 100
## (Intercept) 0.03711720155 0.000973109872 -0.01776733122 -0.000306931148
## tv -0.00109143130 -0.000019982349 -0.00022920641 -0.000098468643
## periodico 0.00056483304 0.000543224560 0.00208921866 -0.000150663149
## radio -0.00139218359 -0.000039807925 -0.00039761456 0.000243435856
## tv:radio 0.00003115892 -0.000002490315 0.00001179114 0.000005777537
## 101 102 103 104
## (Intercept) -0.0225338310 -0.01477602943 0.006267206369 0.007047657756
## tv 0.0014546750 0.00002311872 -0.000136410903 -0.000186875718
## periodico 0.0002772215 0.00191700120 0.000964993794 0.000550081239
## radio 0.0004411293 -0.00122375976 -0.000580253503 -0.000256494905
## tv:radio -0.0000376632 0.00002504818 0.000003477592 0.000002322845
## 105 106 107 108
## (Intercept) -0.015584348939 0.00606196517 0.015896337611 0.021292700367
## tv 0.000008587482 -0.00051842974 0.000293954503 -0.000076845097
## periodico 0.001254319136 0.00016389592 -0.001228173794 -0.000183639474
## radio 0.000591893676 -0.00032987456 0.000150658879 -0.000650144117
## tv:radio -0.000013381956 0.00002443714 -0.000007469916 0.000004553308
## 109 110 111 112
## (Intercept) 0.0259081117932 -0.005860740324 -0.03061336840 -0.018733131894
## tv 0.0000990139053 -0.000268584134 0.00161850702 0.000081711294
## periodico -0.0011455106749 0.001438677829 0.00051166319 0.001331666756
## radio -0.0003291803759 0.000033869706 0.00056348576 0.000439565172
## tv:radio -0.0000002575464 -0.000002297405 -0.00004137368 -0.000009291272
## 113 114 115 116
## (Intercept) 0.02231473162 0.006912749678 0.00350450818 0.0057889554831
## tv -0.00078956888 -0.000340935519 -0.00003798319 0.0001746561394
## periodico 0.00054611658 0.000711647822 -0.00114787787 -0.0010782704889
## radio -0.00063433762 -0.000199410936 0.00087218898 0.0003540656623
## tv:radio 0.00001686494 0.000003836647 -0.00000186540 0.0000009186797
## 117 118 119 120
## (Intercept) 0.007621528529 0.03207311114 0.00348876480 0.023749714363
## tv 0.000084597739 -0.00061246798 0.00002447478 0.000002781209
## periodico -0.000008596092 0.00003257176 -0.00028204609 -0.001641919181
## radio -0.000112009007 -0.00105955721 -0.00037873008 0.000235882012
## tv:radio -0.000003007049 0.00001919778 0.00001640435 -0.000001574375
## 121 122 123 124
## (Intercept) -0.008371458679 0.00579076808 0.0246214091 -0.00227449922
## tv 0.000433253314 0.00077346417 -0.0005429257 0.00007232554
## periodico 0.000163710546 -0.00193425007 0.0005456834 -0.00029428162
## radio 0.000230048851 0.00073326798 -0.0010511815 0.00083239722
## tv:radio -0.000009110407 -0.00001741039 0.0000166963 -0.00001233571
## 125 126 127 128
## (Intercept) -0.016078222443 0.011484301170 0.024292095391 0.04260889135
## tv 0.000403387727 0.000102516148 0.000039115008 -0.00099042499
## periodico 0.000997828788 -0.000313829344 -0.002941077130 -0.00005357818
## radio -0.000299975517 -0.000121803519 0.000741096734 -0.00134594321
## tv:radio 0.000002109399 -0.000002888536 0.000004259761 0.00002996835
## 129 130 131 132
## (Intercept) -0.04184894876 -0.00340210494 0.02635875495 -0.01680636474
## tv 0.00060133947 0.00085499383 0.00030830294 0.00110001765
## periodico 0.00149779495 -0.00064927980 -0.00434810562 0.00066593110
## radio 0.00165705121 0.00038887607 0.00208020059 0.00008961022
## tv:radio -0.00003768111 -0.00002141744 -0.00002206911 -0.00002816220
## 133 134 135 136
## (Intercept) 0.02919608430 -0.012799174302 0.01695419897 -0.00456404857
## tv -0.00022524016 0.000167341441 -0.00004276786 0.00049170331
## periodico -0.00247706077 0.000947591304 -0.00176944627 -0.00200884389
## radio 0.00090330909 0.000048936950 0.00014621124 0.00205837354
## tv:radio -0.00000435789 -0.000002109586 0.00001305045 -0.00003052679
## 137 138 139 140
## (Intercept) 0.00778594521 -0.022596153125 0.014558760948 -0.02706944785
## tv 0.00023732597 0.000501127883 0.000030221347 0.00043062493
## periodico -0.00212530270 0.001495386134 -0.001362977777 0.00069016250
## radio 0.00155063090 -0.000139327979 0.000563911227 0.00150299208
## tv:radio -0.00001977492 -0.000006599755 -0.000005442513 -0.00003119953
## 141 142 143 144
## (Intercept) 0.020217704051 -0.0082726940 -0.014022424717 0.00168474054
## tv -0.000338761324 0.0001422678 0.000137252931 0.00056054532
## periodico -0.000472916959 0.0007072714 0.001111115544 -0.00007458885
## radio -0.000126687558 -0.0004324584 0.000110590807 -0.00006142670
## tv:radio 0.000005712816 0.0000109961 -0.000004031128 -0.00001374589
## 145 146 147 148
## (Intercept) -0.00198829800 0.03763194326 0.02629464036 -0.01543522339
## tv 0.00060063272 -0.00095187216 -0.00081254970 -0.00036412185
## periodico -0.00031062935 0.00025268846 0.00074360036 0.00173390852
## radio 0.00023599016 -0.00130077426 -0.00107440049 -0.00012045260
## tv:radio -0.00001515365 0.00002809328 0.00002194616 0.00001218465
## 149 150 151 152
## (Intercept) 0.00403254897 0.013338479273 -0.01215692775 -0.01696418130
## tv 0.00024965446 0.000027895532 0.00051576139 0.00122943937
## periodico -0.00181365109 -0.001219077194 0.00104733986 0.00004097211
## radio 0.00149714575 0.000528800275 -0.00010471182 0.00041295892
## tv:radio -0.00001947075 -0.000005332999 -0.00001315815 -0.00003126069
## 153 154 155 156
## (Intercept) 0.001417534822 -0.008470399675 0.007947807307 0.04992101740
## tv -0.000186653453 0.000013929832 -0.000357034184 -0.00074406765
## periodico 0.000739771135 0.000456080437 0.000603260782 -0.00226306936
## radio -0.000026698877 0.000326717452 -0.000156384232 -0.00023230271
## tv:radio -0.000000566054 -0.000001153531 0.000003807189 0.00001812879
## 157 158 159 160
## (Intercept) 0.00778287975 0.016044205762 0.0198020921866 -0.00083849089
## tv -0.00021113564 -0.000010068402 0.0001294379729 0.00037522400
## periodico -0.00077083109 0.000141570665 -0.0026390976690 -0.00007396010
## radio 0.00022254888 -0.000633397205 0.0008090011419 0.00015964400
## tv:radio 0.00001127609 0.000002441521 -0.0000006920975 -0.00001006745
## 161 162 163 164
## (Intercept) -0.001790463491 0.0046878588293 0.000641191193 -0.0113002582
## tv 0.000243976491 0.0001313245628 0.000075684680 0.0001188405
## periodico 0.000371379207 -0.0009227419417 0.000527587196 0.0003582426
## radio 0.000004212692 0.0003693030404 -0.000100996155 0.0009174504
## tv:radio -0.000007238300 0.0000007849008 -0.000003640183 -0.0000166088
## 165 166 167 168
## (Intercept) 0.02587324499 -0.07225833746 0.01468504348 0.016653514998
## tv -0.00069122946 0.00355381404 0.00016799347 -0.000278836159
## periodico -0.00003379496 0.00018304394 -0.00246725793 0.000531485509
## radio -0.00050191766 0.00189664307 0.00127380853 -0.000762066393
## tv:radio 0.00001503164 -0.00009459444 -0.00001194893 0.000008659777
## 169 170 171 172
## (Intercept) -0.0202275264 0.02145995752 0.024433979551 -0.01151738746
## tv 0.0007922543 -0.00082193713 -0.000215338084 0.00067982384
## periodico 0.0007450594 0.00106072925 -0.000878903608 0.00019398865
## radio 0.0001652481 -0.00099540350 -0.000200975390 0.00025915032
## tv:radio -0.0000160912 0.00002053574 0.000005026328 -0.00001576635
## 173 174 175 176
## (Intercept) 0.0235756400961 0.02514494140 0.02775837039 -0.021103947078
## tv -0.0001232980468 -0.00059562607 -0.00066961037 -0.000312707438
## periodico -0.0015784204846 0.00031894679 0.00047486399 0.002179846329
## radio 0.0003275548165 -0.00085657634 -0.00108649690 -0.000111475421
## tv:radio -0.0000003579651 0.00001637958 0.00001978627 0.000009318024
## 177 178 179 180
## (Intercept) -0.01266987581 -0.00181364734 0.012597917565 0.015766611551
## tv 0.00001151397 0.00056422607 -0.000052667635 -0.000302612184
## periodico 0.00133066382 0.00016722575 0.000560014361 0.000378948165
## radio 0.00016996396 -0.00005098358 -0.000713286259 -0.000567858400
## tv:radio -0.00000642050 -0.00001404644 0.000003337299 0.000007731302
## 181 182 183 184
## (Intercept) 0.0375360999 0.006249894866 0.01182661003 -0.01048364001
## tv -0.0009729864 0.000166421440 0.00031343971 -0.00052366780
## periodico 0.0002714128 0.000477573166 -0.00052064220 0.00227758259
## radio -0.0012900129 -0.000444787514 -0.00013754438 -0.00108855870
## tv:radio 0.0000283698 -0.000003254349 -0.00000715736 0.00003015597
## 185 186 187 188
## (Intercept) -0.00878542694 -0.02403900996 0.012112048659 -0.004594159372
## tv 0.00018271141 0.00018060293 0.000133139885 -0.000015963520
## periodico 0.00108395802 0.00115362937 0.000138903045 0.000616178833
## radio -0.00007791804 0.00087000126 -0.000505287750 0.000269956649
## tv:radio -0.00000628598 -0.00001563611 -0.000001664489 -0.000005867028
## 189 190 191 192
## (Intercept) 0.01849130882 0.0234202371995 0.00170359876 0.03410061697
## tv -0.00081441237 0.0000187047846 0.00032730615 -0.00080167299
## periodico 0.00112799685 -0.0013998890588 -0.00188675507 -0.00036176657
## radio -0.00084391921 0.0000586567527 0.00172810634 -0.00068749942
## tv:radio 0.00001843973 -0.0000008807575 -0.00002479561 0.00001995866
## 193 194 195 196
## (Intercept) 0.01595370696 -0.02143366996 -0.00858635415 0.03609474671
## tv 0.00044877998 0.00033114925 0.00008016567 -0.00059753273
## periodico -0.00121850007 0.00049794605 0.00023567381 -0.00062652493
## radio 0.00004599813 0.00134310708 0.00088908809 -0.00080528477
## tv:radio -0.00001048754 -0.00002656952 -0.00001558263 0.00001780196
## 197 198 199 200
## (Intercept) 0.03769835435 0.02808613397 -0.01282058462 0.0246735133
## tv -0.00090291732 -0.00085615088 -0.00037778479 -0.0007715129
## periodico -0.00008300596 0.00052741693 0.00213262423 0.0007233832
## radio -0.00107306292 -0.00096939133 -0.00084263918 -0.0009870685
## tv:radio 0.00002548349 0.00002210282 0.00002311277 0.0000202718#Matriz P
matriz_P_1<-matriz_X_1%*%matriz_A_1
print(matriz_A_1)## 1 2 3 4
## (Intercept) -0.01128647020 0.014103779728 0.03506391877 0.00042833810
## tv -0.00006704103 0.000030949142 -0.00061201930 -0.00021206420
## periodico 0.00139818182 -0.001907246902 -0.00254688165 0.00022932432
## radio -0.00058002134 0.000648666541 -0.00010932839 -0.00018997103
## tv:radio 0.00001481990 0.000001945347 0.00003260788 0.00001485919
## 5 6 7 8
## (Intercept) -0.02802264005 0.04966825571 0.006721406058 0.013685771163
## tv 0.00156830834 -0.00096215598 0.000111663780 -0.000326611039
## periodico 0.00023142788 -0.00327078177 -0.001147087178 0.000018738210
## radio 0.00061595623 -0.00033083611 0.000772490837 -0.000100744979
## tv:radio -0.00003937064 0.00004642785 -0.000008511094 0.000004091601
## 9 10 11 12
## (Intercept) 0.06237349027 0.019005753345 0.0190774881507 0.00360558894
## tv -0.00140215257 -0.000194663316 0.0000131809233 -0.00039756395
## periodico -0.00123590258 0.000225171010 -0.0005248899956 0.00094786986
## radio -0.00131344307 -0.000735913630 -0.0003341352189 -0.00006782492
## tv:radio 0.00004044251 0.000007159682 0.0000007922246 0.00000221642
## 13 14 15 16
## (Intercept) 0.014899113818 0.03634139729 -0.011844700457 -0.00371574112
## tv 0.000213286058 -0.00085844430 0.000196044767 -0.00048543648
## periodico -0.002038352833 -0.00023578263 0.000811932802 0.00097903989
## radio 0.000336815220 -0.00089476990 0.000067852348 -0.00029853987
## tv:radio 0.000005251734 0.00002304953 -0.000002557309 0.00002032708
## 17 18 19 20
## (Intercept) 0.01533281090 -0.01698619508 0.013579431193 0.003694513315
## tv 0.00003888414 -0.00013321239 -0.000102411708 -0.000053437811
## periodico -0.00123652539 0.00195099121 -0.000580428810 0.000138506238
## radio -0.00083378557 -0.00044690226 0.000145541883 0.000151937838
## tv:radio 0.00003057413 0.00001102821 -0.000000845325 -0.000003078161
## 21 22 23 24
## (Intercept) -0.015735702231 0.010614027039 0.0039565024 -0.001153283948
## tv 0.000494790312 -0.000052925056 0.0010010281 0.000090362439
## periodico 0.000793792331 0.000598233661 -0.0020168690 0.000758430507
## radio 0.000066848849 -0.000626475545 0.0008051669 -0.000180595683
## tv:radio -0.000008570347 0.000002653757 -0.0000242403 -0.000003687696
## 25 26 27 28
## (Intercept) 0.020932126859 0.017662003516 0.001906479731 0.000459166814
## tv -0.000212519874 -0.000294160295 -0.000060899474 -0.000027522912
## periodico -0.000575150168 0.000571758575 -0.000044632121 0.000892483305
## radio -0.000222315654 -0.000842174697 0.000474802108 -0.000270548547
## tv:radio 0.000004605371 0.000009579634 -0.000006627237 -0.000001017345
## 29 30 31 32
## (Intercept) -0.009229296933 -0.00074255704 -0.020737133237 -0.00169324433
## tv 0.000025600123 0.00064847607 0.000302455839 0.00053023898
## periodico 0.001140034360 -0.00063161328 0.001724818959 -0.00031085451
## radio 0.000072298919 0.00034965909 -0.000097870752 0.00025494172
## tv:radio -0.000005118001 -0.00001618543 -0.000006563958 -0.00001343427
## 33 34 35 36
## (Intercept) 0.008864197705 0.01027581678 0.04236431872 0.03033028408
## tv 0.000349355310 -0.00066871791 -0.00106304463 -0.00092566223
## periodico -0.000041683492 0.00119509517 0.00006961331 0.00082003595
## radio -0.000335141773 -0.00044960780 -0.00136411967 -0.00130298655
## tv:radio -0.000007465636 0.00001074002 0.00003137340 0.00002639173
## 37 38 39 40
## (Intercept) -0.03691371138 0.01189641422 0.008312017572 -0.016421791419
## tv 0.00035119573 -0.00034519692 0.000282477316 0.000055370346
## periodico 0.00199878941 -0.00124947589 -0.001302070784 0.001259943073
## radio 0.00110703560 0.00045954540 0.000562172245 0.000254407845
## tv:radio -0.00002740139 0.00001270793 -0.000008200132 -0.000004675563
## 41 42 43 44
## (Intercept) -0.007098668942 -0.006645942619 -0.007827522911 0.00562418953
## tv 0.000222581337 0.000155364334 -0.000299207325 0.00010847218
## periodico 0.000747334684 0.000268389768 0.001648736739 0.00051002693
## radio 0.000029797138 0.000282618270 0.000056961391 -0.00039368776
## tv:radio -0.000007123256 -0.000004170304 -0.000003001041 -0.00000234203
## 45 46 47 48
## (Intercept) 0.00884005605 -0.002874908479 0.00100977476 -0.02451539730
## tv 0.00046057305 0.000228756456 0.00056234506 0.00014158543
## periodico -0.00171033165 0.000265781998 -0.00025382772 0.00152356607
## radio 0.00063601703 0.000136973068 0.00008691167 0.00066324823
## tv:radio -0.00001057989 -0.000007309017 -0.00001413098 -0.00001407076
## 49 50 51 52
## (Intercept) -0.01878764724 0.00272615988 -0.00178738864 0.03598768964
## tv 0.00095991001 0.00058696509 0.00060244753 -0.00096600682
## periodico 0.00053544810 -0.00060557497 0.00034275541 -0.00005860236
## radio 0.00026402722 0.00021734621 -0.00018745312 -0.00090579432
## tv:radio -0.00002327068 -0.00001482293 -0.00001462919 0.00002476176
## 53 54 55 56
## (Intercept) -0.015681810129 0.00074396845 -0.011168389305 0.00002541883
## tv -0.000067159493 -0.00051349075 -0.000062236867 -0.00074494491
## periodico 0.001345269333 0.00072784252 0.001418155003 0.00125889173
## radio 0.000103698401 -0.00043431222 0.000114865241 -0.00068987402
## tv:radio 0.000001767024 0.00002405979 -0.000005442162 0.00003168547
## 57 58 59 60
## (Intercept) 0.01831083030 0.010653820682 -0.015790159420 -0.005853894390
## tv 0.00037974671 -0.000185452442 -0.000211903555 -0.000115582009
## periodico -0.00266494448 0.000025266826 0.001251591776 0.000891572787
## radio 0.00085114752 -0.000065272194 0.000267070291 0.000307153563
## tv:radio -0.00000890655 0.000001426743 0.000004245123 -0.000006380628
## 61 62 63 64
## (Intercept) 0.024987561800 -0.01302879544 -0.002445181178 0.004949872745
## tv -0.000173057706 -0.00028160660 0.000128664831 -0.000096027177
## periodico -0.000434643049 0.00173643412 0.000887450274 -0.000331339214
## radio -0.000610378781 -0.00043427213 -0.000232545627 0.000581957346
## tv:radio 0.000006783625 0.00001501912 -0.000004268503 -0.000007163908
## 65 66 67 68
## (Intercept) -0.00821712853 0.04273553577 0.0229473149842 0.01820829672
## tv 0.00006580937 -0.00110033999 -0.0003618915503 -0.00052137470
## periodico 0.00001079112 -0.00041586509 -0.0012710528427 0.00037421605
## radio 0.00072776383 -0.00093213639 0.0004476366163 -0.00050459832
## tv:radio -0.00000712428 0.00002837397 0.0000006298279 0.00001117128
## 69 70 71 72
## (Intercept) -0.00587330896 -0.019410327467 -0.010929740347 0.001720050399
## tv -0.00015525538 0.000059280863 0.000208663285 0.000322859874
## periodico 0.00114264022 0.001140588393 0.000751392556 0.000022233140
## radio 0.00012267595 0.000570942606 0.000152655637 0.000009188374
## tv:radio -0.00000395234 -0.000008096259 -0.000005496981 -0.000008621258
## 73 74 75 76
## (Intercept) 0.014155913033 0.004183081585 0.000096498081 0.04133114660
## tv 0.000108055467 0.000387042857 -0.000182478045 -0.00071691915
## periodico -0.001983762315 0.000107969277 0.000778771733 -0.00259939839
## radio 0.001013932711 -0.000213930377 0.000032463947 -0.00064778051
## tv:radio -0.000009976884 -0.000009149524 -0.000001450933 0.00004497655
## 77 78 79 80
## (Intercept) 0.029052219752 0.004019870279 0.02568587940 0.014215520541
## tv -0.000213088043 -0.000054586137 -0.00002185503 -0.000049182026
## periodico -0.000740382123 -0.000223275356 -0.00276030449 0.000039133548
## radio -0.000591894264 0.000460017434 0.00111509664 -0.000423521968
## tv:radio 0.000008015092 -0.000005969221 -0.00000879480 0.000001929415
## 81 82 83 84
## (Intercept) 0.008345104731 -0.0067821677 0.004921867503 0.005786401831
## tv 0.000056112718 0.0007153471 0.000298223994 -0.000016296534
## periodico -0.000802770542 0.0005249677 -0.000613853773 -0.001299317483
## radio 0.000474863427 -0.0001189448 0.000286351652 0.000828294453
## tv:radio -0.000005880627 -0.0000177861 -0.000008593964 -0.000001622635
## 85 86 87 88
## (Intercept) -0.015196826163 -0.02826803967 0.009478256801 0.00697422257
## tv -0.000019107862 0.00135194466 -0.000035948293 -0.00019524819
## periodico 0.001025333548 0.00045642150 -0.000779061930 -0.00042841743
## radio 0.000362986505 0.00046608800 0.000520009708 -0.00015610155
## tv:radio -0.000002475481 -0.00003065541 -0.000005520579 0.00001652708
## 89 90 91 92
## (Intercept) -0.01345947077 0.00933095016 0.0324700557 0.01006201929
## tv 0.00101277100 -0.00042350545 -0.0008446069 0.00055051470
## periodico -0.00055310023 -0.00060621191 0.0003259917 -0.00069641935
## radio 0.00040642574 0.00009972553 -0.0011152082 -0.00006211493
## tv:radio -0.00001719827 0.00001793239 0.0000238975 -0.00001293978
## 93 94 95 96
## (Intercept) -0.012561667584 -0.01359268266 0.0226058265 -0.008511490317
## tv 0.000232457597 0.00002954089 -0.0004968032 0.000309545708
## periodico 0.000875090545 0.00148757677 -0.0001259276 0.000288419189
## radio -0.000128556632 -0.00059940482 -0.0004127950 0.000113139890
## tv:radio 0.000001429484 0.00001306496 0.0000108430 -0.000003180198
## 97 98 99 100
## (Intercept) 0.03711720155 0.000973109872 -0.01776733122 -0.000306931148
## tv -0.00109143130 -0.000019982349 -0.00022920641 -0.000098468643
## periodico 0.00056483304 0.000543224560 0.00208921866 -0.000150663149
## radio -0.00139218359 -0.000039807925 -0.00039761456 0.000243435856
## tv:radio 0.00003115892 -0.000002490315 0.00001179114 0.000005777537
## 101 102 103 104
## (Intercept) -0.0225338310 -0.01477602943 0.006267206369 0.007047657756
## tv 0.0014546750 0.00002311872 -0.000136410903 -0.000186875718
## periodico 0.0002772215 0.00191700120 0.000964993794 0.000550081239
## radio 0.0004411293 -0.00122375976 -0.000580253503 -0.000256494905
## tv:radio -0.0000376632 0.00002504818 0.000003477592 0.000002322845
## 105 106 107 108
## (Intercept) -0.015584348939 0.00606196517 0.015896337611 0.021292700367
## tv 0.000008587482 -0.00051842974 0.000293954503 -0.000076845097
## periodico 0.001254319136 0.00016389592 -0.001228173794 -0.000183639474
## radio 0.000591893676 -0.00032987456 0.000150658879 -0.000650144117
## tv:radio -0.000013381956 0.00002443714 -0.000007469916 0.000004553308
## 109 110 111 112
## (Intercept) 0.0259081117932 -0.005860740324 -0.03061336840 -0.018733131894
## tv 0.0000990139053 -0.000268584134 0.00161850702 0.000081711294
## periodico -0.0011455106749 0.001438677829 0.00051166319 0.001331666756
## radio -0.0003291803759 0.000033869706 0.00056348576 0.000439565172
## tv:radio -0.0000002575464 -0.000002297405 -0.00004137368 -0.000009291272
## 113 114 115 116
## (Intercept) 0.02231473162 0.006912749678 0.00350450818 0.0057889554831
## tv -0.00078956888 -0.000340935519 -0.00003798319 0.0001746561394
## periodico 0.00054611658 0.000711647822 -0.00114787787 -0.0010782704889
## radio -0.00063433762 -0.000199410936 0.00087218898 0.0003540656623
## tv:radio 0.00001686494 0.000003836647 -0.00000186540 0.0000009186797
## 117 118 119 120
## (Intercept) 0.007621528529 0.03207311114 0.00348876480 0.023749714363
## tv 0.000084597739 -0.00061246798 0.00002447478 0.000002781209
## periodico -0.000008596092 0.00003257176 -0.00028204609 -0.001641919181
## radio -0.000112009007 -0.00105955721 -0.00037873008 0.000235882012
## tv:radio -0.000003007049 0.00001919778 0.00001640435 -0.000001574375
## 121 122 123 124
## (Intercept) -0.008371458679 0.00579076808 0.0246214091 -0.00227449922
## tv 0.000433253314 0.00077346417 -0.0005429257 0.00007232554
## periodico 0.000163710546 -0.00193425007 0.0005456834 -0.00029428162
## radio 0.000230048851 0.00073326798 -0.0010511815 0.00083239722
## tv:radio -0.000009110407 -0.00001741039 0.0000166963 -0.00001233571
## 125 126 127 128
## (Intercept) -0.016078222443 0.011484301170 0.024292095391 0.04260889135
## tv 0.000403387727 0.000102516148 0.000039115008 -0.00099042499
## periodico 0.000997828788 -0.000313829344 -0.002941077130 -0.00005357818
## radio -0.000299975517 -0.000121803519 0.000741096734 -0.00134594321
## tv:radio 0.000002109399 -0.000002888536 0.000004259761 0.00002996835
## 129 130 131 132
## (Intercept) -0.04184894876 -0.00340210494 0.02635875495 -0.01680636474
## tv 0.00060133947 0.00085499383 0.00030830294 0.00110001765
## periodico 0.00149779495 -0.00064927980 -0.00434810562 0.00066593110
## radio 0.00165705121 0.00038887607 0.00208020059 0.00008961022
## tv:radio -0.00003768111 -0.00002141744 -0.00002206911 -0.00002816220
## 133 134 135 136
## (Intercept) 0.02919608430 -0.012799174302 0.01695419897 -0.00456404857
## tv -0.00022524016 0.000167341441 -0.00004276786 0.00049170331
## periodico -0.00247706077 0.000947591304 -0.00176944627 -0.00200884389
## radio 0.00090330909 0.000048936950 0.00014621124 0.00205837354
## tv:radio -0.00000435789 -0.000002109586 0.00001305045 -0.00003052679
## 137 138 139 140
## (Intercept) 0.00778594521 -0.022596153125 0.014558760948 -0.02706944785
## tv 0.00023732597 0.000501127883 0.000030221347 0.00043062493
## periodico -0.00212530270 0.001495386134 -0.001362977777 0.00069016250
## radio 0.00155063090 -0.000139327979 0.000563911227 0.00150299208
## tv:radio -0.00001977492 -0.000006599755 -0.000005442513 -0.00003119953
## 141 142 143 144
## (Intercept) 0.020217704051 -0.0082726940 -0.014022424717 0.00168474054
## tv -0.000338761324 0.0001422678 0.000137252931 0.00056054532
## periodico -0.000472916959 0.0007072714 0.001111115544 -0.00007458885
## radio -0.000126687558 -0.0004324584 0.000110590807 -0.00006142670
## tv:radio 0.000005712816 0.0000109961 -0.000004031128 -0.00001374589
## 145 146 147 148
## (Intercept) -0.00198829800 0.03763194326 0.02629464036 -0.01543522339
## tv 0.00060063272 -0.00095187216 -0.00081254970 -0.00036412185
## periodico -0.00031062935 0.00025268846 0.00074360036 0.00173390852
## radio 0.00023599016 -0.00130077426 -0.00107440049 -0.00012045260
## tv:radio -0.00001515365 0.00002809328 0.00002194616 0.00001218465
## 149 150 151 152
## (Intercept) 0.00403254897 0.013338479273 -0.01215692775 -0.01696418130
## tv 0.00024965446 0.000027895532 0.00051576139 0.00122943937
## periodico -0.00181365109 -0.001219077194 0.00104733986 0.00004097211
## radio 0.00149714575 0.000528800275 -0.00010471182 0.00041295892
## tv:radio -0.00001947075 -0.000005332999 -0.00001315815 -0.00003126069
## 153 154 155 156
## (Intercept) 0.001417534822 -0.008470399675 0.007947807307 0.04992101740
## tv -0.000186653453 0.000013929832 -0.000357034184 -0.00074406765
## periodico 0.000739771135 0.000456080437 0.000603260782 -0.00226306936
## radio -0.000026698877 0.000326717452 -0.000156384232 -0.00023230271
## tv:radio -0.000000566054 -0.000001153531 0.000003807189 0.00001812879
## 157 158 159 160
## (Intercept) 0.00778287975 0.016044205762 0.0198020921866 -0.00083849089
## tv -0.00021113564 -0.000010068402 0.0001294379729 0.00037522400
## periodico -0.00077083109 0.000141570665 -0.0026390976690 -0.00007396010
## radio 0.00022254888 -0.000633397205 0.0008090011419 0.00015964400
## tv:radio 0.00001127609 0.000002441521 -0.0000006920975 -0.00001006745
## 161 162 163 164
## (Intercept) -0.001790463491 0.0046878588293 0.000641191193 -0.0113002582
## tv 0.000243976491 0.0001313245628 0.000075684680 0.0001188405
## periodico 0.000371379207 -0.0009227419417 0.000527587196 0.0003582426
## radio 0.000004212692 0.0003693030404 -0.000100996155 0.0009174504
## tv:radio -0.000007238300 0.0000007849008 -0.000003640183 -0.0000166088
## 165 166 167 168
## (Intercept) 0.02587324499 -0.07225833746 0.01468504348 0.016653514998
## tv -0.00069122946 0.00355381404 0.00016799347 -0.000278836159
## periodico -0.00003379496 0.00018304394 -0.00246725793 0.000531485509
## radio -0.00050191766 0.00189664307 0.00127380853 -0.000762066393
## tv:radio 0.00001503164 -0.00009459444 -0.00001194893 0.000008659777
## 169 170 171 172
## (Intercept) -0.0202275264 0.02145995752 0.024433979551 -0.01151738746
## tv 0.0007922543 -0.00082193713 -0.000215338084 0.00067982384
## periodico 0.0007450594 0.00106072925 -0.000878903608 0.00019398865
## radio 0.0001652481 -0.00099540350 -0.000200975390 0.00025915032
## tv:radio -0.0000160912 0.00002053574 0.000005026328 -0.00001576635
## 173 174 175 176
## (Intercept) 0.0235756400961 0.02514494140 0.02775837039 -0.021103947078
## tv -0.0001232980468 -0.00059562607 -0.00066961037 -0.000312707438
## periodico -0.0015784204846 0.00031894679 0.00047486399 0.002179846329
## radio 0.0003275548165 -0.00085657634 -0.00108649690 -0.000111475421
## tv:radio -0.0000003579651 0.00001637958 0.00001978627 0.000009318024
## 177 178 179 180
## (Intercept) -0.01266987581 -0.00181364734 0.012597917565 0.015766611551
## tv 0.00001151397 0.00056422607 -0.000052667635 -0.000302612184
## periodico 0.00133066382 0.00016722575 0.000560014361 0.000378948165
## radio 0.00016996396 -0.00005098358 -0.000713286259 -0.000567858400
## tv:radio -0.00000642050 -0.00001404644 0.000003337299 0.000007731302
## 181 182 183 184
## (Intercept) 0.0375360999 0.006249894866 0.01182661003 -0.01048364001
## tv -0.0009729864 0.000166421440 0.00031343971 -0.00052366780
## periodico 0.0002714128 0.000477573166 -0.00052064220 0.00227758259
## radio -0.0012900129 -0.000444787514 -0.00013754438 -0.00108855870
## tv:radio 0.0000283698 -0.000003254349 -0.00000715736 0.00003015597
## 185 186 187 188
## (Intercept) -0.00878542694 -0.02403900996 0.012112048659 -0.004594159372
## tv 0.00018271141 0.00018060293 0.000133139885 -0.000015963520
## periodico 0.00108395802 0.00115362937 0.000138903045 0.000616178833
## radio -0.00007791804 0.00087000126 -0.000505287750 0.000269956649
## tv:radio -0.00000628598 -0.00001563611 -0.000001664489 -0.000005867028
## 189 190 191 192
## (Intercept) 0.01849130882 0.0234202371995 0.00170359876 0.03410061697
## tv -0.00081441237 0.0000187047846 0.00032730615 -0.00080167299
## periodico 0.00112799685 -0.0013998890588 -0.00188675507 -0.00036176657
## radio -0.00084391921 0.0000586567527 0.00172810634 -0.00068749942
## tv:radio 0.00001843973 -0.0000008807575 -0.00002479561 0.00001995866
## 193 194 195 196
## (Intercept) 0.01595370696 -0.02143366996 -0.00858635415 0.03609474671
## tv 0.00044877998 0.00033114925 0.00008016567 -0.00059753273
## periodico -0.00121850007 0.00049794605 0.00023567381 -0.00062652493
## radio 0.00004599813 0.00134310708 0.00088908809 -0.00080528477
## tv:radio -0.00001048754 -0.00002656952 -0.00001558263 0.00001780196
## 197 198 199 200
## (Intercept) 0.03769835435 0.02808613397 -0.01282058462 0.0246735133
## tv -0.00090291732 -0.00085615088 -0.00037778479 -0.0007715129
## periodico -0.00008300596 0.00052741693 0.00213262423 0.0007233832
## radio -0.00107306292 -0.00096939133 -0.00084263918 -0.0009870685
## tv:radio 0.00002548349 0.00002210282 0.00002311277 0.0000202718n<- nrow(matriz_X_1)
#Matriz M
matriz_M_1<-diag(n)-matriz_P_1#Literal 1/2 Compruebe que los residuos en el objeto “modelo_ventas” son iguales al producto de M*y, donde “y” es la variable endogena en el modelo (“ventas”)
library(magrittr)## Warning: package 'magrittr' was built under R version 4.0.5residuos_modelo_ventas <-modelo_ventas$residuals
datos_modelo<-modelo_ventas$model
residuos_matrices<-matriz_M_1%*%datos_modelo$ventas
cbind(residuos_matrices,residuos_modelo_ventas,residuos_modelo_ventas-residuos_matrices)%>% round(digits = 2) %>% as.data.frame()-> comparacion
names(comparacion)<-c("por_matrices","En_modelo","Diferencia")
comparacion## por_matrices En_modelo Diferencia
## 1 -15.93 -15.93 0
## 2 19.33 19.33 0
## 3 38.02 38.02 0
## 4 -15.43 -15.43 0
## 5 5.16 5.16 0
## 6 80.22 80.22 0
## 7 -16.35 -16.35 0
## 8 -22.89 -22.89 0
## 9 -34.40 -34.40 0
## 10 46.09 46.09 0
## 11 -40.56 -40.56 0
## 12 9.25 9.25 0
## 13 5.82 5.82 0
## 14 -19.64 -19.64 0
## 15 -2.74 -2.74 0
## 16 -20.58 -20.58 0
## 17 -4.89 -4.89 0
## 18 -0.83 -0.83 0
## 19 -36.61 -36.61 0
## 20 -8.11 -8.11 0
## 21 11.00 11.00 0
## 22 56.50 56.50 0
## 23 -7.55 -7.55 0
## 24 31.71 31.71 0
## 25 -40.48 -40.48 0
## 26 86.80 86.80 0
## 27 -1.44 -1.44 0
## 28 35.92 35.92 0
## 29 24.51 24.51 0
## 30 -38.94 -38.94 0
## 31 25.44 25.44 0
## 32 -17.24 -17.24 0
## 33 -44.11 -44.11 0
## 34 47.62 47.62 0
## 35 -38.06 -38.06 0
## 36 104.56 104.56 0
## 37 -22.01 -22.01 0
## 38 6.26 6.26 0
## 39 -20.85 -20.85 0
## 40 -8.14 -8.14 0
## 41 3.52 3.52 0
## 42 7.13 7.13 0
## 43 39.33 39.33 0
## 44 29.22 29.22 0
## 45 -13.40 -13.40 0
## 46 8.21 8.21 0
## 47 -42.24 -42.24 0
## 48 -15.16 -15.16 0
## 49 36.15 36.15 0
## 50 -42.33 -42.33 0
## 51 29.60 29.60 0
## 52 -27.94 -27.94 0
## 53 -25.78 -25.78 0
## 54 -16.15 -16.15 0
## 55 20.52 20.52 0
## 56 -34.10 -34.10 0
## 57 32.84 32.84 0
## 58 -8.84 -8.84 0
## 59 -25.90 -25.90 0
## 60 4.38 4.38 0
## 61 -56.13 -56.13 0
## 62 -5.97 -5.97 0
## 63 34.06 34.06 0
## 64 -22.00 -22.00 0
## 65 -22.92 -22.92 0
## 66 -33.91 -33.91 0
## 67 -25.98 -25.98 0
## 68 -24.10 -24.10 0
## 69 15.18 15.18 0
## 70 -12.96 -12.96 0
## 71 -2.78 -2.78 0
## 72 -39.34 -39.34 0
## 73 9.10 9.10 0
## 74 -23.58 -23.58 0
## 75 16.53 16.53 0
## 76 41.92 41.92 0
## 77 -61.13 -61.13 0
## 78 -11.80 -11.80 0
## 79 42.30 42.30 0
## 80 -28.28 -28.28 0
## 81 -18.10 -18.10 0
## 82 55.90 55.90 0
## 83 -32.84 -32.84 0
## 84 1.99 1.99 0
## 85 -7.85 -7.85 0
## 86 2.06 2.06 0
## 87 -18.80 -18.80 0
## 88 -12.42 -12.42 0
## 89 -34.44 -34.44 0
## 90 -0.71 -0.71 0
## 91 -19.86 -19.86 0
## 92 -71.04 -71.04 0
## 93 5.47 5.47 0
## 94 -1.83 -1.83 0
## 95 -22.68 -22.68 0
## 96 -9.60 -9.60 0
## 97 30.46 30.46 0
## 98 2.80 2.80 0
## 99 -1.39 -1.39 0
## 100 -6.92 -6.92 0
## 101 47.13 47.13 0
## 102 13.67 13.67 0
## 103 74.24 74.24 0
## 104 8.29 8.29 0
## 105 5.20 5.20 0
## 106 -23.16 -23.16 0
## 107 -39.16 -39.16 0
## 108 -36.55 -36.55 0
## 109 -51.39 -51.39 0
## 110 14.93 14.93 0
## 111 31.77 31.77 0
## 112 0.85 0.85 0
## 113 3.55 3.55 0
## 114 20.05 20.05 0
## 115 1.25 1.25 0
## 116 -6.04 -6.04 0
## 117 -5.24 -5.24 0
## 118 -59.49 -59.49 0
## 119 -9.08 -9.08 0
## 120 -15.33 -15.33 0
## 121 -22.74 -22.74 0
## 122 -6.98 -6.98 0
## 123 52.60 52.60 0
## 124 -7.38 -7.38 0
## 125 6.99 6.99 0
## 126 -36.29 -36.29 0
## 127 51.55 51.55 0
## 128 -45.71 -45.71 0
## 129 -38.99 -38.99 0
## 130 -49.74 -49.74 0
## 131 137.77 137.77 0
## 132 68.15 68.15 0
## 133 29.49 29.49 0
## 134 4.04 4.04 0
## 135 2.44 2.44 0
## 136 24.96 24.96 0
## 137 14.83 14.83 0
## 138 18.65 18.65 0
## 139 -13.45 -13.45 0
## 140 -17.13 -17.13 0
## 141 -36.00 -36.00 0
## 142 -7.79 -7.79 0
## 143 -5.46 -5.46 0
## 144 -38.02 -38.02 0
## 145 -33.90 -33.90 0
## 146 -7.02 -7.02 0
## 147 57.02 57.02 0
## 148 -24.42 -24.42 0
## 149 5.76 5.76 0
## 150 -21.33 -21.33 0
## 151 60.96 60.96 0
## 152 -37.40 -37.40 0
## 153 3.70 3.70 0
## 154 -11.53 -11.53 0
## 155 5.57 5.57 0
## 156 20.34 20.34 0
## 157 -6.39 -6.39 0
## 158 0.10 0.10 0
## 159 35.21 35.21 0
## 160 -12.86 -12.86 0
## 161 -0.28 -0.28 0
## 162 -5.49 -5.49 0
## 163 7.07 7.07 0
## 164 -11.48 -11.48 0
## 165 -17.07 -17.07 0
## 166 43.48 43.48 0
## 167 30.49 30.49 0
## 168 32.76 32.76 0
## 169 9.33 9.33 0
## 170 79.25 79.25 0
## 171 -32.15 -32.15 0
## 172 -2.32 -2.32 0
## 173 -19.07 -19.07 0
## 174 11.45 11.45 0
## 175 56.81 56.81 0
## 176 -20.90 -20.90 0
## 177 10.72 10.72 0
## 178 10.72 10.72 0
## 179 99.16 99.16 0
## 180 0.68 0.68 0
## 181 8.07 8.07 0
## 182 43.24 43.24 0
## 183 -53.93 -53.93 0
## 184 -14.41 -14.41 0
## 185 33.06 33.06 0
## 186 -26.80 -26.80 0
## 187 -11.78 -11.78 0
## 188 2.02 2.02 0
## 189 75.55 75.55 0
## 190 -31.29 -31.29 0
## 191 11.46 11.46 0
## 192 -34.64 -34.64 0
## 193 -47.21 -47.21 0
## 194 -20.94 -20.94 0
## 195 -16.29 -16.29 0
## 196 -54.55 -54.55 0
## 197 -31.98 -31.98 0
## 198 8.34 8.34 0
## 199 -9.58 -9.58 0
## 200 49.59 49.59 0#Literal 1/3 Muestre que los autovalores de x’x son positivos (use el comando eigen)
eigen( x= matriz_XX_1,symmetric = TRUE)->descomposicion
auto_valores<- descomposicion$values
print(auto_valores)## [1] 311421698.6388 70252.5341 40973.4590 3714.3627 12.7735print(auto_valores>0)## [1] TRUE TRUE TRUE TRUE TRUE#EJERICIO 2 #Literal 1/4 Estime el modelo propuesto, y coloquele el nombre de “modelo_cajas”
load("C:/Users/HP/Desktop/REBE UNIVERSIDAD/ECONOMETRIA/GUIA ECONOMETRIA/EJERICIO N2/datos_cajas.RData")
modelo_cajas<-lm(formula = Tiempo~Distancia+N_cajas,data = datos_cajas)
Summary=(modelo_cajas)#Literal 2/4 Calcule las matrices A, P, M
options(scipen = 99999999)
matriz_X_2<-model.matrix(modelo_cajas)
matriz_XX_2<-t(matriz_X_2)%*%matriz_X_2
#Matriz A
matriz_A_2<-solve(matriz_XX_2)%*%t(matriz_X_2)
#Matriz P
matriz_P_2<-matriz_X_2%*%matriz_A_2
print(matriz_P_2)## 1 2 3 4 5 6
## 1 0.19781478 0.127154573 0.16766180 0.062524965 -0.03527291 0.057620774
## 2 0.12715457 0.124295239 0.03396629 0.140073563 0.05334477 0.038710181
## 3 0.16766180 0.033966286 0.35585795 -0.137368460 -0.10168744 0.123125512
## 4 0.06252497 0.140073563 -0.13736846 0.257600846 0.15524536 0.006698639
## 5 -0.03527291 0.053344771 -0.10168744 0.155245361 0.18408997 0.046742309
## 6 0.05762077 0.038710181 0.12312551 0.006698639 0.04674231 0.086318088
## 7 0.17558129 0.144648497 0.07654437 0.133523089 0.01345706 0.036955589
## 8 0.11716423 0.050316476 0.21126231 -0.035350897 -0.01751039 0.094896089
## 9 0.09794605 0.077129229 0.10132526 0.055636570 0.03786105 0.067680430
## 10 -0.02906036 -0.056765574 0.20436525 -0.131155907 0.05122193 0.136694350
## 11 -0.01209498 0.081873124 -0.13140718 0.199703669 0.18629079 0.030873007
## 12 0.09285990 0.104513848 0.01812731 0.131114317 0.07550894 0.044246890
## 13 0.15541865 0.125438973 0.08744449 0.109054124 0.01789770 0.046274418
## 14 -0.12402490 -0.005427535 -0.12246527 0.112857904 0.23285894 0.067134558
## 15 -0.05129385 -0.039271650 0.11324781 -0.060157783 0.09995191 0.116029165
## 7 8 9 10 11 12
## 1 0.17558129 0.11716423 0.09794605 -0.02906036 -0.01209498 0.092859897
## 2 0.14464850 0.05031648 0.07712923 -0.05676557 0.08187312 0.104513848
## 3 0.07654437 0.21126231 0.10132526 0.20436525 -0.13140718 0.018127310
## 4 0.13352309 -0.03535090 0.05563657 -0.13115591 0.19970367 0.131114317
## 5 0.01345706 -0.01751039 0.03786105 0.05122193 0.18629079 0.075508940
## 6 0.03695559 0.09489609 0.06768043 0.13669435 0.03087301 0.044246890
## 7 0.18301556 0.07160552 0.08894348 -0.08682757 0.04935470 0.112467995
## 8 0.07160552 0.13896449 0.08399596 0.13551596 -0.03237026 0.042396988
## 9 0.08894348 0.08399596 0.07465547 0.05440619 0.04101064 0.069478345
## 10 -0.08682757 0.13551596 0.05440619 0.34795579 -0.01326471 -0.021162536
## 11 0.04935470 -0.03237026 0.04101064 -0.01326471 0.20329083 0.095597926
## 12 0.11246799 0.04239699 0.06947834 -0.02116254 0.09559793 0.094228911
## 13 0.15702161 0.07705558 0.08545596 -0.04568349 0.04428588 0.099852268
## 14 -0.07689788 -0.02789930 0.01907176 0.16357209 0.20867158 0.042323339
## 15 -0.07939330 0.08995724 0.04540362 0.29018859 0.04818497 -0.001554438
## 13 14 15
## 1 0.15541865 -0.124024902 -0.051293849
## 2 0.12543897 -0.005427535 -0.039271650
## 3 0.08744449 -0.122465266 0.113247813
## 4 0.10905412 0.112857904 -0.060157783
## 5 0.01789770 0.232858944 0.099951911
## 6 0.04627442 0.067134558 0.116029165
## 7 0.15702161 -0.076897883 -0.079393301
## 8 0.07705558 -0.027899299 0.089957240
## 9 0.08545596 0.019071756 0.045403621
## 10 -0.04568349 0.163572088 0.290188586
## 11 0.04428588 0.208671580 0.048184973
## 12 0.09985227 0.042323339 -0.001554438
## 13 0.13743085 -0.052866482 -0.044080529
## 14 -0.05286648 0.352392093 0.210699107
## 15 -0.04408053 0.210699107 0.262089133n<- nrow(matriz_X_2)
#Matriz M
matriz_M_2<-diag(n)-matriz_P_2
print(matriz_M_2)## 1 2 3 4 5 6
## 1 0.80218522 -0.127154573 -0.16766180 -0.062524965 0.03527291 -0.057620774
## 2 -0.12715457 0.875704761 -0.03396629 -0.140073563 -0.05334477 -0.038710181
## 3 -0.16766180 -0.033966286 0.64414205 0.137368460 0.10168744 -0.123125512
## 4 -0.06252497 -0.140073563 0.13736846 0.742399154 -0.15524536 -0.006698639
## 5 0.03527291 -0.053344771 0.10168744 -0.155245361 0.81591003 -0.046742309
## 6 -0.05762077 -0.038710181 -0.12312551 -0.006698639 -0.04674231 0.913681912
## 7 -0.17558129 -0.144648497 -0.07654437 -0.133523089 -0.01345706 -0.036955589
## 8 -0.11716423 -0.050316476 -0.21126231 0.035350897 0.01751039 -0.094896089
## 9 -0.09794605 -0.077129229 -0.10132526 -0.055636570 -0.03786105 -0.067680430
## 10 0.02906036 0.056765574 -0.20436525 0.131155907 -0.05122193 -0.136694350
## 11 0.01209498 -0.081873124 0.13140718 -0.199703669 -0.18629079 -0.030873007
## 12 -0.09285990 -0.104513848 -0.01812731 -0.131114317 -0.07550894 -0.044246890
## 13 -0.15541865 -0.125438973 -0.08744449 -0.109054124 -0.01789770 -0.046274418
## 14 0.12402490 0.005427535 0.12246527 -0.112857904 -0.23285894 -0.067134558
## 15 0.05129385 0.039271650 -0.11324781 0.060157783 -0.09995191 -0.116029165
## 7 8 9 10 11 12
## 1 -0.17558129 -0.11716423 -0.09794605 0.02906036 0.01209498 -0.092859897
## 2 -0.14464850 -0.05031648 -0.07712923 0.05676557 -0.08187312 -0.104513848
## 3 -0.07654437 -0.21126231 -0.10132526 -0.20436525 0.13140718 -0.018127310
## 4 -0.13352309 0.03535090 -0.05563657 0.13115591 -0.19970367 -0.131114317
## 5 -0.01345706 0.01751039 -0.03786105 -0.05122193 -0.18629079 -0.075508940
## 6 -0.03695559 -0.09489609 -0.06768043 -0.13669435 -0.03087301 -0.044246890
## 7 0.81698444 -0.07160552 -0.08894348 0.08682757 -0.04935470 -0.112467995
## 8 -0.07160552 0.86103551 -0.08399596 -0.13551596 0.03237026 -0.042396988
## 9 -0.08894348 -0.08399596 0.92534453 -0.05440619 -0.04101064 -0.069478345
## 10 0.08682757 -0.13551596 -0.05440619 0.65204421 0.01326471 0.021162536
## 11 -0.04935470 0.03237026 -0.04101064 0.01326471 0.79670917 -0.095597926
## 12 -0.11246799 -0.04239699 -0.06947834 0.02116254 -0.09559793 0.905771089
## 13 -0.15702161 -0.07705558 -0.08545596 0.04568349 -0.04428588 -0.099852268
## 14 0.07689788 0.02789930 -0.01907176 -0.16357209 -0.20867158 -0.042323339
## 15 0.07939330 -0.08995724 -0.04540362 -0.29018859 -0.04818497 0.001554438
## 13 14 15
## 1 -0.15541865 0.124024902 0.051293849
## 2 -0.12543897 0.005427535 0.039271650
## 3 -0.08744449 0.122465266 -0.113247813
## 4 -0.10905412 -0.112857904 0.060157783
## 5 -0.01789770 -0.232858944 -0.099951911
## 6 -0.04627442 -0.067134558 -0.116029165
## 7 -0.15702161 0.076897883 0.079393301
## 8 -0.07705558 0.027899299 -0.089957240
## 9 -0.08545596 -0.019071756 -0.045403621
## 10 0.04568349 -0.163572088 -0.290188586
## 11 -0.04428588 -0.208671580 -0.048184973
## 12 -0.09985227 -0.042323339 0.001554438
## 13 0.86256915 0.052866482 0.044080529
## 14 0.05286648 0.647607907 -0.210699107
## 15 0.04408053 -0.210699107 0.737910867#Literal 3/4 Compruebe que los residuos en el objeto “modelo_ventas” son iguales al producto de M*y, donde “y” es la variable endogena en el modelo (“Tiempo”)
library(magrittr)
residuos_modelo<-modelo_cajas$residuals
matriz_y<-datos_cajas$Tiempo
residuos_por_matrices<-matriz_M_2%*%matriz_y
comparativa<-cbind(residuos_modelo,
residuos_por_matrices,
residuos_modelo-residuos_por_matrices) %>% as.data.frame() #%>% round(digits = 2)
names(comparativa)<-c("residuos_modelo","Residuos_matrices","diferencia")
print(comparativa)## residuos_modelo Residuos_matrices diferencia
## 1 -0.7608713 -0.7608713 -0.0000000000002959855
## 2 0.1327095 0.1327095 -0.0000000000003170797
## 3 -0.3200790 -0.3200790 -0.0000000000002652323
## 4 2.9381318 2.9381318 -0.0000000000003237410
## 5 -9.2715743 -9.2715743 -0.0000000000002877698
## 6 0.7655710 0.7655710 -0.0000000000002698952
## 7 1.3084025 1.3084025 -0.0000000000003268497
## 8 -2.0933728 -2.0933728 -0.0000000000002775558
## 9 1.4318218 1.4318218 -0.0000000000002895462
## 10 0.5212280 0.5212280 -0.0000000000002299272
## 11 0.5174718 0.5174718 -0.0000000000002975398
## 12 1.3783003 1.3783003 -0.0000000000003024248
## 13 -1.0247228 -1.0247228 -0.0000000000003108624
## 14 2.8864818 2.8864818 -0.0000000000002602363
## 15 1.5905018 1.5905018 -0.0000000000002442491#Literal 4/4 Muestre que los autovalores de x’x son positivos (use el comando eigen)
eigen(x=matriz_XX_2,symmetric = TRUE)->descomposicion
auto_valores<-descomposicion$values
print(auto_valores)## [1] 16976.7781334 709.9345923 0.2872743print(auto_valores>0)## [1] TRUE TRUE TRUE#EJERCICIO 3
load("C:/Users/HP/Desktop/REBE UNIVERSIDAD/ECONOMETRIA/GUIA ECONOMETRIA/EJERICIO N3/modelo_estimado.RData")
model.matrix(modelo_estimado_1)## (Intercept) poverty single
## 1 1 9.1 14.3
## 2 1 17.4 11.5
## 3 1 20.0 10.7
## 4 1 15.4 12.1
## 5 1 18.2 12.5
## 6 1 9.9 12.1
## 7 1 8.5 10.1
## 8 1 10.2 11.4
## 9 1 17.8 10.6
## 10 1 13.5 13.0
## 11 1 8.0 9.1
## 12 1 10.3 9.0
## 13 1 13.1 9.5
## 14 1 13.6 11.5
## 15 1 12.2 10.8
## 16 1 13.1 9.9
## 17 1 20.4 10.6
## 18 1 26.4 14.9
## 19 1 10.7 10.9
## 20 1 9.7 12.0
## 21 1 10.7 10.6
## 22 1 15.4 13.0
## 23 1 11.6 9.9
## 24 1 16.1 10.9
## 25 1 24.7 14.7
## 26 1 14.9 10.8
## 27 1 14.4 11.1
## 28 1 11.2 8.4
## 29 1 10.3 9.4
## 30 1 9.9 9.2
## 31 1 10.9 9.6
## 32 1 17.4 13.8
## 33 1 9.8 12.4
## 34 1 16.4 12.7
## 35 1 13.0 11.4
## 36 1 19.9 11.1
## 37 1 11.8 11.3
## 38 1 13.2 9.6
## 39 1 11.2 10.8
## 40 1 18.7 12.3
## 41 1 14.2 9.4
## 42 1 19.6 11.2
## 43 1 17.4 11.8
## 44 1 10.7 10.0
## 45 1 9.7 10.3
## 46 1 10.0 11.0
## 47 1 12.1 11.7
## 48 1 12.6 10.4
## 49 1 22.2 9.4
## 50 1 13.3 10.8
## 51 1 26.4 22.1
## attr(,"assign")
## [1] 0 1 2summary(modelo_estimado_1)##
## Call:
## lm(formula = crime ~ poverty + single, data = datos_regresion)
##
## Residuals:
## Min 1Q Median 3Q Max
## -811.14 -114.27 -22.44 121.86 689.82
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -1368.189 187.205 -7.308 0.0000000024786 ***
## poverty 6.787 8.989 0.755 0.454
## single 166.373 19.423 8.566 0.0000000000312 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 243.6 on 48 degrees of freedom
## Multiple R-squared: 0.7072, Adjusted R-squared: 0.695
## F-statistic: 57.96 on 2 and 48 DF, p-value: 0.0000000000001578#Lireral 1/3 Calcule las matrices A, P, M
matriz_X_3<-model.matrix(modelo_estimado_1)
matriz_XX_3<-t(matriz_X_3)%*%matriz_X_3
options(scipen = 99999999)
#Matriz A
matriz_A_3<-solve(matriz_XX_3)%*%t(matriz_X_3)
#Matriz P
matriz_P_3<-matriz_X_3%*%matriz_A_3
print(matriz_P_3)## 1 2 3 4 5
## 1 0.161611077 -0.0127796253 -0.0653080904 0.027207907 0.004995312
## 2 -0.012779625 0.0314650769 0.0450195150 0.021099510 0.030700881
## 3 -0.065308090 0.0450195150 0.0785589509 0.019423663 0.038838489
## 4 0.027207907 0.0210995096 0.0194236628 0.022341212 0.024424173
## 5 0.004995312 0.0307008811 0.0388384889 0.024424173 0.034582636
## 6 0.092237760 -0.0008717607 -0.0291165550 0.020671038 0.005338821
## 7 0.054323132 0.0014555968 -0.0149893805 0.014082445 -0.001730668
## 8 0.069626967 0.0030986911 -0.0171998930 0.018604926 0.005605850
## 9 -0.042019509 0.0366270037 0.0604669926 0.018447422 0.031093776
## 10 0.074183205 0.0099454181 -0.0092620561 0.024537798 0.018826182
## 11 0.033001099 0.0034182300 -0.0061606961 0.010848882 -0.004571400
## 12 0.003083400 0.0130022202 0.0154620904 0.011239145 0.003299178
## 13 -0.016405789 0.0222075842 0.0335527599 0.013630280 0.013568208
## 14 0.032150082 0.0162849295 0.0114826440 0.019945572 0.017514640
## 15 0.029639461 0.0134642545 0.0083959467 0.017363227 0.011882559
## 16 -0.005512265 0.0206235772 0.0282561924 0.014862970 0.014010488
## 17 -0.072760898 0.0470134239 0.0834132808 0.019236959 0.040115944
## 18 -0.026597072 0.0539539932 0.0794281003 0.034310399 0.065690838
## 19 0.050098243 0.0070760900 -0.0061664301 0.017215896 0.006788032
## 20 0.091879085 -0.0012747116 -0.0295575040 0.020302131 0.004534238
## 21 0.041928113 0.0082640934 -0.0021940108 0.016291379 0.006456323
## 22 0.051718351 0.0175354918 0.0075063795 0.025114767 0.025419303
## 23 0.012223150 0.0146314125 0.0150179511 0.014407468 0.008805392
## 24 -0.013749255 0.0286478850 0.0414912433 0.018855703 0.025526380
## 25 -0.011943712 0.0479548814 0.0670730541 0.033177818 0.059570593
## 26 -0.002284283 0.0242501501 0.0322247792 0.018183130 0.021251731
## 27 0.011797679 0.0210647547 0.0238396001 0.018955815 0.019848409
## 28 -0.023898157 0.0189735319 0.0313498957 0.009663407 0.005758814
## 29 0.013976925 0.0114182132 0.0101655228 0.012471835 0.003741458
## 30 0.013259614 0.0106123038 0.0092836039 0.011734023 0.002132290
## 31 0.012329555 0.0130230694 0.0128125179 0.013270384 0.006044635
## 32 0.049858208 0.0223570267 0.0145642200 0.028187487 0.033243991
## 33 0.101590245 -0.0024592395 -0.0339715188 0.021565188 0.005323525
## 34 0.031724585 0.0227182754 0.0203043056 0.024493916 0.028557657
## 35 0.036520858 0.0142840658 0.0075114923 0.019455197 0.015322030
## 36 -0.053232174 0.0430360250 0.0723798179 0.020625989 0.038933762
## 37 0.047985819 0.0098863347 -0.0017549634 0.018782624 0.011047384
## 38 -0.014864749 0.0222110559 0.0331111561 0.013968820 0.014025783
## 39 0.041463070 0.0094694781 -0.0004295475 0.017059559 0.008412495
## 40 -0.006363255 0.0334902728 0.0458995198 0.023959662 0.036096529
## 41 -0.032135147 0.0269978397 0.0445849466 0.013656140 0.017274706
## 42 -0.046961732 0.0414415952 0.0684080408 0.020843060 0.038003315
## 43 -0.004609469 0.0302770697 0.0410470830 0.022024029 0.031032591
## 44 0.025587799 0.0106401078 0.0057508532 0.014442341 0.005792902
## 45 0.045581566 0.0054573242 -0.0070470730 0.015063192 0.002654548
## 46 0.061098162 0.0038837436 -0.0136684226 0.017311502 0.004469558
## 47 0.055332259 0.0095007613 -0.0044038810 0.020106415 0.012530684
## 48 0.014016459 0.0166461782 0.0172227295 0.016252001 0.012828306
## 49 -0.126724035 0.0589560548 0.1151889082 0.016085484 0.045035222
## 50 0.016633486 0.0178585101 0.0181039936 0.017697261 0.015699630
## 51 0.169486584 0.0254418356 -0.0159102166 0.056498849 0.073651878
## 6 7 8 9 10
## 1 0.0922377599 0.054323132 0.069626967 -0.042019509 0.074183205
## 2 -0.0008717607 0.001455597 0.003098691 0.036627004 0.009945418
## 3 -0.0291165550 -0.014989380 -0.017199893 0.060466993 -0.009262056
## 4 0.0206710383 0.014082445 0.018604926 0.018447422 0.024537798
## 5 0.0053388207 -0.001730668 0.005605850 0.031093776 0.018826182
## 6 0.0601833060 0.046325314 0.049657601 -0.014507419 0.045082769
## 7 0.0463253143 0.051525360 0.043514080 -0.002242589 0.026574141
## 8 0.0496576009 0.043514080 0.043047070 -0.005481207 0.035653735
## 9 -0.0145074185 -0.002242589 -0.005481207 0.048317713 -0.002230971
## 10 0.0450827691 0.026574141 0.035653735 -0.002230971 0.042317361
## 11 0.0379595101 0.052952915 0.039313133 0.005088182 0.016572738
## 12 0.0202404094 0.039319171 0.025625073 0.019901969 0.006794286
## 13 0.0061039860 0.023656425 0.013328397 0.031515647 0.002269591
## 14 0.0264276188 0.023732484 0.024553262 0.013858209 0.024140123
## 15 0.0281148133 0.030887233 0.027540841 0.012698365 0.021061366
## 16 0.0108871129 0.024257870 0.016137961 0.027384984 0.007017238
## 17 -0.0331859478 -0.017484675 -0.020160655 0.063896367 -0.011943140
## 18 -0.0248715863 -0.046193174 -0.023833459 0.055442446 0.016681405
## 19 0.0400866623 0.039831103 0.036712140 0.002678021 0.027851446
## 20 0.0604243269 0.047347419 0.050084394 -0.014673105 0.044642941
## 21 0.0364993229 0.039380020 0.034604970 0.005776013 0.024290717
## 22 0.0314330794 0.015435698 0.024926450 0.009153426 0.035220009
## 23 0.0216631859 0.033051380 0.024606872 0.018397301 0.012620412
## 24 0.0012927954 0.008174465 0.006224056 0.035033686 0.007680018
## 25 -0.0150502752 -0.036527926 -0.015640148 0.047321742 0.020657841
## 26 0.0087178833 0.015058917 0.012296802 0.028876195 0.010975654
## 27 0.0158972585 0.018441171 0.017226949 0.022782299 0.016404119
## 28 0.0066000667 0.033140898 0.016329375 0.031490581 -0.003689098
## 29 0.0250235363 0.039920616 0.028434638 0.015771307 0.011541933
## 30 0.0255055964 0.041964833 0.029288235 0.015439919 0.010662291
## 31 0.0231046860 0.036703939 0.026451865 0.017301036 0.011674500
## 32 0.0266312472 0.004913909 0.019253703 0.012875670 0.037244415
## 33 0.0644890461 0.047362628 0.052329362 -0.018204586 0.049017041
## 34 0.0206616799 0.009122273 0.017173332 0.018243217 0.027923819
## 35 0.0295422633 0.027099528 0.027238434 0.011295804 0.025194476
## 36 -0.0236150092 -0.013801698 -0.013825726 0.055737139 -0.004140852
## 37 0.0369673443 0.033983974 0.033311174 0.005138319 0.028490109
## 38 0.0065813727 0.023220557 0.013466200 0.031082152 0.003082966
## 39 0.0352988620 0.036749573 0.033186782 0.006706576 0.024796816
## 40 -0.0006447671 -0.004962560 0.001378097 0.036155002 0.014584634
## 41 -0.0029942508 0.017057493 0.006415471 0.039139282 -0.003026318
## 42 -0.0202640241 -0.011892640 -0.011429560 0.052906946 -0.001833314
## 43 0.0027155901 0.001906681 0.005205868 0.033529002 0.013506159
## 44 0.0293246212 0.038477851 0.030390617 0.011972017 0.017169236
## 45 0.0400960207 0.044791276 0.038143735 0.002882226 0.024465426
## 46 0.0463112824 0.044085102 0.041366692 -0.002548901 0.031653178
## 47 0.0395952552 0.032826716 0.034426956 0.002805195 0.032117120
## 48 0.0204580515 0.027940847 0.022472891 0.019225757 0.014819527
## 49 -0.0604666469 -0.029841231 -0.038752062 0.087073601 -0.032909916
## 50 0.0202123570 0.024438657 0.021330304 0.019289336 0.016952371
## 51 0.0612247882 -0.035367152 0.026738755 -0.018909558 0.102139135
## 11 12 13 14 15
## 1 0.033001099 0.003083400 -0.0164057895 0.032150082 0.0296394607
## 2 0.003418230 0.013002220 0.0222075842 0.016284930 0.0134642545
## 3 -0.006160696 0.015462090 0.0335527599 0.011482644 0.0083959467
## 4 0.010848882 0.011239145 0.0136302797 0.019945572 0.0173632265
## 5 -0.004571400 0.003299178 0.0135682081 0.017514640 0.0118825586
## 6 0.037959510 0.020240409 0.0061039860 0.026427619 0.0281148133
## 7 0.052952915 0.039319171 0.0236564253 0.023732484 0.0308872331
## 8 0.039313133 0.025625073 0.0133283967 0.024553262 0.0275408413
## 9 0.005088182 0.019901969 0.0315156470 0.013858209 0.0126983655
## 10 0.016572738 0.006794286 0.0022695914 0.024140123 0.0210613665
## 11 0.059463778 0.048531233 0.0327063281 0.022149205 0.0318824765
## 12 0.048531233 0.045606447 0.0368271023 0.019221275 0.0273881407
## 13 0.032706328 0.036827102 0.0357916146 0.017007601 0.0219056920
## 14 0.022149205 0.019221275 0.0170076006 0.020763434 0.0208926221
## 15 0.031882477 0.027388141 0.0219056920 0.020892622 0.0236418693
## 16 0.031087826 0.033469599 0.0318979738 0.017876622 0.0218985614
## 17 -0.007727752 0.015646825 0.0350735318 0.010793969 0.0076157965
## 18 -0.054701898 -0.030265921 0.0014273628 0.013064630 -0.0041898611
## 19 0.038871660 0.029003659 0.0188796609 0.022877707 0.0265723376
## 20 0.039349978 0.021407108 0.0068037177 0.026446073 0.0285075624
## 21 0.040085535 0.031521783 0.0217998869 0.022225942 0.0265776855
## 22 0.007207250 0.003684759 0.0048695833 0.021900870 0.0173471827
## 23 0.038481633 0.035924490 0.0298453483 0.019644453 0.0248308123
## 24 0.012253950 0.020166053 0.0262691136 0.016513515 0.0160162331
## 25 -0.045513004 -0.025804963 0.0010478759 0.014633660 -0.0008630804
## 26 0.018573624 0.022969338 0.0256004177 0.017710526 0.0183638179
## 27 0.019824348 0.021269504 0.0219959740 0.018951570 0.0193358870
## 28 0.046522708 0.049169776 0.0438991476 0.016857042 0.0256394868
## 29 0.046912731 0.042248944 0.0329334615 0.020090296 0.0273810101
## 30 0.049693667 0.044582334 0.0343329143 0.020127207 0.0281665101
## 31 0.043145955 0.039588230 0.0318076813 0.019817677 0.0262045455
## 32 -0.005888169 -0.006303442 -0.0001808736 0.021281807 0.0134232536
## 33 0.037238553 0.017885945 0.0030469192 0.027197238 0.0283049477
## 34 0.003491924 0.004566297 0.0091582355 0.020070549 0.0153976967
## 35 0.025511358 0.021042611 0.0171599647 0.021253310 0.0220673058
## 36 -0.007286280 0.012268239 0.0295222676 0.012469523 0.0085843002
## 37 0.031831027 0.023845895 0.0164912701 0.022450321 0.0244148889
## 38 0.031808783 0.035824064 0.0349550407 0.017107002 0.0217084271
## 39 0.036811682 0.029024734 0.0205372750 0.022071176 0.0255967032
## 40 -0.006226751 0.004159633 0.0161992371 0.016490852 0.0109087069
## 41 0.027688833 0.035866230 0.0382702874 0.015493935 0.0197571584
## 42 -0.006212146 0.011919844 0.0281383380 0.013040343 0.0091689663
## 43 0.002204351 0.010484089 0.0192873490 0.016936696 0.0134589066
## 44 0.042513293 0.036558045 0.0276403573 0.020922409 0.0265883814
## 45 0.046228619 0.035676507 0.0233517051 0.022752730 0.0285378674
## 46 0.041917475 0.029309895 0.0169483544 0.023919951 0.0279389383
## 47 0.028733762 0.019997413 0.0130081546 0.022965776 0.0238213077
## 48 0.031529298 0.030091013 0.0263467097 0.019552176 0.0228670652
## 49 -0.011744812 0.022773480 0.0492176249 0.006065502 0.0041184849
## 50 0.026460349 0.025587887 0.0234109512 0.019596212 0.0214915512
## 51 -0.083834974 -0.090701039 -0.0686582458 0.028707029 -0.0043182118
## 16 17 18 19 20
## 1 -0.005512265 -0.072760898 -0.026597072 0.050098243 0.091879085
## 2 0.020623577 0.047013424 0.053953993 0.007076090 -0.001274712
## 3 0.028256192 0.083413281 0.079428100 -0.006166430 -0.029557504
## 4 0.014862970 0.019236959 0.034310399 0.017215896 0.020302131
## 5 0.014010488 0.040115944 0.065690838 0.006788032 0.004534238
## 6 0.010887113 -0.033185948 -0.024871586 0.040086662 0.060424327
## 7 0.024257870 -0.017484675 -0.046193174 0.039831103 0.047347419
## 8 0.016137961 -0.020160655 -0.023833459 0.036712140 0.050084394
## 9 0.027384984 0.063896367 0.055442446 0.002678021 -0.014673105
## 10 0.007017238 -0.011943140 0.016681405 0.027851446 0.044642941
## 11 0.031087826 -0.007727752 -0.054701898 0.038871660 0.039349978
## 12 0.033469599 0.015646825 -0.030265921 0.029003659 0.021407108
## 13 0.031897974 0.035073532 0.001427363 0.018879661 0.006803718
## 14 0.017876622 0.010793969 0.013064630 0.022877707 0.026446073
## 15 0.021898561 0.007615796 -0.004189861 0.026572338 0.028507562
## 16 0.029021418 0.029264481 0.002678779 0.020095101 0.011461679
## 17 0.029264481 0.088678108 0.083419388 -0.008133616 -0.033639964
## 18 0.002678779 0.083419388 0.161434287 -0.020017550 -0.027336512
## 19 0.020095101 -0.008133616 -0.020017550 0.033113678 0.040614464
## 20 0.011461679 -0.033639964 -0.027336512 0.040614464 0.060718819
## 21 0.022252515 -0.003776834 -0.020956111 0.032202099 0.037120999
## 22 0.008390716 0.006166587 0.037126086 0.019950638 0.030782549
## 23 0.027937093 0.014967325 -0.013461762 0.026332583 0.022404096
## 24 0.023998670 0.043336153 0.038088405 0.010658741 0.001221758
## 25 0.002888156 0.070120480 0.142515978 -0.013556129 -0.017264095
## 26 0.023850346 0.033350676 0.024863111 0.015344871 0.008811213
## 27 0.021331484 0.024228162 0.020421494 0.018335613 0.015952162
## 28 0.038435035 0.032938706 -0.022458728 0.023438008 0.007854707
## 29 0.030593043 0.009837774 -0.029014505 0.030219100 0.026065069
## 30 0.031742168 0.008929719 -0.033944364 0.031274710 0.026654070
## 31 0.029588488 0.012652092 -0.021932584 0.028331832 0.024017098
## 32 0.004083363 0.013611347 0.061149643 0.014064879 0.025508592
## 33 0.008657411 -0.038495867 -0.025009056 0.041414071 0.064647282
## 34 0.011271019 0.020054820 0.046947884 0.014880735 0.019994128
## 35 0.018162034 0.006527372 0.006295554 0.025068840 0.029658547
## 36 0.025307341 0.076651069 0.079603479 -0.004535153 -0.024170035
## 37 0.018013710 -0.003458095 -0.006929730 0.029754966 0.037247994
## 38 0.031251119 0.034574400 0.002816248 0.018767693 0.007238723
## 39 0.021175678 -0.001915641 -0.014950222 0.030730659 0.035802507
## 40 0.015810207 0.047786188 0.070445310 0.004101151 -0.001442215
## 41 0.033412287 0.047010377 0.012950898 0.014001648 -0.002385213
## 42 0.024371341 0.072339389 0.076688232 -0.002983801 -0.020817072
## 43 0.018466157 0.042656629 0.054892557 0.007987672 0.002218765
## 44 0.026567356 0.004936756 -0.022833238 0.030378936 0.030134046
## 45 0.023687052 -0.008951477 -0.032655035 0.035448839 0.040922467
## 46 0.018869941 -0.016257890 -0.027236945 0.036328364 0.046885420
## 47 0.015354019 -0.006407713 -0.002450203 0.029722909 0.039717471
## 48 0.025064278 0.017237442 -0.001137130 0.023693564 0.020931608
## 49 0.039195353 0.123261887 0.099033795 -0.019264925 -0.060744778
## 50 0.022693733 0.018100381 0.007646540 0.021998183 0.020483120
## 51 -0.049099283 -0.021143638 0.183959802 0.001860398 0.056506876
## 21 22 23 24 25
## 1 0.0419281127 0.051718351 0.012223150 -0.013749255 -0.0119437117
## 2 0.0082640934 0.017535492 0.014631413 0.028647885 0.0479548814
## 3 -0.0021940108 0.007506379 0.015017951 0.041491243 0.0670730541
## 4 0.0162913794 0.025114767 0.014407468 0.018855703 0.0331778182
## 5 0.0064563227 0.025419303 0.008805392 0.025526380 0.0595705928
## 6 0.0364993229 0.031433079 0.021663186 0.001292795 -0.0150502752
## 7 0.0393800202 0.015435698 0.033051380 0.008174465 -0.0365279257
## 8 0.0346049704 0.024926450 0.024606872 0.006224056 -0.0156401479
## 9 0.0057760130 0.009153426 0.018397301 0.035033686 0.0473217421
## 10 0.0242907171 0.035220009 0.012620412 0.007680018 0.0206578411
## 11 0.0400855353 0.007207250 0.038481633 0.012253950 -0.0455130038
## 12 0.0315217825 0.003684759 0.035924490 0.020166053 -0.0258049627
## 13 0.0217998869 0.004869583 0.029845348 0.026269114 0.0010478759
## 14 0.0222259425 0.021900870 0.019644453 0.016513515 0.0146336599
## 15 0.0265776855 0.017347183 0.024830812 0.016016233 -0.0008630804
## 16 0.0222525147 0.008390716 0.027937093 0.023998670 0.0028881565
## 17 -0.0037768344 0.006166587 0.014967325 0.043336153 0.0701204805
## 18 -0.0209561105 0.037126086 -0.013461762 0.038088405 0.1425159781
## 19 0.0322020992 0.019950638 0.026332583 0.010658741 -0.0135561288
## 20 0.0371209990 0.030782549 0.022404096 0.001221758 -0.0172640947
## 21 0.0318626289 0.017309793 0.027763772 0.012361571 -0.0149363370
## 22 0.0173097926 0.033037319 0.010113892 0.013747203 0.0373184518
## 23 0.0277637720 0.010113892 0.029915925 0.019208787 -0.0102649599
## 24 0.0123615708 0.013747203 0.019208787 0.027902325 0.0337950952
## 25 -0.0149363370 0.037318452 -0.010264960 0.033795095 0.1266889826
## 26 0.0166574230 0.014245466 0.021268915 0.024638023 0.0228125274
## 27 0.0188339802 0.017460712 0.020497332 0.021338560 0.0198083680
## 28 0.0275360867 -0.002630853 0.037599578 0.026445652 -0.0206735215
## 29 0.0319744104 0.007205892 0.034016235 0.017895610 -0.0239646821
## 30 0.0332177672 0.005904840 0.035498051 0.017753527 -0.0283923251
## 31 0.0299962246 0.008277196 0.032270570 0.018676335 -0.0177832958
## 32 0.0108667063 0.037782025 0.003658934 0.015592823 0.0585365059
## 33 0.0372062083 0.034188802 0.020363918 -0.000729358 -0.0145469363
## 34 0.0132961496 0.029247682 0.010225864 0.018643294 0.0447069834
## 35 0.0243172893 0.021709854 0.020913052 0.015165174 0.0089123377
## 36 -0.0013739633 0.011142399 0.013241614 0.038901468 0.0680364614
## 37 0.0286131381 0.022208116 0.022973179 0.011900875 -0.0020702217
## 38 0.0215456293 0.005634993 0.029236361 0.026020824 0.0023848176
## 39 0.0302518570 0.018495967 0.026150034 0.012822977 -0.0096318246
## 40 0.0043929230 0.023084345 0.009099908 0.028258230 0.0630348246
## 41 0.0176451427 0.002725634 0.028871272 0.030349307 0.0102334172
## 42 -0.0001585582 0.012367312 0.013160318 0.037375885 0.0658659127
## 43 0.0086035648 0.020176345 0.013200219 0.026945050 0.0493350940
## 44 0.0311836860 0.012028085 0.030626159 0.015767241 -0.0176967624
## 45 0.0351973290 0.015817723 0.030514187 0.010871150 -0.0250852940
## 46 0.0348871761 0.021635074 0.026778972 0.007855849 -0.0192341757
## 47 0.0279635137 0.025384613 0.020669157 0.010588409 0.0024006838
## 48 0.0246553858 0.013366528 0.026211383 0.019563986 0.0008041373
## 49 -0.0117482331 -0.006464639 0.018317499 0.055895357 0.0803833795
## 50 0.0225360954 0.016083520 0.023379668 0.019528816 0.0087825417
## 51 -0.0128088006 0.100506538 -0.047810391 -0.002779615 0.1756410645
## 26 27 28 29 30
## 1 -0.002284283 0.011797679 -0.023898157 0.013976925 0.013259614
## 2 0.024250150 0.021064755 0.018973532 0.011418213 0.010612304
## 3 0.032224779 0.023839600 0.031349896 0.010165523 0.009283604
## 4 0.018183130 0.018955815 0.009663407 0.012471835 0.011734023
## 5 0.021251731 0.019848409 0.005758814 0.003741458 0.002132290
## 6 0.008717883 0.015897258 0.006600067 0.025023536 0.025505596
## 7 0.015058917 0.018441171 0.033140898 0.039920616 0.041964833
## 8 0.012296802 0.017226949 0.016329375 0.028434638 0.029288235
## 9 0.028876195 0.022782299 0.031490581 0.015771307 0.015439919
## 10 0.010975654 0.016404119 -0.003689098 0.011541933 0.010662291
## 11 0.018573624 0.019824348 0.046522708 0.046912731 0.049693667
## 12 0.022969338 0.021269504 0.049169776 0.042248944 0.044582334
## 13 0.025600418 0.021995974 0.043899148 0.032933462 0.034332914
## 14 0.017710526 0.018951570 0.016857042 0.020090296 0.020127207
## 15 0.018363818 0.019335887 0.025639487 0.027381010 0.028166510
## 16 0.023850346 0.021331484 0.038435035 0.030593043 0.031742168
## 17 0.033350676 0.024228162 0.032938706 0.009837774 0.008929719
## 18 0.024863111 0.020421494 -0.022458728 -0.029014505 -0.033944364
## 19 0.015344871 0.018335613 0.023438008 0.030219100 0.031274710
## 20 0.008811213 0.015952162 0.007854707 0.026065069 0.026654070
## 21 0.016657423 0.018833980 0.027536087 0.031974410 0.033217767
## 22 0.014245466 0.017460712 -0.002630853 0.007205892 0.005904840
## 23 0.021268915 0.020497332 0.037599578 0.034016235 0.035498051
## 24 0.024638023 0.021338560 0.026445652 0.017895610 0.017753527
## 25 0.022812527 0.019808368 -0.020673522 -0.023964682 -0.028392325
## 26 0.023010393 0.020837360 0.027143308 0.021219266 0.021405920
## 27 0.020837360 0.020060941 0.022766732 0.020605014 0.020714818
## 28 0.027143308 0.022766732 0.057867232 0.043705663 0.046214931
## 29 0.021219266 0.020605014 0.043705663 0.039908526 0.041991588
## 30 0.021405920 0.020714818 0.046214931 0.041991588 0.044288531
## 31 0.021376798 0.020606427 0.041307776 0.037369036 0.039193856
## 32 0.014187226 0.017243932 -0.012445150 -0.002039205 -0.004284505
## 33 0.007233237 0.015343282 0.002446292 0.023496437 0.023812930
## 34 0.017278976 0.018515181 0.002024209 0.006679081 0.005343980
## 35 0.017115473 0.018784033 0.017888894 0.022044680 0.022277252
## 36 0.030302607 0.023119498 0.025830072 0.008053313 0.006943244
## 37 0.015487845 0.018282833 0.018586551 0.025368335 0.025929642
## 38 0.025334992 0.021885460 0.042588810 0.032120143 0.033434834
## 39 0.016642864 0.018779786 0.025082516 0.029663138 0.030670433
## 40 0.022987244 0.020458705 0.008769356 0.003770603 0.002175702
## 41 0.027930986 0.022773808 0.045877850 0.031008230 0.032226291
## 42 0.029348806 0.022786546 0.024296960 0.008152847 0.007046736
## 43 0.022937594 0.020566386 0.014875441 0.009662897 0.008669241
## 44 0.019282535 0.019830717 0.035732269 0.035485043 0.037103893
## 45 0.016249025 0.018776247 0.031077207 0.036011854 0.037664753
## 46 0.013702684 0.017780219 0.021682094 0.031231481 0.032379766
## 47 0.014254059 0.017785173 0.013289529 0.022343278 0.022587718
## 48 0.020802277 0.020222820 0.031326401 0.028808582 0.029755693
## 49 0.041698619 0.027222619 0.050333619 0.012751208 0.012194908
## 50 0.020256868 0.019947599 0.026252155 0.024870669 0.025412194
## 51 -0.006638218 0.008460658 -0.120812865 -0.071142077 -0.080577853
## 31 32 33 34 35
## 1 0.0123295548 0.0498582078 0.101590245 0.031724585 0.0365208582
## 2 0.0130230694 0.0223570267 -0.002459239 0.022718275 0.0142840658
## 3 0.0128125179 0.0145642200 -0.033971519 0.020304306 0.0075114923
## 4 0.0132703840 0.0281874867 0.021565188 0.024493916 0.0194551966
## 5 0.0060446349 0.0332439912 0.005323525 0.028557657 0.0153220299
## 6 0.0231046860 0.0266312472 0.064489046 0.020661680 0.0295422633
## 7 0.0367039393 0.0049139095 0.047362628 0.009122273 0.0270995277
## 8 0.0264518654 0.0192537035 0.052329362 0.017173332 0.0272384342
## 9 0.0173010357 0.0128756699 -0.018204586 0.018243217 0.0112958043
## 10 0.0116745000 0.0372444146 0.049017041 0.027923819 0.0251944763
## 11 0.0431459553 -0.0058881689 0.037238553 0.003491924 0.0255113578
## 12 0.0395882297 -0.0063034421 0.017885945 0.004566297 0.0210426110
## 13 0.0318076813 -0.0001808736 0.003046919 0.009158236 0.0171599647
## 14 0.0198176767 0.0212818069 0.027197238 0.020070549 0.0212533105
## 15 0.0262045455 0.0134232536 0.028304948 0.015397697 0.0220673058
## 16 0.0295884880 0.0040833632 0.008657411 0.011271019 0.0181620341
## 17 0.0126520923 0.0136113469 -0.038495867 0.020054820 0.0065273720
## 18 -0.0219325838 0.0611496435 -0.025009056 0.046947884 0.0062955538
## 19 0.0283318322 0.0140648787 0.041414071 0.014880735 0.0250688402
## 20 0.0240170976 0.0255085924 0.064647282 0.019994128 0.0296585468
## 21 0.0299962246 0.0108667063 0.037206208 0.013296150 0.0243172893
## 22 0.0082771964 0.0377820245 0.034188802 0.029247682 0.0217098540
## 23 0.0322705703 0.0036589342 0.020363918 0.010225864 0.0209130523
## 24 0.0186763347 0.0155928232 -0.000729358 0.018643294 0.0151651736
## 25 -0.0177832958 0.0585365059 -0.014546936 0.044706983 0.0089123377
## 26 0.0213767976 0.0141872257 0.007233237 0.017278976 0.0171154734
## 27 0.0206064274 0.0172439320 0.015343282 0.018515181 0.0187840328
## 28 0.0413077762 -0.0124451501 0.002446292 0.002024209 0.0178888943
## 29 0.0373690363 -0.0020392053 0.023496437 0.006679081 0.0220446804
## 30 0.0391938560 -0.0042845049 0.023812930 0.005343980 0.0222772516
## 31 0.0351866025 0.0002626947 0.021619098 0.008153539 0.0214453113
## 32 0.0002626947 0.0468764135 0.029801125 0.034866795 0.0200459709
## 33 0.0216190981 0.0298011246 0.069477338 0.022176589 0.0304772143
## 34 0.0081535391 0.0348667945 0.022176589 0.028359862 0.0191242886
## 35 0.0214453113 0.0200459709 0.030477214 0.019124289 0.0221031999
## 36 0.0107721254 0.0188001716 -0.027580576 0.022347417 0.0086969660
## 37 0.0241457725 0.0186403737 0.038439802 0.017759971 0.0240534980
## 38 0.0310740758 0.0009134858 0.003669120 0.009756111 0.0172270831
## 39 0.0279926005 0.0131403010 0.036109285 0.014700927 0.0239013180
## 40 0.0062602041 0.0312533492 -0.001383890 0.027849650 0.0139039891
## 41 0.0303956229 -0.0009356901 -0.006940478 0.009396484 0.0148920338
## 42 0.0107537449 0.0197813401 -0.023836665 0.022666580 0.0094976844
## 43 0.0113586718 0.0255552094 0.001748636 0.024302866 0.0150356190
## 44 0.0333250198 0.0044703410 0.028790456 0.010126969 0.0228141828
## 45 0.0334486771 0.0073855710 0.040802670 0.011014789 0.0253997482
## 46 0.0290286694 0.0149328762 0.048279736 0.014921194 0.0266031669
## 47 0.0213901624 0.0229894963 0.041708991 0.020081786 0.0245053634
## 48 0.0277085212 0.0092721879 0.019572702 0.013563616 0.0203316281
## 49 0.0160911819 0.0013279314 -0.069375188 0.014970646 0.0002199351
## 50 0.0242376845 0.0137345016 0.019720173 0.016164144 0.0200498918
## 51 -0.0618781056 0.1379059865 0.075979912 0.084978031 0.0243328225
## 36 37 38 39 40
## 1 -0.05323217408 0.04798581916 -0.0148647495 0.0414630703 -0.0063632555
## 2 0.04303602499 0.00988633473 0.0222110559 0.0094694781 0.0334902728
## 3 0.07237981793 -0.00175496339 0.0331111561 -0.0004295475 0.0458995198
## 4 0.02062598918 0.01878262392 0.0139688204 0.0170595586 0.0239596621
## 5 0.03893376181 0.01104738429 0.0140257830 0.0084124945 0.0360965286
## 6 -0.02361500919 0.03696734429 0.0065813727 0.0352988620 -0.0006447671
## 7 -0.01380169777 0.03398397378 0.0232205566 0.0367495729 -0.0049625601
## 8 -0.01382572557 0.03331117441 0.0134662003 0.0331867822 0.0013780971
## 9 0.05573713908 0.00513831910 0.0310821519 0.0067065763 0.0361550017
## 10 -0.00414085174 0.02849010874 0.0030829660 0.0247968158 0.0145846344
## 11 -0.00728628002 0.03183102660 0.0318087828 0.0368116815 -0.0062267510
## 12 0.01226823934 0.02384589473 0.0358240642 0.0290247343 0.0041596327
## 13 0.02952226761 0.01649127006 0.0349550407 0.0205372750 0.0161992371
## 14 0.01246952316 0.02245032083 0.0171070022 0.0220711761 0.0164908525
## 15 0.00858430019 0.02441488885 0.0217084271 0.0255967032 0.0109087069
## 16 0.02530734095 0.01801371004 0.0312511193 0.0211756783 0.0158102073
## 17 0.07665106865 -0.00345809543 0.0345744004 -0.0019156411 0.0477861881
## 18 0.07960347922 -0.00692972971 0.0028162483 -0.0149502220 0.0704453099
## 19 -0.00453515328 0.02975496623 0.0187676927 0.0307306588 0.0041011511
## 20 -0.02417003455 0.03724799441 0.0072387233 0.0358025071 -0.0014422148
## 21 -0.00137396331 0.02861313806 0.0215456293 0.0302518570 0.0043929230
## 22 0.01114239918 0.02220811569 0.0056349928 0.0184959668 0.0230843446
## 23 0.01324161412 0.02297317924 0.0292363608 0.0261500338 0.0090999084
## 24 0.03890146791 0.01190087524 0.0260208239 0.0128229771 0.0282582295
## 25 0.06803646135 -0.00207022173 0.0023848176 -0.0096318246 0.0630348246
## 26 0.03030260695 0.01548784494 0.0253349920 0.0166428640 0.0229872440
## 27 0.02311949799 0.01828283313 0.0218854604 0.0187797857 0.0204587049
## 28 0.02583007240 0.01858655088 0.0425888096 0.0250825157 0.0087693560
## 29 0.00805331268 0.02536833471 0.0321201428 0.0296631376 0.0037706029
## 30 0.00694324425 0.02592964173 0.0334348338 0.0306704326 0.0021757022
## 31 0.01077212543 0.02414577255 0.0310740758 0.0279926005 0.0062602041
## 32 0.01880017158 0.01864037369 0.0009134858 0.0131403010 0.0312533492
## 33 -0.02758057635 0.03843980185 0.0036691196 0.0361092855 -0.0013838897
## 34 0.02234741708 0.01775997109 0.0097561106 0.0147009265 0.0278496499
## 35 0.00869696605 0.02405349796 0.0172270831 0.0239013180 0.0139039891
## 36 0.06736049636 0.00009811276 0.0292729081 0.0005404823 0.0450631341
## 37 0.00009811276 0.02764046450 0.0165412525 0.0277212016 0.0086330079
## 38 0.02927290810 0.01654125249 0.0341633724 0.0203652548 0.0165493298
## 39 0.00054048231 0.02772120165 0.0203652548 0.0289129402 0.0064351744
## 40 0.04506313408 0.00863300790 0.0165493298 0.0064351744 0.0385278098
## 41 0.03942419937 0.01247371607 0.0373585142 0.0167298146 0.0212173782
## 42 0.06389363049 0.00147061226 0.0279439815 0.0016949509 0.0436238207
## 43 0.03987482497 0.01102816653 0.0194331105 0.0099482813 0.0331985000
## 44 0.00494843673 0.02632947446 0.0271015202 0.0292942505 0.0049764686
## 45 -0.00625658117 0.03077761905 0.0229804025 0.0330892908 0.0002111632
## 46 -0.01121956095 0.03244999636 0.0169014875 0.0332116256 0.0008724212
## 47 -0.00170366700 0.02817101001 0.0132402830 0.0273647332 0.0095860387
## 48 0.01601676866 0.02156991823 0.0259496270 0.0236318017 0.0130871533
## 49 0.10377475013 -0.01397678944 0.0481038940 -0.0098000845 0.0570056431
## 50 0.01743250294 0.02077794352 0.0231859171 0.0219488413 0.0158295945
## 51 0.00373471899 0.02047421898 -0.0638544075 -0.0034589501 0.0634427661
## 41 42 43 44 45
## 1 -0.0321351475 -0.04696173199 -0.004609469 0.025587799 0.0455815656
## 2 0.0269978397 0.04144159523 0.030277070 0.010640108 0.0054573242
## 3 0.0445849466 0.06840804084 0.041047083 0.005750853 -0.0070470730
## 4 0.0136561402 0.02084306015 0.022024029 0.014442341 0.0150631921
## 5 0.0172747063 0.03800331459 0.031032591 0.005792902 0.0026545484
## 6 -0.0029942508 -0.02026402406 0.002715590 0.029324621 0.0400960207
## 7 0.0170574929 -0.01189263976 0.001906681 0.038477851 0.0447912755
## 8 0.0064154706 -0.01142955985 0.005205868 0.030390617 0.0381437345
## 9 0.0391392822 0.05290694621 0.033529002 0.011972017 0.0028822256
## 10 -0.0030263180 -0.00183331380 0.013506159 0.017169236 0.0244654258
## 11 0.0276888329 -0.00621214597 0.002204351 0.042513293 0.0462286191
## 12 0.0358662301 0.01191984441 0.010484089 0.036558045 0.0356765073
## 13 0.0382702874 0.02813833798 0.019287349 0.027640357 0.0233517051
## 14 0.0154939355 0.01304034275 0.016936696 0.020922409 0.0227527296
## 15 0.0197571584 0.00916896625 0.013458907 0.026588381 0.0285378674
## 16 0.0334122873 0.02437134100 0.018466157 0.026567356 0.0236870520
## 17 0.0470103772 0.07233938887 0.042656629 0.004936756 -0.0089514769
## 18 0.0129508981 0.07668823171 0.054892557 -0.022833238 -0.0326550350
## 19 0.0140016484 -0.00298380147 0.007987672 0.030378936 0.0354488388
## 20 -0.0023852131 -0.02081707196 0.002218765 0.030134046 0.0409224674
## 21 0.0176451427 -0.00015855822 0.008603565 0.031183686 0.0351973290
## 22 0.0027256341 0.01236731244 0.020176345 0.012028085 0.0158177231
## 23 0.0288712715 0.01316031803 0.013200219 0.030626159 0.0305141871
## 24 0.0303493068 0.03737588549 0.026945050 0.015767241 0.0108711499
## 25 0.0102334172 0.06586591273 0.049335094 -0.017696762 -0.0250852940
## 26 0.0279309862 0.02934880617 0.022937594 0.019282535 0.0162490251
## 27 0.0227738083 0.02278654628 0.020566386 0.019830717 0.0187762474
## 28 0.0458778502 0.02429695980 0.014875441 0.035732269 0.0310772069
## 29 0.0310082300 0.00815284742 0.009662897 0.035485043 0.0360118543
## 30 0.0322262908 0.00704673551 0.008669241 0.037103893 0.0376647528
## 31 0.0303956229 0.01075374486 0.011358672 0.033325020 0.0334486771
## 32 -0.0009356901 0.01978134011 0.025555209 0.004470341 0.0073855710
## 33 -0.0069404777 -0.02383666457 0.001748636 0.028790456 0.0408026702
## 34 0.0093964838 0.02266657998 0.024302866 0.010126969 0.0110147891
## 35 0.0148920338 0.00949768445 0.015035619 0.022814183 0.0253997482
## 36 0.0394241994 0.06389363049 0.039874825 0.004948437 -0.0062565812
## 37 0.0124737161 0.00147061226 0.011028167 0.026329474 0.0307776191
## 38 0.0373585142 0.02794398150 0.019433110 0.027101520 0.0229804025
## 39 0.0167298146 0.00169495094 0.009948281 0.029294251 0.0330892908
## 40 0.0212173782 0.04362382074 0.033198500 0.004976469 0.0002111632
## 41 0.0428148697 0.03730150429 0.023354334 0.024932155 0.0182613048
## 42 0.0373015043 0.06070968685 0.038616343 0.005491946 -0.0048073213
## 43 0.0233543338 0.03861634300 0.029661175 0.009835355 0.0057088348
## 44 0.0249321545 0.00549194624 0.009835355 0.032793191 0.0346943078
## 45 0.0182613048 -0.00480732130 0.005708835 0.034694308 0.0394972418
## 46 0.0106680025 -0.00915736450 0.005324936 0.032004792 0.0387186714
## 47 0.0085239197 -0.00005417871 0.011260163 0.024444711 0.0297475381
## 48 0.0258261094 0.01592558262 0.015684353 0.026579037 0.0263819477
## 49 0.0670336232 0.09709363390 0.051439339 0.003285199 -0.0181500866
## 50 0.0230872378 0.01739038595 0.017320596 0.023611924 0.0235312999
## 51 -0.0744931970 0.00888221413 0.040111081 -0.042147291 -0.0266187841
## 46 47 48 49 50
## 1 0.0610981620 0.05533225892 0.0140164593 -0.1267240353 0.016633486
## 2 0.0038837436 0.00950076134 0.0166461782 0.0589560548 0.017858510
## 3 -0.0136684226 -0.00440388101 0.0172227295 0.1151889082 0.018103994
## 4 0.0173115022 0.02010641476 0.0162520006 0.0160854844 0.017697261
## 5 0.0044695580 0.01253068375 0.0128283056 0.0450352224 0.015699630
## 6 0.0463112824 0.03959525517 0.0204580515 -0.0604666469 0.020212357
## 7 0.0440851020 0.03282671578 0.0279408473 -0.0298412311 0.024438657
## 8 0.0413666923 0.03442695554 0.0224728906 -0.0387520617 0.021330304
## 9 -0.0025489007 0.00280519452 0.0192257569 0.0870736011 0.019289336
## 10 0.0316531775 0.03211711988 0.0148195267 -0.0329099159 0.016952371
## 11 0.0419174754 0.02873376166 0.0315292980 -0.0117448120 0.026460349
## 12 0.0293098947 0.01999741335 0.0300910133 0.0227734801 0.025587887
## 13 0.0169483544 0.01300815463 0.0263467097 0.0492176249 0.023410951
## 14 0.0239199513 0.02296577564 0.0195521764 0.0060655020 0.019596212
## 15 0.0279389383 0.02382130771 0.0228670652 0.0041184849 0.021491551
## 16 0.0188699407 0.01535401915 0.0250642779 0.0391953530 0.022693733
## 17 -0.0162578898 -0.00640771321 0.0172374416 0.1232618867 0.018100381
## 18 -0.0272369450 -0.00245020324 -0.0011371298 0.0990337955 0.007646540
## 19 0.0363283635 0.02972290933 0.0236935635 -0.0192649252 0.021998183
## 20 0.0468854200 0.03971747067 0.0209316081 -0.0607447782 0.020483120
## 21 0.0348871761 0.02796351374 0.0246553858 -0.0117482331 0.022536095
## 22 0.0216350737 0.02538461273 0.0133665276 -0.0064646395 0.016083520
## 23 0.0267789716 0.02066915744 0.0262113827 0.0183174990 0.023379668
## 24 0.0078558491 0.01058840945 0.0195639858 0.0558953572 0.019528816
## 25 -0.0192341757 0.00240068385 0.0008041373 0.0803833795 0.008782542
## 26 0.0137026836 0.01425405946 0.0208022767 0.0416986194 0.020256868
## 27 0.0177802192 0.01778517341 0.0202228196 0.0272226188 0.019947599
## 28 0.0216820941 0.01328952935 0.0313264015 0.0503336190 0.026252155
## 29 0.0312314810 0.02234327786 0.0288085816 0.0127512081 0.024870669
## 30 0.0323797658 0.02258771784 0.0297556925 0.0121949084 0.025412194
## 31 0.0290286694 0.02139016243 0.0277085212 0.0160911819 0.024237684
## 32 0.0149328762 0.02298949630 0.0092721879 0.0013279314 0.013734502
## 33 0.0482797356 0.04170899129 0.0195727024 -0.0693751879 0.019720173
## 34 0.0149211944 0.02008178601 0.0135636164 0.0149706458 0.016164144
## 35 0.0266031669 0.02450536339 0.0203316281 0.0002199351 0.020049892
## 36 -0.0112195610 -0.00170366700 0.0160167687 0.1037747501 0.017432503
## 37 0.0324499964 0.02817101001 0.0215699182 -0.0139767894 0.020777944
## 38 0.0169014875 0.01324028304 0.0259496270 0.0481038940 0.023185917
## 39 0.0332116256 0.02736473324 0.0236318017 -0.0098000845 0.021948841
## 40 0.0008724212 0.00958603873 0.0130871533 0.0570056431 0.015829594
## 41 0.0106680025 0.00852391966 0.0258261094 0.0670336232 0.023087238
## 42 -0.0091573645 -0.00005417871 0.0159255826 0.0970936339 0.017390386
## 43 0.0053249356 0.01126016252 0.0156843529 0.0514393389 0.017320596
## 44 0.0320047921 0.02444471137 0.0265790365 0.0032851987 0.023611924
## 45 0.0387186714 0.02974753809 0.0263819477 -0.0181500866 0.023531300
## 46 0.0404996425 0.03278977545 0.0239082695 -0.0315135010 0.022138980
## 47 0.0327897755 0.02945384619 0.0200580654 -0.0198234879 0.019923538
## 48 0.0239082695 0.02005806536 0.0238436049 0.0197082164 0.022025855
## 49 -0.0315135010 -0.01982348795 0.0197082164 0.1783821911 0.019428916
## 50 0.0221389802 0.01992353828 0.0220258547 0.0194289164 0.020988532
## 51 0.0073516449 0.03977540282 -0.0242209261 -0.0813672909 -0.005263400
## 51
## 1 0.169486584
## 2 0.025441836
## 3 -0.015910217
## 4 0.056498849
## 5 0.073651878
## 6 0.061224788
## 7 -0.035367152
## 8 0.026738755
## 9 -0.018909558
## 10 0.102139135
## 11 -0.083834974
## 12 -0.090701039
## 13 -0.068658246
## 14 0.028707029
## 15 -0.004318212
## 16 -0.049099283
## 17 -0.021143638
## 18 0.183959802
## 19 0.001860398
## 20 0.056506876
## 21 -0.012808801
## 22 0.100506538
## 23 -0.047810391
## 24 -0.002779615
## 25 0.175641065
## 26 -0.006638218
## 27 0.008460658
## 28 -0.120812865
## 29 -0.071142077
## 30 -0.080577853
## 31 -0.061878106
## 32 0.137905986
## 33 0.075979912
## 34 0.084978031
## 35 0.024332822
## 36 0.003734719
## 37 0.020474219
## 38 -0.063854408
## 39 -0.003458950
## 40 0.063442766
## 41 -0.074493197
## 42 0.008882214
## 43 0.040111081
## 44 -0.042147291
## 45 -0.026618784
## 46 0.007351645
## 47 0.039775403
## 48 -0.024220926
## 49 -0.081367291
## 50 -0.005263400
## 51 0.536021500n<- nrow(matriz_X_3)
#Matriz M
matriz_M_3<-diag(n)-matriz_P_3#Literal 2/3Compruebe que los residuos en el objeto “modelo_estimado” son iguales al producto de M*y, donde “y” es la variable endogena en el modelo (“crime”)
library(magrittr)
residuos_modelo<-modelo_estimado_1$residuals
datos_modelo<-modelo_estimado_1$model
residuos_matrices<-matriz_M_3%*%datos_modelo$crime
cbind(residuos_matrices,residuos_modelo,residuos_modelo-residuos_matrices)%>% round(digits = 2) %>% as.data.frame()-> comparacion
names(comparacion)<-c("por matrices", "en modelo", "diferencia")
print(comparacion)## por matrices en modelo diferencia
## 1 -311.71 -311.71 0
## 2 116.80 116.80 0
## 3 45.25 45.25 0
## 4 -34.45 -34.45 0
## 5 243.00 243.00 0
## 6 -145.12 -145.12 0
## 7 86.13 86.13 0
## 8 88.31 88.31 0
## 9 689.82 689.82 0
## 10 -163.29 -163.29 0
## 11 60.90 60.90 0
## 12 126.92 126.92 0
## 13 -19.27 -19.27 0
## 14 322.59 322.59 0
## 15 -22.44 -22.44 0
## 16 128.18 128.18 0
## 17 -70.82 -70.82 0
## 18 -227.95 -227.95 0
## 19 287.10 287.10 0
## 20 303.88 303.88 0
## 21 -341.99 -341.99 0
## 22 -107.18 -107.18 0
## 23 -30.63 -30.63 0
## 24 189.45 189.45 0
## 25 -811.14 -811.14 0
## 26 -351.77 -351.77 0
## 27 102.71 102.71 0
## 28 -23.36 -23.36 0
## 29 73.38 73.38 0
## 30 -91.63 -91.63 0
## 31 324.03 324.03 0
## 32 -115.85 -115.85 0
## 33 113.65 113.65 0
## 34 217.94 217.94 0
## 35 -112.70 -112.70 0
## 36 21.38 21.38 0
## 37 -88.91 -88.91 0
## 38 99.42 99.42 0
## 39 -102.65 -102.65 0
## 40 217.88 217.88 0
## 41 -84.09 -84.09 0
## 42 137.78 137.78 0
## 43 48.89 48.89 0
## 44 -67.16 -67.16 0
## 45 -39.29 -39.29 0
## 46 -415.78 -415.78 0
## 47 -145.50 -145.50 0
## 48 -183.61 -183.61 0
## 49 -138.39 -138.39 0
## 50 -232.91 -232.91 0
## 51 434.17 434.17 0#Literal 3/3 Muestre que los autovalores de x’x son positivos (use el comando eigen)
eigen(x=matriz_XX_3, symmetric = TRUE)->descomposicion
auto_valores<-descomposicion$values
print(auto_valores)## [1] 17956.580914 279.157317 1.681762print(auto_valores>0)## [1] TRUE TRUE TRUE#EJERCICIO 4
#Literal 1/3 Estima la ecuacion de inversion, presenta sus resultados en formato APA
library(readxl)
Investiment_Equation <- read_excel("C:/Users/HP/Desktop/REBE UNIVERSIDAD/ECONOMETRIA/GUIA ECONOMETRIA/EJERICIO N4/Investiment_Equation.xlsx")
ecuacion_inversion<-lm(formula = InvReal~Trend+Inflation+PNBr+Interest,data = Investiment_Equation)
library(stargazer)
stargazer(ecuacion_inversion,title="Ecuacion de Inversion",type="text")##
## Ecuacion de Inversion
## ===============================================
## Dependent variable:
## ---------------------------
## InvReal
## -----------------------------------------------
## Trend -0.016***
## (0.002)
##
## Inflation 0.00002
## (0.001)
##
## PNBr 0.665***
## (0.054)
##
## Interest -0.240*
## (0.120)
##
## Constant -0.503***
## (0.054)
##
## -----------------------------------------------
## Observations 15
## R2 0.973
## Adjusted R2 0.962
## Residual Std. Error 0.007 (df = 10)
## F Statistic 90.089*** (df = 4; 10)
## ===============================================
## Note: *p<0.1; **p<0.05; ***p<0.01#Literal 2/3 Calcule los residuos a traves de la matriz M
model.matrix(ecuacion_inversion)->matriz_X_4
n<-nrow(matriz_X_4)
matriz_M_4<-diag(n)-matriz_X_4%*%solve(t(matriz_X_4)%*%matriz_X_4)%*%t(matriz_X_4)
y<-Investiment_Equation$InvReal
residuos<-matriz_M_4%*%y
print(residuos)## [,1]
## 1 -0.0100602233
## 2 -0.0009290882
## 3 0.0029656679
## 4 0.0078576839
## 5 0.0028109133
## 6 0.0006259732
## 7 0.0075909286
## 8 -0.0055352778
## 9 -0.0037254127
## 10 0.0006953129
## 11 0.0019904770
## 12 -0.0001288433
## 13 -0.0101976729
## 14 0.0068712384
## 15 -0.0008316770#Literal 3/3 Calcule un intervalo de confianza del 93% para el impacto del PNBr en la Inversion, e interpretelo
confint(object = ecuacion_inversion,parm = "PNBr", level = .93)## 3.5 % 96.5 %
## PNBr 0.554777 0.774317Interpretacion:podemos concluir que con el 93% de las ocasiones que estimasemos la ecuacion e esperaria que el impacto de 1 millon de dolares del PNBr se tradujese en el minimo punto de 0.55millones de dolarees en inversion real hasta un maximo de 0.77 millones de dolares en inversion real.
#EJERCICIO 5 #Literal 1/4 Calcule los residuos del modelo
load("C:/Users/HP/Desktop/REBE UNIVERSIDAD/ECONOMETRIA/GUIA ECONOMETRIA/EJERICIO N5/consumption_equation.RData")
n<-nrow(P)
matriz_M_5<-diag(n)-P
residuos<-matriz_M_5%*%C
print(residuos)## [,1]
## 1 -5.859103
## 2 2.605057
## 3 45.765735
## 4 31.102448
## 5 -21.037889
## 6 7.008120
## 7 17.859663
## 8 10.705631
## 9 22.002328
## 10 -2.689665
## 11 7.784083
## 12 -13.127696
## 13 17.521565
## 14 17.304695
## 15 -16.308260
## 16 -5.255508
## 17 2.788211
## 18 -16.379339
## 19 -14.327554
## 20 11.749135
## 21 -31.424669
## 22 -23.329596
## 23 22.171806
## 24 -5.040038
## 25 -36.191398
## 26 -25.211753
## 27 -21.411271
## 28 1.410519
## 29 -24.229564
## 30 20.971808
## 31 43.342653
## 32 36.808458
## 33 17.882297
## 34 -33.100273
## 35 -37.819995
## 36 -49.370820
## 37 23.456143
## 38 -25.510341
## 39 -11.960629
## 40 -9.234201
## 41 21.949616
## 42 3.211123
## 43 -14.511436
## 44 3.197576
## 45 -62.396763
## 46 -66.854500
## 47 8.330745
## 48 91.963380
## 49 61.620735
## 50 48.148861
## 51 -10.717721
## 52 -84.069717
## 53 -56.426627
## 54 125.113605#Literal 2/4 calculo de la varianza del error
#Sigma cuadrado
k<-4
var_error=t(residuos)%*%residuos/(n-k)
print(var_error)## [,1]
## [1,] 1428.746#Literal 3/4 Matriz de Var-Cov del modelo
var_error<-as.vector(var_error)
var.cov<-var_error*solve(XX)
print(var.cov)## (Intercept) Yd W I
## (Intercept) 164.522304918 -0.09333539523 0.009670913575 10.5186890800
## Yd -0.093335395 0.00018911268 -0.000032769561 -0.0072901023
## W 0.009670914 -0.00003276956 0.000006165749 0.0004193421
## I 10.518689080 -0.00729010228 0.000419342092 5.3203789879#Literal 4/4 Calculo de las estimaciones
C_estimada<-P%*%C
cuadro<-as.data.frame(cbind(C,C_estimada,residuos))
names(cuadro)<-c("C","C_estimada","Residuos")
print(head(cuadro))## C C_estimada Residuos
## 1 976.4 982.2591 -5.859103
## 2 998.1 995.4949 2.605057
## 3 1025.3 979.5343 45.765735
## 4 1090.9 1059.7976 31.102448
## 5 1107.1 1128.1379 -21.037889
## 6 1142.4 1135.3919 7.008120#EJERICICIO 6 #Literal 1/3 Estima la ecuacion de ventas, presenta sus resultados en formato APA.
load("C:/Users/HP/Desktop/REBE UNIVERSIDAD/ECONOMETRIA/GUIA ECONOMETRIA/EJERICIO N6/datos_ventas.RData")
ecuacion_ventas<-lm(formula = ventas~tv+radio+periodico, data = datos_ventas)
library(stargazer)
stargazer(ecuacion_ventas, title = "ecuacion de ventas", type = "text")##
## ecuacion de ventas
## ===============================================
## Dependent variable:
## ---------------------------
## ventas
## -----------------------------------------------
## tv 0.045
## (0.118)
##
## radio -3.450***
## (0.206)
##
## periodico 18.485***
## (0.563)
##
## Constant -33.289***
## (7.172)
##
## -----------------------------------------------
## Observations 200
## R2 0.847
## Adjusted R2 0.844
## Residual Std. Error 33.875 (df = 196)
## F Statistic 360.758*** (df = 3; 196)
## ===============================================
## Note: *p<0.1; **p<0.05; ***p<0.01#Literal 2/3 Calcule los residuos a traves de la matriz M
model.matrix(ecuacion_ventas)->matriz_X_6
n<-nrow(matriz_X_6)
matriz_M_6<-diag(n)-matriz_X_6%*%solve(t(matriz_X_6)%*%matriz_X_6)%*%t(matriz_X_6)
y_6<-datos_ventas$ventas
residuos_6<-matriz_M_6%*%y_6
print(residuos_6)## [,1]
## 1 -17.8524638
## 2 19.0821552
## 3 33.7931916
## 4 -17.3508987
## 5 10.2572135
## 6 74.2038531
## 7 -15.2465204
## 8 -23.4242974
## 9 -39.6405207
## 10 45.1613878
## 11 -40.6649446
## 12 8.9590198
## 13 5.1392043
## 14 -22.6236420
## 15 -2.4070766
## 16 -23.2129522
## 17 -8.8531474
## 18 -2.2553254
## 19 -36.4998224
## 20 -7.7078500
## 21 12.1068249
## 22 56.1602254
## 23 -4.4128052
## 24 32.1887265
## 25 -41.0755512
## 26 85.5635292
## 27 -0.5813419
## 28 36.0534051
## 29 25.1740669
## 30 -36.8482228
## 31 26.2868044
## 32 -15.4987679
## 33 -43.1444263
## 34 46.2272722
## 35 -42.1231189
## 36 101.1406281
## 37 -18.4643544
## 38 4.6103727
## 39 -19.7868698
## 40 -7.5311456
## 41 4.4398079
## 42 7.6680340
## 43 39.7209714
## 44 29.5192706
## 45 -12.0296830
## 46 9.1587903
## 47 -40.4095721
## 48 -13.3392963
## 49 39.1663986
## 50 -40.4133097
## 51 31.4947034
## 52 -31.1473697
## 53 -26.0085689
## 54 -19.2664951
## 55 21.2230654
## 56 -38.1989298
## 57 33.9900112
## 58 -9.0293035
## 59 -26.4532617
## 60 5.2080674
## 61 -57.0047861
## 62 -7.9150425
## 63 34.6127519
## 64 -21.0722866
## 65 -22.0009918
## 66 -37.5816067
## 67 -26.0583831
## 68 -25.5511822
## 69 15.6898060
## 70 -11.9155347
## 71 -2.0728694
## 72 -38.2242913
## 73 10.3876466
## 74 -22.3937288
## 75 16.7129825
## 76 36.0899293
## 77 -62.1710571
## 78 -11.0250420
## 79 43.4373665
## 80 -28.5252657
## 81 -17.3349979
## 82 58.2041869
## 83 -31.7292205
## 84 2.1957465
## 85 -7.5263700
## 86 6.0297930
## 87 -18.0886178
## 88 -14.5646838
## 89 -32.2092742
## 90 -3.0357398
## 91 -22.9596676
## 92 -69.3639459
## 93 5.2832261
## 94 -3.5250431
## 95 -24.0859349
## 96 -9.1837060
## 97 26.4214991
## 98 3.1212001
## 99 -2.9193041
## 100 -7.6730899
## 101 52.0041912
## 102 10.4210157
## 103 73.7871802
## 104 7.9854131
## 105 6.9306678
## 106 -26.3200156
## 107 -38.1956964
## 108 -37.1411765
## 109 -51.3551775
## 110 15.2312039
## 111 37.1313178
## 112 2.0545188
## 113 1.3651902
## 114 19.5562302
## 115 1.4943123
## 116 -6.1604300
## 117 -4.8525782
## 118 -61.9776211
## 119 -11.2002151
## 120 -15.1234596
## 121 -21.5609770
## 122 -4.7202419
## 123 50.4386227
## 124 -5.7864987
## 125 6.7136905
## 126 -35.9140143
## 127 50.9977671
## 128 -49.5940852
## 129 -34.1058846
## 130 -46.9621445
## 131 140.6245747
## 132 71.7959550
## 133 30.0581974
## 134 4.3121884
## 135 0.7502335
## 136 28.9102154
## 137 17.3960181
## 138 19.5045244
## 139 -12.7465335
## 140 -13.0911168
## 141 -36.7362018
## 142 -9.2131004
## 143 -4.9409642
## 144 -36.2423122
## 145 -31.9387023
## 146 -10.6584048
## 147 54.1763052
## 148 -25.9962973
## 149 8.2843516
## 150 -20.6385139
## 151 62.6625611
## 152 -33.3544363
## 153 3.7729886
## 154 -11.3760631
## 155 5.0805779
## 156 17.9950120
## 157 -7.8494101
## 158 -0.2201702
## 159 35.3022384
## 160 -11.5540752
## 161 0.6591253
## 162 -5.5871655
## 163 7.5462798
## 164 -9.3303245
## 165 -19.0175482
## 166 55.7263685
## 167 32.0402256
## 168 31.6353338
## 169 11.4108832
## 170 76.5904762
## 171 -32.7991015
## 172 -0.2825814
## 173 -19.0264844
## 174 9.3293014
## 175 54.2493038
## 176 -22.1047290
## 177 11.5543472
## 178 12.5352601
## 179 98.7318828
## 180 -0.3195726
## 181 4.3908269
## 182 43.6649823
## 183 -53.0060931
## 184 -18.3203130
## 185 33.8772333
## 186 -24.7792769
## 187 -11.5610805
## 188 2.7811315
## 189 73.1591927
## 190 -31.1749134
## 191 14.6672356
## 192 -37.2279106
## 193 -45.8529068
## 194 -17.4973688
## 195 -14.2672641
## 196 -56.8557034
## 197 -35.2780738
## 198 5.4730833
## 199 -12.5781921
## 200 46.9637691#Literal 3/3 Calcule un intervalo de confianza del 96.8% para el impacto del gasto de publicidad en TV, en las ventas, e interpretelo
confint(object = ecuacion_ventas,parm = "tv", level = .968)## 1.6 % 98.4 %
## tv -0.2097376 0.2998052Interpretacion:podemos concluir que con el 96.8% de las ocasiones que estimasemos la ecuacion e esperaria que el impacto de 1 millon de dolares en gastos de publicidad tv se tradujese en el minimo punto 0.20 de millones de dolares en las tv hasta un maximo de 0.29 millones en gastos de publicidad en tv.