LISTSERV at the University of Georgia
Menubar Imagemap
Home Browse Manage Request Manuals Register
Previous messageNext messagePrevious in topicNext in topicPrevious by same authorNext by same authorPrevious page (July 2004, week 3)Back to main SAS-L pageJoin or leave SAS-L (or change settings)ReplyPost a new messageSearchProportional fontNon-proportional font
Date:         Mon, 19 Jul 2004 13:43:13 +0200
Reply-To:     Christine Wieck <wieck@AGP.UNI-BONN.DE>
Sender:       "SAS(r) Discussion" <SAS-L@LISTSERV.UGA.EDU>
From:         Christine Wieck <wieck@AGP.UNI-BONN.DE>
Subject:      Better program code
Content-Type: text/plain; charset="ISO-8859-1"


Hello SAS experts,


I want to conduct a correction for heteroscedasticity in a system estimation according to a test proposed in Greene p.224. Basic idea is to identify in a first step by means of the Breusch-Pagan test the variable(s) responsible for the heteroscedasticity. Second, do a regression of the transformed errors from the first step on the variables responsible for heteroscedasticity and third, use the predicted values from that regression as weights in the model from the first step. Then, the heteroscedasticity should be vanished.


I found a suitable way to do it, but there are several shortcomings in the code, as e.g. I wasn't able to address some elements from an ODS data set.

So maybe someone can have a look on the code, and give me a hint, how to do it in a better way.....


Please find below the program code and data. The real model works for a system with 3 equations and 77 variables and some 2000 observations.


Thanks very much for your help,
Christine


P.S. Thanks to all who answered to my question about number of digits and condition numbers.




-----------------------------------------------------------------------


