0% found this document useful (0 votes)
102 views17 pages

Regresi Spasial Stata

The document contains a log of commands executed in Stata for data analysis, including the loading of datasets, setting up panel data structures, and performing statistical summaries and regressions. It details various commands, their outputs, and errors encountered during execution. Key results include a random-effects GLS regression with significant coefficients for certain variables.

Uploaded by

rimamelati606
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
0% found this document useful (0 votes)
102 views17 pages

Regresi Spasial Stata

The document contains a log of commands executed in Stata for data analysis, including the loading of datasets, setting up panel data structures, and performing statistical summaries and regressions. It details various commands, their outputs, and errors encountered during execution. Key results include a random-effects GLS regression with significant coefficients for certain variables.

Uploaded by

rimamelati606
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd

Thursday October 20 [Link] 2022 Page 1

___ ____ ____ ____ ____(R)


/__ / ____/ / ____/
___/ / /___/ / /___/
Statistics/Data Analysis

User: rima

___ ____ ____ ____ ____ (R)


/__ / ____/ / ____/
___/ / /___/ / /___/ 14.2 Copyright 1985-2015 StataCorp LP
Statistics/Data Analysis StataCorp
4905 Lakeway Drive
College Station, Texas 77845 USA
800-STATA-PC [Link]
979-696-4600 stata@[Link]
979-696-4601 (fax)

Single-user Stata perpetual license:


Serial number: 10699393
Licensed to: Andrey

Notes:
1. Unicode is supported; see help unicode_advice.

1 . use "C:\Users\winan\Documents\Analisis Stata\data [Link]", clear

2 . do "C:\Users\winan\AppData\Local\Temp\[Link]"

3 . spmat idistance W_invn LONGITUDE LATITUDE, id(id) normalize(row)

4 .
end of do-file

5 . do "C:\Users\winan\AppData\Local\Temp\[Link]"

6 . smpat summarize W_invn


command smpat is unrecognized
r(199);

end of do-file

r(199);

7 . do "C:\Users\winan\AppData\Local\Temp\[Link]"

8 . spmat summarize W_invn

Summary of spatial-weighting object W_invn

Matrix Description

Dimensions 113 x 113


Stored as 113 x 113
Values
min 0
min>0 .0004474
mean .0088496
max .7718628
Thursday October 20 [Link] 2022 Page 2

9 .
end of do-file

10 . use "C:\Users\winan\Documents\Analisis Stata\data [Link]", clear

11 . do "C:\Users\winan\AppData\Local\Temp\[Link]"

12 . xtset id Tahun
panel variable: id (strongly balanced)
time variable: Tahun, 2015 to 2019
delta: 1 unit

13 .
end of do-file

14 . do "C:\Users\winan\AppData\Local\Temp\[Link]"

15 . global varxpanel PAD DAPER BP BB BM lnPDRBAWAL IPM JALAN

16 .
end of do-file

17 . do "C:\Users\winan\AppData\Local\Temp\[Link]"

18 . xtdescribe

id: 1, 2, ..., 113 n = 113


Tahun: 2015, 2016, ..., 2019 T = 5
Delta(Tahun) = 1 unit
Span(Tahun) = 5 periods
(id*Tahun uniquely identifies each observation)

Distribution of T_i: min 5% 25% 50% 75% 95% max


5 5 5 5 5 5 5

Freq. Percent Cum. Pattern

113 100.00 100.00 11111

113 100.00 XXXXX

19 . xtsum Y

Variable Mean Std. Dev. Min Max Observations

Y overall .0554211 .0850475 -.0358893 .7705331 N = 565


between .0369199 .0047105 .192519 n = 113
within .076679 -.0961691 .6334352 T = 5

20 .
end of do-file
Thursday October 20 [Link] 2022 Page 3

21 . do "C:\Users\winan\AppData\Local\Temp\[Link]"

22 . xtsum Y $varxpanel

Variable Mean Std. Dev. Min Max Observations

Y overall .0554211 .0850475 -.0358893 .7705331 N = 565


between .0369199 .0047105 .192519 n = 113
within .076679 -.0961691 .6334352 T = 5

PAD overall .1874028 .1033401 .0630758 .6140133 N = 565


between .1023878 .0741375 .5973664 n = 113
within .0164396 .1227724 .3086618 T = 5

