2

私はBioconductorを初めて使用し、自分がやりたいことを実行できる適切なパッケージを見つけようとしています。これは、SwissProt IDを取得し、遺伝子記号を取得するか、またはその逆です。

たくさんのパッケージがあり、どれが欲しいのかわかりません。誰かが簡単に答えてくれますか？

4

1 に答える 1

6

あなたが取ることができる1つのアプローチは、生物パッケージを使用することです:

library(org.Hs.eg.db)

私の遺伝子記号がここのキーのようなものだとしましょう:

keys <- c("A1BG","A2M","A2MP1","NAT1","NAT2","AACP")

次に、select() メソッドを使用するだけです (これは R-2.14 以降で機能します)。

select(org.Hs.eg.db, cols=c("SYMBOL", "UNIPROT"), keys= keys, keytype="SYMBOL")

お役に立てれば！

于 2011-11-04T20:29:34.763 に答える