/* --- July 2004 */
/* Correct for heteroscedasticity in model */
/* ------- File names for export + general options --------------- */
*FILENAME in "u:\wieck\coest\data\sas\try.csv";
OPTIONS ps=72 ls=100;
/* ------- Import CSV-data ------------------------------------ */
DATA temp;
* INFILE in dlm=',' dsd truncover LRECL=20000;
LENGTH NOBS $ 8;
INPUT NOBS YFC YFA YFT YFOT YFOA YFMI ZFLA ZFDC PNC PNA PNT;
datalines;
1 2.16E+01 9.94E+01 4.06E+01 -6.31E+00 -4.12E+01 2.94E-01 1.10E+01 6.89E-02 1.32E-02 -1.02E-03 -9.41E-03
2 1.08E+01 1.27E+01 8.72E+00 8.93E+01 -3.85E+01 4.16E-02 -2.66E+00 4.58E-03 -1.39E-02 9.48E-03 6.26E-03
3 -4.05E+00 -1.56E+01 -5.06E+00 8.34E+01 -3.46E+01 2.93E-02 7.36E-01 7.06E-03 -3.14E-02 4.84E-03 2.14E-02
4 -4.13E+00 -3.75E+00 -3.04E+00 -3.78E+01 7.54E+00 4.15E-02 1.85E+00 1.95E-02 8.46E-03 1.78E-03 -7.10E-03
5 7.75E+00 -8.50E+00 2.31E-01 -2.72E+01 1.08E+01 -6.58E-03 1.31E+00 4.12E-03 1.77E-02 -8.57E-03 -9.47E-03
6 9.71E+00 1.57E+01 3.81E+00 2.45E+01 -2.25E+01 9.22E-02 4.51E+00 2.01E-02 6.67E-03 2.83E-03 -6.21E-03
7 -7.10E+00 -1.41E+00 -7.51E+00 1.90E+01 -3.51E+00 -2.28E-02 2.92E-01 -1.05E-02 -1.26E-02 8.73E-03 5.61E-03
8 -1.72E+01 -4.28E+01 -1.86E+01 -1.03E+02 6.99E+01 -2.21E-01 -8.27E+00 -5.33E-02 7.62E-03 -5.15E-03 -3.45E-03
9 -1.74E+01 -5.57E+01 -1.91E+01 -4.23E+01 5.21E+01 -2.48E-01 -8.82E+00 -6.05E-02 4.31E-03 -1.29E-02 2.41E-03
10 1.27E+01 9.54E+01 5.11E-01 -1.16E+01 -9.54E+00 9.27E-02 1.14E+01 4.90E-02 6.10E-03 6.91E-04 -4.86E-03
11 -2.09E+00 1.59E+01 1.08E+01 2.93E+01 -2.28E+01 8.21E-02 3.01E+00 -2.66E-03 -1.32E-02 6.60E-03 6.96E-03
12 -1.06E+01 -1.11E+02 -1.13E+01 -1.77E+01 3.23E+01 -1.75E-01 -1.45E+01 -4.63E-02 7.06E-03 -7.29E-03 -2.10E-03
13 1.41E+00 -4.97E+01 -8.76E+00 -5.25E+00 -9.41E+00 7.63E-02 4.77E+00 -1.55E-02 6.79E-03 -1.92E-03 -4.23E-03
14 -1.68E+00 8.04E+01 -1.48E+01 4.52E+00 -2.08E+01 1.29E-01 6.33E+00 5.01E-02 -2.03E-02 8.58E-03 1.14E-02
15 -1.35E+00 -5.18E+01 -1.56E+01 -1.72E+00 2.92E+01 -1.93E-01 -9.26E+00 -3.99E-02 2.04E-03 8.82E-04 -1.91E-03
16 -2.71E-01 2.38E+01 -7.59E-01 -3.00E+00 1.20E+01 -7.37E-02 -4.98E+00 -1.93E-02 1.12E-02 -9.47E-03 -4.28E-03
17 1.89E+00 -2.71E+00 3.99E+01 5.46E+00 -1.10E+01 6.06E-02 3.13E+00 2.46E-02 2.50E-04 1.93E-03 -1.03E-03
18 8.11E+00 4.62E+01 2.44E+00 0 3.78E+00 -2.55E-02 -1.98E+00 -2.35E-02 2.39E-02 -5.46E-03 -1.55E-02
19 -3.27E+00 -2.23E+01 -1.23E+01 0 1.09E+01 -7.37E-02 -1.03E+00 -2.52E-02 -3.20E-03 5.05E-03 1.93E-04
20 -4.83E+00 -2.39E+01 9.83E+00 0 -1.47E+01 9.92E-02 3.00E+00 4.87E-02 -2.07E-02 4.10E-04 1.53E-02
21 1.52E+01 1.77E+02 4.84E+00 -5.95E+01 0 1.45E-01 2.82E+00 4.03E-02 2.39E-02 -5.46E-03 -1.55E-02
22 -5.02E+00 -6.33E+01 1.15E+00 1.19E+01 0 -2.90E-02 -1.16E-01 -1.55E-02 -3.20E-03 5.05E-03 1.93E-04
23 -1.02E+01 -1.14E+02 -5.99E+00 4.76E+01 0 -1.16E-01 -2.71E+00 -2.48E-02 -2.07E-02 4.10E-04 1.53E-02
24 5.50E+00 2.11E+02 5.39E+01 -6.60E+01 0 1.61E-01 1.51E+01 9.40E-02 2.39E-02 -5.46E-03 -1.55E-02
25 4.38E+00 -8.46E+01 -2.52E+01 2.68E+01 0 -6.56E-02 -2.21E+00 -1.50E-02 -3.20E-03 5.05E-03 1.93E-04
26 -9.88E+00 -1.26E+02 -2.87E+01 3.92E+01 0 -9.59E-02 -1.29E+01 -7.90E-02 -2.07E-02 4.10E-04 1.53E-02
27 -9.95E-01 4.68E+00 -1.55E+00 -1.52E+01 0 3.72E-02 1.88E+00 2.54E-03 8.76E-03 2.32E-03 -7.55E-03
28 9.95E-01 -4.68E+00 1.55E+00 1.52E+01 0 -3.72E-02 -1.88E+00 -2.54E-03 -8.76E-03 -2.32E-03 7.55E-03
29 1.98E+00 5.79E+01 1.92E+01 -1.83E+01 -1.77E+00 5.67E-02 2.43E+00 -1.30E-02 1.36E-02 -5.25E-03 -7.84E-03
30 -1.98E+00 -5.79E+01 -1.92E+01 1.83E+01 1.77E+00 -5.67E-02 -2.43E+00 1.30E-02 -1.36E-02 5.25E-03 7.84E-03
31 1.45E+00 8.13E+01 -1.03E+00 0 0 0 -2.01E+00 -4.98E-03 2.39E-02 -5.46E-03 -1.55E-02
32 -5.35E-01 -1.18E+01 -9.37E+00 0 0 0 5.37E-01 -5.50E-03 -3.20E-03 5.05E-03 1.93E-04
33 -9.13E-01 -6.94E+01 1.04E+01 0 0 0 1.48E+00 1.05E-02 -2.07E-02 4.10E-04 1.53E-02
34 2.43E+01 7.33E+01 2.09E+01 -5.95E+01 0 1.45E-01 5.36E+00 2.57E-02 2.39E-02 -5.46E-03 -1.55E-02
35 -9.20E+00 -1.47E+00 -5.66E+00 -6.15E+00 0 1.50E-02 8.14E-01 4.10E-04 -3.20E-03 5.05E-03 1.93E-04
36 -1.51E+01 -7.18E+01 -1.52E+01 6.56E+01 0 -1.61E-01 -6.17E+00 -2.61E-02 -2.07E-02 4.10E-04 1.53E-02
37 1.27E+01 5.38E+01 1.50E+01 -4.38E+01 0 1.07E-01 5.39E+00 1.17E-02 2.39E-02 -5.46E-03 -1.55E-02
38 -7.78E+00 -2.86E+00 -7.66E+00 -1.10E+01 0 2.69E-02 -2.00E+00 -1.00E-02 -3.20E-03 5.05E-03 1.93E-04
39 -4.91E+00 -5.10E+01 -7.30E+00 5.47E+01 0 -1.34E-01 -3.40E+00 -1.68E-03 -2.07E-02 4.10E-04 1.53E-02
40 2.56E+01 4.74E+01 8.32E+00 -4.73E+01 0 1.16E-01 3.27E+00 1.26E-02 2.39E-02 -5.46E-03 -1.55E-02
41 -1.20E+01 -1.02E+01 -3.54E+00 -7.32E+00 0 1.79E-02 -1.58E+00 -2.15E-03 -3.20E-03 5.05E-03 1.93E-04
42 -1.35E+01 -3.71E+01 -4.78E+00 5.46E+01 0 -1.34E-01 -1.69E+00 -1.05E-02 -2.07E-02 4.10E-04 1.53E-02
43 -5.58E-01 1.50E+02 2.07E+01 -6.17E+01 0 1.51E-01 1.30E+00 1.99E-02 2.39E-02 -5.46E-03 -1.55E-02
44 2.77E+00 -6.87E+01 -7.09E+00 -5.07E+00 0 1.24E-02 1.32E+00 1.47E-03 -3.20E-03 5.05E-03 1.93E-04
45 -2.22E+00 -8.16E+01 -1.36E+01 6.68E+01 0 -1.63E-01 -2.62E+00 -2.14E-02 -2.07E-02 4.10E-04 1.53E-02
46 -1.95E-01 6.95E+01 1.05E+01 -4.93E+01 0 1.21E-01 7.87E+00 4.98E-03 2.39E-02 -5.46E-03 -1.55E-02
47 8.69E-01 -1.84E+01 -7.30E+00 1.20E+01 0 -2.93E-02 -3.51E+00 3.77E-03 -3.20E-03 5.05E-03 1.93E-04
48 -6.74E-01 -5.11E+01 -3.16E+00 3.73E+01 0 -9.13E-02 -4.35E+00 -8.75E-03 -2.07E-02 4.10E-04 1.53E-02
49 -3.31E+00 1.27E+02 1.40E+01 -5.04E+01 -5.36E+00 1.59E-01 1.10E+01 2.74E-02 2.39E-02 -5.46E-03 -1.55E-02
50 6.40E+00 -6.14E+01 2.91E+00 3.28E+01 -1.02E+01 -1.14E-02 -2.78E+00 2.16E-03 -3.20E-03 5.05E-03 1.93E-04
51 -3.08E+00 -6.57E+01 -1.69E+01 1.76E+01 1.56E+01 -1.48E-01 -8.23E+00 -2.95E-02 -2.07E-02 4.10E-04 1.53E-02
52 -3.92E+00 9.84E+01 1.61E+01 -4.35E+01 -1.11E+01 1.82E-01 8.23E+00 2.32E-02 2.39E-02 -5.46E-03 -1.55E-02
53 5.74E+00 -3.95E+01 -4.04E+00 2.84E+01 -7.17E+00 -2.12E-02 -1.88E+00 3.11E-03 -3.20E-03 5.05E-03 1.93E-04
54 -1.82E+00 -5.89E+01 -1.21E+01 1.51E+01 1.83E+01 -1.60E-01 -6.35E+00 -2.63E-02 -2.07E-02 4.10E-04 1.53E-02
55 3.20E+00 1.50E+02 3.39E+01 -4.31E+01 7.65E+00 5.39E-02 1.03E+01 4.35E-02 2.39E-02 -5.46E-03 -1.55E-02
56 2.86E+00 -7.54E+01 -1.40E+01 1.64E+01 1.12E+00 -4.78E-02 -2.20E+00 -9.72E-03 -3.20E-03 5.05E-03 1.93E-04
57 -6.06E+00 -7.49E+01 -1.99E+01 2.67E+01 -8.78E+00 -6.08E-03 -8.12E+00 -3.38E-02 -2.07E-02 4.10E-04 1.53E-02
58 1.74E+00 4.78E+01 5.52E+00 -4.44E+01 2.29E+00 9.31E-02 -1.79E-01 -5.06E-02 2.39E-02 -5.46E-03 -1.55E-02
59 -1.03E+00 -1.79E+01 -8.25E+00 3.11E+01 1.01E+01 -1.45E-01 -1.90E+00 -5.34E-03 -3.20E-03 5.05E-03 1.93E-04
60 -7.03E-01 -2.99E+01 2.74E+00 1.32E+01 -1.24E+01 5.15E-02 2.08E+00 5.59E-02 -2.07E-02 4.10E-04 1.53E-02
61 -1.07E+00 4.17E+01 1.23E+01 -2.71E+01 -5.78E+00 1.05E-01 1.88E+00 2.32E-02 2.39E-02 -5.46E-03 -1.55E-02
62 2.99E+00 1.56E+01 -2.63E+00 2.36E+01 -2.61E+01 1.18E-01 1.74E+00 1.73E-02 -3.20E-03 5.05E-03 1.93E-04
63 -1.92E+00 -5.73E+01 -9.68E+00 3.48E+00 3.19E+01 -2.24E-01 -3.62E+00 -4.06E-02 -2.07E-02 4.10E-04 1.53E-02
64 -4.44E+00 3.45E+01 4.00E+00 -2.00E+01 8.85E+00 -1.09E-02 4.39E-01 1.22E-02 1.36E-02 -5.25E-03 -7.84E-03
65 4.44E+00 -3.45E+01 -4.00E+00 2.00E+01 -8.85E+00 1.09E-02 -4.39E-01 -1.22E-02 -1.36E-02 5.25E-03 7.84E-03
66 -6.34E+00 3.08E+01 2.17E+01 -5.06E+01 -3.62E+00 1.48E-01 1.33E+01 3.02E-02 2.39E-02 -5.46E-03 -1.55E-02
67 9.92E+00 1.04E+01 1.12E+00 3.28E+01 -1.51E+01 2.13E-02 3.04E+00 -2.53E-03 -3.20E-03 5.05E-03 1.93E-04
68 -3.59E+00 -4.12E+01 -2.29E+01 1.77E+01 1.87E+01 -1.69E-01 -1.63E+01 -2.77E-02 -2.07E-02 4.10E-04 1.53E-02
69 4.86E+00 1.85E+00 1.11E+01 1.92E+01 -2.21E+01 1.02E-01 6.27E+00 2.09E-02 8.76E-03 2.32E-03 -7.55E-03
70 -4.86E+00 -1.85E+00 -1.11E+01 -1.92E+01 2.21E+01 -1.02E-01 -6.27E+00 -2.09E-02 -8.76E-03 -2.32E-03 7.55E-03
71 -2.48E+00 1.51E+02 1.93E+01 -4.39E+01 -8.14E+00 1.62E-01 1.22E+01 2.25E-02 2.39E-02 -5.46E-03 -1.55E-02
72 5.32E+00 -6.18E+01 -5.65E+00 3.45E+01 -1.68E+01 2.88E-02 -7.50E-01 4.51E-03 -3.20E-03 5.05E-03 1.93E-04
73 -2.84E+00 -8.91E+01 -1.36E+01 9.46E+00 2.49E+01 -1.91E-01 -1.15E+01 -2.70E-02 -2.07E-02 4.10E-04 1.53E-02
74 1.07E+01 1.11E+02 9.79E+00 -4.25E+01 -9.65E+00 1.69E-01 8.25E-01 2.71E-02 2.39E-02 -5.46E-03 -1.55E-02
75 5.23E-01 -1.82E+01 3.87E+00 3.87E+01 -1.79E+01 2.60E-02 3.08E+00 6.17E-03 -3.20E-03 5.05E-03 1.93E-04
76 -1.12E+01 -9.24E+01 -1.37E+01 3.88E+00 2.75E+01 -1.95E-01 -3.90E+00 -3.33E-02 -2.07E-02 4.10E-04 1.53E-02
77 -3.97E-01 9.53E+01 1.16E+01 -3.19E+01 -9.80E+00 1.44E-01 -4.00E-01 3.30E-02 2.39E-02 -5.46E-03 -1.55E-02
78 2.67E+00 -2.14E+01 -3.80E+00 2.61E+01 -1.62E+01 4.58E-02 2.75E+00 4.08E-03 -3.20E-03 5.05E-03 1.93E-04
79 -2.28E+00 -7.39E+01 -7.77E+00 5.84E+00 2.60E+01 -1.90E-01 -2.35E+00 -3.70E-02 -2.07E-02 4.10E-04 1.53E-02
80 6.17E+00 1.22E+02 3.14E+01 -3.95E+01 -1.02E+01 1.66E-01 -3.63E-01 3.14E-02 2.39E-02 -5.46E-03 -1.55E-02
81 8.26E-01 -3.51E+01 -8.32E+00 2.74E+01 -1.41E+01 2.81E-02 3.52E+00 2.96E-03 -3.20E-03 5.05E-03 1.93E-04
82 -6.99E+00 -8.69E+01 -2.31E+01 1.21E+01 2.43E+01 -1.94E-01 -3.16E+00 -3.44E-02 -2.07E-02 4.10E-04 1.53E-02
83 1.38E+00 5.17E+01 1.66E+01 -4.68E+01 -6.92E+00 1.61E-01 1.51E+00 1.41E-02 2.39E-02 -5.46E-03 -1.55E-02
84 2.84E+00 2.46E+01 -1.22E+00 2.92E+01 -2.37E+01 8.82E-02 -8.57E-02 1.11E-02 -3.20E-03 5.05E-03 1.93E-04
85 -4.22E+00 -7.63E+01 -1.54E+01 1.75E+01 3.06E+01 -2.49E-01 -1.42E+00 -2.52E-02 -2.07E-02 4.10E-04 1.53E-02
86 -2.11E+00 7.73E+01 3.65E+01 -4.22E+01 -1.23E+01 1.86E-01 4.90E+00 3.39E-02 2.39E-02 -5.46E-03 -1.55E-02
87 4.90E+00 -1.76E+01 -1.38E+01 2.62E+01 -1.21E+01 1.74E-02 1.95E-02 -4.71E-04 -3.20E-03 5.05E-03 1.93E-04
88 -2.79E+00 -5.96E+01 -2.27E+01 1.60E+01 2.43E+01 -2.04E-01 -4.92E+00 -3.35E-02 -2.07E-02 4.10E-04 1.53E-02
89 -1.29E-01 1.86E+02 1.68E+01 -4.97E+01 -6.01E+00 1.62E-01 6.26E+00 2.57E-02 2.39E-02 -5.46E-03 -1.55E-02
90 3.53E+00 -8.48E+01 -2.90E+00 2.92E+01 -1.20E+01 9.74E-03 1.78E-01 2.51E-04 -3.20E-03 5.05E-03 1.93E-04
91 -3.40E+00 -1.01E+02 -1.39E+01 2.05E+01 1.80E+01 -1.72E-01 -6.43E+00 -2.60E-02 -2.07E-02 4.10E-04 1.53E-02
92 1.16E+01 8.30E+00 1.44E+01 -2.01E+01 -1.34E+01 1.39E-01 2.68E+00 2.25E-02 2.39E-02 -5.46E-03 -1.55E-02
93 -5.66E+00 1.34E+01 -3.30E+00 5.53E+00 -1.15E+01 6.39E-02 1.80E+00 1.02E-02 -3.20E-03 5.05E-03 1.93E-04
94 -5.96E+00 -2.17E+01 -1.11E+01 1.46E+01 2.49E+01 -2.03E-01 -4.49E+00 -3.27E-02 -2.07E-02 4.10E-04 1.53E-02
95 -3.60E+00 1.70E+02 -6.30E-01 -4.28E+01 -1.10E+01 1.79E-01 3.10E-01 2.17E-02 2.39E-02 -5.46E-03 -1.55E-02
96 4.87E+00 -6.20E+01 6.01E+00 3.20E+01 -1.34E+01 1.20E-02 6.17E+00 8.62E-03 -3.20E-03 5.05E-03 1.93E-04
97 -1.27E+00 -1.08E+02 -5.38E+00 1.08E+01 2.44E+01 -1.91E-01 -6.48E+00 -3.03E-02 -2.07E-02 4.10E-04 1.53E-02
98 -4.98E+00 5.28E+01 7.16E+00 -2.90E+01 -2.71E+00 8.93E-02 1.06E+00 1.41E-02 1.36E-02 -5.25E-03 -7.84E-03
99 4.98E+00 -5.28E+01 -7.16E+00 2.90E+01 2.71E+00 -8.93E-02 -1.06E+00 -1.41E-02 -1.36E-02 5.25E-03 7.84E-03
100 1.01E+00 1.33E+02 1.84E+01 -3.17E+01 -1.18E+01 1.58E-01 2.15E+00 2.59E-02 2.39E-02 -5.46E-03 -1.55E-02
101 1.70E+00 -3.55E+01 -4.67E+00 1.14E+01 -1.06E+01 4.33E-02 1.43E+00 5.97E-03 -3.20E-03 5.05E-03 1.93E-04
102 -2.71E+00 -9.72E+01 -1.37E+01 2.03E+01 2.24E+01 -2.01E-01 -3.59E+00 -3.19E-02 -2.07E-02 4.10E-04 1.53E-02
103 -4.28E+00 1.56E+02 1.90E+01 -4.09E+01 4.17E+00 7.20E-02 9.36E-01 1.52E-02 1.36E-02 -5.25E-03 -7.84E-03
104 4.28E+00 -1.56E+02 -1.90E+01 4.09E+01 -4.17E+00 -7.20E-02 -9.36E-01 -1.52E-02 -1.36E-02 5.25E-03 7.84E-03
105 1.36E+00 5.09E+01 1.13E+01 -1.61E+01 -1.59E+01 1.47E-01 3.64E+00 3.52E-02 2.39E-02 -5.46E-03 -1.55E-02
106 1.10E+00 -1.13E+01 -7.79E+00 1.81E+01 -1.75E+01 7.37E-02 2.81E-01 1.56E-02 -3.20E-03 5.05E-03 1.93E-04
107 -2.47E+00 -3.96E+01 -3.48E+00 -2.06E+00 3.34E+01 -2.21E-01 -3.92E+00 -5.08E-02 -2.07E-02 4.10E-04 1.53E-02
108 -6.59E-01 2.06E+01 -2.83E+00 0 -2.30E+01 1.55E-01 8.61E+00 2.81E-02 8.76E-03 2.32E-03 -7.55E-03
109 6.59E-01 -2.06E+01 2.83E+00 0 2.30E+01 -1.55E-01 -8.61E+00 -2.81E-02 -8.76E-03 -2.32E-03 7.55E-03
110 2.20E+00 3.99E+01 9.17E+00 -6.29E+00 -3.58E+00 3.95E-02 -2.02E+00 9.09E-04 1.36E-02 -5.25E-03 -7.84E-03
111 -2.20E+00 -3.99E+01 -9.17E+00 6.29E+00 3.58E+00 -3.95E-02 2.02E+00 -9.09E-04 -1.36E-02 5.25E-03 7.84E-03
112 -2.08E+00 6.78E+01 1.79E+01 -3.94E+01 -9.57E+00 1.61E-01 5.27E-01 2.87E-02 2.39E-02 -5.46E-03 -1.55E-02
113 3.48E+00 -1.81E+01 -9.39E+00 1.27E+01 -1.50E+01 7.06E-02 3.97E-01 8.88E-03 -3.20E-03 5.05E-03 1.93E-04
114 -1.41E+00 -4.97E+01 -8.51E+00 2.67E+01 2.46E+01 -2.32E-01 -9.24E-01 -3.76E-02 -2.07E-02 4.10E-04 1.53E-02
115 3.14E-01 5.67E+01 6.55E+00 -1.66E+01 6.52E-01 3.62E-02 -6.31E-01 6.28E-03 1.36E-02 -5.25E-03 -7.84E-03
116 -3.14E-01 -5.67E+01 -6.55E+00 1.66E+01 -6.52E-01 -3.62E-02 6.31E-01 -6.28E-03 -1.36E-02 5.25E-03 7.84E-03
117 -4.02E+00 5.47E+01 8.97E+00 -1.65E+01 1.38E+00 3.10E-02 -1.36E+00 5.34E-03 1.36E-02 -5.25E-03 -7.84E-03
118 4.02E+00 -5.47E+01 -8.97E+00 1.65E+01 -1.38E+00 -3.10E-02 1.36E+00 -5.34E-03 -1.36E-02 5.25E-03 7.84E-03
119 3.53E+01 1.53E+02 3.25E+01 -1.95E+01 0 4.78E-02 4.37E+00 7.89E-02 2.39E-02 -5.46E-03 -1.55E-02
120 -1.50E+01 -6.41E+01 -1.51E+01 2.22E+01 0 -5.44E-02 5.98E+00 1.39E-02 -3.20E-03 5.05E-03 1.93E-04
121 -2.03E+01 -8.92E+01 -1.74E+01 -2.70E+00 0 6.59E-03 -1.03E+01 -9.28E-02 -2.07E-02 4.10E-04 1.53E-02
122 -1.48E+00 1.94E+01 6.76E+00 -2.43E+01 0 5.95E-02 -3.08E+00 7.90E-03 2.39E-02 -5.46E-03 -1.55E-02
123 -6.49E-01 -9.29E+00 -7.87E+00 -1.92E+00 0 4.69E-03 1.28E+00 -8.90E-03 -3.20E-03 5.05E-03 1.93E-04
124 2.13E+00 -1.02E+01 1.10E+00 2.62E+01 0 -6.42E-02 1.80E+00 9.94E-04 -2.07E-02 4.10E-04 1.53E-02
125 -1.77E+00 1.14E+02 2.14E+01 -4.61E+01 -9.65E+00 1.78E-01 2.48E+00 3.33E-02 2.39E-02 -5.46E-03 -1.55E-02
126 3.55E-01 -2.79E+01 -1.50E+01 1.66E+01 -1.18E+01 3.89E-02 2.31E+00 4.41E-03 -3.20E-03 5.05E-03 1.93E-04
127 1.41E+00 -8.65E+01 -6.40E+00 2.95E+01 2.14E+01 -2.17E-01 -4.79E+00 -3.78E-02 -2.07E-02 4.10E-04 1.53E-02
128 5.64E+00 1.44E+02 5.24E+01 -5.41E+01 0 1.32E-01 1.81E+00 1.28E-02 2.39E-02 -5.46E-03 -1.55E-02
129 -6.70E+00 -3.60E+01 -2.68E+01 -1.87E+00 0 4.57E-03 5.45E-01 1.24E-03 -3.20E-03 5.05E-03 1.93E-04
130 1.06E+00 -1.08E+02 -2.56E+01 5.60E+01 0 -1.37E-01 -2.35E+00 -1.40E-02 -2.07E-02 4.10E-04 1.53E-02
131 -4.37E-01 2.99E+00 -1.31E+00 5.25E+00 -1.65E+00 -1.74E-03 -3.03E+00 -4.33E-03 9.63E-03 -2.95E-03 -5.91E-03
132 4.37E-01 -2.99E+00 1.31E+00 -5.25E+00 1.65E+00 1.74E-03 3.03E+00 4.33E-03 -9.63E-03 2.95E-03 5.91E-03
133 2.72E+01 9.92E+01 3.43E+01 -4.65E+01 -2.74E+01 2.99E-01 4.60E+00 6.78E-02 1.32E-02 -1.02E-03 -9.41E-03
134 1.03E+01 9.19E+01 2.96E+00 2.71E+01 -3.28E+01 1.55E-01 1.13E+01 6.83E-02 -1.39E-02 9.48E-03 6.26E-03
135 -9.19E+00 -1.44E+01 -1.11E+01 -6.56E+00 1.88E+01 -1.11E-01 -2.64E+00 -1.53E-02 -3.14E-02 4.84E-03 2.14E-02
136 3.41E+00 -2.82E+01 -1.42E+01 -1.06E+00 -9.64E-01 9.11E-03 4.97E+00 -1.17E-03 8.46E-03 1.78E-03 -7.10E-03
137 2.23E+00 -2.13E+01 -2.28E+00 -6.94E+00 3.67E+00 -7.82E-03 -1.19E+00 1.35E-03 1.77E-02 -8.57E-03 -9.47E-03
138 1.66E+00 -2.26E+01 2.02E+00 8.26E+01 -3.22E+01 1.52E-02 1.73E+00 -2.73E-02 6.67E-03 2.83E-03 -6.21E-03
139 3.89E+00 -1.17E+00 1.22E+01 7.13E+01 -2.40E+01 -1.22E-02 -1.70E-01 -1.05E-02 -1.26E-02 8.73E-03 5.61E-03
140 -1.56E+01 -4.24E+01 4.43E+00 -2.34E+01 1.96E+01 -7.48E-02 -2.03E+00 -2.93E-02 7.62E-03 -5.15E-03 -3.45E-03
141 -2.40E+01 -6.09E+01 -2.83E+01 -9.66E+01 7.54E+01 -2.73E-01 -1.66E+01 -5.39E-02 4.31E-03 -1.29E-02 2.41E-03
142 4.36E+01 4.18E+01 4.29E+01 -3.61E+01 1.09E+00 8.10E-02 9.86E+00 2.84E-02 1.36E-02 -5.25E-03 -7.84E-03
143 -4.36E+01 -4.18E+01 -4.29E+01 3.61E+01 -1.09E+00 -8.10E-02 -9.86E+00 -2.84E-02 -1.36E-02 5.25E-03 7.84E-03
144 5.69E-01 6.95E+01 1.29E+00 4.33E+00 -3.24E+01 2.08E-01 1.58E+00 5.23E-02 5.17E-03 4.46E-03 -5.80E-03
145 9.74E-01 6.61E-01 2.53E+00 5.25E+00 -5.04E+01 3.27E-01 4.53E-01 4.08E-02 -1.41E-02 1.04E-02 6.02E-03
146 -5.76E-01 -1.01E+01 -9.99E-01 -3.04E+00 1.84E+01 -1.17E-01 6.10E-01 -1.81E-02 6.12E-03 -3.52E-03 -3.04E-03
147 -9.68E-01 -6.01E+01 -2.82E+00 -6.55E+00 6.43E+01 -4.18E-01 -2.64E+00 -7.51E-02 2.81E-03 -1.13E-02 2.82E-03
148 -4.44E-01 5.64E+01 3.03E+00 -2.40E+01 -2.29E+01 2.14E-01 4.41E+00 4.42E-02 2.39E-02 -5.46E-03 -1.55E-02
149 1.17E+00 2.39E+01 -4.30E+00 -5.12E+00 -4.87E+00 4.54E-02 2.07E+00 7.92E-03 -3.20E-03 5.05E-03 1.93E-04
150 -7.30E-01 -8.03E+01 1.27E+00 2.92E+01 2.78E+01 -2.59E-01 -6.48E+00 -5.21E-02 -2.07E-02 4.10E-04 1.53E-02
151 1.80E+01 1.17E+02 2.79E+01 -2.43E+01 -2.83E+01 2.50E-01 7.15E+00 5.49E-02 1.32E-02 -1.02E-03 -9.41E-03
152 -4.52E+00 9.49E+01 2.70E+01 8.01E+00 -2.29E+01 1.35E-01 8.70E+00 4.23E-02 -1.39E-02 9.48E-03 6.26E-03
153 -3.50E+00 -1.83E+01 -4.79E+00 4.98E+01 1.29E+00 -1.31E-01 -4.37E+00 -3.34E-02 -3.14E-02 4.84E-03 2.14E-02
154 -1.17E+00 -8.23E+00 2.23E+00 1.26E+01 -1.15E+01 4.68E-02 -3.39E-01 1.11E-02 8.46E-03 1.78E-03 -7.10E-03
155 1.77E+00 -1.40E+01 -9.32E+00 -4.40E+01 7.82E+00 5.49E-02 -1.86E-01 6.24E-03 1.77E-02 -8.57E-03 -9.47E-03
156 5.61E+00 1.42E+01 2.82E-01 4.86E+01 -3.27E+01 1.02E-01 7.29E+00 1.80E-02 6.67E-03 2.83E-03 -6.21E-03
157 1.05E+00 -1.90E+01 -1.04E+00 2.17E+01 -1.65E+01 5.86E-02 -4.86E-01 7.53E-03 -1.26E-02 8.73E-03 5.61E-03
158 -4.81E+00 -6.03E+01 -1.61E+01 3.23E+00 1.85E+01 -1.33E-01 -3.13E+00 -2.68E-02 7.62E-03 -5.15E-03 -3.45E-03
159 -1.25E+01 -1.06E+02 -2.61E+01 -7.56E+01 8.43E+01 -3.84E-01 -1.46E+01 -7.99E-02 4.31E-03 -1.29E-02 2.41E-03
160 1.86E+01 1.56E+02 8.02E+01 -5.48E+01 -8.47E-01 1.40E-01 1.30E+01 4.57E-02 1.32E-02 -1.02E-03 -9.41E-03
161 1.15E+01 2.60E+01 1.96E+01 3.40E+01 -2.21E+01 6.60E-02 3.95E+00 1.89E-02 -1.39E-02 9.48E-03 6.26E-03
162 -9.09E+00 -2.78E+01 -1.55E+01 2.97E+01 -2.83E+00 -5.35E-02 -3.22E+00 -2.67E-02 -3.14E-02 4.84E-03 2.14E-02
163 -3.77E+00 -2.06E+01 -9.98E+00 -2.80E+00 7.18E+00 -4.16E-02 -3.57E+00 -1.64E-02 8.46E-03 1.78E-03 -7.10E-03
164 -7.21E+00 -3.88E+01 -1.79E+01 -7.13E+01 1.70E+01 5.93E-02 -1.30E+00 1.23E-02 1.77E-02 -8.57E-03 -9.47E-03
165 5.75E+00 3.90E+01 -1.01E+01 3.52E+01 -3.34E+01 1.39E-01 2.31E+00 3.70E-02 6.67E-03 2.83E-03 -6.21E-03
166 9.44E+00 -1.53E+01 2.14E+00 5.60E+01 -2.35E+01 2.16E-02 -1.13E+00 9.19E-03 -1.26E-02 8.73E-03 5.61E-03
167 -7.03E+00 -4.01E+01 -1.18E+01 1.88E+01 2.23E+00 -6.10E-02 -1.18E+00 -1.56E-02 7.62E-03 -5.15E-03 -3.45E-03
168 -1.82E+01 -7.88E+01 -3.67E+01 -4.47E+01 5.62E+01 -2.70E-01 -8.86E+00 -6.44E-02 4.31E-03 -1.29E-02 2.41E-03
169 3.35E+01 1.52E+02 4.28E+01 -2.80E+01 -3.47E+01 3.03E-01 1.13E+01 1.00E-01 1.32E-02 -1.02E-03 -9.41E-03
170 2.05E+01 7.06E+01 8.32E+00 5.60E+01 -1.17E+01 -5.83E-02 8.43E+00 4.67E-03 -1.39E-02 9.48E-03 6.26E-03
171 -9.32E+00 -5.02E+01 -5.80E+00 1.37E+02 -2.20E+01 -1.86E-01 -7.25E+00 -3.93E-02 -3.14E-02 4.84E-03 2.14E-02
172 -5.56E+00 -4.05E+01 8.93E-01 2.47E+01 -3.43E+00 -3.74E-02 -3.34E+00 -1.95E-02 8.46E-03 1.78E-03 -7.10E-03
173 -1.07E+00 -3.38E+01 -9.79E+00 -6.88E+01 4.04E+01 -1.05E-01 -3.53E+00 -6.77E-03 1.77E-02 -8.57E-03 -9.47E-03
174 6.93E-01 2.26E+01 1.03E+01 3.05E+01 -1.37E+01 1.79E-02 1.45E+01 1.24E-02 6.67E-03 2.83E-03 -6.21E-03
175 -4.21E+00 1.87E+01 5.82E+00 -3.14E+01 -1.38E+01 1.70E-01 1.54E+00 1.22E-02 -1.26E-02 8.73E-03 5.61E-03
176 -1.56E+01 -5.85E+01 -2.34E+01 -5.31E+01 7.43E+00 7.98E-02 -6.99E+00 -2.03E-02 7.62E-03 -5.15E-03 -3.45E-03
177 -1.89E+01 -8.04E+01 -2.91E+01 -6.68E+01 5.14E+01 -1.84E-01 -1.47E+01 -4.35E-02 4.31E-03 -1.29E-02 2.41E-03
178 1.20E+00 1.86E+02 4.85E+01 -4.13E+01 -1.46E+01 1.99E-01 1.31E+01 6.00E-02 1.32E-02 -1.02E-03 -9.41E-03
179 9.93E+00 9.57E+01 2.66E+01 -1.80E+01 -4.78E-01 4.71E-02 1.50E+01 6.45E-02 -1.39E-02 9.48E-03 6.26E-03
180 -2.22E+00 -4.53E+01 -3.70E+00 5.73E+01 -1.33E+01 -5.03E-02 -2.72E+00 -2.40E-02 -3.14E-02 4.84E-03 2.14E-02
181 -4.66E+00 -3.70E+01 1.90E+01 5.67E+01 -1.62E+01 -2.93E-02 -3.95E+00 -1.69E-02 8.46E-03 1.78E-03 -7.10E-03
182 7.15E+00 -3.20E+01 -7.74E+00 -5.72E+01 8.47E+00 8.28E-02 7.42E-01 1.23E-02 1.77E-02 -8.57E-03 -9.47E-03
183 6.70E+00 -1.27E+00 -6.76E+00 5.82E+01 -2.88E+01 5.24E-02 -3.09E+00 1.04E-02 6.67E-03 2.83E-03 -6.21E-03
184 8.73E+00 -7.97E+00 -1.81E+01 8.48E+01 -3.25E+01 1.19E-02 2.16E+00 1.40E-02 -1.26E-02 8.73E-03 5.61E-03
185 -1.34E+01 -7.42E+01 -2.66E+01 -8.03E+01 2.93E+01 -9.96E-04 -8.38E+00 -4.80E-02 7.62E-03 -5.15E-03 -3.45E-03
186 -1.34E+01 -8.34E+01 -3.11E+01 -6.03E+01 6.82E+01 -3.13E-01 -1.29E+01 -7.23E-02 4.31E-03 -1.29E-02 2.41E-03
187 5.11E+01 1.51E+02 4.95E+01 2.02E+01 -3.17E+01 1.64E-01 1.50E+01 6.80E-02 1.32E-02 -1.02E-03 -9.41E-03
188 2.09E+01 1.21E+00 2.60E+01 8.93E+01 -2.25E+01 -6.63E-02 1.03E+00 -2.65E-03 -1.39E-02 9.48E-03 6.26E-03
189 -6.97E+00 -1.81E+01 -1.95E+01 9.08E+01 -2.03E+01 -8.52E-02 -3.38E+00 -1.37E-02 -3.14E-02 4.84E-03 2.14E-02
190 -4.54E+00 -2.73E+01 -1.25E+01 4.47E+00 -3.10E+00 1.00E-02 2.71E+00 6.27E-03 8.46E-03 1.78E-03 -7.10E-03
191 -9.39E+00 -4.58E+00 -1.67E+01 -7.32E+01 1.39E+01 8.54E-02 -1.32E+00 -1.23E-02 1.77E-02 -8.57E-03 -9.47E-03
192 -4.78E-01 9.68E-01 -8.56E+00 5.31E+01 -2.95E+01 6.94E-02 2.26E-01 -2.74E-03 6.67E-03 2.83E-03 -6.21E-03
193 -8.45E+00 1.94E+01 2.19E+01 -5.10E+01 -8.28E+00 1.81E-01 -9.84E-01 2.52E-02 -1.26E-02 8.73E-03 5.61E-03
194 -1.76E+01 -4.49E+01 -1.08E+01 -4.34E+01 3.40E+01 -1.23E-01 -2.93E+00 -1.37E-02 7.62E-03 -5.15E-03 -3.45E-03
195 -2.47E+01 -7.78E+01 -2.94E+01 -9.02E+01 6.76E+01 -2.35E-01 -1.04E+01 -5.43E-02 4.31E-03 -1.29E-02 2.41E-03
196 -3.61E+00 1.90E+02 6.15E+01 -5.82E+01 -3.19E+01 3.58E-01 1.92E+01 1.07E-01 1.32E-02 -1.02E-03 -9.41E-03
197 2.64E+01 7.46E+01 1.30E+01 8.56E+00 -2.77E+01 1.66E-01 8.58E+00 4.97E-02 -1.39E-02 9.48E-03 6.26E-03
198 -4.84E+00 -4.92E+01 -1.80E+01 9.66E+01 -1.51E+01 -1.34E-01 -3.91E+00 -3.15E-02 -3.14E-02 4.84E-03 2.14E-02
199 1.63E+00 -4.49E+01 -1.11E+01 3.21E+01 -9.32E-02 -7.80E-02 -1.08E+00 -2.34E-02 8.46E-03 1.78E-03 -7.10E-03
200 -8.63E+00 -4.42E+01 -6.73E+00 -7.17E+01 3.23E+01 -4.25E-02 -3.21E+00 -2.04E-02 1.77E-02 -8.57E-03 -9.47E-03
201 7.43E+00 2.04E+01 9.91E+00 6.59E+01 -3.35E+01 6.52E-02 1.12E+00 1.78E-02 6.67E-03 2.83E-03 -6.21E-03
202 1.30E+00 4.82E-01 1.63E+00 1.74E+01 -1.53E+01 6.06E-02 -3.47E+00 -4.44E-03 -1.26E-02 8.73E-03 5.61E-03
203 -8.64E+00 -6.44E+01 -1.71E+01 -2.46E+01 3.39E+01 -1.69E-01 -5.42E+00 -3.36E-02 7.62E-03 -5.15E-03 -3.45E-03
204 -1.10E+01 -8.32E+01 -3.31E+01 -6.60E+01 5.74E+01 -2.26E-01 -1.18E+01 -6.08E-02 4.31E-03 -1.29E-02 2.41E-03
205 3.08E+00 9.19E+01 1.94E+01 -8.50E+01 -6.04E+00 2.49E-01 7.14E+00 6.37E-02 2.23E-02 -2.93E-03 -1.54E-02
206 -3.08E+00 -9.19E+01 -1.94E+01 8.50E+01 6.04E+00 -2.49E-01 -7.14E+00 -6.37E-02 -2.23E-02 2.93E-03 1.54E-02
207 1.78E+01 1.29E+02 1.11E+01 -5.02E+01 -1.73E+01 2.39E-01 9.79E+00 7.07E-02 1.32E-02 -1.02E-03 -9.41E-03
208 7.90E+00 9.33E+01 1.56E+00 -5.41E+01 -1.97E+00 1.46E-01 1.33E+01 6.89E-02 -1.39E-02 9.48E-03 6.26E-03
209 -1.40E+01 -1.57E+01 1.45E+00 2.94E+01 -1.08E+01 6.67E-04 -8.03E-01 -4.16E-03 -3.14E-02 4.84E-03 2.14E-02
210 1.57E+01 -5.33E+01 1.48E+01 1.46E+02 -3.71E+01 -1.07E-01 -5.12E+00 -3.88E-02 8.46E-03 1.78E-03 -7.10E-03
211 7.16E+00 -9.99E+00 -2.36E+00 -2.67E+01 4.85E+00 3.27E-02 3.21E+00 -1.12E-02 1.77E-02 -8.57E-03 -9.47E-03
212 1.00E+01 -2.37E+01 1.01E+01 1.07E+02 -3.02E+01 -5.89E-02 1.42E+00 -1.92E-02 6.67E-03 2.83E-03 -6.21E-03
213 -2.49E+00 -7.11E+00 3.00E+00 4.60E+01 -3.53E+01 1.26E-01 -1.30E+00 1.36E-02 -1.26E-02 8.73E-03 5.61E-03
214 -1.73E+01 -3.63E+01 -8.38E+00 -5.90E+01 2.78E+01 -4.34E-02 -1.71E+00 6.41E-03 7.62E-03 -5.15E-03 -3.45E-03
215 -2.48E+01 -7.62E+01 -3.13E+01 -1.39E+02 9.99E+01 -3.34E-01 -1.88E+01 -8.62E-02 4.31E-03 -1.29E-02 2.41E-03
216 1.01E+01 7.73E+01 4.65E+01 -3.98E+00 -1.74E+01 1.27E-01 1.32E+01 5.47E-02 1.32E-02 -1.02E-03 -9.41E-03
217 2.18E+01 1.39E+01 9.70E+00 1.19E+02 -2.87E+01 -9.72E-02 8.91E+00 3.45E-03 -1.39E-02 9.48E-03 6.26E-03
218 1.70E+00 -1.21E+01 5.46E+00 1.02E+02 -2.73E+01 -6.46E-02 -3.28E+00 -3.28E-02 -3.14E-02 4.84E-03 2.14E-02
219 1.45E+01 -2.56E+01 -7.28E+00 -1.25E+01 -6.73E-01 3.50E-02 6.34E-01 1.59E-02 8.46E-03 1.78E-03 -7.10E-03
220 -7.54E+00 1.30E+01 -7.69E+00 -9.36E+01 1.43E+01 1.32E-01 -5.42E+00 2.07E-03 1.77E-02 -8.57E-03 -9.47E-03
221 -5.19E+00 1.63E+01 -5.94E+00 -7.33E+00 -1.71E+01 1.34E-01 1.03E+01 1.72E-02 6.67E-03 2.83E-03 -6.21E-03
222 -4.93E+00 1.23E+01 -9.39E+00 1.67E+01 -1.45E+01 5.67E-02 -7.07E+00 8.74E-04 -1.26E-02 8.73E-03 5.61E-03
223 -1.25E+01 -3.95E+01 -1.08E+01 -5.50E+01 2.90E+01 -6.13E-02 -5.98E+00 -1.42E-02 7.62E-03 -5.15E-03 -3.45E-03
224 -1.79E+01 -5.57E+01 -2.06E+01 -6.49E+01 6.23E+01 -2.62E-01 -1.12E+01 -4.71E-02 4.31E-03 -1.29E-02 2.41E-03
225 4.79E+01 4.35E+01 1.78E+02 8.36E+01 -3.25E+01 1.48E-02 3.80E+01 2.54E-02 1.32E-02 -1.02E-03 -9.41E-03
226 2.70E+01 1.33E+01 5.03E+01 1.29E+02 -3.25E+01 -9.55E-02 1.22E+01 -1.19E-02 -1.39E-02 9.48E-03 6.26E-03
227 -7.83E+00 -2.27E+01 -1.74E+01 5.28E+01 8.74E+00 -1.88E-01 -4.93E+00 -1.72E-02 -3.14E-02 4.84E-03 2.14E-02
228 -6.85E+00 1.21E+01 -3.28E+01 3.79E+01 -4.74E+00 -6.06E-02 -5.91E+00 -1.00E-02 8.46E-03 1.78E-03 -7.10E-03
229 -1.24E+01 7.12E+00 -2.47E+01 -7.94E+01 5.75E+00 1.55E-01 3.20E+00 6.61E-03 1.77E-02 -8.57E-03 -9.47E-03
230 2.51E+00 1.25E+01 -2.76E+01 -1.70E+01 -2.37E+01 2.02E-01 -5.98E+00 2.09E-02 6.67E-03 2.83E-03 -6.21E-03
231 -1.18E+01 -9.36E+00 -3.08E+01 -6.52E+00 -2.06E+01 1.55E-01 -7.66E+00 2.05E-02 -1.26E-02 8.73E-03 5.61E-03
232 -1.77E+01 -1.98E+01 -4.57E+01 -8.08E+01 2.86E+01 4.81E-03 -1.18E+01 -2.41E-03 7.62E-03 -5.15E-03 -3.45E-03
233 -2.08E+01 -3.67E+01 -4.94E+01 -1.19E+02 7.10E+01 -1.87E-01 -1.72E+01 -3.19E-02 4.31E-03 -1.29E-02 2.41E-03
234 1.19E+01 1.75E+01 -3.15E+00 6.34E-01 -2.00E+00 1.19E-02 2.57E+01 1.23E-01 1.37E-02 -2.64E-03 -9.11E-03
235 8.96E+00 3.83E+00 -4.91E+00 1.43E+01 -8.33E+00 2.12E-02 1.45E+01 5.77E-02 -1.34E-02 7.87E-03 6.56E-03
236 -2.00E+00 -2.89E+01 1.51E+01 6.09E+01 1.22E+01 -2.31E-01 6.70E-02 1.53E-02 -3.09E-02 3.23E-03 2.17E-02
237 -4.04E+00 4.34E+01 -6.34E-01 -6.29E+00 5.54E+00 -2.20E-02 -6.40E+00 -3.29E-03 9.00E-03 1.66E-04 -6.79E-03
238 -6.39E+00 -2.23E+01 -4.31E+00 -4.43E+01 1.11E+01 3.33E-02 -1.23E+01 -5.66E-02 1.82E-02 -1.02E-02 -9.17E-03
239 2.65E+00 -3.44E+01 1.82E+01 6.13E+01 -4.07E+01 1.25E-01 -1.11E+01 -5.00E-02 7.21E-03 1.21E-03 -5.91E-03
240 -5.16E+00 4.56E+01 -9.09E+00 -3.19E+01 -1.22E+01 1.60E-01 -6.93E+00 -2.73E-02 -1.21E-02 7.12E-03 5.91E-03
241 -5.90E+00 -2.46E+01 -1.12E+01 -5.47E+01 3.44E+01 -9.87E-02 -3.55E+00 -5.89E-02 8.16E-03 -6.76E-03 -3.15E-03
242 1.85E+01 -7.90E-01 3.40E+01 -1.94E+01 -4.56E+00 7.82E-02 1.10E+01 7.45E-02 2.39E-02 -5.46E-03 -1.55E-02
243 -3.53E+00 -5.33E+00 -1.29E+01 5.47E+01 -1.65E+01 -2.24E-02 -2.82E+00 -2.19E-02 -3.20E-03 5.05E-03 1.93E-04
244 -1.50E+01 6.13E+00 -2.12E+01 -3.53E+01 2.11E+01 -5.58E-02 -8.22E+00 -5.26E-02 -2.07E-02 4.10E-04 1.53E-02
245 3.73E+01 1.66E+01 1.24E+02 -3.39E+01 7.30E+00 3.37E-02 2.39E+01 3.99E-02 1.32E-02 -1.02E-03 -9.41E-03
246 1.22E+01 1.72E+01 4.43E+01 3.17E+01 -9.69E+00 -1.22E-02 8.31E+00 2.10E-02 -1.39E-02 9.48E-03 6.26E-03
247 3.87E-01 1.34E+00 -1.64E+01 1.14E+02 -3.17E+01 -6.46E-02 3.21E+00 -1.19E-02 -3.14E-02 4.84E-03 2.14E-02
248 6.52E+00 2.63E+01 1.39E+01 9.04E+01 -3.17E+01 -7.30E-03 5.68E+00 -6.14E-03 8.46E-03 1.78E-03 -7.10E-03
249 -9.11E+00 2.12E+01 -1.90E+01 -9.62E+01 2.33E+01 7.80E-02 3.39E+00 -4.64E-03 1.77E-02 -8.57E-03 -9.47E-03
250 -6.21E+00 3.34E+00 -1.93E+01 -1.32E+01 -2.03E+01 1.69E-01 -5.92E+00 1.02E-02 6.67E-03 2.83E-03 -6.21E-03
251 -7.83E+00 -1.44E+01 -3.50E+01 3.90E+01 -2.30E+01 5.97E-02 -9.46E+00 -2.94E-04 -1.26E-02 8.73E-03 5.61E-03
252 -1.31E+01 -1.68E+01 -4.12E+01 -4.08E+01 1.57E+01 -5.90E-03 -1.07E+01 -3.55E-03 7.62E-03 -5.15E-03 -3.45E-03
253 -2.02E+01 -5.48E+01 -5.18E+01 -9.07E+01 7.00E+01 -2.51E-01 -1.85E+01 -4.45E-02 4.31E-03 -1.29E-02 2.41E-03
254 1.83E-01 1.62E+01 1.40E+00 0 1.92E+01 -1.30E-01 -1.22E+00 1.12E-02 -2.97E-02 5.49E-03 1.98E-02
255 2.17E+00 2.92E+01 9.57E+00 0 -9.94E+00 6.71E-02 3.77E+00 1.56E-02 1.02E-02 2.43E-03 -8.70E-03
256 -2.35E+00 -4.55E+01 -1.10E+01 0 -9.28E+00 6.26E-02 -2.55E+00 -2.68E-02 1.94E-02 -7.92E-03 -1.11E-02
257 5.26E+00 3.82E+01 2.74E+01 -6.31E+01 0 1.54E-01 2.19E+00 1.61E-02 1.36E-02 -5.25E-03 -7.84E-03
RUN;
ods trace off;
ODS output ResidSummary=resi
(keep= SS);
/* ------- Define model --------------------------------------------------- */
title 'SUR regression with normal variables ';
PROC MODEL data=temp;
yfC= C1*YFOT +C2*YFMI +C3*ZFDC +C4*ZFLA +C5*PNC;
yfA= C1*YFOT +C2*YFMI +C3*ZFDC +C4*ZFLA +C6*PNA;
yfT= C1*YFOT +C2*YFMI +C3*ZFDC +C4*ZFLA +C7*PNT;
Parms C C1-C7;
FIT YFC YFA YFT /ITSUR METHOD=Marquardt CONVERGE=1e-6 MAXITER=1000 BREUSCH=(ZFLA)
OUTEST=EXP OUTRESID OUT=RESID1 ;
RUN;
/*
title 'print resid1';
proc print data=resid1;
run;
*/
/* Create new left hand side variable for each of the equations above
sigmaC_i^2 = errorC_i^2/(sseC/N)
sigmaA_i^2 = errorA_i^2/(sseA/N)
sigmaT_i^2 = errorT_i^2/(sseT/N) */
/* here I would like to use the SSE from the residual summary in order to calculate new variable
from the resi ods data set */
*%let ssC=ssC;
*%let ssA=ssA;
*%let ssT=ssT;
/* thats the way I do it now....I copy SSE and number of obsverations manually from the listing ... */
data sser;
set resid1;
siC = (yfc*yfc)/(16684.2/257);
siA = (yfa*yfa)/(1005031/257);
siT = (yft*yft)/(83052.6/257);
yfce = yfc;
run;
* merge the two set in order to get FYI from original data in set, new name: sser1;
data sser1;
set sser;
set temp;
merge;
run;
/*
title 'print merged data set sser1';
proc print data=sser1;
run;
*/
/* Run a regression of ZFLA on the sigma variables */
title 'Regression on error terms';
PROC MODEL data=sser1;
PARMS A A1-A4;
siC= A1 + A2*ZFLA + A3*ZFLA*ZFLA;
siA= A1 + A2*ZFLA + A3*ZFLA*ZFLA;
siT= A1 + A2*ZFLA + A3*ZFLA*ZFLA;
FIT siC siA siT /out=pre outpredict;
RUN;
/*
title 'Print predicted values of lhs';
proc print data=pre;
run;
*/
/* Now I have to transfor each equation with the weights from above */
/* Calculate the weights and get all data on one data set */
data tempW;
set pre;
set temp;
wC=1/sqrt(siC);
wA=1/sqrt(siA);
wT=1/sqrt(siT);
run;
/*
title 'Print weights for each equation';
proc print data = tempW;
run;
*/
/* Transform orignal model variables for each of the equations */
title 'SUR regression with transformed variables ';
PROC Model data=tempW;
endogenous yfc yfa yft;
yfC= C1*YFOT +C2*YFMI +C3*ZFDC +C4*ZFLA +C5*PNC;
yfA= C1*YFOT +C2*YFMI +C3*ZFDC +C4*ZFLA +C6*PNA;
yfT= C1*YFOT +C2*YFMI +C3*ZFDC +C4*ZFLA +C7*PNT;
resid.yfC=resid.yfc*wC;
resid.yfA=resid.yfa*wA;
resid.yfT=resid.yfT*wT;
Parms C C1-C7;
FIT yfc yfa yft /ITSUR METHOD=Marquardt CONVERGE=1e-6 MAXITER=1000 BREUSCH=(ZFLA);
RUN;
ods output close;
quit;
Christine Wieck
Institute for Agricultural Policy, University of Bonn
Back to: Top of message | Previous page | Main SAS-L page
listserv.uga.edu
Enterprise Information Technology Services
The University of Georgia
listhelp@uga.edu