home *** CD-ROM | disk | FTP | other *** search
/ DOS Wares / doswares.zip / doswares / PROGS / SPSS4 / SPSSNTA9.ZIP / ADVTEST.LIS < prev    next >
Encoding:
File List  |  1992-03-12  |  34.1 KB  |  779 lines
  1.            SPSS/PC+ The Statistical Package for IBM PC                 10/2/90 
  2. INC '\spss\advtest.inc'                                                        
  3. title 'Test Job for SPSS/PC+ ADVANCED STATS'.
  4. SET LOG OFF VIEW 24.
  5. -------------------------------------------------------------------------------
  6. Page   2   Test Job for SPSS/PC+ ADVANCED STATS                        10/2/90 
  7. DATA LIST  /CASEID 1-2 SURVIVAL 4 TREATMNT 6 TIME 8-10(1)
  8.   WEIGHT 12-15(3) APGAR 17-18 SEX 20 AGE 22-23
  9.   PO2 25-27 PCO2 29-31 PH 33-35(2) RESP 37.
  10. VAR LABELS  SURVIVAL 'Infant survival'
  11.   TREATMNT 'Treatment administered'
  12.   TIME 'Time to spontaneous respiration'
  13.   WEIGHT 'Birthweight in kilograms'
  14.   APGAR 'APGAR score'
  15.   SEX 'Sex of respondent'
  16.   AGE 'Gestation age'
  17.   PO2 'O2 level'
  18.   PCO2 'CO2 level'
  19.   PH 'PH level'
  20.   RESP 'Respiratory level'.
  21. RECODE SURVIVAL (0=1)(1=2).
  22. VALUE LABELS SURVIVAL 1 'Die' 2 'Survive'/
  23.   TREATMNT 1 'THAM' 0 'Sodium bicarbonate'/
  24.   SEX 1 'Male' 0 'Female'/
  25.   RESP 1 'Yes' 0 'No' 9 'No answer'/.
  26. MISSING VALUES  PCO2(-1)/RESP(9).
  27. BEGIN DATA
  28. END DATA.
  29.      50 cases are written to the compressed active file.
  30.  
  31. This procedure was completed at 14:29:09
  32. -------------------------------------------------------------------------------
  33. Page   3   Test Job for SPSS/PC+ ADVANCED STATS                        10/2/90 
  34.  
  35. * next procedure requires that M module be installed.
  36.  
  37. DSCRIMINANT GROUPS=SURVIVAL(1,2)
  38.   /VARIABLES=TREATMNT TO AGE PH RESP
  39.   /METHOD=DIRECT
  40.  
  41. Since ANALYSIS= was omitted for the first analysis all variables
  42. on the VARIABLES= list will be entered at level 1.
  43.   /STATISTICS=13 15
  44.   /OPTIONS= 6 7.
  45.  
  46.  
  47. This Discriminant Analysis requires      2240 (       2.2K) BYTES of workspace.
  48. -------------------------------------------------------------------------------
  49. Page   4   Test Job for SPSS/PC+ ADVANCED STATS                        10/2/90 
  50.  
  51. - - - - - - - -   D I S C R I M I N A N T   A N A L Y S I S   - - - - - - - - 
  52.  
  53. On groups defined by SURVIVAL  Infant survival                         
  54.  
  55.             50 (unweighted) cases were processed.
  56.              2 of these were excluded from the analysis.
  57.                 0 had missing or out-of-range group codes.
  58.                 2 had at least one missing discriminating variable.
  59.             48 (unweighted) cases will be used in the analysis.
  60.  
  61. Number of Cases by Group
  62.  
  63.                Number of Cases
  64.   SURVIVAL  Unweighted     Weighted  Label
  65.          1          26         26.0  Die                 
  66.          2          22         22.0  Survive             
  67.      Total          48         48.0
  68. -------------------------------------------------------------------------------
  69. Page   5   Test Job for SPSS/PC+ ADVANCED STATS                        10/2/90 
  70.  
  71. - - - - - - - -   D I S C R I M I N A N T   A N A L Y S I S   - - - - - - - - 
  72.  
  73. On groups defined by SURVIVAL  Infant survival                         
  74.  
  75. Analysis number       1
  76. Direct method:  All variables passing the tolerance test are entered.
  77.      Minimum Tolerance Level..................  .00100
  78.  
  79. Canonical Discriminant Functions
  80.  
  81.      Maximum number of functions..............       1
  82.      Minimum cumulative percent of variance...  100.00
  83.      Maximum significance of Wilks' Lambda....  1.0000
  84.  
  85. Prior probability for each group is  .50000
  86. -------------------------------------------------------------------------------
  87. Page   6   Test Job for SPSS/PC+ ADVANCED STATS                        10/2/90 
  88.  
  89.                       Canonical Discriminant Functions
  90.  
  91.                 Pct of   Cum  Canonical  After  Wilks'
  92.  Fcn Eigenvalue Variance  Pct     Corr      Fcn  Lambda  Chisquare    DF  Sig
  93.                                          :    0   .6150     20.419     8  .0089
  94.    1*     .6261  100.00  100.00    .6205 : 
  95.  
  96.    * marks the   1 canonical discriminant functions remaining in the analysis.
  97.  
  98. Standardized Canonical Discriminant Function Coefficients
  99.  
  100.              FUNC  1
  101. TREATMNT      .21531 
  102. TIME          .13170 
  103. WEIGHT       1.16789 
  104. APGAR         .34638 
  105. SEX           .00337 
  106. AGE          -.64084 
  107. PH            .16862 
  108. RESP         -.58743 
  109. -------------------------------------------------------------------------------
  110. Page   7   Test Job for SPSS/PC+ ADVANCED STATS                        10/2/90 
  111.  
  112. Structure Matrix:
  113.  
  114. Pooled-within-groups correlations between discriminating variables
  115.                                   and canonical discriminant functions
  116. (Variables ordered by size of correlation within function)
  117.  
  118.              FUNC  1
  119. WEIGHT        .73338 
  120. RESP         -.51940 
  121. AGE           .49290 
  122. TREATMNT      .26572 
  123. PH            .25946 
  124. APGAR         .19210 
  125. TIME         -.10157 
  126. SEX           .03738 
  127. -------------------------------------------------------------------------------
  128. Page   8   Test Job for SPSS/PC+ ADVANCED STATS                        10/2/90 
  129.  
  130. Canonical Discriminant Functions evaluated at Group Means (Group Centroids)
  131.  
  132.   Group      FUNC   1
  133.        1       -.71252
  134.        2        .84206
  135.  
  136. Symbols used in Plots
  137.  
  138. Symbol  Group  Label
  139. ------  -----  --------------------
  140.   1         1  Die                 
  141.   2         2  Survive             
  142. -------------------------------------------------------------------------------
  143. Page   9   Test Job for SPSS/PC+ ADVANCED STATS                        10/2/90 
  144.                 All-groups stacked Histogram
  145.  
  146.               Canonical Discriminant Function 1
  147.        4 ┼              2   1 2                    ┼
  148.          │              2   1 2                    │
  149.          │              2   1 2                    │
  150.          │              2   1 2                    │
  151.        3 ┼              1   11222 2  2             ┼
  152.          │              1   11222 2  2             │
  153.          │              1   11222 2  2             │
  154.          │              1   11222 2  2             │
  155.        2 ┼         11 211  11111222222             ┼
  156.          │         11 211  11111222222             │
  157.          │         11 211  11111222222             │
  158.          │         11 211  11111222222             │
  159.        1 ┼         11 1111 111112211222   22  2    ┼
  160.          │         11 1111 111112211222   22  2    │
  161.          │         11 1111 111112211222   22  2    │
  162.          │         11 1111 111112211222   22  2    │
  163.           X────┼────┼────┼────┼────┼────┼────┼────X
  164.         Out -3.0 -2.0 -1.0   .0  1.0  2.0  3.0    Out
  165.   Class   11111111111111111111122222222222222222222
  166. Centroids                  1       2
  167. -------------------------------------------------------------------------------
  168. Page  10   Test Job for SPSS/PC+ ADVANCED STATS                        10/2/90 
  169.  
  170. Classification Results -
  171.  
  172.                       No. of    Predicted Group Membership
  173.    Actual Group        Cases          1          2
  174. --------------------  ------   --------   --------
  175. Group       1             26         22          4
  176. Die                                84.6%      15.4%
  177.  
  178. Group       2             22          4         18
  179. Survive                            18.2%      81.8%
  180.  
  181. Percent of "grouped" cases correctly classified:  83.33%
  182.  
  183. Classification Processing Summary
  184.         50 Cases were processed.
  185.          0 Cases were excluded for missing or out-of-range group codes.
  186.          2 Cases had at least one missing discriminating variable.
  187.         48 Cases were used for printed output.
  188. -------------------------------------------------------------------------------
  189. Page  11   Test Job for SPSS/PC+ ADVANCED STATS                        10/2/90 
  190.  
  191. This procedure was completed at 14:29:23
  192.  
  193. * next procedure requires that D module be installed.
  194.  
  195. LOGISTIC REGRESSION SURVIVAL WITH TREATMNT AGE SEX WEIGHT APGAR.
  196.  
  197.  
  198. -------------------------------------------------------------------------------
  199. Page  12   Test Job for SPSS/PC+ ADVANCED STATS                        10/2/90 
  200.       Total number of cases:      50 (Unweighted)
  201.       Number of selected cases:   50
  202.       Number of unselected cases: 0
  203.  
  204.       Number of selected cases:                 50
  205.       Number rejected because of missing data:  0
  206.       Number of cases included in the analysis: 50
  207.  
  208.  
  209.  
  210. Dependent Variable Encoding:
  211.  
  212. Original       Internal
  213. Value          Value
  214. 1              0
  215. 2              1
  216. -------------------------------------------------------------------------------
  217. Page  13   Test Job for SPSS/PC+ ADVANCED STATS                        10/2/90 
  218. Dependent Variable..   SURVIVAL   Infant survival
  219.  
  220. Beginning Block Number  0.  Initial Log Likelihood Function
  221.  
  222. -2 Log Likelihood   68.994376
  223.  
  224. * Constant is included in the model.
  225.  
  226.  
  227. Beginning Block Number  1.  Method: Enter
  228.  
  229. Variable(s) Entered on Step Number
  230. 1..       TREATMNT  Treatment administered
  231.           AGE       Gestation age
  232.           SEX       Sex of respondent
  233.           WEIGHT    Birthweight in kilograms
  234.           APGAR     APGAR score
  235.  
  236. Estimation terminated at iteration number 4 because
  237. Log Likelihood decreased by less than .01 percent.
  238.  
  239. -------------------------------------------------------------------------------
  240. Page  14   Test Job for SPSS/PC+ ADVANCED STATS                        10/2/90 
  241.                      Chi-Square    df Significance
  242.  -2 Log Likelihood       53.751    44        .1489
  243.  Model Chi-Square        15.243     5        .0094
  244.  Improvement             15.243     5        .0094
  245.  Goodness of Fit         44.541    44        .4489
  246.  
  247. Classification Table for SURVIVAL
  248.                         Predicted
  249.                       Die    Survive    Percent Correct
  250.                         D   │    S
  251. Observed           ┼────────┼────────┼
  252.    Die         D   │   22   │    5   │   81.48%
  253.                    ┼────────┼────────┼
  254.    Survive     S   │    9   │   14   │   60.87%
  255.                    ┼────────┼────────┼
  256.                                 Overall  72.00%
  257.  
  258. -------------------------------------------------------------------------------
  259. Page  15   Test Job for SPSS/PC+ ADVANCED STATS                        10/2/90 
  260. ---------------------- Variables in the Equation -----------------------
  261.  
  262. Variable           B      S.E.     Wald    df      Sig       R   Exp(B)
  263.  
  264. TREATMNT       .8414     .6804   1.5295     1    .2162   .0000   2.3197
  265. AGE           -.2907     .2419   1.4447     1    .2294   .0000    .7477
  266. SEX           -.4179     .7677    .2963     1    .5862   .0000    .6584
  267. WEIGHT        3.1243    1.3639   5.2472     1    .0220   .2169  22.7438
  268. APGAR          .0704     .1352    .2708     1    .6028   .0000   1.0729
  269. Constant      2.7913    5.8809    .2253     1    .6350
  270. -------------------------------------------------------------------------------
  271. Page  16   Test Job for SPSS/PC+ ADVANCED STATS                        10/2/90 
  272.  
  273. This procedure was completed at 14:29:38
  274.  
  275. * next procedure requires that E module be installed.
  276. TITLE 'PROBIT Analysis'.
  277. DATA LIST / Root 1 X 3-7 N 9-11 R 13-15.
  278. BEGIN DATA
  279. END DATA.
  280.      10 cases are written to the compressed active file.
  281.  
  282. This procedure was completed at 14:29:42
  283.  
  284. PROBIT R of N by ROOT(1,2) WITH X
  285.   /MODEL = BOTH
  286.   /NATRES
  287.  /PRINT ALL.
  288. PROBIT requires 776 BYTES of workspace for execution.   
  289. -------------------------------------------------------------------------------
  290. Page  17   PROBIT Analysis                                             10/2/90 
  291.     
  292. * * * * * * * * * * *  P R O B I T    A N A L Y S I S  * * * * * * * * * * *
  293.     
  294. DATA  Information   
  295.     
  296.           9 unweighted cases accepted.  
  297.           0 cases rejected because of out-of-range group values.
  298.           0 cases rejected because of missing data. 
  299.           0 Cases rejected because LOG-transform can't be done. 
  300.     
  301.     
  302. Group Information   
  303.     
  304.     ROOT      Level  N of Cases    Label
  305.                   1           4        1        
  306.                   2           5        2        
  307.     
  308.     
  309. MODEL Information   
  310.     
  311.       BOTH Probit and Logit models are requested.   
  312.     
  313.     
  314. Natural Response rate to be estimated   
  315.     
  316.     The number of  subjects in the CONTROL group   129.0
  317.     The number of responses in the CONTROL group    21.0
  318. -------------------------------------------------------------------------------
  319. Page  18   PROBIT Analysis                                             10/2/90 
  320.     
  321. - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
  322.     
  323. -------------------------------------------------------------------------------
  324. Page  19   PROBIT Analysis                                             10/2/90 
  325.  
  326.                                   LOG(X) VS. LOGIT                             
  327.                                                                                
  328.           -+----+----+----+----+----+----+----+----+----+----+----+----+--     
  329.           :                                                              :     
  330.           :                                                    1         :     
  331.  L        :                                        2                     :     
  332.  O        :                                                              :     
  333.  G      6 +                                   2     1                    +     
  334.  I        :                        2                                     :     
  335.  T        :                                                              :     
  336.           :                     1                                        :     
  337.           :                                                              :     
  338.         4 +    2                                                         +     
  339.           :                                                              :     
  340.           -+----+----+----+----+----+----+----+----+----+----+----+----+--     
  341.          .45       .75      1.05      1.35      1.65      1.95      2.25       
  342.                .6        .9       1.2       1.5       1.8       2.1            
  343.                                                                                
  344.                                        LOG(X)                                  
  345. -------------------------------------------------------------------------------
  346. Page  20   PROBIT Analysis                                             10/2/90 
  347.     
  348.  ML converged at iteration  6.  The converge criterion =   .00095   
  349.     
  350.  Parameter Estimates (LOGIT model:  (LOG(p/(1-p))/2 + 5) = Intercept + BX ):
  351.                       Note 5 added to intercept and logit divided by 2. 
  352.     
  353.            Regression Coeff.  Standard Error     Coeff./S.E.
  354.     
  355.    X                 2.45569          .19850        12.37135
  356.     
  357.     
  358.                    Intercept  Standard Error  Intercept/S.E.  ROOT      
  359.     
  360.                      2.02382          .30064         6.73182         1         
  361.                      2.59473          .28483         9.10988         2         
  362.     
  363.   Estimate of Natural Response Rate = .172976  with  S.E. =    .03147   
  364.     
  365.   Pearson  Goodness-of-Fit  Chi Square =      8.966    DF = 6   P =  .176   
  366.            PARALLELISM TEST CHI SQUARE =       .048    DF = 1   P =  .826   
  367.     
  368.   Since Goodness-of-Fit Chi square is NOT significant, no heterogeneity 
  369.   factor is used in the calculation of confidence limits.   
  370. -------------------------------------------------------------------------------
  371. Page  21   PROBIT Analysis                                             10/2/90 
  372. * * * * * * * * * * *  P R O B I T    A N A L Y S I S  * * * * * * * * * * *
  373.  Covariance(below) and Correlation(above) Matrices of Parameter Estimates   
  374.     
  375.                      X   NAT RESP   
  376.     
  377.  X              .03940     .42344   
  378.  NAT RESP       .00264     .00099   
  379. -------------------------------------------------------------------------------
  380. Page  22   PROBIT Analysis                                             10/2/90 
  381.     
  382.  Observed and Expected Frequencies  
  383.     
  384.     
  385.                        Number of    Observed    Expected
  386.   ROOT              X   Subjects   Responses   Responses    Residual     Prob  
  387.     
  388.       1          2.17      142.0       142.0     140.728       1.272   .99105  
  389.       1          2.00      127.0       126.0     124.406       1.594   .97958  
  390.       1          1.68      128.0       115.0     116.388      -1.388   .90928  
  391.       1          1.08      126.0        58.0      43.158      14.842   .34252  
  392.     
  393.       2          1.79      125.0       125.0     122.735       2.265   .98188  
  394.       2          1.66      117.0       115.0     113.057       1.943   .96630  
  395.       2          1.49      127.0       114.0     117.491      -3.491   .92512  
  396.       2          1.17       51.0        40.0      36.644       3.356   .71850  
  397.       2           .58      132.0        37.0      16.255      20.745   .12314  
  398.     
  399. -------------------------------------------------------------------------------
  400. Page  23   PROBIT Analysis                                             10/2/90 
  401.     
  402.  Confidence Limits for Effective X          
  403.     
  404.     
  405.   ROOT        1         1           
  406.     
  407.                             95% Confidence Limits   
  408.   Prob             X         Lower         Upper
  409.     
  410.    .01       1.88949       1.09495       2.83761
  411.    .02       2.62750       1.61246       3.78991
  412.    .03       3.19291       2.02582       4.49950
  413.    .04       3.67174       2.38515       5.09032
  414.    .05       4.09673       2.71020       5.60842
  415.    .06       4.48451       3.01123       6.07679
  416.    .07       4.84483       3.29438       6.50879
  417.    .08       5.18404       3.56370       6.91298
  418.    .09       5.50649       3.82200       7.29523
  419.    .10       5.81537       4.07136       7.65976
  420.    .15       7.22389       5.22909       9.30609
  421.    .20       8.50530       6.30599      10.78770
  422.    .25       9.73339       7.35390      12.19915
  423.    .30      10.95051       8.40392      13.59378
  424.    .35      12.18736       9.47954      15.00993
  425. -------------------------------------------------------------------------------
  426. Page  24   PROBIT Analysis                                             10/2/90 
  427.    .40      13.47086      10.60208      16.48107
  428.    .45      14.82836      11.79372      18.04134
  429.    .50      16.29114      13.08011      19.73008
  430.    .55      17.89822      14.49337      21.59690
  431.    .60      19.70188      16.07639      23.70910
  432.    .65      21.77676      17.88994      26.16406
  433.    .70      24.23644      20.02552      29.11216
  434.    .75      27.26708      22.63118      32.80390
  435.    .80      31.20422      25.96981      37.69896
  436.    .85      36.73936      30.57370      44.76500
  437.    .90      45.63792      37.77256      56.53912
  438.    .91      48.19789      39.80239      60.01242
  439.    .92      51.19587      42.15862      64.12498
  440.    .93      54.78029      44.94814      69.10254
  441.    .94      59.18183      48.33569      75.29989
  442.    .95      64.78363      52.59221      83.31382
  443.    .96      72.28220      58.20465      94.24584
  444.    .97      83.12198      66.16875     110.42335
  445.    .98     101.00910      78.99264     137.96693
  446.    .99     140.46189     106.22774     201.78711
  447. -------------------------------------------------------------------------------
  448. Page  25   PROBIT Analysis                                             10/2/90 
  449.     
  450.  Confidence Limits for Effective X          
  451.     
  452.     
  453.   ROOT        2         2           
  454.     
  455.                             95% Confidence Limits   
  456.   Prob             X         Lower         Upper
  457.     
  458.    .01       1.10627        .60305       1.73822
  459.    .02       1.53837        .88826       2.32107
  460.    .03       1.86941       1.11620       2.75507
  461.    .04       2.14976       1.31446       3.11617
  462.    .05       2.39859       1.49390       3.43263
  463.    .06       2.62562       1.66017       3.71854
  464.    .07       2.83659       1.81665       3.98208
  465.    .08       3.03519       1.96556       4.22851
  466.    .09       3.22398       2.10846       4.46142
  467.    .10       3.40483       2.24647       4.68340
  468.    .15       4.22950       2.88824       5.68417
  469.    .20       4.97975       3.48680       6.58206
  470.    .25       5.69878       4.07090       7.43472
  471.    .30       6.41139       4.65794       8.27438
  472.    .35       7.13555       5.26127       9.12395
  473. -------------------------------------------------------------------------------
  474. Page  26   PROBIT Analysis                                             10/2/90 
  475.    .40       7.88702       5.89319      10.00308
  476.    .45       8.68182       6.56666      10.93157
  477.    .50       9.53826       7.29687      11.93193
  478.    .55      10.47919       8.10303      13.03228
  479.    .60      11.53521       9.01097      14.27051
  480.    .65      12.75003      10.05751      15.70113
  481.    .70      14.19014      11.29835      17.40806
  482.    .75      15.96454      12.82395      19.53073
  483.    .80      18.26969      14.79521      22.32464
  484.    .85      21.51044      17.53817      26.32751
  485.    .90      26.72044      21.86628      32.95019
  486.    .91      28.21927      23.09232      34.89717
  487.    .92      29.97455      24.51759      37.20002
  488.    .93      32.07319      26.20723      39.98453
  489.    .94      34.65024      28.26152      43.44843
  490.    .95      37.93002      30.84530      47.92442
  491.    .96      42.32034      34.25466      54.02659
  492.    .97      48.66690      39.09473      63.05261
  493.    .98      59.13959      46.88897      78.41486
  494.    .99      82.23871      63.43687     113.99793
  495. -------------------------------------------------------------------------------
  496. Page  27   PROBIT Analysis                                             10/2/90 
  497.     
  498.  Estimates of Relative Median Potency   
  499.     
  500.     
  501.                             95% Confidence Limits   
  502.   ROOT        Estimate       Lower         Upper
  503.     
  504.   1 VS.  2      1.7080     1.31475       2.30747
  505. -------------------------------------------------------------------------------
  506. Page  28   PROBIT Analysis                                             10/2/90 
  507.  
  508.                                   LOG(X) VS. PROBIT                            
  509.                                                                                
  510.           -+----+----+----+----+----+----+----+----+----+----+----+----+--     
  511.           :                                                              :     
  512.           :                                        2           1         :     
  513.  P        :                                                              :     
  514.  R        :                                                              :     
  515.  O      6 +                                   2     1                    +     
  516.  B        :                        2                                     :     
  517.  I        :                                                              :     
  518.  T        :                     1                                        :     
  519.           :                                                              :     
  520.         4 +    2                                                         +     
  521.           :                                                              :     
  522.           -+----+----+----+----+----+----+----+----+----+----+----+----+--     
  523.          .45       .75      1.05      1.35      1.65      1.95      2.25       
  524.                .6        .9       1.2       1.5       1.8       2.1            
  525.                                                                                
  526.                                        LOG(X)                                  
  527. -------------------------------------------------------------------------------
  528. Page  29   PROBIT Analysis                                             10/2/90 
  529.     
  530.  ML converged at iteration  6.  The converge criterion =   .00063   
  531.     
  532.  Parameter Estimates (PROBIT model:  (PROBIT(p) + 5) = Intercept + BX ):
  533.                       Note 5 added to intercept.
  534.     
  535.            Regression Coeff.  Standard Error     Coeff./S.E.
  536.     
  537.    X                 2.78382          .20654        13.47809
  538.     
  539.     
  540.                    Intercept  Standard Error  Intercept/S.E.  ROOT      
  541.     
  542.                      1.59842          .33212         4.81282         1         
  543.                      2.26158          .30366         7.44765         2         
  544.     
  545.   Estimate of Natural Response Rate = .168731  with  S.E. =    .03177   
  546.     
  547.   Pearson  Goodness-of-Fit  Chi Square =      5.933    DF = 6   P =  .431   
  548.            PARALLELISM TEST CHI SQUARE =       .050    DF = 1   P =  .823   
  549.     
  550.   Since Goodness-of-Fit Chi square is NOT significant, no heterogeneity 
  551.   factor is used in the calculation of confidence limits.   
  552. -------------------------------------------------------------------------------
  553. Page  30   PROBIT Analysis                                             10/2/90 
  554. * * * * * * * * * * *  P R O B I T    A N A L Y S I S  * * * * * * * * * * *
  555.  Covariance(below) and Correlation(above) Matrices of Parameter Estimates   
  556.     
  557.                      X   NAT RESP   
  558.     
  559.  X              .04266     .46719   
  560.  NAT RESP       .00307     .00101   
  561. -------------------------------------------------------------------------------
  562. Page  31   PROBIT Analysis                                             10/2/90 
  563.     
  564.  Observed and Expected Frequencies  
  565.     
  566.     
  567.                        Number of    Observed    Expected
  568.   ROOT              X   Subjects   Responses   Responses    Residual     Prob  
  569.     
  570.       1          2.17      142.0       142.0     141.411        .589   .99586  
  571.       1          2.00      127.0       126.0     125.075        .925   .98485  
  572.       1          1.68      128.0       115.0     115.136       -.136   .89950  
  573.       1          1.08      126.0        58.0      43.540      14.460   .34556  
  574.     
  575.       2          1.79      125.0       125.0     123.477       1.523   .98782  
  576.       2          1.66      117.0       115.0     113.566       1.434   .97065  
  577.       2          1.49      127.0       114.0     116.994      -2.994   .92121  
  578.       2          1.17       51.0        40.0      35.611       4.389   .69825  
  579.       2           .58      132.0        37.0      17.215      19.785   .13042  
  580.     
  581. -------------------------------------------------------------------------------
  582. Page  32   PROBIT Analysis                                             10/2/90 
  583.     
  584.  Confidence Limits for Effective X          
  585.     
  586.     
  587.   ROOT        1         1           
  588.     
  589.                             95% Confidence Limits   
  590.   Prob             X         Lower         Upper
  591.     
  592.    .01       2.43358       1.49073       3.54170
  593.    .02       3.04907       1.93470       4.32588
  594.    .03       3.51801       2.28205       4.91252
  595.    .04       3.91772       2.58346       5.40658
  596.    .05       4.27613       2.85742       5.84556
  597.    .06       4.60690       3.11310       6.24774
  598.    .07       4.91792       3.35578       6.62356
  599.    .08       5.21417       3.58886       6.97966
  600.    .09       5.49906       3.81465       7.32052
  601.    .10       5.77505       4.03482       7.64936
  602.    .15       7.07305       5.08714       9.18084
  603.    .20       8.30973       6.11123      10.62212
  604.    .25       9.54166       7.14783      12.04574
  605.    .30      10.80286       8.22256      13.49459
  606.    .35      12.11989       9.35637      15.00163
  607. -------------------------------------------------------------------------------
  608. Page  33   PROBIT Analysis                                             10/2/90 
  609.    .40      13.51780      10.56970      16.59764
  610.    .45      15.02357      11.88506      18.31571
  611.    .50      16.66914      13.32937      20.19508
  612.    .55      18.49495      14.93677      22.28584
  613.    .60      20.55513      16.75274      24.65563
  614.    .65      22.92596      18.84086      27.40086
  615.    .70      25.72100      21.29482      30.66690
  616.    .75      29.12076      24.26177      34.68812
  617.    .80      33.43794      27.99319      39.87663
  618.    .85      39.28436      32.97398      47.05260
  619.    .90      48.11392      40.33718      58.20447
  620.    .91      50.52866      42.31938      61.31623
  621.    .92      53.28948      44.57022      64.90463
  622.    .93      56.49958      47.16756      69.11657
  623.    .94      60.31388      50.22764      74.17419
  624.    .95      64.97940      53.93463      80.43461
  625.    .96      70.92392      58.60509      88.52216
  626.    .97      78.98231      64.85142      99.66931
  627.    .98      91.12942      74.10611     116.83393
  628.    .99     114.17771      91.23085     150.43809
  629. -------------------------------------------------------------------------------
  630. Page  34   PROBIT Analysis                                             10/2/90 
  631.     
  632.  Confidence Limits for Effective X          
  633.     
  634.     
  635.   ROOT        2         2           
  636.     
  637.                             95% Confidence Limits   
  638.   Prob             X         Lower         Upper
  639.     
  640.    .01       1.40612        .83171       2.10009
  641.    .02       1.76176       1.07972       2.56433
  642.    .03       2.03271       1.27386       2.91142
  643.    .04       2.26367       1.44239       3.20359
  644.    .05       2.47075       1.59564       3.46308
  645.    .06       2.66188       1.73871       3.70072
  646.    .07       2.84158       1.87455       3.92270
  647.    .08       3.01275       2.00506       4.13295
  648.    .09       3.17737       2.13153       4.33413
  649.    .10       3.33683       2.25489       4.52814
  650.    .15       4.08682       2.84506       5.43077
  651.    .20       4.80137       3.42031       6.27872
  652.    .25       5.51318       4.00355       7.11474
  653.    .30       6.24190       4.60933       7.96390
  654.    .35       7.00289       5.24966       8.84530
  655. -------------------------------------------------------------------------------
  656. Page  35   PROBIT Analysis                                             10/2/90 
  657.    .40       7.81060       5.93636       9.77656
  658.    .45       8.68064       6.68260      10.77649
  659.    .50       9.63145       7.50419      11.86722
  660.    .55      10.68641       8.42131      13.07685
  661.    .60      11.87678       9.46096      14.44321
  662.    .65      13.24665      10.66101      16.02005
  663.    .70      14.86163      12.07738      17.88828
  664.    .75      16.82601      13.79805      20.17822
  665.    .80      19.32048      15.97345      23.11903
  666.    .85      22.69855      18.89332      27.16715
  667.    .90      27.80029      23.23350      33.43063
  668.    .91      29.19553      24.40516      35.17470
  669.    .92      30.79073      25.73678      37.18464
  670.    .93      32.64553      27.27463      39.54249
  671.    .94      34.84944      29.08777      42.37235
  672.    .95      37.54519      31.28554      45.87376
  673.    .96      40.97994      34.05583      50.39566
  674.    .97      45.63608      37.76201      56.62695
  675.    .98      52.65470      43.25369      66.22123
  676.    .99      65.97204      53.41329      85.00561
  677. -------------------------------------------------------------------------------
  678. Page  36   PROBIT Analysis                                             10/2/90 
  679.     
  680.  Estimates of Relative Median Potency   
  681.     
  682.     
  683.                             95% Confidence Limits   
  684.   ROOT        Estimate       Lower         Upper
  685.     
  686.   1 VS.  2      1.7307     1.34843       2.30803
  687. -------------------------------------------------------------------------------
  688. Page  37   PROBIT Analysis                                             10/2/90 
  689.  
  690. This procedure was completed at 14:30:04
  691.  
  692.  
  693.  
  694. Title Non-linear model for windspeed measurement.
  695.  
  696. DATA LIST / X 1-3 Y 5-9.
  697. VAR LABELS X 'anomometer height' Y 'windspeed'.
  698. BEGIN DATA
  699. END DATA.
  700.       5 cases are written to the compressed active file.
  701.  
  702. This procedure was completed at 14:30:07
  703.  
  704. * next procedure requires that E module be installed.
  705.  
  706. MODEL PROGRAM
  707.   A=100 B=7.06 C=-223.26.
  708. -------------------------------------------------------------------------------
  709. Page  38   NON-LINEAR MODEL FOR WINDSPEED MEASUREMENT                  10/2/90 
  710.  
  711. This procedure was completed at 14:30:11
  712. COMPUTE PRED=A*LN(B*X+C).
  713. NLR Y WITH X/PRED=PRED/SAVE PRED.
  714.  
  715. All the derivatives will be calculated numerically.
  716. -------------------------------------------------------------------------------
  717. Page  39   NON-LINEAR MODEL FOR WINDSPEED MEASUREMENT                  10/2/90 
  718.  
  719. There are 5 cases.  There is enough memory for them all.
  720.  
  721.  Iteration  Residual SS           A           B           C
  722.  
  723.      1      6822.703884  100.000000  7.06000000  -223.26000
  724.      1.1    3899.279813  107.235597  2.71055158  -18.043743
  725.      2      3899.279813  107.235597  2.71055158  -18.043743
  726.      2.1    107.8417730  115.003329  2.27328009  -24.273351
  727.      3      107.8417730  115.003329  2.27328009  -24.273351
  728.      3.1    7.040578209  115.150921  2.30999733  -22.101959
  729.      4      7.040578209  115.150921  2.30999733  -22.101959
  730.      4.1    7.013266157  115.147444  2.31056594  -22.026921
  731.      5      7.013266157  115.147444  2.31056594  -22.026921
  732.      5.1    7.013265848  115.146862  2.31064751  -22.028813
  733.      6      7.013265848  115.146862  2.31064751  -22.028813
  734.      6.1    7.013265848  115.146870  2.31064641  -22.028787
  735.  
  736. Run stopped after 12 model evaluations and 6 derivative evaluations.
  737. Iterations have been stopped because the relative reduction between successive
  738. residual sums of squares is at most SSCON = 1.000E-08
  739. -------------------------------------------------------------------------------
  740. Page  40   NON-LINEAR MODEL FOR WINDSPEED MEASUREMENT                  10/2/90 
  741.  
  742.  
  743. Nonlinear Regression Summary Statistics     Dependent Variable Y
  744.  
  745.   Source                 DF  Sum of Squares  Mean Square
  746.  
  747.   Regression              3  2314527.69673   771509.23224
  748.   Residual                2        7.01327        3.50663
  749.   Uncorrected Total       5  2314534.71000
  750.  
  751.   (Corrected Total)       4    75525.34800
  752.  
  753.   R squared = 1 - Residual SS / Corrected SS =     .99991
  754.  
  755.                                            Asymptotic 95 %
  756.                           Asymptotic     Confidence Interval
  757.   Parameter   Estimate    Std. Error     Lower         Upper
  758.  
  759.   A         115.14686999  2.040556318 106.36706478 123.92667520
  760.   B          2.310646412   .280311752  1.104562288  3.516730537
  761.   C         -22.02878660  6.409417923 -49.60628612  5.548712919
  762. -------------------------------------------------------------------------------
  763. Page  41   NON-LINEAR MODEL FOR WINDSPEED MEASUREMENT                  10/2/90 
  764.  
  765.   Asymptotic Correlation Matrix of the Parameter Estimates
  766.  
  767.                    A         B         C
  768.  
  769.   A           1.0000    -.9963     .9666
  770.   B           -.9963    1.0000    -.9802
  771.   C            .9666    -.9802    1.0000
  772. -------------------------------------------------------------------------------
  773. Page  42   NON-LINEAR MODEL FOR WINDSPEED MEASUREMENT                  10/2/90 
  774.  
  775. This procedure was completed at 14:30:20
  776. STOP.
  777.  
  778. End of Include file.
  779.