home *** CD-ROM | disk | FTP | other *** search
/ Liren Large Software Subsidy 7 / 07.iso / c / c129 / 1.ddi / REGDEM.C < prev    next >
Encoding:
C/C++ Source or Header  |  1989-08-18  |  2.6 KB  |  96 lines
  1. # include <stdio.h>
  2. # include <stdlib.h>
  3. # include <math.h>
  4. # include "gj.h"
  5. # include "mulreg.h"
  6. # include "miscio.h"
  7.  
  8. float storedata[10][2];
  9. float sales[10];
  10. float eqncoef[3];
  11. float coefsig[3];
  12. int numobs;
  13. int numiv;
  14. int i;
  15. char error;
  16. float yest[10];
  17. float resid[10];
  18. float rsq;
  19. float r;
  20. float see;
  21. char varstring[8][32];
  22.  
  23.  
  24. void DisplayRegData()
  25. {
  26.    int i;
  27.    int j;
  28.  
  29.    ClrScr();
  30.    printf("multiple regression - store sales vs district income and population\n");
  31.    printf("\n");
  32.    printf("            actual     district      district   estimated            \n");
  33.    printf("index       sales       income      population    sales     residuals\n");
  34.    printf("=====       =====      ========     ==========  =========   =========\n");
  35.    for ( i = 0; i <= numobs-1; ++i ) {
  36.       printf("%2d         %8.2f    %8.2f    %8.2f    %8.2f    %8.2f\n", i,
  37.       sales[i],  storedata[i][0], storedata[i][1],
  38.       yest[i], resid[i]);
  39.    }
  40.    printf("\n");
  41.    printf("estimated regression equation is: \n");
  42.    printf("%s %s\n", varstring[0] ," (y) =    variable   x   coefficient         error" );
  43.    printf( "                ========       ===========         =====\n");
  44.    for ( i = 0; i <= 2; ++i ) {
  45.       printf(" %s %2d %s   %10s      %8.2f         %8.2f\n", " x( ",i ,
  46.        " ) ",(varstring[i]),  eqncoef[i],  coefsig[i]);
  47.    }
  48.    printf(" correlation coef (r ) = %4.3f  coefficient OF determination rsqr = %4.3f\n",  r,  rsq);
  49.    printf("standard error of estimate (see) = %10.2f\n", see);
  50.    printf( "press carriage return to end program");
  51. }
  52.  
  53.  
  54.  
  55. void main()
  56. {
  57.    storedata[0][0] = 305.0;
  58.    storedata[0][1] = 35.0;
  59.    storedata[1][0] = 130.0;
  60.    storedata[1][1] = 98.0;
  61.    storedata[2][0] = 189.0;
  62.    storedata[2][1] = 83.0;
  63.    storedata[3][0] = 175.0;
  64.    storedata[3][1] = 76.0;
  65.    storedata[4][0]= 101.0;
  66.    storedata[4][1]= 93.0;
  67.    storedata[5][0] = 269.0;
  68.    storedata[5][1] = 77.0;
  69.    storedata[6][0] = 421.0;
  70.    storedata[6][1] = 44.0;
  71.    storedata[7][0] = 195.0;
  72.    storedata[7][1] = 57.0;
  73.    storedata[8][0] = 282.0;
  74.    storedata[8][1] = 31.0;
  75.    storedata[9][0] = 203.0;
  76.    storedata[9][1] = 92.0;
  77.    sales[0]= 20.0;
  78.    sales[1]= 15.0;
  79.    sales[2] = 17.0;
  80.    sales[3] = 9.0;
  81.    sales[4] = 16.0;
  82.    sales[5] = 27.0;
  83.    sales[6] = 35.0;
  84.    sales[7] = 7.0;
  85.    sales[8] = 22.0;
  86.    sales[9] = 23.0;
  87.    strcpy(varstring[0], "sales");
  88.    strcpy(varstring[1], "constant");
  89.    strcpy(varstring[2], "income");
  90.    strcpy(varstring[3], "population");
  91.    numobs = 10;
  92.    numiv = 2;
  93.    MultipleReg(&storedata[0][0],sales,numiv,numobs,eqncoef,yest,resid,&see,coefsig,&rsq,&r,&error);
  94.    DisplayRegData();
  95. }
  96.