ここから入手できるデータを使用してヒストグラムを作成しようとしています。このデータベースのCSVバージョンを使用して、1年に発見された外植片の数を表示しています。簡単なスクリプトは
bulkdata <- read.csv('file.csv',head=1,sep=',')
pdf(file="yearcount.pdf",family="Times")
bins <- seq(min(bulkdata$discovered,na.rm=T),max(bulkdata$discovered,na.rm=T),by=1)
hist(bulkdata$discovered,breaks=bins,col='gray',ylab="Discovered",xlab="Year",main="",ylim=c(0,100),axes=FALSE)
axis(1, at=seq(1989,2012,by=1))
axis(2, at=seq(0,100,by=10))
grid(nx=10)
hist(bulkdata$discovered,breaks=bins,col='gray',ylab="Discovered",xlab="Year",main="", add=TRUE)
dev.off()
問題は、がのポイントと整列してxaxis
いないことです。によって描かれた線は目盛りと整列していないため、何の意味も持たないため、これは問題です。修正するオプションを追加しようとしましたが、この方法では軸は正しく描画されますが、グリッドは軸の下から始まります。たぶん非標準のパッケージが必要ですか?0
yaxis
grid()
axis(1, at=seq(1989,2012,by=1))
line=-1