DAPER overall .6058808 .1040745 .2092107 .7979695 N = 565


between .092763 .2553034 .7075371 n = 113
within .0478284 .5036778 .7421331 T = 5

BP overall .435803 .0656658 .2256855 .6403283 N = 565


between .0471806 .2925538 .5403104 n = 113
within .0458452 .312509 .5965983 T = 5

BB overall .2284096 .0738348 0 .463549 N = 565


between .0644476 .1133992 .3967522 n = 113
within .0364355 -.008458 .3465105 T = 5

BM overall .1921365 .0533169 .0429147 .4126614 N = 565


between .0420703 .1184212 .3589152 n = 113
within .0329444 .0929669 .3236293 T = 5

lnPDRB~L overall 10.13769 .629289 9.266579 12.71835 N = 565


between .6220136 9.348675 12.51616 n = 113
within .1088468 9.464232 10.38127 T = 5

IPM overall 70.97543 5.29034 58.18 86.65 N = 565


between 5.235828 60.022 85.626 n = 113
within .8764859 69.04543 72.89343 T = 5

JALAN overall .8050443 .2120316 0 1 N = 565


between .153333 .3416725 .9996946 n = 113
within .147014 .0676735 1.38599 T = 5

23 .
end of do-file

24 . do "C:\Users\winan\AppData\Local\Temp\[Link]"

25 . ssc install xtcsd


checking xtcsd consistency and verifying not already installed...
all files already exist and are up to date.
Thursday October 20 [Link] 2022 Page 4

26 .
27 .
end of do-file

28 . do "C:\Users\winan\AppData\Local\Temp\[Link]"

29 . xtreg Y $varxpanel

Random-effects GLS regression Number of obs = 565


Group variable: id Number of groups = 113

R-sq: Obs per group:


within = 0.5866 min = 5
between = 0.0572 avg = 5.0
overall = 0.0047 max = 5

Wald chi2(8) = 133.45


corr(u_i, X) = 0 (assumed) Prob > chi2 = 0.0000

Y Coef. Std. Err. z P>|z| [95% Conf. Interval]

PAD .4109484 .1244908 3.30 0.001 .1669509 .6549458


DAPER -.1211693 .0644538 -1.88 0.060 -.2474965 .0051578
BP .0852955 .0976374 0.87 0.382 -.1060702 .2766613
BB .1535818 .1315783 1.17 0.243 -.1043069 .4114705
BM .1756479 .1085453 1.62 0.106 -.0370969 .3883927
lnPDRBAWAL -.178378 .017503 -10.19 0.000 -.2126832 -.1440728
IPM .0031466 .0022565 1.39 0.163 -.001276 .0075692
JALAN -.0247307 .0193393 -1.28 0.201 -.062635 .0131736
_cons 1.55074 .1994708 7.77 0.000 1.159784 1.941695

sigma_u .03153183
sigma_e .02377301
rho .63758405 (fraction of variance due to u_i)

30 .
end of do-file

31 . do "C:\Users\winan\AppData\Local\Temp\[Link]"

32 . xtcsd pes abs


varlist not allowed
r(101);

end of do-file

r(101);

33 . do "C:\Users\winan\AppData\Local\Temp\[Link]"
Thursday October 20 [Link] 2022 Page 5

34 . xtcsd pes abs


varlist not allowed
r(101);

end of do-file

r(101);

35 . do "C:\Users\winan\AppData\Local\Temp\[Link]"

36 . ssc install xtcsd


checking xtcsd consistency and verifying not already installed...
all files already exist and are up to date.

37 .
end of do-file

38 . do "C:\Users\winan\AppData\Local\Temp\[Link]"

39 . xtcsd pes abs


varlist not allowed
r(101);

end of do-file

r(101);

40 . do "C:\Users\winan\AppData\Local\Temp\[Link]"

41 . clear

42 . #INPUT DATA LONGLAT


Unknown #command
43 . import excel
using required
r(100);

end of do-file

r(100);

44 . do "C:\Users\winan\AppData\Local\Temp\[Link]"

45 . xtcsd, pes abs


time variable not set, use tsset varname ...
r(111);

end of do-file

r(111);
Thursday October 20 [Link] 2022 Page 6

46 . do "C:\Users\winan\AppData\Local\Temp\[Link]"

