================================================================================================================================================================================ ./plink2 --vcf /home/peixott/beegfs/TopMed/WithoutHWE/chr22.dose.vcf.gz --max-alleles 2 --rm-dup exclude-all --snps-only --make-pgen --maf 0.005 --out ./PFILES/LARGE_chr22_Try2 PLINK v2.00a3LM AVX2 Intel (11 Oct 2021) www.cog-genomics.org/plink/2.0/ (C) 2005-2021 Shaun Purcell, Christopher Chang GNU General Public License v3 Logging to ./PFILES/LARGE_chr22_Try2.log. Options in effect: --maf 0.005 --make-pgen --max-alleles 2 --out ./PFILES/LARGE_chr22_Try2 --rm-dup exclude-all --snps-only --vcf /home/peixott/beegfs/TopMed/WithoutHWE/chr22.dose.vcf.gz Start time: Wed Jan 12 13:09:03 2022 192026 MiB RAM detected; reserving 96013 MiB for main workspace. Using up to 40 threads (change this with --threads). --vcf: 596489 variants scanned. --vcf: ./PFILES/LARGE_chr22_Try2-temporary.pgen + ./PFILES/LARGE_chr22_Try2-temporary.pvar.zst + ./PFILES/LARGE_chr22_Try2-temporary.psam written. 1498 samples (0 females, 0 males, 1498 ambiguous; 1498 founders) loaded from ./PFILES/LARGE_chr22_Try2-temporary.psam. 553786 out of 596489 variants loaded from ./PFILES/LARGE_chr22_Try2-temporary.pvar.zst. Note: No phenotype data present. Note: Skipping --rm-dup since no duplicate IDs are present. Calculating allele frequencies... done. 413393 variants removed due to allele frequency threshold(s) (--maf/--max-maf/--mac/--max-mac). 140393 variants remaining after main filters. Writing ./PFILES/LARGE_chr22_Try2.psam ... done. Writing ./PFILES/LARGE_chr22_Try2.pvar ... done. Writing ./PFILES/LARGE_chr22_Try2.pgen ... done. End time: Wed Jan 12 13:09:15 2022 ================================================================================================================================================================================ admix assoc --pfile ./PFILES/LARGE_chr22_Biallelic_COVAR --pheno ./PFILES/LARGE_chr22_Biallelic_COVAR_pheno.txt --pheno-col DISEASE --method ATT,TRACTOR,SNP1 --out AssocTest22 --family binary --covar ./PFILES/LARGE_chr22_Biallelic_COVAR_covarNonNA.txt 2022-01-12 13:21.31 [info ] Received parameters: assoc --pfile=./PFILES/LARGE_chr22_Biallelic_COVAR --pheno=./PFILES/LARGE_chr22_Biallelic_COVAR_pheno.txt --pheno_col=DISEASE --out=AssocTest22 --covar=./PFILES/LARGE_chr22_Biallelic_COVAR_covarNonNA.txt --method=('ATT', 'TRACTOR', 'SNP1') --family=binary --fast=True 2022-01-12 13:21.31 [info ] admix.Dataset: read local ancestry from ./PFILES/LARGE_chr22_Biallelic_COVAR.lanc 2022-01-12 13:21.31 [info ] admix.Dataset: `n_anc` is not provided, infered n_anc from the first 1,000 SNPs is 3. If this is not correct, provide `n_anc` when constructing admix.Dataset 2022-01-12 13:21.31 [info ] 140393 SNPs and 1480 individuals are loaded 2022-01-12 13:21.31 [warning ] admix.dataset.append_indiv_info: 1/1481 individuals in the new dataframe not in the dataset; These individuals will be ignored. 2022-01-12 13:21.31 [warning ] admix.dataset.append_indiv_info: 1/1481 individuals in the new dataframe not in the dataset; These individuals will be ignored. 2022-01-12 13:21.31 [info ] 12 covariates are loaded: AGE,SEX,PC1,PC2,PC3,PC4,PC5,PC6,PC7,PC8,PC9,PC10 2022-01-12 13:21.31 [info ] 140393 SNPs and 1480 individuals left after filtering for missing phenotype, or completely missing covariate 2022-01-12 13:21.31 [info ] Detected categorical columns: 2022-01-12 13:21.31 [info ] Added dummy variables: 2022-01-12 13:21.31 [info ] Performing association analysis with method ATT admix.assoc.marginal: 5%|█████████ | 7/138 [00:24<07:30, 3.44s/it]/cm/shared/apps/python/3.8.6/lib/python3.8/site-packages/scipy/stats/_distn_infrastructure.py:903: RuntimeWarning: invalid value encountered in greater return (a < x) & (x < b) /cm/shared/apps/python/3.8.6/lib/python3.8/site-packages/scipy/stats/_distn_infrastructure.py:903: RuntimeWarning: invalid value encountered in less return (a < x) & (x < b) /cm/shared/apps/python/3.8.6/lib/python3.8/site-packages/scipy/stats/_distn_infrastructure.py:1912: RuntimeWarning: invalid value encountered in less_equal cond2 = cond0 & (x <= _a) admix.assoc.marginal: 7%|███████████▌ | 9/138 [00:36<10:56, 5.09s/it]Warning: [tinygwas.logistic_reg] A fraction of 6% of the observations can be perfectly predicted (abs(y - g) < 1e-4). This might indicate complete quasi-separation. Exiting iterations early and the corresponding loglik is set to NaN Warning: [tinygwas.logistic_reg] A fraction of 31% of the observations can be perfectly predicted (abs(y - g) < 1e-4). This might indicate complete quasi-separation. Exiting iterations early and the corresponding loglik is set to NaN Warning: [tinygwas.logistic_reg] A fraction of 31% of the observations can be perfectly predicted (abs(y - g) < 1e-4). This might indicate complete quasi-separation. Exiting iterations early and the corresponding loglik is set to NaN Warning: [tinygwas.logistic_reg] A fraction of 22% of the observations can be perfectly predicted (abs(y - g) < 1e-4). This might indicate complete quasi-separation. Exiting iterations early and the corresponding loglik is set to NaN Warning: [tinygwas.logistic_reg] A fraction of 22% of the observations can be perfectly predicted (abs(y - g) < 1e-4). This might indicate complete quasi-separation. Exiting iterations early and the corresponding loglik is set to NaN Warning: [tinygwas.logistic_reg] A fraction of 22% of the observations can be perfectly predicted (abs(y - g) < 1e-4). This might indicate complete quasi-separation. Exiting iterations early and the corresponding loglik is set to NaN Warning: [tinygwas.logistic_reg] A fraction of 22% of the observations can be perfectly predicted (abs(y - g) < 1e-4). This might indicate complete quasi-separation. Exiting iterations early and the corresponding loglik is set to NaN Warning: [tinygwas.logistic_reg] A fraction of 31% of the observations can be perfectly predicted (abs(y - g) < 1e-4). This might indicate complete quasi-separation. Exiting iterations early and the corresponding loglik is set to NaN Warning: [tinygwas.logistic_reg] A fraction of 22% of the observations can be perfectly predicted (abs(y - g) < 1e-4). This might indicate complete quasi-separation. Exiting iterations early and the corresponding loglik is set to NaN admix.assoc.marginal: 10%|█████████████████▉ | 14/138 [00:53<07:39, 3.71s/it]Traceback (most recent call last): File "/home/peixott/.local/bin/admix", line 33, in sys.exit(load_entry_point('admix-kit', 'console_scripts', 'admix')()) File "/mnt/beegfs/peixott/admix-kit/admix/cli.py", line 542, in cli fire.Fire() File "/home/peixott/.local/lib/python3.8/site-packages/fire/core.py", line 141, in Fire component_trace = _Fire(component, args, parsed_flag_args, context, name) File "/home/peixott/.local/lib/python3.8/site-packages/fire/core.py", line 466, in _Fire component, remaining_args = _CallAndUpdateTrace( File "/home/peixott/.local/lib/python3.8/site-packages/fire/core.py", line 681, in _CallAndUpdateTrace component = fn(*varargs, **kwargs) File "/mnt/beegfs/peixott/admix-kit/admix/cli.py", line 526, in assoc dict_rls[m] = admix.assoc.marginal( File "/mnt/beegfs/peixott/admix-kit/admix/assoc/__init__.py", line 333, in marginal _block_test( File "/mnt/beegfs/peixott/admix-kit/admix/assoc/__init__.py", line 71, in _block_test var = var.compute() File "/home/peixott/.local/lib/python3.8/site-packages/dask/base.py", line 288, in compute (result,) = compute(self, traverse=False, **kwargs) File "/home/peixott/.local/lib/python3.8/site-packages/dask/base.py", line 571, in compute results = schedule(dsk, keys, **kwargs) File "/home/peixott/.local/lib/python3.8/site-packages/dask/threaded.py", line 79, in get results = get_async( File "/home/peixott/.local/lib/python3.8/site-packages/dask/local.py", line 507, in get_async raise_exception(exc, tb) File "/home/peixott/.local/lib/python3.8/site-packages/dask/local.py", line 315, in reraise raise exc File "/home/peixott/.local/lib/python3.8/site-packages/dask/local.py", line 220, in execute_task result = _execute_task(task, data) File "/home/peixott/.local/lib/python3.8/site-packages/dask/core.py", line 119, in _execute_task return func(*(_execute_task(a, cache) for a in args)) File "/home/peixott/.local/lib/python3.8/site-packages/dask/core.py", line 119, in return func(*(_execute_task(a, cache) for a in args)) File "/home/peixott/.local/lib/python3.8/site-packages/dask/core.py", line 119, in _execute_task return func(*(_execute_task(a, cache) for a in args)) File "/home/peixott/.local/lib/python3.8/site-packages/dask/core.py", line 119, in return func(*(_execute_task(a, cache) for a in args)) File "/home/peixott/.local/lib/python3.8/site-packages/dask/core.py", line 119, in _execute_task return func(*(_execute_task(a, cache) for a in args)) File "/home/peixott/.local/lib/python3.8/site-packages/dask/core.py", line 119, in return func(*(_execute_task(a, cache) for a in args)) File "/home/peixott/.local/lib/python3.8/site-packages/dask/core.py", line 119, in _execute_task return func(*(_execute_task(a, cache) for a in args)) File "/home/peixott/.local/lib/python3.8/site-packages/dask/utils.py", line 37, in apply return func(*args, **kwargs) File "/home/peixott/.local/lib/python3.8/site-packages/dapgen/_read.py", line 350, in _read_pgen_chunk pgen.read_alleles_range(snp_start, snp_stop, geno) File "pgenlib.pyx", line 894, in pgenlib.PgenReader.read_alleles_range File "pgenlib.pyx", line 922, in pgenlib.PgenReader.read_alleles_range RuntimeError: variant_idx 14875 read_alleles_range() error 6 admix.assoc.marginal: 10%|█████████████████▉