@@ -1270,14 +1270,14 @@ static void print_stats(args_t *args)
1270
1270
for (id = 0 ; id < args -> nstats ; id ++ )
1271
1271
{
1272
1272
stats_t * stats = & args -> stats [id ];
1273
- printf ("SN\t%d\tnumber of records:\t%" PRId64 "\n" , id , stats -> n_records );
1274
- printf ("SN\t%d\tnumber of no-ALTs:\t%" PRId64 "\n" , id , stats -> n_noalts );
1275
- printf ("SN\t%d\tnumber of SNPs:\t%" PRId64 "\n" , id , stats -> n_snps );
1276
- printf ("SN\t%d\tnumber of MNPs:\t%" PRId64 "\n" , id , stats -> n_mnps );
1277
- printf ("SN\t%d\tnumber of indels:\t%" PRId64 "\n" , id , stats -> n_indels );
1278
- printf ("SN\t%d\tnumber of others:\t%" PRId64 "\n" , id , stats -> n_others );
1279
- printf ("SN\t%d\tnumber of multiallelic sites:\t%" PRId64 "\n" , id , stats -> n_mals );
1280
- printf ("SN\t%d\tnumber of multiallelic SNP sites:\t%" PRId64 "\n" , id , stats -> n_snp_mals );
1273
+ printf ("SN\t%d\tnumber of records:\t%" PRIu64 "\n" , id , stats -> n_records );
1274
+ printf ("SN\t%d\tnumber of no-ALTs:\t%" PRIu64 "\n" , id , stats -> n_noalts );
1275
+ printf ("SN\t%d\tnumber of SNPs:\t%" PRIu64 "\n" , id , stats -> n_snps );
1276
+ printf ("SN\t%d\tnumber of MNPs:\t%" PRIu64 "\n" , id , stats -> n_mnps );
1277
+ printf ("SN\t%d\tnumber of indels:\t%" PRIu64 "\n" , id , stats -> n_indels );
1278
+ printf ("SN\t%d\tnumber of others:\t%" PRIu64 "\n" , id , stats -> n_others );
1279
+ printf ("SN\t%d\tnumber of multiallelic sites:\t%" PRIu64 "\n" , id , stats -> n_mals );
1280
+ printf ("SN\t%d\tnumber of multiallelic SNP sites:\t%" PRIu64 "\n" , id , stats -> n_snp_mals );
1281
1281
}
1282
1282
printf ("# TSTV, transitions/transversions:\n# TSTV\t[2]id\t[3]ts\t[4]tv\t[5]ts/tv\t[6]ts (1st ALT)\t[7]tv (1st ALT)\t[8]ts/tv (1st ALT)\n" );
1283
1283
for (id = 0 ; id < args -> nstats ; id ++ )
@@ -1419,7 +1419,7 @@ static void print_stats(args_t *args)
1419
1419
{
1420
1420
if ( usr -> vals_ts [j ]+ usr -> vals_tv [j ] == 0 ) continue ; // skip empty bins
1421
1421
float val = usr -> min + (usr -> max - usr -> min )* j /(usr -> nbins - 1 );
1422
- const char * fmt = usr -> type == BCF_HT_REAL ? "USR:%s/%d\t%d\t%e\t%" PRId64 "\t%" PRId64 "\t%" PRId64 "\n" : "USR:%s/%d\t%d\t%.0f\t%" PRId64 "\t%" PRId64 "\t%" PRId64 "\n" ;
1422
+ const char * fmt = usr -> type == BCF_HT_REAL ? "USR:%s/%d\t%d\t%e\t%" PRIu64 "\t%" PRIu64 "\t%" PRIu64 "\n" : "USR:%s/%d\t%d\t%.0f\t%" PRIu64 "\t%" PRIu64 "\t%" PRIu64 "\n" ;
1423
1423
printf (fmt ,usr -> tag ,usr -> idx ,id ,val ,usr -> vals_ts [j ]+ usr -> vals_tv [j ],usr -> vals_ts [j ],usr -> vals_tv [j ]);
1424
1424
}
1425
1425
}
@@ -1512,8 +1512,8 @@ static void print_stats(args_t *args)
1512
1512
}
1513
1513
double af = args -> af_bins ? (bin_get_value (args -> af_bins ,i )+ bin_get_value (args -> af_bins ,i - 1 ))* 0.5 : (double )(i - 1 )/(args -> m_af - 1 );
1514
1514
printf ("GC%cAF\t2\t%f" , x == 0 ? 's' : 'i' , af );
1515
- printf ("\t%" PRId64 "\t%" PRId64 "\t%" PRId64 "" , m [T2S (GT_HOM_RR )],m [T2S (GT_HET_RA )],m [T2S (GT_HOM_AA )]);
1516
- printf ("\t%" PRId64 "\t%" PRId64 "\t%" PRId64 "" , mm [T2S (GT_HOM_RR )],mm [T2S (GT_HET_RA )],mm [T2S (GT_HOM_AA )]);
1515
+ printf ("\t%" PRIu64 "\t%" PRIu64 "\t%" PRIu64 "" , m [T2S (GT_HOM_RR )],m [T2S (GT_HET_RA )],m [T2S (GT_HOM_AA )]);
1516
+ printf ("\t%" PRIu64 "\t%" PRIu64 "\t%" PRIu64 "" , mm [T2S (GT_HOM_RR )],mm [T2S (GT_HET_RA )],mm [T2S (GT_HOM_AA )]);
1517
1517
if ( stats [i ].n && !isnan (r2 ) ) printf ("\t%f" , r2 );
1518
1518
else printf ("\t" NA_STRING );
1519
1519
printf ("\t%.0f\n" , stats [i ].n );
@@ -1571,11 +1571,11 @@ static void print_stats(args_t *args)
1571
1571
r2 *= r2 ;
1572
1572
}
1573
1573
printf ("GC%cS\t2\t%s\t%.3f" , x == 0 ? 's' : 'i' , args -> files -> samples [i ], m + mm ? mm * 100.0 /(m + mm ) : 0 );
1574
- printf ("\t%" PRId64 "\t%" PRId64 "\t%" PRId64 "" ,
1574
+ printf ("\t%" PRIu64 "\t%" PRIu64 "\t%" PRIu64 "" ,
1575
1575
stats [i ].gt2gt [T2S (GT_HOM_RR )][T2S (GT_HOM_RR )],
1576
1576
stats [i ].gt2gt [T2S (GT_HET_RA )][T2S (GT_HET_RA )],
1577
1577
stats [i ].gt2gt [T2S (GT_HOM_AA )][T2S (GT_HOM_AA )]);
1578
- printf ("\t%" PRId64 "\t%" PRId64 "\t%" PRId64 "" ,
1578
+ printf ("\t%" PRIu64 "\t%" PRIu64 "\t%" PRIu64 "" ,
1579
1579
stats [i ].gt2gt [T2S (GT_HOM_RR )][T2S (GT_HET_RA )] + stats [i ].gt2gt [T2S (GT_HOM_RR )][T2S (GT_HOM_AA )],
1580
1580
stats [i ].gt2gt [T2S (GT_HET_RA )][T2S (GT_HOM_RR )] + stats [i ].gt2gt [T2S (GT_HET_RA )][T2S (GT_HOM_AA )],
1581
1581
stats [i ].gt2gt [T2S (GT_HOM_AA )][T2S (GT_HOM_RR )] + stats [i ].gt2gt [T2S (GT_HOM_AA )][T2S (GT_HET_RA )]);
@@ -1631,7 +1631,7 @@ static void print_stats(args_t *args)
1631
1631
printf ("GCT%c\t%s" , x == 0 ? 's' : 'i' , args -> files -> samples [i ]);
1632
1632
for (j = 0 ; j < 5 ; j ++ )
1633
1633
for (k = 0 ; k < 5 ; k ++ )
1634
- printf ("\t%" PRId64 , stats [i ].gt2gt [j ][k ]);
1634
+ printf ("\t%" PRIu64 , stats [i ].gt2gt [j ][k ]);
1635
1635
printf ("\n" );
1636
1636
}
1637
1637
}
@@ -1650,8 +1650,8 @@ static void print_stats(args_t *args)
1650
1650
if ( i == 0 ) printf ("<%d" , stats -> dp .min );
1651
1651
else if ( i + 1 == stats -> dp .m_vals ) printf (">%d" , stats -> dp .max );
1652
1652
else printf ("%d" , idist_i2bin (& stats -> dp ,i ));
1653
- printf ("\t%" PRId64 "\t%f" , stats -> dp .vals [i ], sum ? stats -> dp .vals [i ]* 100. /sum : 0 );
1654
- printf ("\t%" PRId64 "\t%f\n" , stats -> dp_sites .vals [i ], sum_sites ? stats -> dp_sites .vals [i ]* 100. /sum_sites : 0 );
1653
+ printf ("\t%" PRIu64 "\t%f" , stats -> dp .vals [i ], sum ? stats -> dp .vals [i ]* 100. /sum : 0 );
1654
+ printf ("\t%" PRIu64 "\t%f\n" , stats -> dp_sites .vals [i ], sum_sites ? stats -> dp_sites .vals [i ]* 100. /sum_sites : 0 );
1655
1655
}
1656
1656
}
1657
1657
0 commit comments