BeginPackage["psefalse`"]





FalsePostive::usage="FalsePositive[{pF,pS,wF,wS},k]"



Clear[FalsePositive]

  

     

Begin["Private`"] 







FalsePositive[{pF_,pS_,wF_,wS_},k_:9]:=

Module[{est,var,se,lower,upper,res,tab},

 est=pF+k pS;

 var= pF (1-pF)/wF +  k^2 pS(1-pS)/wS;

 se=Sqrt[var]; 

 lower=est-1.96 se;

 upper=est+1.96 se;

 res={pF,pS,est,lower,upper}//N;

 resr=Round[res 1000]/1000//N;

 tab=TableForm[{resr},TableHeadings->{Automatic,{"pF","pS","est","low","up"}}];

 Print @ tab;

Return[Null]] 





End[]

EndPackage[]



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

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

  

                

                

                

                

                

                

                     

