stata 编程分析数据问题

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

用编程分析数据，出表和图时无法显示，具体程序如下，请问各位专家，能给看出其中的问题不？
4. Table

1. quietly {
   
2.     log using patent_survey_5.log
   
3.    
   
4. forvalues i = 2013(1)2015 {
   
5.         noisily tabulate rdex`i'interval rdex`i'check
   
6.   noisily list RDexpense`i'_thRMB rdex`i'check rdex`i'interval if rdex`i'check > rdex`i'interval & rdex`i'check !=. & rdex`i'interval !=.
   
7.     }
   
8. 
   
9. global var_people year employee HE_ratio HEemployee RD_ratio RDemployee
   
10. global var_patent Nb_patent effectivepatent_per_em effectivepatent_per_rdem invention utilitymodel design invention_proportion
    
11.     global var_expense RDexpense2015_thRMB RDexpense2015_per_em RDexpense2015_per_rdem RDexpense2015_per_patent growth_rate_2015 patentexpense_thRMB
    
12. global var_interval rdratio_code pexratio_code govpatentratio_code patentenforceratio_code inventionenforceratio_code UMenforceratio_code ///
    
13.                      designenforceratio_code patentproduct_code inventionproduct_code umproduct_code designproduct_code patentlicensing_code ///
    
14.       inventionlicensing_code umlicensing_code designlicensing_code tranferratio_code inventiontransfer_code umtransfer_code designtransfer_code ///
    
15.       patentincome_code rdperiod_code
    
16. global var_group1 province_code region
    
17. global var_group2 industry_regroup2
    
18. global var_group3 year_interval HE_level
    
19. global var_group4 scale2 type1 ownership_regroup
    
20.     global var_group5 patentenforceratio_code patentproduct_code patentlicensing_code tranferratio_code patentincome_code rdperiod_code
    
21.    
    
22. * mytabstat
    
23.     * syntax varlist, i(varlist) [j(varlist)]
    
24. noisily: mytabstat $var_people , i($var_group1 $var_group2 $var_group3 $var_group4 patent_interval)
    
25.     noisily: mytabstat $var_patent , i($var_group1 $var_group2 $var_group3 $var_group4 patent_interval)
    
26.     noisily: mytabstat $var_expense , i($var_group1 $var_group2 $var_group3 $var_group4 patent_interval)
    
27. noisily: mytabstat $var_people , i($var_group1) j($var_group2 $var_group3 $var_group4 patent_interval)   
    
28. noisily: mytabstat $var_patent , i($var_group1) j($var_group2 $var_group3 $var_group4 patent_interval)
    
29.     noisily: mytabstat $var_expense , i($var_group1) j($var_group2 $var_group3 $var_group4 patent_interval)
    
30. noisily: mytabstat $var_people , i($var_group2) j($var_group3 $var_group4 patentNb_if)
    
31.     noisily: mytabstat $var_patent , i($var_group2) j($var_group3 $var_group4 patentNb_if)
    
32.     noisily: mytabstat $var_expense , i($var_group2) j($var_group3 $var_group4 patentNb_if)
    
33.     noisily: mytabstat $var_people , i($var_group5)
    
34.     noisily: mytabstat $var_patent , i($var_group5)
    
35.     noisily: mytabstat $var_expense , i($var_group5)
    
36.     *mytabulate
    
37.     *syntax varlist, i(varlist) [j(varlist)]
    
38.     noisily: mytabulate $var_group1 $var_group2 $var_group3 $var_group4 , i($var_group1 $var_group2 $var_group3 $var_group4)
    
39. noisily: mytabulate $var_interval , i($var_group1 $var_group2 $var_group3 $var_group4)
    
40.     * myHHI
    
41.     * syntax varlist [, i(varlist)]
    
42.     noisily: myHHI Nb_patent invention utilitymodel design employee HEemployee RDemployee
    
43. noisily: myHHI Nb_patent invention utilitymodel design employee HEemployee RDemployee, i($var_group1 $var_group2)
    
44.     * myHHI2
    
45.     * syntax varlist, i(varlist)
    
46.     noisily: myHHI2 Nb_patent invention utilitymodel design employee HEemployee RDemployee, i(industry_code city province_code)
    
47.     * mybig
    
48.     * syntax varlist, n(integer) [i(varlist)] [j(varlist)] [k(varlist)]
    
49.     noisily: mybig Nb_patent invention employee RDemployee RDexpense2015_per_patent RDexpense2015_thRMB, n(25) ///
    
50.           j(Nb_patent invention employee RDexpense2015_thRMB year) ///
    
51.     k(industry_code region province_code city type1 ownership_regroup HElevel)
    
52.    
    
53.     * mybig2
    
54.     * syntax varlist, n(integer) i(varlist)
    
55.     gen Nb_firm = 1
    
56. noisily: mybig2 Nb_firm Nb_patent employee HEemployee RDemployee, n(25) i(province_code industry_code)
    
57.     noisily: mybig2 Nb_firm Nb_patent employee HEemployee RDemployee, n(10) i(city)
    
58. * mybig2_share
    
59.     * syntax varlist, n(integer) i(varlist)
    
60. noisily: mybig2_share HEhigh HElow IPRreduce IPRimprove gf patentinfringe cord IPmanagement, n(25) i(city industry_code)
    
61. noisily: mybig2_share HEhigh HElow IPRreduce IPRimprove gf patentinfringe cord IPmanagement, n(10) i(province_code)
    
62. 
    
63. * patent>100
    
64. preserve
    
65. keep if patent_interval == 10
    
66. noisily display "In the firms holding more than 100 effective patents"
    
67. noisily: mytabstat $var_people , i($var_group1 $var_group2 $var_group3 $var_group4)
    
68. noisily display "In the firms holding more than 100 effective patents"
    
69.     noisily: mytabstat $var_patent , i($var_group1 $var_group2 $var_group3 $var_group4)
    
70. noisily display "In the firms holding more than 100 effective patents"
    
71.     noisily: mytabstat $var_expense , i($var_group1 $var_group2 $var_group3 $var_group4)
    
72. restore
    
73.     log close
    
74. }
    
75. * myscatter
    
76. * syntax varlist, y(varlist) [i(varlist)]
    
77. myscatter RDexpense2015_thRMB y(Nb_patent effectivepatent_per_em effectivepatent_per_rdem invention_proportion)
    
78. myscatter RDemployee RD_ratio employee year, y(Nb_patent invention_proportion)
    
79. myscatter year, y(employee)
    
80. myscatter rdratio_code, y(Nb_patent effectivepatent_per_em effectivepatent_per_rdem invention_proportion)
    
81. * myscatter_province
    
82. * syntax varlist, y(varlist)
    
83. myscatter_province RDexpense2015_thRMB y(Nb_patent effectivepatent_per_em effectivepatent_per_rdem invention_proportion)
    
84. myscatter_province RDemployee RD_ratio employee year, y(Nb_patent invention_proportion)
    
85. myscatter_province year, y(employee)
    
86. myscatter_province rdratio_code, y(Nb_patent effectivepatent_per_em effectivepatent_per_rdem invention_proportion)
    
87. * myscatter_sector
    
88. * syntax varlist, y(varlist)
    
89. myscatter_sector RDexpense2015_thRMB y(Nb_patent effectivepatent_per_em effectivepatent_per_rdem invention_proportion)
    
90. myscatter_sector RDemployee RD_ratio employee year, y(Nb_patent invention_proportion)
    
91. myscatter_sector year, y(employee)
    
92. myscatter_sector rdratio_code, y(Nb_patent effectivepatent_per_em effectivepatent_per_rdem invention_proportion)
    
93. * myscatter_weight
    
94. * syntax varlist, y(varlist) [i(varlist)]
    
95. myscatter_weight rdratio_code pexratio_code HElevel, y($var_group5)
    
96. * mybar
    
97. * syntax varlist, type(string) i(varlist) [j(varlist)]
    
98. mybar employee HEemployee RDemployee Nb_patent invention utilitymodel design RDexpense2015_thRMB patentexpense_thRMB, ///
    
99.       type(sum)　i(scale2 year_interval ownership_regroup)
    
100. mybar HElow HEhigh HE_ratio RD_ratio effectivepatent_per_em effectivepatent_per_rdem invention_proportion RDexpense2015_per_em RDexpense2015_per_rdem RDexpense2015_per_patent growth_rate_2015, ///
     
101.       type(mean)　i(scale2 year_interval ownership_regroup)
     
102. mybar patentNb_if patentadvantage_if IPRreduce IPRimprove gf patentinfringe noproductpatent nocraftpatent cord IPmanagement, ///
     
103.       type(mean)　i(scale2 year_interval ownership_regroup)
     
104. * mybar_province
     
105. * syntax varlist, type(string) i(varlist)
     
106. mybar_province employee HEemployee RDemployee Nb_patent invention utilitymodel design RDexpense2015_thRMB patentexpense_thRMB, ///
     
107.                type(sum)　i(scale2 year_interval ownership_regroup)
     
108. mybar_province HElow HE high HE_ratio RD_ratio effectivepatent_per_em effectivepatent_per_rdem invention_proportion RDexpense2015_per_em RDexpense2015_per_rdem RDexpense2015_per_patent growth_rate_2015, ///
     
109.                type(mean)　i(scale2 year_interval ownership_regroup)
     
110. mybar_province patentNb_if patentadvantage_if IPRreduce IPRimprove gf patentinfringe noproductpatent nocraftpatent cord IPmanagement, ///
     
111.                type(mean)　i(scale2 year_interval ownership_regroup)
     
112. * mybar_sector
     
113. * syntax varlist, type(string) i(varlist)
     
114. mybar_sector employee HEemployee RDemployee Nb_patent invention utilitymodel design RDexpense2015_thRMB patentexpense_thRMB, ///
     
115.              type(sum)　i(scale2 year_interval ownership_regroup)
     
116. mybar_sector HElow HE high HE_ratio RD_ratio effectivepatent_per_em effectivepatent_per_rdem invention_proportion RDexpense2015_per_em RDexpense2015_per_rdem RDexpense2015_per_patent growth_rate_2015, ///
     
117.              type(mean)　i(scale2 year_interval ownership_regroup)
     
118. mybar_sector patentNb_if patentadvantage_if IPRreduce IPRimprove gf patentinfringe noproductpatent nocraftpatent cord IPmanagement, ///
     
119.              type(mean)　i(scale2 year_interval ownership_regroup)
     
120. 
     
121. 表格部分除了list可以显示之外，其他都无法显示。以下是关于图的显示结果：
     
122. myscatter RDexpense2015_thRMB y(Nb_patent effectivepatent_per_em effectivepatent_per_rdem invention_proportion)
     
123. factor variables and time-series operators not allowed
     
124. r(101);
     
125. .
     
126. . myscatter RDemployee RD_ratio employee year, y(Nb_patent invention_proportion)
     
127. command xtitle is unrecognized
     
128. r(199);
     
129. .
     
130. . myscatter year, y(employee)
     
131. command xtitle is unrecognized
     
132. r(199);
     
133. .
     
134. . myscatter rdratio_code, y(Nb_patent effectivepatent_per_em effectivepatent_per_rdem invention_proportion)
     
135. command xtitle is unrecognized
     
136. r(199);
     
137. .

复制代码

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：Stata tata unrecognized time-series Proportion

大段代码用代码格式发就好看多了。
另外，最好能提供一些数据作为测试用。

而且，你这里头这么多自己写的程序，你都不发出来，就给个syntax，然后报错，那怎么可能知道具体是哪里错呢？自己定义的mytabstat这些程序，涉及到的最好是发出来。比如

1. myscatter RDexpense2015_thRMB y(Nb_patent effectivepatent_per_em effectivepatent_per_rdem invention_proportion)
   
2. factor variables and time-series operators not allowed
   
3. r(101);

复制代码

看得到的变量列表里没有factor variables and time-series operators， 那就只可能是你的myscatter程序里有。你不发程序出来，自然就看不到。

另外有个不解的地方就是：你都自己写了这么多程序了，怎么factor variables and time-series operators not allowed和command xtitle is unrecognized这种错误都找不出来呢？程序是你自己写的吗还是别处找来的？