理学療法士がまとめたノート

統計学備忘録 since2016

Rを使って統計学を勉強するブログです

x軸操作いろいろ

覚えきれないことは素直にかくことにしました

サンプル
標準正規分布N(0,1)
x<-seq(-3,3,0.01)
y<-dnorm(x)
plot(x,y,type="l",ylab ="確率密度",xlab="")

f:id:yoshida931:20170609162253p:plain

 

目盛り操作でN(3,1)に変更する方法
plot(x,y,type="l",ylab ="確率密度",xlab="",xaxt="n")  #xaxt="n"で軸を消去
name<-c(0:6)
axis(side=1,-3:3,labels=name)  #side=1:x軸

f:id:yoshida931:20170609162347p:plain



目盛りと目盛りラベルを消去してテキストに変更する方法

plot(y,type="l")
plot(y,xaxp=c(0,600,4),type="l",xaxt="n")   #xaxt="n"で軸を消去
name<-c("あ","い","う")
axis(side=1,at=c(0,400,600),labels=name)  #指定した場所に挿入f:id:yoshida931:20170609162510p:plain


消す
枠を消す+x軸も消す
plot(y,type="l",bty="n",xaxt="n")  #グラフを重ねるときに使うかも

f:id:yoshida931:20170609162732p:plain

ついでにy軸も消す
plot(y,type="l",bty="n",xaxt="n",yaxt="n")

f:id:yoshida931:20170609162832p:plain

目盛りの向きを変更
plot(y,type="l",las=2)

f:id:yoshida931:20170609162927p:plain

 

目盛りを自由に入れる
(r-de-r様の助言をもとに修正 2017.6.12)
N(300,100^2)のグラフに目盛りを追加してみます
curve(dnorm(x,300,100),0,600)
axis(side=1,at=c(270,350))                 #x軸に270と350を追加しています

f:id:yoshida931:20170612084533p:plain

当分割
plot(x,y,xaxp=c(-3,3,20),type="l")          # axpを使って20分割

f:id:yoshida931:20170609163113p:plain

N(300,100^2)のグラフ
curve(dnorm(x,300,100),0,600,xaxp=c(0,600,8))      #axpを使って8分割

f:id:yoshida931:20170612085145p:plain


以下のサイトがとても参考になりました
バイオインフォ 道場 [bioinfo-Dojo]

R 使い方 軸・ラベルの調整(向き・サイズ・色など) グラフの描き方 – バイオインフォ 道場 [bioinfo-Dojo]