47 . ssc install xtcsd


checking xtcsd consistency and verifying not already installed...
all files already exist and are up to date.

48 .
end of do-file

49 . do "C:\Users\winan\AppData\Local\Temp\[Link]"

50 . xtreg Y $varxpanel
variable Y not found
r(111);

end of do-file

r(111);

51 . import excel "C:\Users\winan\Documents\KOTAKAB\[Link]", sheet("Sheet1") firstrow

52 . do "C:\Users\winan\AppData\Local\Temp\[Link]"

53 . xtcsd, pes abs


time variable not set, use tsset varname ...
r(111);

end of do-file

r(111);

54 . do "C:\Users\winan\AppData\Local\Temp\[Link]"

55 . xtset id Tahun
panel variable: id (strongly balanced)
time variable: Tahun, 2015 to 2019
delta: 1 unit

56 .
end of do-file

57 . do "C:\Users\winan\AppData\Local\Temp\[Link]"

58 . global varxpanel PAD DAPER BP BB BM lnPDRBAWAL IPM JALAN

59 .
end of do-file

60 . do "C:\Users\winan\AppData\Local\Temp\[Link]"

61 . xtdescribe

id: 1, 2, ..., 113 n = 113


Tahun: 2015, 2016, ..., 2019 T = 5
Delta(Tahun) = 1 unit
Span(Tahun) = 5 periods
(id*Tahun uniquely identifies each observation)

Distribution of T_i: min 5% 25% 50% 75% 95% max


5 5 5 5 5 5 5
Thursday October 20 [Link] 2022 Page 7

Freq. Percent Cum. Pattern

113 100.00 100.00 11111

113 100.00 XXXXX

62 .
end of do-file

63 . do "C:\Users\winan\AppData\Local\Temp\[Link]"

64 . xtsum Y

Variable Mean Std. Dev. Min Max Observations

Y overall .0554211 .0850475 -.0358893 .7705331 N = 565


between .0369199 .0047105 .192519 n = 113
within .076679 -.0961691 .6334352 T = 5

65 .
end of do-file

66 . do "C:\Users\winan\AppData\Local\Temp\[Link]"

67 . ssc install xtcsd


checking xtcsd consistency and verifying not already installed...
all files already exist and are up to date.

68 .
end of do-file

69 . do "C:\Users\winan\AppData\Local\Temp\[Link]"

70 . xtreg Y $varxpanel

Random-effects GLS regression Number of obs = 565


Group variable: id Number of groups = 113

R-sq: Obs per group:


within = 0.5866 min = 5
between = 0.0572 avg = 5.0
overall = 0.0047 max = 5

Wald chi2(8) = 133.45


corr(u_i, X) = 0 (assumed) Prob > chi2 = 0.0000

Y Coef. Std. Err. z P>|z| [95% Conf. Interval]

PAD .4109484 .1244908 3.30 0.001 .1669509 .6549458


DAPER -.1211693 .0644538 -1.88 0.060 -.2474965 .0051578
BP .0852955 .0976374 0.87 0.382 -.1060702 .2766613
BB .1535818 .1315783 1.17 0.243 -.1043069 .4114705
BM .1756479 .1085453 1.62 0.106 -.0370969 .3883927
lnPDRBAWAL -.178378 .017503 -10.19 0.000 -.2126832 -.1440728
IPM .0031466 .0022565 1.39 0.163 -.001276 .0075692
JALAN -.0247307 .0193393 -1.28 0.201 -.062635 .0131736
_cons 1.55074 .1994708 7.77 0.000 1.159784 1.941695

sigma_u .03153183
sigma_e .02377301
rho .63758405 (fraction of variance due to u_i)
Thursday October 20 [Link] 2022 Page 8

71 .
end of do-file

72 . do "C:\Users\winan\AppData\Local\Temp\[Link]"

73 . xtcsd, pes abs

Pesaran's test of cross sectional independence = 70.641, Pr = 0.0000

Average absolute value of the off-diagonal elements = 0.639

74 .
end of do-file

75 . do "C:\Users\winan\AppData\Local\Temp\[Link]"

76 . xsmle Y $varxpanel, wmat(W_invn) mod(sar) hausman nolog


convergence not achieved
... estimating fixed-effects model to perform Hausman test

