nrm=1;w1=5;fa0=.1;pic=0;shp=1;
[A nnn nnp npn npp]=roc_ratios(nrm,w1,fa0,pic,shp);


hit=npp./(npp+nnp)*100; fa=npn./(nnn+npn)*100; 
ppv=npp./(npp+npn)*100; npv=nnn./(nnn+nnp)*100;
round([hit;fa;ppv;npv])

fa0=.1

%    48    37    52    38    30    64    34    37
%    10    10    10    10    10    10    10    10
%    30    24     6     5    10    56    13     1
%    95    94    99    99    97    93    97   100

fa0=.01;

%    17    12     8    12     6    28     5     0
%     1     1     1     1     1     1     1     1
%    61    51     8    13    18    85    16     0
%    93    93    99    99    97    87    96   100

fa0=.2
%    60    51    70    54    51    77    51    51
%    20    20    20    20    20    20    20    20
%    21    18     4     3     9    43    10     1
%    96    95   100    99    98    94    97   100

round(npv+ppv)
% 117   113   103   103   106   138   107   100 fa0=.2
% 125   118   105   104   107   149   110   101 fa0=.1;
% 154   144   107   112   115   172   112   100 fa0=.01;
% 171   147   107   113    97   180    96   100 fa0=.001;

round(npv.*ppv)
%2039  1678   376   339   843  4102   964    55 fa0=.2
%2850  2256   594   495   970  5208  1261   100 fa0=.1
%5673  4743   792  1287  1746  7395  1536    50 fa0=.01
%2016  1710   400   297   882  4042   970   100 fa0=.001

%hit*(1-fa)
%4806   4083  5634  4311  4047  6124  4115  4098
%4320   3330  4680  3420  2700  5760  3060  3330
%1683   1188   792  1188   594  2772   495     0

fa0=0;



