PubmedをRで検索するライブラリRISmedを試してしてみます。
こちらのサイトを参考にさせていただきました。
Musings from a PhD candidate
まずはRISmedをインストール。
install.packages("RISmed")
レトロトランスポゾンで検索した結果を、barplot()を使ってグラフを書いてみます。
1970年から2013年の件数を抽出します。
library(RISmed)
tally <- array()
x <- 1
for (i in 1970:2013){
Sys.sleep(1)
r <- EUtilsSummary('retrotransposon', type='esearch', db='pubmed', mindate=i, maxdate=i)
tally[x] <- QueryCount(r)
x <- x + 1
}
names(tally) <- 1970:2013 #列に名前を付ける
max(tally)
barplot(tally, las=2, ylim=c(0,600), main="Number of PubMed articles containing retrotransposon")
実行結果は以下のようになります。
