picme є Python пакет, який містить програми для оцінки та побудови філогенетичного інформативність для великих наборів даних.
<Сильний> Установка
На даний момент, найпростіший спосіб встановити програму:
Git клон Git: //github.com/faircloth-lab/picme.git / шлях / до / picme
Щоб запустити тести:
кд / шлях / до / picme /
Тест Python / test_townsend_code.py
<Сильний> Використовувати
Код estimate_p_i.py викликає пакетний файл для Hyphy, який знаходиться в шаблони /. Цей файл повинен бути в тій же позиції по відношенню до скрізь, де ви поклали estimate_p_i.py. Якщо ви встановите рідшає, як зазначено вище, все буде в порядку, на даний момент.
Бігти:
кд / шлях / до / picme /
пітон picme_compute.py Input_Folder_of_Nexus_Files / Input.tree
& NBSP; - вихід Output_Directory
& NBSP; - епохи = 32-42,88-98,95-105,164-174
& NBSP; - раз = 37,93,100,170
& NBSP; - многопроцессорная
--multiprocessing є обов'язковим, без нього, кожен локус буде працювати послідовно.
Якщо ви вже запустили вище і збережені показники в папці виведення (див нижче), ви можете використовувати вже існуючі записи на сайті-ставки, а не оцінки тих знову:
пітон picme_compute.py Input_Folder_of_Site_Rate_JSON_Files / Input.tree
& NBSP; - вихід Output_Directory
& NBSP; - епохи = 32-42,88-98,95-105,164-174
& NBSP; - раз = 37,93,100,170
& NBSP; - багатопроцесорних
& NBSP; - сайт-ставки
Результати
picme пише результати в базі даних SQLite у вихідний каталог за вашим вибором. Цей каталог також зберігає файли швидкості сайту в форматі JSON для кожного локусу, що пройшов через picme_compute.py.
Ви можете отримати доступ до результатів в базі даних таким чином. Додаткові приклади, в тому числі креслення, наведені в документації
- Провернути SQLite:
& NBSP; sqlite3 Output_Directory / філіпченкове-informativeness.sqlite
- Отримати інтегральні дані для всіх епох:
& NBSP; виберіть локус, інтервал, пі від локусів, інтервал, де loci.id = interval.id
- Отримати інтегральні дані для конкретного епохи:
& NBSP; виберіть локус, інтервал, пі від локусів, інтервал
& NBSP; де інтервал = '95 -105 'і loci.id = interval.id;
- Отримати кількість локусів, що мають макс (PI) при різних епох:
& NBSP; створити тимчасовий стіл макс також вибрати ідентифікатор, макс (PI), як максимум з групи інтервалу за ID;
& NBSP; створити тимчасовий таблиці Т як виберіть interval.id, інтервал, макс з інтервалу, макс
& NBSP; де interval.pi = max.max;
& NBSP; виберіть інтервал, COUNT (*) FROM T групи по інтервалу;
<Сильний> Посилаючись picme
При використанні picme, будь ласка, привести:
- Ферклот до н.е., Чанг J, Альфаро Я: picme забезпечує високу пропускну здатність аналіз філогенетичного інформативності.
- Таунсенд JP: Профілізація філіпченкове інформативність. Систематичне Biol. 2007, 56: 222-231.
- Ставок SLK, Мороз ТПВ, Муза SV: Hyphy: гіпотеза Тестування з використанням філогенез. Біоінформатика 2005, 21: 676-679.
Вимоги
- Python
- hyphy2
- NumPy
- SciPy
- DendroPy
Коментар не знайдено