A. Jack - thank you for adding the extra functions I will test them
after I have installed the newly compiled version of gretl.
B. Allin - the crash in gretl for Windows occurred when I going
through the process of storing the package locally and then upload it
to the server. It happened only once, which makes the problem
difficult to find because it is clear that it depends on other things
in my environment and perhaps in the outside world. However, as I
noted the packaging process managed to attach the wrong subsidiary
functions on at least one occasion. Subsequently, I was careful to
ensure that I had only one function or linked set of functions loaded
when creating the package - when things went wrong I had everything
loaded by reading and executing the function script files.
C. The first attachment (sfa_panel1.txt) is the full contents of the
script output for the script that I have been testing. I have added
various "print x -o" commands as it progresses. The first one
printing panel_unit seems fine, but the second one printing the first
dummy variables shows all missing values. As you will see my dataset
is quite large, but gretl recognises that it is a panel dataset and
has set the observation range correctly.
As a further thought I have just tested the script without the smpl
command -see the second attachment (sfa_panel2.txt). Then the dummy
variables are created properly and the script works as intended, so
there is a bug linked to the way in which the restricted sample is handled.
Gordon
gretl version 1.7.5
Current session: 2008/07/02 18:56
? set shelldir "F:\Work_Files\CCPR_Paper\Gretl\"
? open LEC_data_gretl.gdt
Read datafile F:\Work_Files\CCPR_Paper\Gretl\LEC_data_gretl.gdt
periodicity: 12, maxobs: 816,
observations range: 1:01-68:12
Listing 113 variables:
0) const 1) year 2) cref
3) cosa 4) company 5) parent
6) state 7) vdummy 8) nyne
9) midwest 10) south 11) rboc
12) wacc 13) coc 14) switched_lines
15) leased_lines 16) total_sheath 17) aerial_sheath
18) non_aerial_shea 19) fibre_prop 20) duct
21) local_switch_mi 22) main_switch_min 23) swminutes
24) central_switche 25) remote_switches 26) busres
27) opc_total 28) opc_network 29) opc_network_xpp
30) tc_total 31) tc_network 32) tc_network_xpp
33) tz_total 34) tz_network 35) tz_network_xpp
36) sname 37) popden 38) population_dist
39) ilec_total_07 40) clec_resale_07 41) clec_une_07
42) clec_own_07 43) clec_total_07 44) lec_total_07
45) wireless_subs_0 46) pop_06 47) hhinc_06
48) popden_00 49) popden_06 50) pct_rur
51) pop_fb_05 52) pct_pop_fb 53) pct_cbsa_06
54) pct_nomove_05 55) pct_union_all 56) pct_union_priv
57) pct_forest 58) price_elect 59) clec_zip_0
60) clec_zip_1to4 61) clec_zip_5plus 62) clec_pct_lines
63) clec_own_lines 64) clec_une_lines 65) clec_resale_lin
66) clec_total_line 67) time 68) avpay
69) cost_aerial 70) cost_under 71) cost_duct
72) ln_gdpdef 73) ln_ncost 74) ln_opcost
75) ln_sl 76) ln_ll 77) ln_swmins
78) ln_ts 79) ln_duct 80) ln_popden
81) ln_fibreprop 82) ln_busres 83) ln_pay
84) ln_caerial 85) ln_cunder 86) ln_cduct
87) ln_ncostpsl 88) ln_opcostpsl 89) ln_llpsl
90) ln_swminpsl 91) ln_tspsl 92) ln_ductpsl
93) d9699 94) ln_sl_d9699 95) ln_ll_d9699
96) ln_swmins_d9699 97) ln_ts_d9699 98) ln_duct_d9699
99) ln_llpsl_d9699 100) ln_swminpsl_d96 101) ln_tspsl_d9699
102) ln_ductpsl_d969 103) time_d9699 104) scale
105) scalesq 106) pct_clec_total 107) pct_clec_own
108) mobilepsl 109) peak_lines 110) strand
111) ln_strand 112) ln_strand_lag
? smpl year > 2002 --restrict
Full data set: 816 observations
Current sample: 1:1 - 68:5 (n = 340)
? set debug 2
? set bfgs_toler 0.000001
? series panel_unit=$unit
Generated series panel_unit (ID 113)
# series panel_time=$period
? print panel_unit -o
Obs panel_unit
1:1 1
1:2 1
1:3 1
1:4 1
1:5 1
2:1 2
2:2 2
2:3 2
2:4 2
2:5 2
3:1 3
3:2 3
3:3 3
3:4 3
3:5 3
4:1 4
4:2 4
4:3 4
4:4 4
4:5 4
5:1 5
5:2 5
5:3 5
5:4 5
5:5 5
6:1 6
6:2 6
6:3 6
6:4 6
6:5 6
7:1 7
7:2 7
7:3 7
7:4 7
7:5 7
8:1 8
8:2 8
8:3 8
8:4 8
8:5 8
9:1 9
9:2 9
9:3 9
9:4 9
9:5 9
10:1 10
10:2 10
10:3 10
10:4 10
10:5 10
11:1 11
11:2 11
11:3 11
11:4 11
11:5 11
12:1 12
12:2 12
12:3 12
12:4 12
12:5 12
13:1 13
13:2 13
13:3 13
13:4 13
13:5 13
14:1 14
14:2 14
14:3 14
14:4 14
14:5 14
15:1 15
15:2 15
15:3 15
15:4 15
15:5 15
16:1 16
16:2 16
16:3 16
16:4 16
16:5 16
17:1 17
17:2 17
17:3 17
17:4 17
17:5 17
18:1 18
18:2 18
18:3 18
18:4 18
18:5 18
19:1 19
19:2 19
19:3 19
19:4 19
19:5 19
20:1 20
20:2 20
20:3 20
20:4 20
20:5 20
21:1 21
21:2 21
21:3 21
21:4 21
21:5 21
22:1 22
22:2 22
22:3 22
22:4 22
22:5 22
23:1 23
23:2 23
23:3 23
23:4 23
23:5 23
24:1 24
24:2 24
24:3 24
24:4 24
24:5 24
25:1 25
25:2 25
25:3 25
25:4 25
25:5 25
26:1 26
26:2 26
26:3 26
26:4 26
26:5 26
27:1 27
27:2 27
27:3 27
27:4 27
27:5 27
28:1 28
28:2 28
28:3 28
28:4 28
28:5 28
29:1 29
29:2 29
29:3 29
29:4 29
29:5 29
30:1 30
30:2 30
30:3 30
30:4 30
30:5 30
31:1 31
31:2 31
31:3 31
31:4 31
31:5 31
32:1 32
32:2 32
32:3 32
32:4 32
32:5 32
33:1 33
33:2 33
33:3 33
33:4 33
33:5 33
34:1 34
34:2 34
34:3 34
34:4 34
34:5 34
35:1 35
35:2 35
35:3 35
35:4 35
35:5 35
36:1 36
36:2 36
36:3 36
36:4 36
36:5 36
37:1 37
37:2 37
37:3 37
37:4 37
37:5 37
38:1 38
38:2 38
38:3 38
38:4 38
38:5 38
39:1 39
39:2 39
39:3 39
39:4 39
39:5 39
40:1 40
40:2 40
40:3 40
40:4 40
40:5 40
41:1 41
41:2 41
41:3 41
41:4 41
41:5 41
42:1 42
42:2 42
42:3 42
42:4 42
42:5 42
43:1 43
43:2 43
43:3 43
43:4 43
43:5 43
44:1 44
44:2 44
44:3 44
44:4 44
44:5 44
45:1 45
45:2 45
45:3 45
45:4 45
45:5 45
46:1 46
46:2 46
46:3 46
46:4 46
46:5 46
47:1 47
47:2 47
47:3 47
47:4 47
47:5 47
48:1 48
48:2 48
48:3 48
48:4 48
48:5 48
49:1 49
49:2 49
49:3 49
49:4 49
49:5 49
50:1 50
50:2 50
50:3 50
50:4 50
50:5 50
51:1 51
51:2 51
51:3 51
51:4 51
51:5 51
52:1 52
52:2 52
52:3 52
52:4 52
52:5 52
53:1 53
53:2 53
53:3 53
53:4 53
53:5 53
54:1 54
54:2 54
54:3 54
54:4 54
54:5 54
55:1 55
55:2 55
55:3 55
55:4 55
55:5 55
56:1 56
56:2 56
56:3 56
56:4 56
56:5 56
57:1 57
57:2 57
57:3 57
57:4 57
57:5 57
58:1 58
58:2 58
58:3 58
58:4 58
58:5 58
59:1 59
59:2 59
59:3 59
59:4 59
59:5 59
60:1 60
60:2 60
60:3 60
60:4 60
60:5 60
61:1 61
61:2 61
61:3 61
61:4 61
61:5 61
62:1 62
62:2 62
62:3 62
62:4 62
62:5 62
63:1 63
63:2 63
63:3 63
63:4 63
63:5 63
64:1 64
64:2 64
64:3 64
64:4 64
64:5 64
65:1 65
65:2 65
65:3 65
65:4 65
65:5 65
66:1 66
66:2 66
66:3 66
66:4 66
66:5 66
67:1 67
67:2 67
67:3 67
67:4 67
67:5 67
68:1 68
68:2 68
68:3 68
68:4 68
68:5 68
? genr unitdum
Panel dummy variables generated.
? list dlist = du_*
Generated list dlist
? print du_1 -o
Obs du_1
1:1
1:2
1:3
1:4
1:5
2:1
2:2
2:3
2:4
2:5
3:1
3:2
3:3
3:4
3:5
4:1
4:2
4:3
4:4
4:5
5:1
5:2
5:3
5:4
5:5
6:1
6:2
6:3
6:4
6:5
7:1
7:2
7:3
7:4
7:5
8:1
8:2
8:3
8:4
8:5
9:1
9:2
9:3
9:4
9:5
10:1
10:2
10:3
10:4
10:5
11:1
11:2
11:3
11:4
11:5
12:1
12:2
12:3
12:4
12:5
13:1
13:2
13:3
13:4
13:5
14:1
14:2
14:3
14:4
14:5
15:1
15:2
15:3
15:4
15:5
16:1
16:2
16:3
16:4
16:5
17:1
17:2
17:3
17:4
17:5
18:1
18:2
18:3
18:4
18:5
19:1
19:2
19:3
19:4
19:5
20:1
20:2
20:3
20:4
20:5
21:1
21:2
21:3
21:4
21:5
22:1
22:2
22:3
22:4
22:5
23:1
23:2
23:3
23:4
23:5
24:1
24:2
24:3
24:4
24:5
25:1
25:2
25:3
25:4
25:5
26:1
26:2
26:3
26:4
26:5
27:1
27:2
27:3
27:4
27:5
28:1
28:2
28:3
28:4
28:5
29:1
29:2
29:3
29:4
29:5
30:1
30:2
30:3
30:4
30:5
31:1
31:2
31:3
31:4
31:5
32:1
32:2
32:3
32:4
32:5
33:1
33:2
33:3
33:4
33:5
34:1
34:2
34:3
34:4
34:5
35:1
35:2
35:3
35:4
35:5
36:1
36:2
36:3
36:4
36:5
37:1
37:2
37:3
37:4
37:5
38:1
38:2
38:3
38:4
38:5
39:1
39:2
39:3
39:4
39:5
40:1
40:2
40:3
40:4
40:5
41:1
41:2
41:3
41:4
41:5
42:1
42:2
42:3
42:4
42:5
43:1
43:2
43:3
43:4
43:5
44:1
44:2
44:3
44:4
44:5
45:1
45:2
45:3
45:4
45:5
46:1
46:2
46:3
46:4
46:5
47:1
47:2
47:3
47:4
47:5
48:1
48:2
48:3
48:4
48:5
49:1
49:2
49:3
49:4
49:5
50:1
50:2
50:3
50:4
50:5
51:1
51:2
51:3
51:4
51:5
52:1
52:2
52:3
52:4
52:5
53:1
53:2
53:3
53:4
53:5
54:1
54:2
54:3
54:4
54:5
55:1
55:2
55:3
55:4
55:5
56:1
56:2
56:3
56:4
56:5
57:1
57:2
57:3
57:4
57:5
58:1
58:2
58:3
58:4
58:5
59:1
59:2
59:3
59:4
59:5
60:1
60:2
60:3
60:4
60:5
61:1
61:2
61:3
61:4
61:5
62:1
62:2
62:3
62:4
62:5
63:1
63:2
63:3
63:4
63:5
64:1
64:2
64:3
64:4
64:5
65:1
65:2
65:3
65:4
65:5
66:1
66:2
66:3
66:4
66:5
67:1
67:2
67:3
67:4
67:5
68:1
68:2
68:3
68:4
68:5
? matrix dunit={dlist}
Data error
Error executing script: halting
matrix dunit={dlist}
gretl version 1.7.5
Current session: 2008/07/02 18:59
? set shelldir "F:\Work_Files\CCPR_Paper\Gretl\"
? open LEC_data_gretl.gdt
Read datafile F:\Work_Files\CCPR_Paper\Gretl\LEC_data_gretl.gdt
periodicity: 12, maxobs: 816,
observations range: 1:01-68:12
Listing 113 variables:
0) const 1) year 2) cref
3) cosa 4) company 5) parent
6) state 7) vdummy 8) nyne
9) midwest 10) south 11) rboc
12) wacc 13) coc 14) switched_lines
15) leased_lines 16) total_sheath 17) aerial_sheath
18) non_aerial_shea 19) fibre_prop 20) duct
21) local_switch_mi 22) main_switch_min 23) swminutes
24) central_switche 25) remote_switches 26) busres
27) opc_total 28) opc_network 29) opc_network_xpp
30) tc_total 31) tc_network 32) tc_network_xpp
33) tz_total 34) tz_network 35) tz_network_xpp
36) sname 37) popden 38) population_dist
39) ilec_total_07 40) clec_resale_07 41) clec_une_07
42) clec_own_07 43) clec_total_07 44) lec_total_07
45) wireless_subs_0 46) pop_06 47) hhinc_06
48) popden_00 49) popden_06 50) pct_rur
51) pop_fb_05 52) pct_pop_fb 53) pct_cbsa_06
54) pct_nomove_05 55) pct_union_all 56) pct_union_priv
57) pct_forest 58) price_elect 59) clec_zip_0
60) clec_zip_1to4 61) clec_zip_5plus 62) clec_pct_lines
63) clec_own_lines 64) clec_une_lines 65) clec_resale_lin
66) clec_total_line 67) time 68) avpay
69) cost_aerial 70) cost_under 71) cost_duct
72) ln_gdpdef 73) ln_ncost 74) ln_opcost
75) ln_sl 76) ln_ll 77) ln_swmins
78) ln_ts 79) ln_duct 80) ln_popden
81) ln_fibreprop 82) ln_busres 83) ln_pay
84) ln_caerial 85) ln_cunder 86) ln_cduct
87) ln_ncostpsl 88) ln_opcostpsl 89) ln_llpsl
90) ln_swminpsl 91) ln_tspsl 92) ln_ductpsl
93) d9699 94) ln_sl_d9699 95) ln_ll_d9699
96) ln_swmins_d9699 97) ln_ts_d9699 98) ln_duct_d9699
99) ln_llpsl_d9699 100) ln_swminpsl_d96 101) ln_tspsl_d9699
102) ln_ductpsl_d969 103) time_d9699 104) scale
105) scalesq 106) pct_clec_total 107) pct_clec_own
108) mobilepsl 109) peak_lines 110) strand
111) ln_strand 112) ln_strand_lag
# smpl year > 2002 --restrict
# set debug 2
# set bfgs_toler 0.000001
? series panel_unit=$unit
Generated series panel_unit (ID 113)
# series panel_time=$period
? print panel_unit -o
Obs panel_unit
1:01 1
1:02 1
1:03 1
1:04 1
1:05 1
1:06 1
1:07 1
1:08 1
1:09 1
1:10 1
1:11 1
1:12 1
2:01 2
2:02 2
2:03 2
2:04 2
2:05 2
2:06 2
2:07 2
2:08 2
2:09 2
2:10 2
2:11 2
2:12 2
3:01 3
3:02 3
3:03 3
3:04 3
3:05 3
3:06 3
3:07 3
3:08 3
3:09 3
3:10 3
3:11 3
3:12 3
4:01 4
4:02 4
4:03 4
4:04 4
4:05 4
4:06 4
4:07 4
4:08 4
4:09 4
4:10 4
4:11 4
4:12 4
5:01 5
5:02 5
5:03 5
5:04 5
5:05 5
5:06 5
5:07 5
5:08 5
5:09 5
5:10 5
5:11 5
5:12 5
6:01 6
6:02 6
6:03 6
6:04 6
6:05 6
6:06 6
6:07 6
6:08 6
6:09 6
6:10 6
6:11 6
6:12 6
7:01 7
7:02 7
7:03 7
7:04 7
7:05 7
7:06 7
7:07 7
7:08 7
7:09 7
7:10 7
7:11 7
7:12 7
8:01 8
8:02 8
8:03 8
8:04 8
8:05 8
8:06 8
8:07 8
8:08 8
8:09 8
8:10 8
8:11 8
8:12 8
9:01 9
9:02 9
9:03 9
9:04 9
9:05 9
9:06 9
9:07 9
9:08 9
9:09 9
9:10 9
9:11 9
9:12 9
10:01 10
10:02 10
10:03 10
10:04 10
10:05 10
10:06 10
10:07 10
10:08 10
10:09 10
10:10 10
10:11 10
10:12 10
11:01 11
11:02 11
11:03 11
11:04 11
11:05 11
11:06 11
11:07 11
11:08 11
11:09 11
11:10 11
11:11 11
11:12 11
12:01 12
12:02 12
12:03 12
12:04 12
12:05 12
12:06 12
12:07 12
12:08 12
12:09 12
12:10 12
12:11 12
12:12 12
13:01 13
13:02 13
13:03 13
13:04 13
13:05 13
13:06 13
13:07 13
13:08 13
13:09 13
13:10 13
13:11 13
13:12 13
14:01 14
14:02 14
14:03 14
14:04 14
14:05 14
14:06 14
14:07 14
14:08 14
14:09 14
14:10 14
14:11 14
14:12 14
15:01 15
15:02 15
15:03 15
15:04 15
15:05 15
15:06 15
15:07 15
15:08 15
15:09 15
15:10 15
15:11 15
15:12 15
16:01 16
16:02 16
16:03 16
16:04 16
16:05 16
16:06 16
16:07 16
16:08 16
16:09 16
16:10 16
16:11 16
16:12 16
17:01 17
17:02 17
17:03 17
17:04 17
17:05 17
17:06 17
17:07 17
17:08 17
17:09 17
17:10 17
17:11 17
17:12 17
18:01 18
18:02 18
18:03 18
18:04 18
18:05 18
18:06 18
18:07 18
18:08 18
18:09 18
18:10 18
18:11 18
18:12 18
19:01 19
19:02 19
19:03 19
19:04 19
19:05 19
19:06 19
19:07 19
19:08 19
19:09 19
19:10 19
19:11 19
19:12 19
20:01 20
20:02 20
20:03 20
20:04 20
20:05 20
20:06 20
20:07 20
20:08 20
20:09 20
20:10 20
20:11 20
20:12 20
21:01 21
21:02 21
21:03 21
21:04 21
21:05 21
21:06 21
21:07 21
21:08 21
21:09 21
21:10 21
21:11 21
21:12 21
22:01 22
22:02 22
22:03 22
22:04 22
22:05 22
22:06 22
22:07 22
22:08 22
22:09 22
22:10 22
22:11 22
22:12 22
23:01 23
23:02 23
23:03 23
23:04 23
23:05 23
23:06 23
23:07 23
23:08 23
23:09 23
23:10 23
23:11 23
23:12 23
24:01 24
24:02 24
24:03 24
24:04 24
24:05 24
24:06 24
24:07 24
24:08 24
24:09 24
24:10 24
24:11 24
24:12 24
25:01 25
25:02 25
25:03 25
25:04 25
25:05 25
25:06 25
25:07 25
25:08 25
25:09 25
25:10 25
25:11 25
25:12 25
26:01 26
26:02 26
26:03 26
26:04 26
26:05 26
26:06 26
26:07 26
26:08 26
26:09 26
26:10 26
26:11 26
26:12 26
27:01 27
27:02 27
27:03 27
27:04 27
27:05 27
27:06 27
27:07 27
27:08 27
27:09 27
27:10 27
27:11 27
27:12 27
28:01 28
28:02 28
28:03 28
28:04 28
28:05 28
28:06 28
28:07 28
28:08 28
28:09 28
28:10 28
28:11 28
28:12 28
29:01 29
29:02 29
29:03 29
29:04 29
29:05 29
29:06 29
29:07 29
29:08 29
29:09 29
29:10 29
29:11 29
29:12 29
30:01 30
30:02 30
30:03 30
30:04 30
30:05 30
30:06 30
30:07 30
30:08 30
30:09 30
30:10 30
30:11 30
30:12 30
31:01 31
31:02 31
31:03 31
31:04 31
31:05 31
31:06 31
31:07 31
31:08 31
31:09 31
31:10 31
31:11 31
31:12 31
32:01 32
32:02 32
32:03 32
32:04 32
32:05 32
32:06 32
32:07 32
32:08 32
32:09 32
32:10 32
32:11 32
32:12 32
33:01 33
33:02 33
33:03 33
33:04 33
33:05 33
33:06 33
33:07 33
33:08 33
33:09 33
33:10 33
33:11 33
33:12 33
34:01 34
34:02 34
34:03 34
34:04 34
34:05 34
34:06 34
34:07 34
34:08 34
34:09 34
34:10 34
34:11 34
34:12 34
35:01 35
35:02 35
35:03 35
35:04 35
35:05 35
35:06 35
35:07 35
35:08 35
35:09 35
35:10 35
35:11 35
35:12 35
36:01 36
36:02 36
36:03 36
36:04 36
36:05 36
36:06 36
36:07 36
36:08 36
36:09 36
36:10 36
36:11 36
36:12 36
37:01 37
37:02 37
37:03 37
37:04 37
37:05 37
37:06 37
37:07 37
37:08 37
37:09 37
37:10 37
37:11 37
37:12 37
38:01 38
38:02 38
38:03 38
38:04 38
38:05 38
38:06 38
38:07 38
38:08 38
38:09 38
38:10 38
38:11 38
38:12 38
39:01 39
39:02 39
39:03 39
39:04 39
39:05 39
39:06 39
39:07 39
39:08 39
39:09 39
39:10 39
39:11 39
39:12 39
40:01 40
40:02 40
40:03 40
40:04 40
40:05 40
40:06 40
40:07 40
40:08 40
40:09 40
40:10 40
40:11 40
40:12 40
41:01 41
41:02 41
41:03 41
41:04 41
41:05 41
41:06 41
41:07 41
41:08 41
41:09 41
41:10 41
41:11 41
41:12 41
42:01 42
42:02 42
42:03 42
42:04 42
42:05 42
42:06 42
42:07 42
42:08 42
42:09 42
42:10 42
42:11 42
42:12 42
43:01 43
43:02 43
43:03 43
43:04 43
43:05 43
43:06 43
43:07 43
43:08 43
43:09 43
43:10 43
43:11 43
43:12 43
44:01 44
44:02 44
44:03 44
44:04 44
44:05 44
44:06 44
44:07 44
44:08 44
44:09 44
44:10 44
44:11 44
44:12 44
45:01 45
45:02 45
45:03 45
45:04 45
45:05 45
45:06 45
45:07 45
45:08 45
45:09 45
45:10 45
45:11 45
45:12 45
46:01 46
46:02 46
46:03 46
46:04 46
46:05 46
46:06 46
46:07 46
46:08 46
46:09 46
46:10 46
46:11 46
46:12 46
47:01 47
47:02 47
47:03 47
47:04 47
47:05 47
47:06 47
47:07 47
47:08 47
47:09 47
47:10 47
47:11 47
47:12 47
48:01 48
48:02 48
48:03 48
48:04 48
48:05 48
48:06 48
48:07 48
48:08 48
48:09 48
48:10 48
48:11 48
48:12 48
49:01 49
49:02 49
49:03 49
49:04 49
49:05 49
49:06 49
49:07 49
49:08 49
49:09 49
49:10 49
49:11 49
49:12 49
50:01 50
50:02 50
50:03 50
50:04 50
50:05 50
50:06 50
50:07 50
50:08 50
50:09 50
50:10 50
50:11 50
50:12 50
51:01 51
51:02 51
51:03 51
51:04 51
51:05 51
51:06 51
51:07 51
51:08 51
51:09 51
51:10 51
51:11 51
51:12 51
52:01 52
52:02 52
52:03 52
52:04 52
52:05 52
52:06 52
52:07 52
52:08 52
52:09 52
52:10 52
52:11 52
52:12 52
53:01 53
53:02 53
53:03 53
53:04 53
53:05 53
53:06 53
53:07 53
53:08 53
53:09 53
53:10 53
53:11 53
53:12 53
54:01 54
54:02 54
54:03 54
54:04 54
54:05 54
54:06 54
54:07 54
54:08 54
54:09 54
54:10 54
54:11 54
54:12 54
55:01 55
55:02 55
55:03 55
55:04 55
55:05 55
55:06 55
55:07 55
55:08 55
55:09 55
55:10 55
55:11 55
55:12 55
56:01 56
56:02 56
56:03 56
56:04 56
56:05 56
56:06 56
56:07 56
56:08 56
56:09 56
56:10 56
56:11 56
56:12 56
57:01 57
57:02 57
57:03 57
57:04 57
57:05 57
57:06 57
57:07 57
57:08 57
57:09 57
57:10 57
57:11 57
57:12 57
58:01 58
58:02 58
58:03 58
58:04 58
58:05 58
58:06 58
58:07 58
58:08 58
58:09 58
58:10 58
58:11 58
58:12 58
59:01 59
59:02 59
59:03 59
59:04 59
59:05 59
59:06 59
59:07 59
59:08 59
59:09 59
59:10 59
59:11 59
59:12 59
60:01 60
60:02 60
60:03 60
60:04 60
60:05 60
60:06 60
60:07 60
60:08 60
60:09 60
60:10 60
60:11 60
60:12 60
61:01 61
61:02 61
61:03 61
61:04 61
61:05 61
61:06 61
61:07 61
61:08 61
61:09 61
61:10 61
61:11 61
61:12 61
62:01 62
62:02 62
62:03 62
62:04 62
62:05 62
62:06 62
62:07 62
62:08 62
62:09 62
62:10 62
62:11 62
62:12 62
63:01 63
63:02 63
63:03 63
63:04 63
63:05 63
63:06 63
63:07 63
63:08 63
63:09 63
63:10 63
63:11 63
63:12 63
64:01 64
64:02 64
64:03 64
64:04 64
64:05 64
64:06 64
64:07 64
64:08 64
64:09 64
64:10 64
64:11 64
64:12 64
65:01 65
65:02 65
65:03 65
65:04 65
65:05 65
65:06 65
65:07 65
65:08 65
65:09 65
65:10 65
65:11 65
65:12 65
66:01 66
66:02 66
66:03 66
66:04 66
66:05 66
66:06 66
66:07 66
66:08 66
66:09 66
66:10 66
66:11 66
66:12 66
67:01 67
67:02 67
67:03 67
67:04 67
67:05 67
67:06 67
67:07 67
67:08 67
67:09 67
67:10 67
67:11 67
67:12 67
68:01 68
68:02 68
68:03 68
68:04 68
68:05 68
68:06 68
68:07 68
68:08 68
68:09 68
68:10 68
68:11 68
68:12 68
? genr unitdum
Panel dummy variables generated.
? list dlist = du_*
Generated list dlist
? print du_1 -o
Obs du_1
1:01 1
1:02 1
1:03 1
1:04 1
1:05 1
1:06 1
1:07 1
1:08 1
1:09 1
1:10 1
1:11 1
1:12 1
2:01 0
2:02 0
2:03 0
2:04 0
2:05 0
2:06 0
2:07 0
2:08 0
2:09 0
2:10 0
2:11 0
2:12 0
3:01 0
3:02 0
3:03 0
3:04 0
3:05 0
3:06 0
3:07 0
3:08 0
3:09 0
3:10 0
3:11 0
3:12 0
4:01 0
4:02 0
4:03 0
4:04 0
4:05 0
4:06 0
4:07 0
4:08 0
4:09 0
4:10 0
4:11 0
4:12 0
5:01 0
5:02 0
5:03 0
5:04 0
5:05 0
5:06 0
5:07 0
5:08 0
5:09 0
5:10 0
5:11 0
5:12 0
6:01 0
6:02 0
6:03 0
6:04 0
6:05 0
6:06 0
6:07 0
6:08 0
6:09 0
6:10 0
6:11 0
6:12 0
7:01 0
7:02 0
7:03 0
7:04 0
7:05 0
7:06 0
7:07 0
7:08 0
7:09 0
7:10 0
7:11 0
7:12 0
8:01 0
8:02 0
8:03 0
8:04 0
8:05 0
8:06 0
8:07 0
8:08 0
8:09 0
8:10 0
8:11 0
8:12 0
9:01 0
9:02 0
9:03 0
9:04 0
9:05 0
9:06 0
9:07 0
9:08 0
9:09 0
9:10 0
9:11 0
9:12 0
10:01 0
10:02 0
10:03 0
10:04 0
10:05 0
10:06 0
10:07 0
10:08 0
10:09 0
10:10 0
10:11 0
10:12 0
11:01 0
11:02 0
11:03 0
11:04 0
11:05 0
11:06 0
11:07 0
11:08 0
11:09 0
11:10 0
11:11 0
11:12 0
12:01 0
12:02 0
12:03 0
12:04 0
12:05 0
12:06 0
12:07 0
12:08 0
12:09 0
12:10 0
12:11 0
12:12 0
13:01 0
13:02 0
13:03 0
13:04 0
13:05 0
13:06 0
13:07 0
13:08 0
13:09 0
13:10 0
13:11 0
13:12 0
14:01 0
14:02 0
14:03 0
14:04 0
14:05 0
14:06 0
14:07 0
14:08 0
14:09 0
14:10 0
14:11 0
14:12 0
15:01 0
15:02 0
15:03 0
15:04 0
15:05 0
15:06 0
15:07 0
15:08 0
15:09 0
15:10 0
15:11 0
15:12 0
16:01 0
16:02 0
16:03 0
16:04 0
16:05 0
16:06 0
16:07 0
16:08 0
16:09 0
16:10 0
16:11 0
16:12 0
17:01 0
17:02 0
17:03 0
17:04 0
17:05 0
17:06 0
17:07 0
17:08 0
17:09 0
17:10 0
17:11 0
17:12 0
18:01 0
18:02 0
18:03 0
18:04 0
18:05 0
18:06 0
18:07 0
18:08 0
18:09 0
18:10 0
18:11 0
18:12 0
19:01 0
19:02 0
19:03 0
19:04 0
19:05 0
19:06 0
19:07 0
19:08 0
19:09 0
19:10 0
19:11 0
19:12 0
20:01 0
20:02 0
20:03 0
20:04 0
20:05 0
20:06 0
20:07 0
20:08 0
20:09 0
20:10 0
20:11 0
20:12 0
21:01 0
21:02 0
21:03 0
21:04 0
21:05 0
21:06 0
21:07 0
21:08 0
21:09 0
21:10 0
21:11 0
21:12 0
22:01 0
22:02 0
22:03 0
22:04 0
22:05 0
22:06 0
22:07 0
22:08 0
22:09 0
22:10 0
22:11 0
22:12 0
23:01 0
23:02 0
23:03 0
23:04 0
23:05 0
23:06 0
23:07 0
23:08 0
23:09 0
23:10 0
23:11 0
23:12 0
24:01 0
24:02 0
24:03 0
24:04 0
24:05 0
24:06 0
24:07 0
24:08 0
24:09 0
24:10 0
24:11 0
24:12 0
25:01 0
25:02 0
25:03 0
25:04 0
25:05 0
25:06 0
25:07 0
25:08 0
25:09 0
25:10 0
25:11 0
25:12 0
26:01 0
26:02 0
26:03 0
26:04 0
26:05 0
26:06 0
26:07 0
26:08 0
26:09 0
26:10 0
26:11 0
26:12 0
27:01 0
27:02 0
27:03 0
27:04 0
27:05 0
27:06 0
27:07 0
27:08 0
27:09 0
27:10 0
27:11 0
27:12 0
28:01 0
28:02 0
28:03 0
28:04 0
28:05 0
28:06 0
28:07 0
28:08 0
28:09 0
28:10 0
28:11 0
28:12 0
29:01 0
29:02 0
29:03 0
29:04 0
29:05 0
29:06 0
29:07 0
29:08 0
29:09 0
29:10 0
29:11 0
29:12 0
30:01 0
30:02 0
30:03 0
30:04 0
30:05 0
30:06 0
30:07 0
30:08 0
30:09 0
30:10 0
30:11 0
30:12 0
31:01 0
31:02 0
31:03 0
31:04 0
31:05 0
31:06 0
31:07 0
31:08 0
31:09 0
31:10 0
31:11 0
31:12 0
32:01 0
32:02 0
32:03 0
32:04 0
32:05 0
32:06 0
32:07 0
32:08 0
32:09 0
32:10 0
32:11 0
32:12 0
33:01 0
33:02 0
33:03 0
33:04 0
33:05 0
33:06 0
33:07 0
33:08 0
33:09 0
33:10 0
33:11 0
33:12 0
34:01 0
34:02 0
34:03 0
34:04 0
34:05 0
34:06 0
34:07 0
34:08 0
34:09 0
34:10 0
34:11 0
34:12 0
35:01 0
35:02 0
35:03 0
35:04 0
35:05 0
35:06 0
35:07 0
35:08 0
35:09 0
35:10 0
35:11 0
35:12 0
36:01 0
36:02 0
36:03 0
36:04 0
36:05 0
36:06 0
36:07 0
36:08 0
36:09 0
36:10 0
36:11 0
36:12 0
37:01 0
37:02 0
37:03 0
37:04 0
37:05 0
37:06 0
37:07 0
37:08 0
37:09 0
37:10 0
37:11 0
37:12 0
38:01 0
38:02 0
38:03 0
38:04 0
38:05 0
38:06 0
38:07 0
38:08 0
38:09 0
38:10 0
38:11 0
38:12 0
39:01 0
39:02 0
39:03 0
39:04 0
39:05 0
39:06 0
39:07 0
39:08 0
39:09 0
39:10 0
39:11 0
39:12 0
40:01 0
40:02 0
40:03 0
40:04 0
40:05 0
40:06 0
40:07 0
40:08 0
40:09 0
40:10 0
40:11 0
40:12 0
41:01 0
41:02 0
41:03 0
41:04 0
41:05 0
41:06 0
41:07 0
41:08 0
41:09 0
41:10 0
41:11 0
41:12 0
42:01 0
42:02 0
42:03 0
42:04 0
42:05 0
42:06 0
42:07 0
42:08 0
42:09 0
42:10 0
42:11 0
42:12 0
43:01 0
43:02 0
43:03 0
43:04 0
43:05 0
43:06 0
43:07 0
43:08 0
43:09 0
43:10 0
43:11 0
43:12 0
44:01 0
44:02 0
44:03 0
44:04 0
44:05 0
44:06 0
44:07 0
44:08 0
44:09 0
44:10 0
44:11 0
44:12 0
45:01 0
45:02 0
45:03 0
45:04 0
45:05 0
45:06 0
45:07 0
45:08 0
45:09 0
45:10 0
45:11 0
45:12 0
46:01 0
46:02 0
46:03 0
46:04 0
46:05 0
46:06 0
46:07 0
46:08 0
46:09 0
46:10 0
46:11 0
46:12 0
47:01 0
47:02 0
47:03 0
47:04 0
47:05 0
47:06 0
47:07 0
47:08 0
47:09 0
47:10 0
47:11 0
47:12 0
48:01 0
48:02 0
48:03 0
48:04 0
48:05 0
48:06 0
48:07 0
48:08 0
48:09 0
48:10 0
48:11 0
48:12 0
49:01 0
49:02 0
49:03 0
49:04 0
49:05 0
49:06 0
49:07 0
49:08 0
49:09 0
49:10 0
49:11 0
49:12 0
50:01 0
50:02 0
50:03 0
50:04 0
50:05 0
50:06 0
50:07 0
50:08 0
50:09 0
50:10 0
50:11 0
50:12 0
51:01 0
51:02 0
51:03 0
51:04 0
51:05 0
51:06 0
51:07 0
51:08 0
51:09 0
51:10 0
51:11 0
51:12 0
52:01 0
52:02 0
52:03 0
52:04 0
52:05 0
52:06 0
52:07 0
52:08 0
52:09 0
52:10 0
52:11 0
52:12 0
53:01 0
53:02 0
53:03 0
53:04 0
53:05 0
53:06 0
53:07 0
53:08 0
53:09 0
53:10 0
53:11 0
53:12 0
54:01 0
54:02 0
54:03 0
54:04 0
54:05 0
54:06 0
54:07 0
54:08 0
54:09 0
54:10 0
54:11 0
54:12 0
55:01 0
55:02 0
55:03 0
55:04 0
55:05 0
55:06 0
55:07 0
55:08 0
55:09 0
55:10 0
55:11 0
55:12 0
56:01 0
56:02 0
56:03 0
56:04 0
56:05 0
56:06 0
56:07 0
56:08 0
56:09 0
56:10 0
56:11 0
56:12 0
57:01 0
57:02 0
57:03 0
57:04 0
57:05 0
57:06 0
57:07 0
57:08 0
57:09 0
57:10 0
57:11 0
57:12 0
58:01 0
58:02 0
58:03 0
58:04 0
58:05 0
58:06 0
58:07 0
58:08 0
58:09 0
58:10 0
58:11 0
58:12 0
59:01 0
59:02 0
59:03 0
59:04 0
59:05 0
59:06 0
59:07 0
59:08 0
59:09 0
59:10 0
59:11 0
59:12 0
60:01 0
60:02 0
60:03 0
60:04 0
60:05 0
60:06 0
60:07 0
60:08 0
60:09 0
60:10 0
60:11 0
60:12 0
61:01 0
61:02 0
61:03 0
61:04 0
61:05 0
61:06 0
61:07 0
61:08 0
61:09 0
61:10 0
61:11 0
61:12 0
62:01 0
62:02 0
62:03 0
62:04 0
62:05 0
62:06 0
62:07 0
62:08 0
62:09 0
62:10 0
62:11 0
62:12 0
63:01 0
63:02 0
63:03 0
63:04 0
63:05 0
63:06 0
63:07 0
63:08 0
63:09 0
63:10 0
63:11 0
63:12 0
64:01 0
64:02 0
64:03 0
64:04 0
64:05 0
64:06 0
64:07 0
64:08 0
64:09 0
64:10 0
64:11 0
64:12 0
65:01 0
65:02 0
65:03 0
65:04 0
65:05 0
65:06 0
65:07 0
65:08 0
65:09 0
65:10 0
65:11 0
65:12 0
66:01 0
66:02 0
66:03 0
66:04 0
66:05 0
66:06 0
66:07 0
66:08 0
66:09 0
66:10 0
66:11 0
66:12 0
67:01 0
67:02 0
67:03 0
67:04 0
67:05 0
67:06 0
67:07 0
67:08 0
67:09 0
67:10 0
67:11 0
67:12 0
68:01 0
68:02 0
68:03 0
68:04 0
68:05 0
68:06 0
68:07 0
68:08 0
68:09 0
68:10 0
68:11 0
68:12 0
? matrix dunit={dlist}
Generated matrix dunit
? matrix pnumb=sumc(dunit)
Generated matrix pnumb
? print pnumb
pnumb (1 x 68)
12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12
12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12
12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12
12 12 12 12 12 12 12 12 12 12 12 12 12 12
? matrix dnumb=dunit*transp(pnumb)
Generated matrix dnumb
? matrix dnper=sumr(dnumb)
Generated matrix dnper
? print dnper
dnper (816 x 1)
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
12
# series y = ln_ncostpsl
# list xlist = const ln_strand_lag ln_tspsl ln_cunder time
# matrix sfa_result=sfa_mod(y, xlist, "t", "C")
# list efflist=sfa_eff(y, xlist, sfa_result, "t", "C")
# print sfa_result
# print y efflist --byobs