



(*:Mathematica:: V2.1 *)



(*:Context: "psedata`" *)          



(*Title: periodic screeening *)



(*References:   Baker, S.G. and Chu, K.C.*)



(*Keywords:  periodic screening evaluation: incidence data *)



(*Requirements:  matrixlx *)



(*History: 1994 Stuart G. Baker *)









dataOBS={detectOBS,mortOBS,nodeOBS};    

   detectOBS={dataF,dataI,dataS,datanumFIS,datacensFIS,0,0,0,0};

   mortOBS={dF,nF,dI,nI,dS,nS,dR,nR,0,0};

   nodeOBS={xdF,xnF,xdI,xnI,xdS,xnS,xdR,xnR,0,0};







dataRAN={detectRAN,mortRAN,nodeRAN};

   detectRAN={dataF,dataI,dataS,datanumFIS,datacensFIS,

              dataR,dataC,datanumR,datanumC};

   mortRAN={dF,nF,dI,nI,dS,nS,dR,nR,dC,nC};

   nodeRAN={xdF,xnF,xdI,xnI,xdS,xnS,xdR,xnR,xdC,xnC};









(*-----------------------------------------------------------------*)

(*---------NONRANDOMIZED-------------------------------------------*)

(*-----------------------------------------------------------------*)





(*---------age-specific incidence of detection------------*)







(*number of detections by age*)

(*each row is an age 40 to 64*)

(*each column refers to a screen*)

dataF ={0,1,4,1,0,3,1,3,2,1,1,2,3,3,4,3,3,4,3,3,3,4,0,2,0};

dataI={{0,0,0,0},{0,1,0,0},{2,0,1,2},{0,0,1,2},{2,0,1,1},{0,0,1,1},{0,0,2,1},

       {0,1,0,0},{2,0,2,0},{0,1,0,0},{0,0,0,1},{3,0,0,0},{3,0,1,0},{0,2,1,2},

       {1,1,0,1},{0,0,0,1},{1,1,0,3},{0,1,2,0},{1,1,0,0},{0,1,0,0},{0,0,0,0},

       {0,0,0,0},{0,0,0,0},{0,0,0,0},{1,0,0,0}};

dataS={{0,0,0},{0,0,0},{1,0,1},{2,0,0},{1,0,1},{2,0,0},{1,2,3},{2,1,1},

       {2,0,1},{0,0,2},{0,1,2},{1,2,1},{2,0,2},{1,0,1},{3,1,0},{2,1,1},

       {1,1,2},{2,1,0},{3,1,0},{1,1,2},{0,1,2},{1,0,0},{0,3,0},{2,0,0},

       {0,0,0}};

                                                        





        

(*number receiving first screen by age*)

datanumFIS={308,839,1047,1087,1057,1044,981,1018,934,980,912,958,909,880,927,874,876,760,

            718,632,548,501,538,429,376};

     



(*number who drop-out after each screen by age*)

datacensFIS={{29,18,31},{94,55,87},{123,70,107},  {126,68,90},{130,68,100},{131,58,111},

        {132,68,93},{120,55,88},{118,60,85}, {126,72,94},{94,52,88},{105,55,111},

        {103,76,87},{110,60,85},{121,60,112},{106,68,103},{111,57,104},{100,58,69},

        {91,53,71},{89,42,73},{66,31,67},{73,38,52},{84,45,0},{70,0,0},{0,0,0}};    





datanumR={196,396,505,502,470,457,477,484,395,455,493,470,450,440,465,434,432,416,393,

      338,326,314,243,218,213};  





(*------mortality----------*)

 



(*deaths each year after detection*)

  dF={2, 2, 0, 1, 3, 3, 0, 1, 0, 3, 1, 0, 1, 0, 0};

  dI={2, 7, 4, 1, 2, 3, 0, 2, 2, 1, 0, 1, 0, 0, 0};

  dS={1, 3, 0, 1, 3, 2, 1, 3, 1, 3, 0, 2, 1, 1, 0};

  dR={43, 18, 14, 11, 10,  8,  7,  4,  1,  4, 4, 0, 0, 0, 0};



(*at risk at each year after detection*) 

 nF={55, 53, 51, 51, 49, 45, 41, 38, 36, 35, 31, 28, 26, 25, 24};

 nI={48, 46, 39, 34, 33, 30, 26, 25, 23, 20, 19, 19, 18, 18, 16};

 nS={77, 76, 73, 73, 72, 68, 65, 63, 59, 58, 54, 54, 52, 49, 45};

 nR={278, 233, 210, 191, 173, 155, 139, 118, 104, 87, 71, 54, 47, 41, 29};





 



(*-----------node---------------*)

  

(*positive*)

  xdF=13;

  xdI=20;

  xdS=17;

  xdR=22;

  

 (*total*)

  xnF=49;

  xnI=48;

  xnS=74;

  xnR=46;





(*========================================================================*)

(*=======================RANDOMIZED====================================*)

(*========================================================================*)





        

(*number of detections in controls*)

(*rows are age*)

(*columns are intervals after randomization*)



dataC={{1,0,0,1},{2,1,4,1},{2,2,0,5},{2,3,4,4},{2,3,0,4},{2,1,4,3},{2,4,0,1},

       {1,4,2,7},{2,6,3,1},{5,5,1,4},{7,5,2,5},{1,4,3,2},{3,2,3,5},{1,1,2,5},

       {1,3,0,5},{1,1,1,3},{3,6,4,4},{5,3,3,9},{5,1,1,1},{0,4,2,3},{4,2,3,3},

       {3,4,0,2},{3,3,3,1},{2,3,4,1},{4,2,0,1}};





(*number of controls at randomization by age*)



datanumC={552,1362,1530,1578,1496,1550,1455,1423,1432,1428,1413,1416,1392,1350,1389,

        1343,1350,1207,1116,993,909,825,785,666,609};



(*number of refusers at time of randomization by age*)



dataR={{0,0,0,0},{0,0,1,0},{1,0,0,0},{1,1,2,0},{1,0,1,4},{0,2,2,2},{1,1,1,1},

       {1,0,0,0},{0,1,2,1},{0,0,1,0},{1,0,0,0},{0,3,3,0},{0,1,1,1},{1,0,2,1},

       {0,2,1,0},{1,0,1,2},{0,2,0,1},{0,3,1,0},{1,1,0,1},{1,0,1,1},{2,0,1,0},

       {2,0,0,0},{0,2,0,0},{0,0,0,0},{1,0,0,0}};





datanumR={196,396,505,502,470,457,477,484,395,455,493,470,450,440,465,434,432,416,393,

      338,326,314,243,218,213};  

      



(*-------------mortality------------------*)





 (*deaths each year after detection*)

  dC={81, 62, 66, 46, 42, 26, 28, 15, 13, 10, 3, 7, 6, 2, 3};



  (*at risk each year after detection*)

  nC={940, 851, 782, 702, 636, 548, 479, 414, 354, 296, 253, 209, 176, 135, 99};



(*-----------node---------------*)

   

(*positive*)

  xdC=88;

 (*total*)

  xnC=181;



(*----------------false positives------------*)



 parSickles={.014,.006,29694,42451}//N;

 parHIP={.0125,.0051,20211,44056}//N;

