Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

error in building model for likelihood (LH) #6

Open
LimKaiShi opened this issue Oct 10, 2024 · 0 comments
Open

error in building model for likelihood (LH) #6

LimKaiShi opened this issue Oct 10, 2024 · 0 comments

Comments

@LimKaiShi
Copy link

LimKaiShi commented Oct 10, 2024

Hi encountered this error while building the model for LH, and when I try to use the output model for creating the universe it failed as well

the code I am running:

%%bash

# Extract unique categories from the metadata file
categories=($(awk -F"\t" '$6 != "disease_adult" && $6 != "tumor" && FNR > 1 {print $6}' "$HOME/scratch/retina_files/master_metadata.txt" | sort | uniq))

for categ in "${categories[@]}"; do
    echo "Processing category: $categ"

    file_no=$(wc -l < "$HOME/scratch/retina_files/analysis/universe/file_list_${categ}.txt")

    geniml lh build_model \
        --model-file "$HOME/scratch/retina_files/analysis/universe/coverage_${categ}/model.tar" \
        --coverage-folder "$HOME/scratch/retina_files/analysis/universe/coverage_${categ}/" \
        --file-no $file_no &
        
done
wait

The output:

Command: lh
/home/users/astar/gis/stukaishi/miniforge3/lib/python3.12/site-packages/geniml/likelihood/build_model.py:27: RuntimeWarning: invalid value encountered in log10
  np.log10((file_no - i) / (no_possible - no_cov) + 1e-10),
Function 'main' executed in 6.2146min
Processing category: healthy_fetus
Command: lh
/home/users/astar/gis/stukaishi/miniforge3/lib/python3.12/site-packages/geniml/likelihood/build_model.py:27: RuntimeWarning: invalid value encountered in log10
  np.log10((file_no - i) / (no_possible - no_cov) + 1e-10),
Function 'main' executed in 6.0920min

I have checked the file_no there is no problem and the coverage folder also works well with the other type of universe
Hope you can check it out, thanks for creating this amazing tools!

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant