基本的なことをメモ
> df <- read.table("aaa.dat")ファイルを読み込む
> df
表示
> df[c(4)]
4列目を表示
> x <- df[,4] > hist(x)
4列目の値のヒストグラムを作成
> hist(x, xlim=c(1500,4000), breaks=100)
ヒストグラムのx軸を1500-4000に、binを100刻みにして表示
> table(x)
度数分布表を作成
> y = table(x)
> source("http://aoki2.si.gunma-u.ac.jp/R/src/fit_normal.R", encoding="euc-jp")http://aoki2.si.gunma-u.ac.jp/R/fit_normal.html
上記URLからソースを取ってくる。
> fit.normal(y, 1800, 100, method="density")
正規分布に当てはめ、度数分布を計算する
> source("aaa.R")上記の作業手順をaaa.Rに保存したのちに、sourceで呼び出して同じ操作をする
> y <- hist(x) > y$count
ヒストグラムの各binのカウント値を出す
> echo source\(\"aaa.R\"\) | R --vanilla
シェルからRコマンドをRに送る。今回は画像の保存もやってくれる。
par(new=T)
グラフを重ねて表示
lines(c(X1,X2),c(Y1,Y2))
(X1,Y1)から(X2,Y2)に線を引く
