Skip to content

Installing the pipeline and question on halPhyloPTrain.py #1

@SimonaSecomandi

Description

@SimonaSecomandi

Hello,

I'm trying to train a neutral model that I will then use to estimate PhyloP scores on my model species (Hirundo rustica), which is one of the species included in my HAL alignment generated with Cactus 1.2.3.

For the neutral model I would like to reproduce this analysis with your pipeline, but I don't know how to install it! However, I tried using halPhyloPTrain.py on ancestral repeats.

This is the command I'm using:
halPhyloPTrain.py --numProc 32 --no4d --substMod REV ../../9_genomes.hal Anc0 Anc0_REPEATS_NEUTRAL_MODEL.bed Repeats_neutral_model.mod

For this command I'm using the root ancestral genome Anc0 as a reference, while for PhyloP scores I'm planning to use halPhyloPMP.py with my species (Hirundo_rustica) as a reference.

Does it sound good to you? It is correct to use the ancestral genome (from which I extracted the repeats) as a reference for training the neutral model, even if I'm interested in another species?

Many Thanks

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Type

    No type

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions