Найдите одноэкзонные модели, которые являются частью многоэкзонного гена

#r #model #statistics #subset #bioinformatics

#r #Модель #Статистика #подмножество #биоинформатика

Вопрос:

У меня есть наборы моделей, которые включают в себя модели с одним экзоном и модели с несколькими экзонами (.psl). Я должен найти модели с одним экзоном, которые являются частью моделей с несколькими экзонами

 names(a) [1] "qseqid" "a_stitle" "a_qlen" "a_slen" "a_length" "a_qstart" "a_qend" "a_sstart" "a_send" "a_evalue" "a_pident" "a_qcovs" "a_qcovhsp"  

Как написать код для сравнения sstart и отправки моделей?

Ожидаемый результат (фрейм данных):

 Mult_model Sing_models m_qseqid_1 s_qseqid_1, s_qseqid_2 ... m_qseqid_2 s_qseqid_3 m_qseqid_3 s_qseqid_4, s_qseqid_5 ...  

Комментарии:

1. Вопрос не ясен. Каков ожидаемый результат?

2. Извините 🙁 Ожидаемый результат — это кадр со следующими параметрами: 1 — qseqid модели с несколькими экзонами, 2-qseqid модели(ов) с одним экзоном, которые являются частью модели с несколькими экзонами

3. Потому qsqid что не имеет «_»? Попробуй a[, grep("_", a, fixed = TRUE, value = TRUE, invert = TRUE)] .