1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
|
Author: Andreas Tille <tille@debian.org>
Last-Update: 2023-11-28
Description: Suppress warnings to clean up build log
--- a/emmax.c
+++ b/emmax.c
@@ -101,7 +101,7 @@ double tcdf(double t, double nu);
int main(int argc, char** argv) {
int i, j, k, l, n, nf, q0, ngrids, ndigits, istart, iend, nelems, nmiss, *wids;
char *kinf, *phenof, *tpedf, *covf, *outf, *inf, *delims, *lbuf;
- int mphenoflag, write_eig_flag, gz_flag, tped_nheadercols, tfam_nheadercols, zero_miss_flag, isnpid, gen_kin_method, dosage_flag, gls_flag, c;
+ int mphenoflag, write_eig_flag, gz_flag, tped_nheadercols, tfam_nheadercols, zero_miss_flag, isnpid, gen_kin_method, /*dosage_flag,*/ gls_flag, c;
double *phenos, *covs, *kins, *snps;
double sum, llim, ulim, stat, p;
clock_t cstart, cend, sum0, sum1, sum2, clapstart, clapend;
@@ -120,7 +120,7 @@ int main(int argc, char** argv) {
tfam_nheadercols = DEFAULT_TFAM_NUM_HEADER_COLS;
zero_miss_flag = 1;
isnpid = DEFAULT_TPED_SNPID_INDEX;
- dosage_flag = 0;
+ /*dosage_flag = 0;*/
gls_flag = 1;
// Read optional parameters
@@ -188,7 +188,8 @@ int main(int argc, char** argv) {
zero_miss_flag = 0;
break;
case 'O' :
- dosage_flag = 1;
+ /* dosage_flag = 1; */
+ /* This flag is not used */
break;
case 'N' :
gls_flag = 0;
@@ -515,12 +516,12 @@ int main(int argc, char** argv) {
eLvecs = (double*)malloc(sizeof(double)*nf*nf);
remlh = open_file_with_suffix(inf,"reml",0,0);
- fscanf(remlh.fp,"%lf",&REML);
- fscanf(remlh.fp,"%lf",&REML0);
- fscanf(remlh.fp,"%lf",&optdelta);
- fscanf(remlh.fp,"%lf",&optvg);
- fscanf(remlh.fp,"%lf",&optve);
- fscanf(remlh.fp,"%lf",&hg);
+ if ( fscanf(remlh.fp,"%lf",&REML) != 1 ) { emmax_error("Parser error"); }
+ if ( fscanf(remlh.fp,"%lf",&REML0) != 1 ) { emmax_error("Parser error"); }
+ if ( fscanf(remlh.fp,"%lf",&optdelta) != 1 ) { emmax_error("Parser error"); }
+ if ( fscanf(remlh.fp,"%lf",&optvg) != 1 ) { emmax_error("Parser error"); }
+ if ( fscanf(remlh.fp,"%lf",&optve) != 1 ) { emmax_error("Parser error"); }
+ if ( fscanf(remlh.fp,"%lf",&hg) != 1 ) { emmax_error("Parser error"); }
fclose(remlh.fp);
eLvalsh = open_file_with_suffix(inf,"eLvals",0,0);
--- a/emmax-kin.c
+++ b/emmax-kin.c
@@ -59,7 +59,7 @@ void print_help(void) {
int main(int argc, char** argv) {
int i, j, k, n, ac0, ac1, nmiss, nelems, nex, *nexes;
- int verbose, ndigits, tped_nheadercols, tfam_nheadercols, flag_autosomal, rand_fill_flag, ibs_flag, gz_flag, hetero_division_flag;
+ int verbose, ndigits, tped_nheadercols, tfam_nheadercols, flag_autosomal, rand_fill_flag, ibs_flag, hetero_division_flag;
unsigned char *snprow;
char buf[SZBUF];
int c;
@@ -82,7 +82,6 @@ int main(int argc, char** argv) {
rand_fill_flag = 0;
hetero_division_flag = 0;
ibs_flag = 0;
- gz_flag = 0;
verbose = 0;
ndigits = DEFAULT_NDIGITS;
|
