home *** CD-ROM | disk | FTP | other *** search
/ Liren Large Software Subsidy 10 / 10.iso / l / l292 / 1.ddi / STATSDEM.FOR < prev    next >
Encoding:
Text File  |  1989-10-05  |  1.9 KB  |  60 lines
  1.       PROGRAM STATSDEM
  2.       INCLUDE 'STDHDR.FOR'
  3.       REAL storedata(0:maxr, 0:maxc)
  4.       REAL minima(0:maxc), maxima(0:maxc)
  5.       REAL range(0:maxc), sumxx(0:maxc)
  6.       REAL mean(0:maxc), variance(0:maxc), stddev(0:maxc)
  7.       REAL semean(0:maxc), mode(0:maxc)
  8.       INTEGER numobs, numiv
  9.  
  10.       storedata(0, 0) = 305.0
  11.       storedata(0, 1) = 35.0
  12.       storedata(1, 0) = 130.0
  13.       storedata(1, 1) = 98.0
  14.       storedata(2, 0) = 189.0
  15.       storedata(2, 1) = 83.0
  16.       storedata(3, 0) = 175.0
  17.       storedata(3, 1) = 76.0
  18.       storedata(4, 0) = 101.0
  19.       storedata(4, 1) = 93.0
  20.       storedata(5, 0) = 269.0
  21.       storedata(5, 1) = 77.0
  22.       storedata(6, 0) = 421.0
  23.       storedata(6, 1) = 44.0
  24.       storedata(7, 0) = 195.0
  25.       storedata(7, 1) = 57.0
  26.       storedata(8, 0) = 282.0
  27.       storedata(8, 1) = 31.0
  28.       storedata(9, 0) = 203.0
  29.       storedata(9, 1) = 92.0
  30.  
  31.       numobs = 10
  32.       numiv = 2
  33.  
  34.       CALL SummaryStats(storedata, numobs, numiv, minima,
  35.      + maxima, range, sumxx, mean, variance, stddev, semean, mode)
  36.  
  37.       WRITE (*,*) 'Summary statistics for storedata'
  38.       DO i = 0, numiv - 1
  39.         WRITE (*,*) '        Variable #', i + 1
  40.         WRITE (*,*) 'Number of observations =', numobs
  41.         WRITE (*,10)  'minima =                  ', minima(i)
  42.         WRITE (*,10)  'maxima =                  ', maxima(i)
  43.         WRITE (*,10)  'range =                   ', range(i)
  44.         WRITE (*,10)  'mean =                    ', mean(i)
  45.         WRITE (*,10)  'sum =                     ', sumxx(i)
  46.         WRITE (*,10)  'variance =                ', variance(i)
  47.         WRITE (*,10)  'standard deviation =      ', stddev(i)
  48.         WRITE (*,10)  'standard error of mean =  ', semean(i)
  49.         WRITE (*,10)  'mode =                    ', mode(i)
  50. 10      FORMAT (A28,f9.2)
  51.       END DO
  52.       WRITE (*,*) 'Press carriage return to end program. '
  53.       READ (*,*)
  54.  
  55.  
  56.       END
  57.  
  58.  
  59.  
  60.