home *** CD-ROM | disk | FTP | other *** search
/ Monster Media 1993 #2 / Image.iso / math / ksprob21.zip / KSDOCS.EXE / KSPDNS.EXE / NOT.DOC < prev    next >
Text File  |  1990-05-28  |  7KB  |  127 lines
  1.        Creating a series of graphics screens. Joseph C. Hudson
  2.        
  3.        The screens in not1.com - not6.com illustrate an idea in 
  4.        elementary statistics. I will briefly explain the idea, how 
  5.        the screens were created, and how they can be used.
  6.        
  7.        1. The idea
  8.        
  9.        Tables of Student's t percentage points often stop at 30 or so 
  10.        df, with a last line giving standard normal percentiles. It is 
  11.        implied that when df is above 30, normal percentiles and p values 
  12.        are adequate approximations to those of the Student's t. The 
  13.        graphs in not1.scn through not6.scn compare the normal and t 
  14.        distributions, with not4.scn, not5.scn and not6.scn showing 
  15.        details of the 32 df t and normal relevant to the impliction 
  16.        mentioned above.
  17.        
  18.        The series of graphs in not1.com through not6.com explore the
  19.        relationship between the t and normal distributions, with a 
  20.        detailed look at the right tail areas of the normal and 32 df t. 
  21.        
  22.        Since programs are available to compute Student's t percentage 
  23.        points and p values for any df, Why use an approximation that, 
  24.        for "close calls", can lead to the wrong conclusion? The normal 
  25.        approximation to the t is not conservative. Better rational 
  26.        function approximations can be had - see ksprbas.bas for an 
  27.        example of one.
  28.        
  29.        2. how these files were created.
  30.        
  31.        Once the idea is clear the screen contents can be decided on: 
  32.        
  33.        screen contents
  34.         not1  std normal pdf and Student's t pdf with  2 df
  35.         not2  std normal pdf and Student's t pdf with  8 df
  36.         not3  std normal pdf and Student's t pdf with 32 df
  37.         not4  std normal and 32 df t reliab fcns (1-cdf) -5   < x < 5
  38.         not5  std normal and 32 df t reliab fcns (1-cdf)  1.5 < x < 3
  39.         not6  % diff between the fcns shown in not5, 100 * (t - n) / t
  40.        
  41.        The functions plotted on these screens cannot be represented by 
  42.        simple formulas that plotting programs can take as input. It was 
  43.        necessary, then, to use kspdat to generate data sets containing 
  44.        values of the functions and then plot the data. 
  45.        
  46.        I first laid out the relationships between the graph files and 
  47.        data files:
  48.        
  49.        graph file   data files needed          The *.c## are produced   
  50.           not1      not.d01, not.d02           by kspdat - they give
  51.           not2      not.d01, not.d03           column names and missing
  52.           not3      not.d01, not.d04           value info for the cor-
  53.           not4      not.d05, not.d06           esponding .d## files.
  54.           not5      not.d07, not.d08           The .c## files are not 
  55.           not6      not.d09                    needed here. 
  56.        
  57.        You might wonder why not.d01 - not.d04 are separate files. The 
  58.        kspdat program can easily create a 5 column data file with the 
  59.        contents of these files as columns. I would have created just one 
  60.        file insetad of four if the graphing program I used could handle 
  61.        it that way. It can't. Hence the separate files. (The program - 
  62.        gnuplot - is really great. occasional quirks like this are not 
  63.        enough to stop me from using it.)
  64.        
  65.        There are three .spc files that can be read in by kspdat to 
  66.        produce these output files. notp.spc has the specs for not.d01 - 
  67.        not.do4. notr1.spc has the specs for not.d05 and not.d06. 
  68.        notr2.spc has the specs for not.d07 and not.d08. 
  69.        
  70.        not.d09 has the common column 1 from not.d07 and not.d08 and a 
  71.        second column created by combining the second columns in these 
  72.        two files. Since most spreadsheets and general purpose statistics 
  73.        programs can do this kind of thing and output not.d09, I won't 
  74.        describe the (hopefully) obvious here, except to say that if you 
  75.        output one table instead of two using notr2.spc in kspdat, you 
  76.        will have a three column data file with the contents of not.d07 
  77.        and not.d08 that can be read by most stat programs.
  78.        
  79.        With the data files in hand, it only remains to produce the 
  80.        graphs and the executable files that contain them. For this, I 
  81.        use gnuplot, version 2.0, and grabber. grabber produces the .com 
  82.        files from the graphics screens. It is shareware and the current 
  83.        version costs about $25. gnuplot is copyrighted but free. I don't 
  84.        see it distributed in the usual shareware channels, but it seems 
  85.        to be available on academic bbs's. If you can't find it anywhere, 
  86.        send me a postage-paid return mailer and a blank formatted disk 
  87.        and I'll be glad to send you the current version I have. (Please 
  88.        don't send money - I don't want to get involved in any copyright 
  89.        infringements, no matter how innocent.)
  90.        
  91.        The file not.gnu contains gnuplot commands to produce the 6 
  92.        graphs. To produce the .com files, I
  93.        
  94.        a. set my vga into cga emulation.
  95.        
  96.        b. loaded grabber.
  97.        
  98.        c. ran gnuplot, used the command load"b:not.gnu to execute the 
  99.        commands in not.gnu (carriage returning past a meaningless grerr 
  100.        message at first).
  101.        
  102.        d. as each screen appeared, used grabber's hot key to capture the 
  103.        screen in a .com file.
  104.        
  105.        e. hit a carriage return to go on to the next screen.
  106.        
  107.        The current version of grabber can do vga screen captures (The 
  108.        copy I have (and registered) does not) and gnuplot 2.0 does a 
  109.        gorgeous job in vga mode, so If you have vga, use it for this. 
  110.        The quality is much much better. 
  111.        
  112.        3. How to use the results
  113.        
  114.        You can, of course, simply type the name of the .com files to 
  115.        display the screens. I find that using dirmagic as a controller 
  116.        program makes a nice presentation manager. Dirmagic presents a 
  117.        menu of the files and allows them to be executed by pressing 
  118.        ctrl-enter.
  119.        
  120.        You can use a cga screen capture program to extract the screens 
  121.        from the .com files, and then use the screens in a slide show 
  122.        presentation program, like present or flicker. 
  123.        
  124.        It might be nice to automate this process so that a program could 
  125.        take the commands of kspdat and quickly produce the graphs 
  126.        directly. Maybe some day.
  127.