SAR with random-effects Number of obs = 565

Group variable: id Number of groups = 113


Time variable: Tahun Panel length = 5

R-sq: within = 0.0272


between = 0.0001
overall = 0.0185

Log-likelihood = 760.8760

Y Coef. Std. Err. z P>|z| [95% Conf. Interval]

Main
PAD .1922362 .0700122 2.75 0.006 .0550148 .3294576
DAPER .0063101 .0463518 0.14 0.892 -.0845378 .097158
BP -.0318979 .0600582 -0.53 0.595 -.1496097 .085814
BB .1071028 .0702699 1.52 0.127 -.0306236 .2448293
BM .051306 .0680503 0.75 0.451 -.0820702 .1846821
lnPDRBAWAL -.0457992 .0077857 -5.88 0.000 -.0610589 -.0305395
IPM -.0009958 .0009718 -1.02 0.306 -.0029005 .000909
JALAN -.0110201 .0133398 -0.83 0.409 -.0371656 .0151255
_cons .4861848 .0896197 5.42 0.000 .3105333 .6618362

Spatial
rho .9691737 .013626 71.13 0.000 .9424673 .9958802

Variance
lgt_theta .6368401 . . . . .
sigma2_e .0030927 .0001891 16.35 0.000 .002722 .0034633

Ho: difference in coeffs not systematic chi2(9) = 249.83 Prob>=chi2 = 0.0000


Thursday October 20 [Link] 2022 Page 9

77 .
end of do-file

78 . do "C:\Users\winan\AppData\Local\Temp\[Link]"

79 . estat ic

Akaike's information criterion and Bayesian information criterion

Model Obs ll(null) ll(model) df AIC BIC

. 565 . 760.876 11 -1499.752 -1452.047

Note: N=Obs used in calculating BIC; see [R] BIC note.

80 .
end of do-file

81 . do "C:\Users\winan\AppData\Local\Temp\[Link]"

82 . estimates store SARpnel

83 .
end of do-file

84 . do "C:\Users\winan\AppData\Local\Temp\[Link]"

85 . xsmle Y $varxpanel, emat(W_invn) mod(sem) hausman nolog


... estimating fixed-effects model to perform Hausman test

SEM with random-effects Number of obs = 565

Group variable: id Number of groups = 113


Time variable: Tahun Panel length = 5

R-sq: within = 0.0030


between = 0.1263
overall = 0.0258

Log-likelihood = 770.2903

Y Coef. Std. Err. z P>|z| [95% Conf. Interval]

Main
PAD .0232489 .0805439 0.29 0.773 -.1346141 .181112
DAPER -.0134794 .0698687 -0.19 0.847 -.1504197 .1234608
BP -.0038933 .0779874 -0.05 0.960 -.1567457 .1489592
BB .0870093 .064819 1.34 0.179 -.0400336 .2140521
BM .0524397 .0706202 0.74 0.458 -.0859733 .1908526
lnPDRBAWAL -.0016876 .0059646 -0.28 0.777 -.0133779 .0100027
IPM -.0009241 .0008941 -1.03 0.301 -.0026766 .0008283
JALAN -.0047498 .0120016 -0.40 0.692 -.0282725 .0187728
_cons .140556 .109103 1.29 0.198 -.0732819 .354394

Spatial
lambda .9704617 .0130536 74.34 0.000 .9448771 .9960462

Variance
ln_phi -16.25903 440.2437 -0.04 0.971 -879.1209 846.6028
sigma2_e .0035424 .000213 16.63 0.000 .0031249 .00396

Ho: difference in coeffs not systematic chi2(9) = 893.12 Prob>=chi2 = 0.0000


Thursday October 20 [Link] 2022 Page 10

86 .
87 .
end of do-file

88 . do "C:\Users\winan\AppData\Local\Temp\[Link]"

89 . xsmle Y $varxpanel, fe wmat(W_invn) emat(W_invn) mod(sac)


Iteration 0: Log-likelihood = 1392.6681 (not concave)
Iteration 1: Log-likelihood = 1411.2807 (not concave)
Iteration 2: Log-likelihood = 1466.2702 (not concave)
Iteration 3: Log-likelihood = 1509.396
Iteration 4: Log-likelihood = 1523.1896
Iteration 5: Log-likelihood = 1531.1716
Iteration 6: Log-likelihood = 1536.2189
Iteration 7: Log-likelihood = 1541.9404
Iteration 8: Log-likelihood = 1545.7335
Iteration 9: Log-likelihood = 1546.5963
Iteration 10: Log-likelihood = 1546.6008
Iteration 11: Log-likelihood = 1546.6008

