meineko’s blog

元つくばの某独立行政法人勤務の植物屋です。最近は、ほぼ、突発天体の話題です。

遺伝子listを元にデータを抜き出し

みたいなことを、Rでやろうとしました。
全データをalldataに読み込んでおいてから、先頭のalldata$GeneNameというのをラベルにして、

alldata[alldata$GeneName=="AK072482", ]

みたにすると、お目当てのデータが抜き出せました。
気を良くして、遺伝子のlistを読み込んで、抜き出そうとすると

オブジェクトの長さが短いオブジェクトの長さの倍数になっていません

と、怒られます。
alldata[alldata$GeneName==c("AK072482", "AK102606", ]みたいなもので、すでに駄目でした。
なら、1つずつ呼び出そうとしたら、どうも、値でなくて、ポインタで渡すようでうまくいきません。うーむ。


Rのデータの型って、普段は、上手く隠蔽されていて、それが便利なことも多いのですが、こういう時は使いにくいです。
あと、Rを使っていて、普段なら、for loopを回さないいけないようなところを、1行でかけたりして便利なことが多いのですが、でも、これって、"奥が深い症候群”?