#!/usr/bin/perl
print "Content-type: text/html\n\n";
### (C) COPYRIGHT G. WILLIAM MOORE, MD, PhD, BALTIMORE, MD, 2003.
### toknswap.cgi: PERL script to perform the TOKEN SWAP TEST.
### Last modified: 9/21/2003, G. William Moore, MD, PhD.
###
### PRINT HEADER.
print qq| STATISTICS CALCULATOR. |;
print qq|\n|;
print qq|\n STATISTICS CALCULATOR. |;
print qq|\n
|;
print qq|\n http://www.netautopsy.org/stdnfezf.htm |;
print qq|\n
© 2003, G. William Moore, MD, PhD,|;
print qq|\n
Baltimore, MD, 2003. |;
print qq|\n
All Rights Reserved.
|;
###
### CAPTURE INPUT DATA FROM ENVIRONMENTAL VARIABLES.
$inpdt=lc($ENV{'QUERY_STRING'}); chop; $srqln=$amp.$inpdt.$amp;
###
### SPLIT THE INPUT QUERY BY AMPERSAND TO OBTAIN INPUT ARRAY.
($scr,$xhi)=split(/aaaa=/,$srqln,2); ($aaaa,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/bbbb=/,$srqln,2); ($bbbb,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/cccc=/,$srqln,2); ($cccc,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/dddd=/,$srqln,2); ($dddd,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/titl=/,$srqln,2); ($titl,$scr)=split(/\&/,$xhi,2);
$uctitl=uc($titl); print qq|\n
Title: $uctitl. |;
($scr,$xhi)=split(/x000=/,$srqln,2); ($x000,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/y000=/,$srqln,2); ($y000,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/x001=/,$srqln,2); ($x001,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/y001=/,$srqln,2); ($y001,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/x002=/,$srqln,2); ($x002,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/y002=/,$srqln,2); ($y002,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/x003=/,$srqln,2); ($x003,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/y003=/,$srqln,2); ($y003,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/x004=/,$srqln,2); ($x004,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/y004=/,$srqln,2); ($y004,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/x005=/,$srqln,2); ($x005,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/y005=/,$srqln,2); ($y005,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/x006=/,$srqln,2); ($x006,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/y006=/,$srqln,2); ($y006,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/x007=/,$srqln,2); ($x007,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/y007=/,$srqln,2); ($y007,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/x008=/,$srqln,2); ($x008,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/y008=/,$srqln,2); ($y008,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/x009=/,$srqln,2); ($x009,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/y009=/,$srqln,2); ($y009,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/x010=/,$srqln,2); ($x010,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/y010=/,$srqln,2); ($y010,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/x011=/,$srqln,2); ($x011,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/y011=/,$srqln,2); ($y011,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/x012=/,$srqln,2); ($x012,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/y012=/,$srqln,2); ($y012,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/x013=/,$srqln,2); ($x013,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/y013=/,$srqln,2); ($y013,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/x014=/,$srqln,2); ($x014,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/y014=/,$srqln,2); ($y014,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/x015=/,$srqln,2); ($x015,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/y015=/,$srqln,2); ($y015,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/x016=/,$srqln,2); ($x016,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/y016=/,$srqln,2); ($y016,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/x017=/,$srqln,2); ($x017,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/y017=/,$srqln,2); ($y017,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/x018=/,$srqln,2); ($x018,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/y018=/,$srqln,2); ($y018,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/x019=/,$srqln,2); ($x019,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/y019=/,$srqln,2); ($y019,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/x020=/,$srqln,2); ($x020,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/y020=/,$srqln,2); ($y020,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/x021=/,$srqln,2); ($x021,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/y021=/,$srqln,2); ($y021,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/x022=/,$srqln,2); ($x022,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/y022=/,$srqln,2); ($y022,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/x023=/,$srqln,2); ($x023,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/y023=/,$srqln,2); ($y023,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/x024=/,$srqln,2); ($x024,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/y024=/,$srqln,2); ($y024,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/x025=/,$srqln,2); ($x025,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/y025=/,$srqln,2); ($y025,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/x026=/,$srqln,2); ($x026,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/y026=/,$srqln,2); ($y026,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/x027=/,$srqln,2); ($x027,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/y027=/,$srqln,2); ($y027,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/x028=/,$srqln,2); ($x028,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/y028=/,$srqln,2); ($y028,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/x029=/,$srqln,2); ($x029,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/y029=/,$srqln,2); ($y029,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/x030=/,$srqln,2); ($x030,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/y030=/,$srqln,2); ($y030,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/x031=/,$srqln,2); ($x031,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/y031=/,$srqln,2); ($y031,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/x032=/,$srqln,2); ($x032,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/y032=/,$srqln,2); ($y032,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/x033=/,$srqln,2); ($x033,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/y033=/,$srqln,2); ($y033,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/x034=/,$srqln,2); ($x034,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/y034=/,$srqln,2); ($y034,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/x035=/,$srqln,2); ($x035,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/y035=/,$srqln,2); ($y035,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/x036=/,$srqln,2); ($x036,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/y036=/,$srqln,2); ($y036,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/x037=/,$srqln,2); ($x037,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/y037=/,$srqln,2); ($y037,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/x038=/,$srqln,2); ($x038,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/y038=/,$srqln,2); ($y038,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/x039=/,$srqln,2); ($x039,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/y039=/,$srqln,2); ($y039,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/x040=/,$srqln,2); ($x040,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/y040=/,$srqln,2); ($y040,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/x041=/,$srqln,2); ($x041,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/y041=/,$srqln,2); ($y041,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/x042=/,$srqln,2); ($x042,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/y042=/,$srqln,2); ($y042,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/x043=/,$srqln,2); ($x043,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/y043=/,$srqln,2); ($y043,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/x044=/,$srqln,2); ($x044,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/y044=/,$srqln,2); ($y044,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/x045=/,$srqln,2); ($x045,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/y045=/,$srqln,2); ($y045,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/x046=/,$srqln,2); ($x046,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/y046=/,$srqln,2); ($y046,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/x047=/,$srqln,2); ($x047,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/y047=/,$srqln,2); ($y047,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/x048=/,$srqln,2); ($x048,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/y048=/,$srqln,2); ($y048,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/x049=/,$srqln,2); ($x049,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/y049=/,$srqln,2); ($y049,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/x050=/,$srqln,2); ($x050,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/y050=/,$srqln,2); ($y050,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/x051=/,$srqln,2); ($x051,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/y051=/,$srqln,2); ($y051,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/x052=/,$srqln,2); ($x052,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/y052=/,$srqln,2); ($y052,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/x053=/,$srqln,2); ($x053,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/y053=/,$srqln,2); ($y053,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/x054=/,$srqln,2); ($x054,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/y054=/,$srqln,2); ($y054,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/x055=/,$srqln,2); ($x055,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/y055=/,$srqln,2); ($y055,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/x056=/,$srqln,2); ($x056,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/y056=/,$srqln,2); ($y056,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/x057=/,$srqln,2); ($x057,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/y057=/,$srqln,2); ($y057,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/x058=/,$srqln,2); ($x058,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/y058=/,$srqln,2); ($y058,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/x059=/,$srqln,2); ($x059,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/y059=/,$srqln,2); ($y059,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/x060=/,$srqln,2); ($x060,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/y060=/,$srqln,2); ($y060,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/x061=/,$srqln,2); ($x061,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/y061=/,$srqln,2); ($y061,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/x062=/,$srqln,2); ($x062,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/y062=/,$srqln,2); ($y062,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/x063=/,$srqln,2); ($x063,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/y063=/,$srqln,2); ($y063,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/x064=/,$srqln,2); ($x064,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/y064=/,$srqln,2); ($y064,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/x065=/,$srqln,2); ($x065,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/y065=/,$srqln,2); ($y065,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/x066=/,$srqln,2); ($x066,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/y066=/,$srqln,2); ($y066,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/x067=/,$srqln,2); ($x067,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/y067=/,$srqln,2); ($y067,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/x068=/,$srqln,2); ($x068,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/y068=/,$srqln,2); ($y068,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/x069=/,$srqln,2); ($x069,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/y069=/,$srqln,2); ($y069,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/x070=/,$srqln,2); ($x070,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/y070=/,$srqln,2); ($y070,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/x071=/,$srqln,2); ($x071,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/y071=/,$srqln,2); ($y071,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/x072=/,$srqln,2); ($x072,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/y072=/,$srqln,2); ($y072,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/x073=/,$srqln,2); ($x073,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/y073=/,$srqln,2); ($y073,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/x074=/,$srqln,2); ($x074,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/y074=/,$srqln,2); ($y074,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/x075=/,$srqln,2); ($x075,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/y075=/,$srqln,2); ($y075,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/x076=/,$srqln,2); ($x076,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/y076=/,$srqln,2); ($y076,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/x077=/,$srqln,2); ($x077,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/y077=/,$srqln,2); ($y077,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/x078=/,$srqln,2); ($x078,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/y078=/,$srqln,2); ($y078,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/x079=/,$srqln,2); ($x079,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/y079=/,$srqln,2); ($y079,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/x080=/,$srqln,2); ($x080,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/y080=/,$srqln,2); ($y080,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/x081=/,$srqln,2); ($x081,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/y081=/,$srqln,2); ($y081,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/x082=/,$srqln,2); ($x082,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/y082=/,$srqln,2); ($y082,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/x083=/,$srqln,2); ($x083,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/y083=/,$srqln,2); ($y083,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/x084=/,$srqln,2); ($x084,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/y084=/,$srqln,2); ($y084,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/x085=/,$srqln,2); ($x085,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/y085=/,$srqln,2); ($y085,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/x086=/,$srqln,2); ($x086,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/y086=/,$srqln,2); ($y086,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/x087=/,$srqln,2); ($x087,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/y087=/,$srqln,2); ($y087,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/x088=/,$srqln,2); ($x088,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/y088=/,$srqln,2); ($y088,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/x089=/,$srqln,2); ($x089,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/y089=/,$srqln,2); ($y089,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/x090=/,$srqln,2); ($x090,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/y090=/,$srqln,2); ($y090,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/x091=/,$srqln,2); ($x091,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/y091=/,$srqln,2); ($y091,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/x092=/,$srqln,2); ($x092,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/y092=/,$srqln,2); ($y092,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/x093=/,$srqln,2); ($x093,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/y093=/,$srqln,2); ($y093,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/x094=/,$srqln,2); ($x094,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/y094=/,$srqln,2); ($y094,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/x095=/,$srqln,2); ($x095,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/y095=/,$srqln,2); ($y095,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/x096=/,$srqln,2); ($x096,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/y096=/,$srqln,2); ($y096,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/x097=/,$srqln,2); ($x097,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/y097=/,$srqln,2); ($y097,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/x098=/,$srqln,2); ($x098,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/y098=/,$srqln,2); ($y098,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/x099=/,$srqln,2); ($x099,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/y099=/,$srqln,2); ($y099,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/x100=/,$srqln,2); ($x100,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/y100=/,$srqln,2); ($y100,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/x101=/,$srqln,2); ($x101,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/y101=/,$srqln,2); ($y101,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/x102=/,$srqln,2); ($x102,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/y102=/,$srqln,2); ($y102,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/x103=/,$srqln,2); ($x103,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/y103=/,$srqln,2); ($y103,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/x104=/,$srqln,2); ($x104,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/y104=/,$srqln,2); ($y104,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/x105=/,$srqln,2); ($x105,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/y105=/,$srqln,2); ($y105,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/x106=/,$srqln,2); ($x106,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/y106=/,$srqln,2); ($y106,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/x107=/,$srqln,2); ($x107,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/y107=/,$srqln,2); ($y107,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/x108=/,$srqln,2); ($x108,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/y108=/,$srqln,2); ($y108,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/x109=/,$srqln,2); ($x109,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/y109=/,$srqln,2); ($y109,$scr)=split(/\&/,$yhi,2);
($scr,$xhi)=split(/x110=/,$srqln,2); ($x110,$scr)=split(/\&/,$xhi,2);
($scr,$yhi)=split(/y110=/,$srqln,2); ($y110,$scr)=split(/\&/,$yhi,2);
$iii=0; while($iii<199){$iii++; $p[$iii]=0;};
$iii=-1; while($iii<6){$iii++; $jjj=-1; while($jjj<6){$jjj++;
$ar[$iii][$jjj]=0;};};
# $aaaa="13";
# $bbbb="26";
# $cccc="7";
# $dddd="57";
# $aaaa="9";
# $bbbb="4";
# $cccc="11";
# $dddd="47";
# $x001="1"; $y001="1.1";
# $x002="2"; $y002="4.1";
# $x003="3"; $y003="6.1";
# $x004="4"; $y004="7.1";
# $x005="5"; $y005="9.1";
if(($x000 ne "")&&($y000 ne "")){$p[0]=1;};
if(($x001 ne "")&&($y001 ne "")){$p[1]=1;};
if(($x002 ne "")&&($y002 ne "")){$p[2]=1;};
if(($x003 ne "")&&($y003 ne "")){$p[3]=1;};
if(($x004 ne "")&&($y004 ne "")){$p[4]=1;};
if(($x005 ne "")&&($y005 ne "")){$p[5]=1;};
if(($x006 ne "")&&($y006 ne "")){$p[6]=1;};
if(($x007 ne "")&&($y007 ne "")){$p[7]=1;};
if(($x008 ne "")&&($y008 ne "")){$p[8]=1;};
if(($x009 ne "")&&($y009 ne "")){$p[9]=1;};
if(($x010 ne "")&&($y010 ne "")){$p[10]=1;};
if(($x011 ne "")&&($y011 ne "")){$p[11]=1;};
if(($x012 ne "")&&($y012 ne "")){$p[12]=1;};
if(($x013 ne "")&&($y013 ne "")){$p[13]=1;};
if(($x014 ne "")&&($y014 ne "")){$p[14]=1;};
if(($x015 ne "")&&($y015 ne "")){$p[15]=1;};
if(($x016 ne "")&&($y016 ne "")){$p[16]=1;};
if(($x017 ne "")&&($y017 ne "")){$p[17]=1;};
if(($x018 ne "")&&($y018 ne "")){$p[18]=1;};
if(($x019 ne "")&&($y019 ne "")){$p[19]=1;};
if(($x020 ne "")&&($y020 ne "")){$p[20]=1;};
if(($x021 ne "")&&($y021 ne "")){$p[21]=1;};
if(($x022 ne "")&&($y022 ne "")){$p[22]=1;};
if(($x023 ne "")&&($y023 ne "")){$p[23]=1;};
if(($x024 ne "")&&($y024 ne "")){$p[24]=1;};
if(($x025 ne "")&&($y025 ne "")){$p[25]=1;};
if(($x026 ne "")&&($y026 ne "")){$p[26]=1;};
if(($x027 ne "")&&($y027 ne "")){$p[27]=1;};
if(($x028 ne "")&&($y028 ne "")){$p[28]=1;};
if(($x029 ne "")&&($y029 ne "")){$p[29]=1;};
if(($x030 ne "")&&($y030 ne "")){$p[30]=1;};
if(($x031 ne "")&&($y031 ne "")){$p[31]=1;};
if(($x032 ne "")&&($y032 ne "")){$p[32]=1;};
if(($x033 ne "")&&($y033 ne "")){$p[33]=1;};
if(($x034 ne "")&&($y034 ne "")){$p[34]=1;};
if(($x035 ne "")&&($y035 ne "")){$p[35]=1;};
if(($x036 ne "")&&($y036 ne "")){$p[36]=1;};
if(($x037 ne "")&&($y037 ne "")){$p[37]=1;};
if(($x038 ne "")&&($y038 ne "")){$p[38]=1;};
if(($x039 ne "")&&($y039 ne "")){$p[39]=1;};
if(($x040 ne "")&&($y040 ne "")){$p[40]=1;};
if(($x041 ne "")&&($y041 ne "")){$p[41]=1;};
if(($x042 ne "")&&($y042 ne "")){$p[42]=1;};
if(($x043 ne "")&&($y043 ne "")){$p[43]=1;};
if(($x044 ne "")&&($y044 ne "")){$p[44]=1;};
if(($x045 ne "")&&($y045 ne "")){$p[45]=1;};
if(($x046 ne "")&&($y046 ne "")){$p[46]=1;};
if(($x047 ne "")&&($y047 ne "")){$p[47]=1;};
if(($x048 ne "")&&($y048 ne "")){$p[48]=1;};
if(($x049 ne "")&&($y049 ne "")){$p[49]=1;};
if(($x050 ne "")&&($y050 ne "")){$p[50]=1;};
if(($x051 ne "")&&($y051 ne "")){$p[51]=1;};
if(($x052 ne "")&&($y052 ne "")){$p[52]=1;};
if(($x053 ne "")&&($y053 ne "")){$p[53]=1;};
if(($x054 ne "")&&($y054 ne "")){$p[54]=1;};
if(($x055 ne "")&&($y055 ne "")){$p[55]=1;};
if(($x056 ne "")&&($y056 ne "")){$p[56]=1;};
if(($x057 ne "")&&($y057 ne "")){$p[57]=1;};
if(($x058 ne "")&&($y058 ne "")){$p[58]=1;};
if(($x059 ne "")&&($y059 ne "")){$p[59]=1;};
if(($x060 ne "")&&($y060 ne "")){$p[60]=1;};
if(($x061 ne "")&&($y061 ne "")){$p[61]=1;};
if(($x062 ne "")&&($y062 ne "")){$p[62]=1;};
if(($x063 ne "")&&($y063 ne "")){$p[63]=1;};
if(($x064 ne "")&&($y064 ne "")){$p[64]=1;};
if(($x065 ne "")&&($y065 ne "")){$p[65]=1;};
if(($x066 ne "")&&($y066 ne "")){$p[66]=1;};
if(($x067 ne "")&&($y067 ne "")){$p[67]=1;};
if(($x068 ne "")&&($y068 ne "")){$p[68]=1;};
if(($x069 ne "")&&($y069 ne "")){$p[69]=1;};
if(($x070 ne "")&&($y070 ne "")){$p[70]=1;};
if(($x071 ne "")&&($y071 ne "")){$p[71]=1;};
if(($x072 ne "")&&($y072 ne "")){$p[72]=1;};
if(($x073 ne "")&&($y073 ne "")){$p[73]=1;};
if(($x074 ne "")&&($y074 ne "")){$p[74]=1;};
if(($x075 ne "")&&($y075 ne "")){$p[75]=1;};
if(($x076 ne "")&&($y076 ne "")){$p[76]=1;};
if(($x077 ne "")&&($y077 ne "")){$p[77]=1;};
if(($x078 ne "")&&($y078 ne "")){$p[78]=1;};
if(($x079 ne "")&&($y079 ne "")){$p[79]=1;};
if(($x080 ne "")&&($y080 ne "")){$p[80]=1;};
if(($x081 ne "")&&($y081 ne "")){$p[81]=1;};
if(($x082 ne "")&&($y082 ne "")){$p[82]=1;};
if(($x083 ne "")&&($y083 ne "")){$p[83]=1;};
if(($x084 ne "")&&($y084 ne "")){$p[84]=1;};
if(($x085 ne "")&&($y085 ne "")){$p[85]=1;};
if(($x086 ne "")&&($y086 ne "")){$p[86]=1;};
if(($x087 ne "")&&($y087 ne "")){$p[87]=1;};
if(($x088 ne "")&&($y088 ne "")){$p[88]=1;};
if(($x089 ne "")&&($y089 ne "")){$p[89]=1;};
if(($x090 ne "")&&($y090 ne "")){$p[90]=1;};
if(($x091 ne "")&&($y091 ne "")){$p[91]=1;};
if(($x092 ne "")&&($y092 ne "")){$p[92]=1;};
if(($x093 ne "")&&($y093 ne "")){$p[93]=1;};
if(($x094 ne "")&&($y094 ne "")){$p[94]=1;};
if(($x095 ne "")&&($y095 ne "")){$p[95]=1;};
if(($x096 ne "")&&($y096 ne "")){$p[96]=1;};
if(($x097 ne "")&&($y097 ne "")){$p[97]=1;};
if(($x098 ne "")&&($y098 ne "")){$p[98]=1;};
if(($x099 ne "")&&($y099 ne "")){$p[99]=1;};
if(($x100 ne "")&&($y100 ne "")){$p[100]=1;};
if(($x101 ne "")&&($y101 ne "")){$p[101]=1;};
if(($x102 ne "")&&($y102 ne "")){$p[102]=1;};
if(($x103 ne "")&&($y103 ne "")){$p[103]=1;};
if(($x104 ne "")&&($y104 ne "")){$p[104]=1;};
if(($x105 ne "")&&($y105 ne "")){$p[105]=1;};
if(($x106 ne "")&&($y106 ne "")){$p[106]=1;};
if(($x107 ne "")&&($y107 ne "")){$p[107]=1;};
if(($x108 ne "")&&($y108 ne "")){$p[108]=1;};
if(($x109 ne "")&&($y109 ne "")){$p[109]=1;};
if(($x110 ne "")&&($y110 ne "")){$p[110]=1;};
$aobs=$aaaa-0; $bobs=$bbbb-0; $cobs=$cccc-0; $dobs=$dddd-0;
$x[0]=$x000-0; $y[0]=$y000-0;
$x[1]=$x001-0; $y[1]=$y001-0; $x[2]=$x002-0; $y[2]=$y002-0;
$x[3]=$x003-0; $y[3]=$y003-0; $x[4]=$x004-0; $y[4]=$y004-0;
$x[5]=$x005-0; $y[5]=$y005-0; $x[6]=$x006-0; $y[6]=$y006-0;
$x[7]=$x007-0; $y[7]=$y007-0; $x[8]=$x008-0; $y[8]=$y008-0;
$x[9]=$x009-0; $y[9]=$y009-0; $x[10]=$x010-0; $y[10]=$y010-0;
$x[11]=$x011-0; $y[11]=$y011-0; $x[12]=$x012-0; $y[12]=$y012-0;
$x[13]=$x013-0; $y[13]=$y013-0; $x[14]=$x014-0; $y[14]=$y014-0;
$x[15]=$x015-0; $y[15]=$y015-0; $x[16]=$x016-0; $y[16]=$y016-0;
$x[17]=$x017-0; $y[17]=$y017-0; $x[18]=$x018-0; $y[18]=$y018-0;
$x[19]=$x019-0; $y[19]=$y019-0; $x[20]=$x020-0; $y[20]=$y020-0;
$x[21]=$x021-0; $y[21]=$y021-0; $x[22]=$x022-0; $y[22]=$y022-0;
$x[23]=$x023-0; $y[23]=$y023-0; $x[24]=$x024-0; $y[24]=$y024-0;
$x[25]=$x025-0; $y[25]=$y025-0; $x[26]=$x026-0; $y[26]=$y026-0;
$x[27]=$x027-0; $y[27]=$y027-0; $x[28]=$x028-0; $y[28]=$y028-0;
$x[29]=$x029-0; $y[29]=$y029-0; $x[30]=$x030-0; $y[30]=$y030-0;
$x[31]=$x031-0; $y[31]=$y031-0; $x[32]=$x032-0; $y[32]=$y032-0;
$x[33]=$x033-0; $y[33]=$y033-0; $x[34]=$x034-0; $y[34]=$y034-0;
$x[35]=$x035-0; $y[35]=$y035-0; $x[36]=$x036-0; $y[36]=$y036-0;
$x[37]=$x037-0; $y[37]=$y037-0; $x[38]=$x038-0; $y[38]=$y038-0;
$x[39]=$x039-0; $y[39]=$y039-0; $x[40]=$x040-0; $y[40]=$y040-0;
$x[41]=$x041-0; $y[41]=$y041-0; $x[42]=$x042-0; $y[42]=$y042-0;
$x[43]=$x043-0; $y[43]=$y043-0; $x[44]=$x044-0; $y[44]=$y044-0;
$x[45]=$x045-0; $y[45]=$y045-0; $x[46]=$x046-0; $y[46]=$y046-0;
$x[47]=$x047-0; $y[47]=$y047-0; $x[48]=$x048-0; $y[48]=$y048-0;
$x[49]=$x049-0; $y[49]=$y049-0; $x[50]=$x050-0; $y[50]=$y050-0;
$x[51]=$x051-0; $y[51]=$y051-0; $x[52]=$x052-0; $y[52]=$y052-0;
$x[53]=$x053-0; $y[53]=$y053-0; $x[54]=$x054-0; $y[54]=$y054-0;
$x[55]=$x055-0; $y[55]=$y055-0; $x[56]=$x056-0; $y[56]=$y056-0;
$x[57]=$x057-0; $y[57]=$y057-0; $x[58]=$x058-0; $y[58]=$y058-0;
$x[59]=$x059-0; $y[59]=$y059-0; $x[60]=$x060-0; $y[60]=$y060-0;
$x[61]=$x061-0; $y[61]=$y061-0; $x[62]=$x062-0; $y[62]=$y062-0;
$x[63]=$x063-0; $y[63]=$y063-0; $x[64]=$x064-0; $y[64]=$y064-0;
$x[65]=$x065-0; $y[65]=$y065-0; $x[66]=$x066-0; $y[66]=$y066-0;
$x[67]=$x067-0; $y[67]=$y067-0; $x[68]=$x068-0; $y[68]=$y068-0;
$x[69]=$x069-0; $y[69]=$y069-0; $x[70]=$x070-0; $y[70]=$y070-0;
$x[71]=$x071-0; $y[71]=$y071-0; $x[72]=$x072-0; $y[72]=$y072-0;
$x[73]=$x073-0; $y[73]=$y073-0; $x[74]=$x074-0; $y[74]=$y074-0;
$x[75]=$x075-0; $y[75]=$y075-0; $x[76]=$x076-0; $y[76]=$y076-0;
$x[77]=$x077-0; $y[77]=$y077-0; $x[78]=$x078-0; $y[78]=$y078-0;
$x[79]=$x079-0; $y[79]=$y079-0; $x[80]=$x080-0; $y[80]=$y080-0;
$x[81]=$x081-0; $y[81]=$y081-0; $x[82]=$x082-0; $y[82]=$y082-0;
$x[83]=$x083-0; $y[83]=$y083-0; $x[84]=$x084-0; $y[84]=$y084-0;
$x[85]=$x085-0; $y[85]=$y085-0; $x[86]=$x086-0; $y[86]=$y086-0;
$x[87]=$x087-0; $y[87]=$y087-0; $x[88]=$x088-0; $y[88]=$y088-0;
$x[89]=$x089-0; $y[89]=$y089-0; $x[90]=$x090-0; $y[90]=$y090-0;
$x[91]=$x091-0; $y[91]=$y091-0; $x[92]=$x092-0; $y[92]=$y092-0;
$x[93]=$x093-0; $y[93]=$y093-0; $x[94]=$x094-0; $y[94]=$y094-0;
$x[95]=$x095-0; $y[95]=$y095-0; $x[96]=$x096-0; $y[96]=$y096-0;
$x[97]=$x097-0; $y[97]=$y097-0; $x[98]=$x098-0; $y[98]=$y098-0;
$x[99]=$x099-0; $y[99]=$y099-0; $x[100]=$x100-0; $y[100]=$y100-0;
$x[101]=$x101-0; $y[101]=$y101-0; $x[102]=$x102-0; $y[102]=$y102-0;
$x[103]=$x103-0; $y[103]=$y103-0; $x[104]=$x104-0; $y[104]=$y104-0;
$x[105]=$x105-0; $y[105]=$y105-0; $x[106]=$x106-0; $y[106]=$y106-0;
$x[107]=$x107-0; $y[107]=$y107-0; $x[108]=$x108-0; $y[108]=$y108-0;
$x[109]=$x109-0; $y[109]=$y109-0; $x[110]=$x110-0; $y[110]=$y110-0;
###
### REORDER CONTINGENCY TABLE.
$afnl=$aobs; $bfnl=$bobs; $cfnl=$cobs; $dfnl=$dobs;
$ao=$afnl; $bo=$bfnl; $co=$cfnl; $do=$dfnl;
$ncel=0; if($ao>0){$ncel++;}; if($bo>0){$ncel++;}; if($co>0){$ncel++;};
if($do>0){$ncel++;}; $vvvv=$ao+$bo; $wwww=$co+$do; $xxxx=$ao+$co;
$yyyy=$bo+$do; $zzzz=$vvvv+$wwww;
if($ncel>2){$vv=$ao+$bo; $ww=$co+$do; $xx=$ao+$co; $yy=$bo+$do;
if($vv>$ww){$afnl=$co; $bfnl=$do; $cfnl=$ao; $dfnl=$bo;
$ao=$afnl; $bo=$bfnl; $co=$cfnl; $do=$dfnl; $xxx=$vv; $vv=$ww; $ww=$xxx;
$vv=$ao+$bo; $ww=$co+$do; $xx=$ao+$co; $yy=$bo+$do;};
if($xx>$yy){$afnl=$bo; $bfnl=$ao; $cfnl=$do; $dfnl=$co;
$ao=$afnl; $bo=$bfnl; $co=$cfnl; $do=$dfnl; $xxx=$xx; $xx=$yy; $yy=$xxx;
$vv=$ao+$bo; $ww=$co+$do; $xx=$ao+$co; $yy=$bo+$do;}; $zz=$vv+$ww;
$ax=($vv*$xx)/$zz; $bx=($vv*$yy)/$zz; $cx=($ww*$xx)/$zz; $dx=($ww*$yy)/$zz;
$aw=int($ax); $bw=int($bx); $cw=int($cx); $dw=int($dx);
print qq|\n
|;
print qq|\n $titl |;
print qq|\n ORIGINAL OBSERVED VALUES |;
print qq|\n\n ENTITY LOW ENTITY HIGH TOTALS |;
print qq|\n -------------------------------------------- |;
print qq|\n DATA LOW : $aaaa : $bbbb : $vvvv |;
print qq|\n --------------------------------------------- |;
print qq|\n DATA HIGH : $cccc : $dddd : $wwww |;
print qq|\n --------------------------------------------- |;
print qq|\n TOTALS : $xxxx : $yyyy : $zzzz |;
print qq|\n
|;
if($ncel>2){print qq|\n
|;
print qq|\n RECONFIGURED OBSERVED VALUES |;
print qq|\n\n ENTITY LOW ENTITY HIGH TOTALS |;
print qq|\n -------------------------------------------- |;
print qq|\n DATA LOW : $ao : $bo : $vv |;
print qq|\n --------------------------------------------- |;
print qq|\n DATA HIGH : $co : $do : $ww |;
print qq|\n --------------------------------------------- |;
print qq|\n TOTALS : $xx : $yy : $zz |;
print qq|\n
|;
print qq|\n
|;
print qq|\n RECONFIGURED EXPECTED VALUES |;
print qq|\n\n ENTITY LOW ENTITY HIGH TOTALS |;
print qq|\n -------------------------------------------- |;
print qq|\n DATA LOW : $aw : $bw : $vv |;
print qq|\n --------------------------------------------- |;
print qq|\n DATA HIGH : $cw : $dw : $ww |;
print qq|\n --------------------------------------------- |;
print qq|\n TOTALS : $xx : $yy : $zz |;
print qq|\n
|;
$sq=($ax-$ao)*($ax-$ao); $aq=$sq/$ax; $bq=$sq/$bx; $cq=$sq/$cx; $dq=$sq/$dx;
$chsq=$aq+$bq+$cq+$dq; $chsqw=int($chsq*1000000)/1000000;
print qq|\n
Chisquare: $chsqw |;};
print qq|\n Chisquare Mini-table, 1 d.f. |;
print qq|\n P: : 0.05 : 0.025 : 0.01 : 0.005 |;
print qq|\n X2: : 3.841 : 5.024 : 6.635 : 7.879 |;
print qq|\n
|;};
###
### CALCULATE SUMS OF SQUARES.
$ntot=0; $sx=0; $sy=0; $sxx=0; $syy=0; $sxy=0;
$iii=0; while($iii<110){$iii++; $xi=$x[$iii]; $yi=$y[$iii];
if($p[$iii]>0){$ntot++; $sx=$sx+$xi; $sy=$sy+$yi; $xint=int($xi);
$yint=int($yi); if($xint>5){$xint=5;}; if($yint>5){$yint=5;};
$ar[$xint][$yint]++;
$sxx=$sxx+($xi*$xi); $syy=$syy+($yi*$yi); $sxy=$sxy+($xi*$yi);};};
if(($ncel<3)&&($ntot<2)){
###
### ERROR MESSAGE FOR INSUFFICIENT DATA.
print qq|\n
THERE ARE INSUFFICIENT DATA FOR ANALYSIS. |;
print qq|\n
JOB TERMINATED. ntot=$ntot ncel=$ncel |;
print qq|\n
Last modified: 9/21/2003, G. William Moore, MD, PhD. |;
print qq|\n
\n\n |; exit;};
###
### CALCULATE VARIANCE, STANDARD DEVIATION, R-SQUARE.
if($ntot>1){$mnx=$sx/$ntot; $mny=$sy/$ntot;
$ssxx=$sxx-($sx*$mnx); $ssyy=$syy-($sy*$mny); $ssxy=$sxy-($sx*$mny);
$vrx=$ssxx/($ntot-1); $vry=$ssyy/($ntot-1); $sdx=sqrt($vrx); $sdy=sqrt($vry);
$sex=sqrt($vrx/$ntot); $sey=sqrt($vry/$ntot);
$nztx=$mnx/$sex; $nzty=$mny/$sey;
$bet1=$ssxy/$ssxx; $bet0=$mny-($bet1*$mnx);
$sse=$ssyy-($bet1*$ssxy); $var=$sse/($ntot-2);
$teed=$var/$ssxx; $tees=sqrt($teed); $tee=$bet1/$tees;
$df2=$ntot-2; $df1=$ntot-1;
$are=$ssxy/sqrt($ssxx*$ssyy); $rsq=$are*$are;
$teer=0; if($rsq>0){$teer=sqrt(($rsq*$df2)/(1-$rsq));};
$mnxw=int($mnx*1000000)/1000000; $mnyw=int($mny*1000000)/1000000;
$sdxw=int($sdx*1000000)/1000000; $sdyw=int($sdy*1000000)/1000000;
$teerw=int($teer*1000000)/1000000; $teew=int($tee*1000000)/1000000;
$arew=int($are*1000000)/1000000; $rsqw=int($rsq*1000000)/1000000;
print qq|\n
Sample size: $ntot |;
print qq| Mean X: $mnxw Mean Y: $mnyw |;
print qq|\n
Stdev X: $sdxw Stdev Y: $sdyw |;
print qq|\n
r $arew r2 $rsqw |;
print qq|\n
Correlation-t: $teerw $df2 d.f. |;
print qq|\n
Paired-t: $teew $df2 d.f. |;
print qq|\n Student-t Two-tail Mini-table, 1 d.f. : |;
print qq|\n P : 0.10 : 0.05 : 0.02 : 0.01 : 0.001 |;
print qq|\n --------------------------------------------------------|;
print qq|\n 1 : 6.314 : 12.706 : 31.821 : 63.657 : 636.619 |;
print qq|\n 5 : 2.015 : 2.571 : 3.365 : 4.032 : 6.859 |;
print qq|\n 10 : 1.812 : 2.228 : 2.764 : 3.169 : 4.587 |;
print qq|\n 20 : 1.725 : 2.086 : 2.528 : 2.845 : 3.850 |;
print qq|\n 40 : 1.684 : 2.021 : 2.423 : 2.704 : 3.551 |;
print qq|\n 120 : 1.658 : 1.980 : 2.358 : 2.617 : 3.373 |;
print qq|\n inf : 1.645 : 1.960 : 2.326 : 2.576 : 3.291
|;
print qq|\n
|;
print qq|\n
Non-zero-t-X: $nztx $df1 d.f. |;
print qq|\n
Non-zero-t-Y: $nzty $df1 d.f. |;
print qq|\n RAW CONTINGENCY TABLE |;
print qq|\n 0 1 2 3 4 5 |;
print qq|\n 0: $ar[0][0] $ar[0][1] $ar[0][2] $ar[0][3] $ar[0][4] $ar[0][5]|;
print qq|\n 1: $ar[1][0] $ar[1][1] $ar[1][2] $ar[1][3] $ar[1][4] $ar[1][5]|;
print qq|\n 2: $ar[2][0] $ar[2][1] $ar[2][2] $ar[2][3] $ar[2][4] $ar[2][5]|;
print qq|\n 3: $ar[3][0] $ar[3][1] $ar[3][2] $ar[3][3] $ar[3][4] $ar[3][5]|;
print qq|\n 4: $ar[4][0] $ar[4][1] $ar[4][2] $ar[4][3] $ar[4][4] $ar[4][5]|;
print qq|\n 5: $ar[5][0] $ar[5][1] $ar[5][2] $ar[5][3] $ar[5][4] $ar[5][5]|;
print qq|\n
|;};
###
### CALCULATE ONE-CELL FISHER EXACT TEST.
if($ncel>2){
$fisher=1; $iii=$bo; while($iii<$vv){$iii++; $fisher=$fisher*$iii;};
$iii=$xx; while($iii<$zz){$iii++; $fisher=$fisher/$iii;};
$iii=$co; while($iii<$ww){$iii++; $fisher=$fisher*$iii;};
$iii=$do; while($iii<$yy){$iii++; $fisher=$fisher*$iii;};
$iii=0; while($iii<$ao){$iii++; $fisher=$fisher/$iii;}; $sumfisher=$fisher;
$fisherw=int($fisher*1000000)/1000000;
$sumfisherw=int($sumfisher*1000000)/1000000;
print qq| FISHER EXACT TEST |;
print qq|\n d-cell : value : total : |;
print qq|\n One-Cell: $do : $fisherw : $sumfisherw : |;
###
### CALCULATE ONE-TAIL FISHER EXACT TEST.
$as=$ao; $bs=$bo; $cs=$co; $ds=$do; $mn=$bo; if($co<$mn){$mn=$co;};
$iwhile=0; while($iwhile<$mn){$iwhile++;
$as=$as+1; $bs=$bs-1; $cs=$cs-1; $ds=$ds+1; $fisher=1;
if($bs>=0){if($cs>=0){
$iii=$bs; while($iii<$vv){$iii++; $fisher=$fisher*$iii;};
$iii=$xx; while($iii<$zz){$iii++; $fisher=$fisher/$iii;};
$iii=$cs; while($iii<$ww){$iii++; $fisher=$fisher*$iii;};
$iii=$ds; while($iii<$yy){$iii++; $fisher=$fisher*$iii;};
$iii=0; while($iii<$as){$iii++; $fisher=$fisher/$iii;};
$sumfisher=$sumfisher+$fisher;
$fisherw=int($fisher*1000000)/1000000;
if($fisher<0.000005){$fisherw="0.000000";};
$sumfisherw=int($sumfisher*1000000)/1000000;
print qq|\n High-Tail: $ds : $fisherw : $sumfisherw : |;};};};
$fisherone=$sumfisherw;
###
### CALCULATE TWO-TAIL FISHER EXACT TEST.
$dfx=$do-int($dx); $dfn=$dfx; if($dfn<0){$dfn=-$dfx;};
$tpd=$do-(2*$dfn); $bta=$ao-(2*$dfn);
if($bta<0){$bta=0;}; if($tpd<0){$tpd=0;}; $btd=$tpd-$bta;
if($btd<0){$btd=0;};
$ds=$btd-1; $cs=$ww-$ds; $bs=$yy-$ds; $as=$vv-$bs;
$iwhile=$btd-1; while($iwhile<$tpd){$iwhile++;
$ds=$iwhile; $cs=$ww-$ds; $bs=$yy-$ds; $as=$vv-$bs;
if($as>-1){$iii=$bs; while($iii<$vv){$iii++; $fisher=$fisher*$iii;};
$iii=$xx; while($iii<$zz){$iii++; $fisher=$fisher/$iii;};
$iii=$cs; while($iii<$ww){$iii++; $fisher=$fisher*$iii;};
$iii=$ds; while($iii<$yy){$iii++; $fisher=$fisher*$iii;};
$iii=0; while($iii<$as){$iii++; $fisher=$fisher/$iii;};
$sumfisher=$sumfisher+$fisher;};
$fisherw=int($fisher*1000000)/1000000;
if($fisher<0.000001){$fisherw="0.000000";};
$sumfisherw=int($sumfisher*1000000)/1000000;
print qq|\n Low-Tail: $ds : $fisherw : $sumfisherw : |;};
print qq|
|;
$fishertwo=$sumfisherw;
print qq|\n
FISHER ONE-TAIL: $fisherone |;
print qq|\n
FISHER TWO-TAIL: $fishertwo |;
###
### UNARY PROBABILITY DISTRIBUTION FOR TOKEN SWAP.
$gtot=$zz; $iii=-1; while($iii<1000){$iii++; $pdsn[$iii]=0;};
$pdsn[$ax]=1; $itera=0; while($itera<$gtot){$itera++;
$iii=0; while($iii<$gtot){$iii++; $ai=$iii; $bi=$vv-$ai; $ci=$xx-$ai;
$di=$yy-$bi; $ah=$iii-1; $bh=$vv-$ah; $ch=$xx-$ah; $dh=$yy-$bh;
$qh=($bh*$ch)/(($bh*$ch)+($ah*$dh)); $aj=$iii+1; $bj=$vv-$aj;
$cj=$xx-$aj; $dj=$yy-$bj; $qj=($bj*$cj)/(($aj*$dj)+($aj*$dj));
$pdsx[$iii]=($qh*$pdsn[$ah])+($qj*$pdsn[$aj]);};
$iii=0; while($iii<$gtot){$iii++;$pdsn[$iii]=$pdsx[$iii];};
if($pdsn[$ao]>0){$itera=$zz;};}; $tknsw=$pdsn[$ao];
$itknswd=int($tknsw*10000000); $tknswd=$itknswd/10000000;
if($itknswd<1){$tknswd="<0.00000001";};
###
### PRINT OUTPUT VALUES: CHISQUARE, TOKENSWAP.
print qq|\n
CHISQUARE: $chsqw, 1 d.f. |;
print qq|\n
TOKENSWAP, p: $tknswd |;};
###
### REPRINT RAW DATA.
if($ntot>2){print qq|\n RAW DATA |;
print qq|\n : n : i : x : y : |;
$iii=0; $jjj=0; while($iii<110){$iii++; $xi=$x[$iii]; $yi=$y[$iii];
if($p[$iii]>0){$jjj++;
print qq|\n : $iii : $jjj : $xi : $yi : |;};};};
###
### END JOB.
print qq|\n
Last modified: 9/21/2003, G. William Moore, MD, PhD. |;
print qq|\n