#r #expression #bioconductor
#r #выражение #биопроводник
Вопрос:
У меня есть ExpressionSet
объект, который я хочу подмножествовать. Например,
> str(ESet)
Formal class 'ExpressionSet' [package "Biobase"] ..
..@ assayData :..
..@ phenoData :
.. .. .. ..$ STATUS : num [1:210] 1 1 1 1 1 1 1 1 1 1 ...
....
Я хочу извлечь подмножество, где STATUS==0
. Я пытался:
exprs(ESet@phenoData$STATUS==0)
но это не работает.
Ответ №1:
Вы почти на месте. Исходя из вашей структуры данных, я думаю, что должно сработать следующее:
exprs(ESet)[ESet@phenoData$STATUS==0,]
Если вы посмотрите на эту статью, то увидите множество примеров подстановки данных в зависимости от нагрузки и времени. Я должен признать, что я автор этой статьи.
Комментарии:
1. @Joris: Меня всегда немного смущает самореклама — в конце концов, я британец, — но когда я проанализировал свой первый набор данных микрочипов, я пожалел, что у меня нет похожей статьи, на которую можно было бы взглянуть.
2. Использование
ESet$STATUS == 0
; смотрите целевую страницу Biobase и раздел «Основы набора выражений» «Введение в Biobase и ExpressionSet» в дополнение к?ExpressionSet
.3. @csgillespie : Не могли бы вы, пожалуйста, обновить ссылку, старая, похоже, мертва