SAC with spatial fixed-effects Number of obs = 565

Group variable: id Number of groups = 113


Time variable: Tahun Panel length = 5

R-sq: within = 0.4324


between = 0.1170
overall = 0.0000

Mean of fixed-effects = 8.8900

Log-likelihood = 1546.6008

Y Coef. Std. Err. z P>|z| [95% Conf. Interval]

Main
PAD .0383919 .0453388 0.85 0.397 -.0504705 .1272543
DAPER .0301 .0285555 1.05 0.292 -.0258678 .0860678
BP .0253361 .0372245 0.68 0.496 -.0476226 .0982949
BB .0448943 .0365927 1.23 0.220 -.0268262 .1166148
BM .1005354 .0348344 2.89 0.004 .0322611 .1688096
lnPDRBAWAL -.8503843 .0150628 -56.46 0.000 -.8799068 -.8208618
IPM -.0042451 .0030074 -1.41 0.158 -.0101396 .0016494
JALAN -.0026049 .0043641 -0.60 0.551 -.0111583 .0059485

Spatial
rho .4345515 .0559341 7.77 0.000 .3249227 .5441804
lambda .9757173 .0075755 128.80 0.000 .9608695 .9905651

Variance
sigma2_e .0002804 .0000134 20.88 0.000 .000254 .0003067
Thursday October 20 [Link] 2022 Page 11

90 .
end of do-file

91 . do "C:\Users\winan\AppData\Local\Temp\[Link]"

92 . estat ic

Akaike's information criterion and Bayesian information criterion

Model Obs ll(null) ll(model) df AIC BIC

. 565 . 1546.601 11 -3071.202 -3023.497

Note: N=Obs used in calculating BIC; see [R] BIC note.

93 .
end of do-file

94 . do "C:\Users\winan\AppData\Local\Temp\[Link]"

95 . estat ic

Akaike's information criterion and Bayesian information criterion

Model Obs ll(null) ll(model) df AIC BIC

. 565 . 1546.601 11 -3071.202 -3023.497

Note: N=Obs used in calculating BIC; see [R] BIC note.

96 .
end of do-file

97 . do "C:\Users\winan\AppData\Local\Temp\[Link]"

98 . estat ic

Akaike's information criterion and Bayesian information criterion

Model Obs ll(null) ll(model) df AIC BIC

. 565 . 1546.601 11 -3071.202 -3023.497

Note: N=Obs used in calculating BIC; see [R] BIC note.

99 .
end of do-file
Thursday October 20 [Link] 2022 Page 12

100 . do "C:\Users\winan\AppData\Local\Temp\[Link]"

101 . estimates store GSMpnel

102 .
end of do-file

103 . do "C:\Users\winan\AppData\Local\Temp\[Link]"

104 . estimates tabele SARpnel SEMpnel, b


subcommand estimates tabele is unrecognized
r(198);

end of do-file

r(198);

105 . do "C:\Users\winan\AppData\Local\Temp\[Link]"

106 . estimates tabele SARpnel SEMpnel GSMpnel, b


subcommand estimates tabele is unrecognized
r(198);

end of do-file

r(198);

107 . do "C:\Users\winan\AppData\Local\Temp\[Link]"

108 . estimates table SARpnel SEMpnel, b


estimation result SEMpnel not found
r(111);

end of do-file

r(111);

109 . do "C:\Users\winan\AppData\Local\Temp\[Link]"

110 . xsmle Y $varxpanel, wmat(W_invn) mod(sar) hausman nolog


convergence not achieved
... estimating fixed-effects model to perform Hausman test

SAR with random-effects Number of obs = 565

Group variable: id Number of groups = 113


Time variable: Tahun Panel length = 5

R-sq: within = 0.0272


between = 0.0001
overall = 0.0185
Thursday October 20 [Link] 2022 Page 13

Log-likelihood = 760.8760

Y Coef. Std. Err. z P>|z| [95% Conf. Interval]

