KangchengHou / admix-kit

Toolkit for analyzing genetics data from admixed populations
https://kangchenghou.github.io/admix-kit
22 stars 5 forks source link

Index of df_indiv does not match the columns of lanc #27

Open michaelofrancis opened 1 month ago

michaelofrancis commented 1 month ago

Describe the bug Hi, I am running admix lanc-convert with pfile and msp, the msp spans the entire chr1 and the pfile is a shard of chr1, not sure if that matters, but I get the following error after running for ~1 hour :

df_snp["POS"][-1] + 1, df_rfmix_info.loc[len(df_rfmix_info) - 1, "epos"] Traceback (most recent call last): File "/group/tools/Anaconda/Anaconda3-2019.10/envs/admix-kit/bin/admix", line 8, in sys.exit(cli()) File "/group/tools/Anaconda/Anaconda3-2019.10/envs/admix-kit/lib/python3.9/site-packages/admix/cli/init.py", line 39, in cli fire.Fire() File "/group/tools/Anaconda/Anaconda3-2019.10/envs/admix-kit/lib/python3.9/site-packages/fire/core.py", line 143, in Fire component_trace = _Fire(component, args, parsed_flag_args, context, name) File "/group/tools/Anaconda/Anaconda3-2019.10/envs/admix-kit/lib/python3.9/site-packages/fire/core.py", line 477, in _Fire component, remaining_args = _CallAndUpdateTrace( File "/group/tools/Anaconda/Anaconda3-2019.10/envs/admix-kit/lib/python3.9/site-packages/fire/core.py", line 693, in _CallAndUpdateTrace component = fn(*varargs, **kwargs) File "/group/tools/Anaconda/Anaconda3-2019.10/envs/admix-kit/lib/python3.9/site-packages/admix/cli/lanc.py", line 135, in lanc_convert lanc = admix.io.read_rfmix( File "/group/tools/Anaconda/Anaconda3-2019.10/envs/admix-kit/lib/python3.9/site-packages/admix/io/_read.py", line 320, in read_rfmix assert np.all(df_indiv.index == lanc.columns) AssertionError

Admix kit function:

admix lanc-convert \ --pfile $pfilestem \ --rfmix $mspchr1 \ --out $out/test-lanc

Thanks, Mike