せっかくなのだからと、今年は授業も受けてみることに。
#って、いっても他学部のですが。
#断っておきますが、理学研究科のじゃないよ(w
授業なんて何年ぶりだかわからないほどなので、1時間半の講義、寝ずに聞けるか心配だったのですが、Rを使った結構ハードな授業だったので、あっと言う間に終わりました。
(追記)
メモ
授業の応用で、アレイの結果から抜き出した遺伝子のリストgene_list1とgene_list2があり、遺伝子名はaccessionsというラベルが付いている時に、共通の遺伝子listを得るには、
data_1 <- read.csv("gene_list1.txt", sep="\t", header=TRUE)
data_2 <- read.csv("gene_list2.txt", sep="\t", header=TRUE)
res <- data_1[is.element(data_1$accessions,intersect(data_1$accessions, data_2$accessions))==T,]
おぉ、できた。
これまでは、sort+uniqでやっていました。
これで、Rでできます。
#来週の授業あたりで、もっと、簡単な方法を教わったりして？？