Main
PAD .1922362 .0700122 2.75 0.006 .0550148 .3294576
DAPER .0063101 .0463518 0.14 0.892 -.0845378 .097158
BP -.0318979 .0600582 -0.53 0.595 -.1496097 .085814
BB .1071028 .0702699 1.52 0.127 -.0306236 .2448293
BM .051306 .0680503 0.75 0.451 -.0820702 .1846821
lnPDRBAWAL -.0457992 .0077857 -5.88 0.000 -.0610589 -.0305395
IPM -.0009958 .0009718 -1.02 0.306 -.0029005 .000909
JALAN -.0110201 .0133398 -0.83 0.409 -.0371656 .0151255
_cons .4861848 .0896197 5.42 0.000 .3105333 .6618362

Spatial
rho .9691737 .013626 71.13 0.000 .9424673 .9958802

Variance
lgt_theta .6368401 . . . . .
sigma2_e .0030927 .0001891 16.35 0.000 .002722 .0034633

Ho: difference in coeffs not systematic chi2(9) = 249.83 Prob>=chi2 = 0.0000

111 . estat ic

Akaike's information criterion and Bayesian information criterion

Model Obs ll(null) ll(model) df AIC BIC

. 565 . 760.876 11 -1499.752 -1452.047

Note: N=Obs used in calculating BIC; see [R] BIC note.

112 . estimates store SARpnel

113 .
end of do-file

114 . do "C:\Users\winan\AppData\Local\Temp\[Link]"

115 . estat ic

Akaike's information criterion and Bayesian information criterion

Model Obs ll(null) ll(model) df AIC BIC

SARpnel 565 . 760.876 11 -1499.752 -1452.047

Note: N=Obs used in calculating BIC; see [R] BIC note.


Thursday October 20 [Link] 2022 Page 14

116 .
end of do-file

117 . do "C:\Users\winan\AppData\Local\Temp\[Link]"

118 . estimates store SARpnel

119 .
end of do-file

120 . do "C:\Users\winan\AppData\Local\Temp\[Link]"

121 . xsmle Y $varxpanel, emat(W_invn) mod(sem) hausman nolog


... estimating fixed-effects model to perform Hausman test

SEM with random-effects Number of obs = 565

Group variable: id Number of groups = 113


Time variable: Tahun Panel length = 5

R-sq: within = 0.0030


between = 0.1263
overall = 0.0258

Log-likelihood = 770.2903

Y Coef. Std. Err. z P>|z| [95% Conf. Interval]

Main
PAD .0232489 .0805439 0.29 0.773 -.1346141 .181112
DAPER -.0134794 .0698687 -0.19 0.847 -.1504197 .1234608
BP -.0038933 .0779874 -0.05 0.960 -.1567457 .1489592
BB .0870093 .064819 1.34 0.179 -.0400336 .2140521
BM .0524397 .0706202 0.74 0.458 -.0859733 .1908526
lnPDRBAWAL -.0016876 .0059646 -0.28 0.777 -.0133779 .0100027
IPM -.0009241 .0008941 -1.03 0.301 -.0026766 .0008283
JALAN -.0047498 .0120016 -0.40 0.692 -.0282725 .0187728
_cons .140556 .109103 1.29 0.198 -.0732819 .354394

Spatial
lambda .9704617 .0130536 74.34 0.000 .9448771 .9960462

Variance
ln_phi -16.25903 440.2437 -0.04 0.971 -879.1209 846.6028
sigma2_e .0035424 .000213 16.63 0.000 .0031249 .00396

Ho: difference in coeffs not systematic chi2(9) = 893.12 Prob>=chi2 = 0.0000

122 .
end of do-file
Thursday October 20 [Link] 2022 Page 15

123 . do "C:\Users\winan\AppData\Local\Temp\[Link]"

124 . estimates store SEMpnel

125 .
end of do-file

126 . do "C:\Users\winan\AppData\Local\Temp\[Link]"

127 . estimates store SEMpnel

128 .
end of do-file

129 . do "C:\Users\winan\AppData\Local\Temp\[Link]"

130 . xsmle Y $varxpanel, fe wmat(W_invn) emat(W_invn) mod(sac)


