#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)]
.