求助：关于双对数坐标下如何作一元线性回归拟合

我是把横纵坐标轴给调成了双对数坐标轴，想研究幂律分布问题的，然后想做一个一元线性回归来拟合，但是做出来图成了这样，该怎么样才能正常画出拟合的图呢？小白一个，刚刚学习stata。。。

这个是我的程序代码：
clear
input ///
x        y
1        6547
2        4074
3        2851
4        2289
5        1656
6        1438
7        1319
8        1133
9        991
10        758
11        581
12        474
13        365
14        381
15        298
16        278
17        248
18        221
19        166
20        246
21        177
22        139
23        128
24        128
25        119
26        105
27        107
28        96
29        81
30        110
31        78
32        76
33        50
34        63
35        57
36        70
37        63
38        51
39        43
40        61
41        44
42        53
43        50
44        43
45        30
46        40
47        41
48        31
49        21
50        38
51        21
52        26
53        26
54        17
55        25
56        24
57        18
58        23
59        15
60        19
61        23
62        24
63        17
64        15
65        15
66        9
67        19
68        15
69        11
70        19
71        9
72        14
73        9
74        12
75        9
76        6
77        10
78        6
79        11
80        16
81        12
82        7
83        7
84        3
85        11
86        8
87        4
88        4
89        5
90        12
91        9
92        9
93        12
94        7
95        10
96        5
97        7
98        7
99        3
100        12
101        3
102        5
103        7
104        6
105        5
106        9
107        7
108        7
109        4
110        7
111        4
112        4
113        8
114        3
115        2
116        5
118        2
119        4
120        8
121        5
122        3
123        7
124        2
127        9
128        5
129        3
130        4
132        2
133        5
134        2
135        3
136        2
137        2
138        3
141        4
145        3
147        2
149        3
150        5
151        3
153        2
158        4
159        2
163        4
166        3
167        2
171        2
173        3
174        3
177        2
181        2
183        3
185        2
190        2
192        3
195        2
200        6
204        2
205        2
206        2
212        3
217        3
222        3
230        2
241        3
248        2
252        2
260        2
264        2
280        3
295        3
300        3
315        2
330        2
400        2
429        2
438        2
1002        2
end
twoway(scatter y x)(lfit y x), xsca(log) ysca(log) ///
xlabel(1 10 100 1000) ylabel( 10 100 1000 10000)

没人知道吗。。。