Iteration 0: Log-likelihood = 1392.6681 (not concave)
Iteration 1: Log-likelihood = 1411.2807 (not concave)
Iteration 2: Log-likelihood = 1466.2702 (not concave)
Iteration 3: Log-likelihood = 1509.396
Iteration 4: Log-likelihood = 1523.1896
Iteration 5: Log-likelihood = 1531.1716
Iteration 6: Log-likelihood = 1536.2189
Iteration 7: Log-likelihood = 1541.9404
Iteration 8: Log-likelihood = 1545.7335
Iteration 9: Log-likelihood = 1546.5963
Iteration 10: Log-likelihood = 1546.6008
Iteration 11: Log-likelihood = 1546.6008

SAC with spatial fixed-effects Number of obs = 565

Group variable: id Number of groups = 113


Time variable: Tahun Panel length = 5

R-sq: within = 0.4324


between = 0.1170
overall = 0.0000

Mean of fixed-effects = 8.8900

Log-likelihood = 1546.6008

Y Coef. Std. Err. z P>|z| [95% Conf. Interval]

Main
PAD .0383919 .0453388 0.85 0.397 -.0504705 .1272543
DAPER .0301 .0285555 1.05 0.292 -.0258678 .0860678
BP .0253361 .0372245 0.68 0.496 -.0476226 .0982949
BB .0448943 .0365927 1.23 0.220 -.0268262 .1166148
BM .1005354 .0348344 2.89 0.004 .0322611 .1688096
lnPDRBAWAL -.8503843 .0150628 -56.46 0.000 -.8799068 -.8208618
IPM -.0042451 .0030074 -1.41 0.158 -.0101396 .0016494
JALAN -.0026049 .0043641 -0.60 0.551 -.0111583 .0059485

Spatial
rho .4345515 .0559341 7.77 0.000 .3249227 .5441804
lambda .9757173 .0075755 128.80 0.000 .9608695 .9905651

Variance
sigma2_e .0002804 .0000134 20.88 0.000 .000254 .0003067
Thursday October 20 [Link] 2022 Page 16

131 .
end of do-file

132 . do "C:\Users\winan\AppData\Local\Temp\[Link]"

133 . estat ic

Akaike's information criterion and Bayesian information criterion

Model Obs ll(null) ll(model) df AIC BIC

. 565 . 1546.601 11 -3071.202 -3023.497

Note: N=Obs used in calculating BIC; see [R] BIC note.

134 .
end of do-file

135 . do "C:\Users\winan\AppData\Local\Temp\[Link]"

136 . estimates store GSMpnel

137 .
end of do-file

138 . do "C:\Users\winan\AppData\Local\Temp\[Link]"

139 . estimates table SARpnel SEMpnel, b

Variable SARpnel SEMpnel

Main
PAD .19223621 .02324893
DAPER .00631006 -.01347943
BP -.03189786 -.00389326
BB .10710283 .08700926
BM .05130598 .05243967
lnPDRBAWAL -.0457992 -.00168758
IPM -.00099576 -.00092414
JALAN -.01102006 -.00474984
_cons .48618477 .14055604

Spatial
rho .96917371
lambda .97046167

Variance
lgt_theta .63684013
sigma2_e .00309269 .00354245
ln_phi -16.259028
Thursday October 20 [Link] 2022 Page 17

140 .
end of do-file

141 . do "C:\Users\winan\AppData\Local\Temp\[Link]"

142 . estimates table SARpnel SEMpnel GSMpnel, b

Variable SARpnel SEMpnel GSMpnel

Main
PAD .19223621 .02324893 .03839192
DAPER .00631006 -.01347943 .03009998
BP -.03189786 -.00389326 .02533614
BB .10710283 .08700926 .0448943
BM .05130598 .05243967 .10053536
lnPDRBAWAL -.0457992 -.00168758 -.85038431
IPM -.00099576 -.00092414 -.00424512
JALAN -.01102006 -.00474984 -.00260492
_cons .48618477 .14055604

Spatial
rho .96917371 .43455153
lambda .97046167 .97571734

Variance
lgt_theta .63684013
sigma2_e .00309269 .00354245 .00028036
ln_phi -16.259028

143 .
end of do-file

144 . save "C:\Users\winan\Documents\Analisis Stata\percobaan [Link]"


file C:\Users\winan\Documents\Analisis Stata\percobaan [Link] saved

145 .

You might also like