png("hoge.png", width=800, height=800, res=120) bins = 20 h <- hist(data3, plot = FALSE, breaks=bins) x <- seq(max(min(h$breaks), (a[1] - a[2]/a[3] + 0.001)), max(h$breaks), length = 100) # lengthは曲線の目の細かさ y <- gev.dens(a, x) hist…

dir_r_result = sprintf("r_%s", header) dir.create(dir_r_result) library(fs) dir_create(dir_r_result) みたいな感じ2種類あるけど、どちらもディレクトリを作れる dir.create()はすでにディレクトリがあると、エラーで止まるfsというパッケージの方のdi…

あるベクトルに他のベクトルを追記したいとき x <- 1:10="" br="">x = 1:10 y = 11:20 x [1] 1 2 3 4 5 6 7 8 9 10 y [1] 11 12 13 14 15 16 17 18 19 20 x = c(x, y) x [1] 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 for文で回した計算をベクト…

#fitdiftrを使うためのパッケージ library(MASS) # tableでデータ読み込み、1列目のみ抜き出しているが必要かはわからん data1 <- read.table("hoge.txt") data2 <- data1[, 1] # chi square分布でフィットする # chi_kには推定された自由度を代入しておく c…

2つのデータに対してKS検定をする場合 % x=rnorm(10, mean = 1) % y=rnorm(10, mean = 1) % x [1] 0.91430534 1.59262261 1.52565219 1.68065990 2.35755910 1.91386531 [7] -2.26359817 0.08465784 0.06163455 -0.10943752 % y [1] 2.2200235 2.2241460 1.0…

rm(list=ls()) qiita.com

% .libPaths() でわかる自分の場合は /Library/Frameworks/R.framework/Versions/3.5/Resources/libraryだった 自分で追加したい場合は % .libPaths("./hoge") のようにする これを ~/.Rprofile に書いておくと勝手に起動時に読み込まれる (パッケージのイン…

plotの基本 plot(x, y) www.okadajp.org 対数表示、常用対数のみ使える plot(x, y, log="x") plot(x, y, log="y") plot(x, y, log="xy") 範囲指定 plot(x, y, xlim=c(0, 1)) plot(x, y, ylim=c(0, 1)) タイトルとかラベル plot(x, y, main="上につくタイトル…

ここがわかりやすかった といってもヒストグラムを1から書いたわけではなくて、どういうオプションが使えるのかをちょっと知りたかっただけ data-science.gr.jp グラフの中に最尤推定法で得られたパラメーターを書き込みたい Rの数式描画の中で変数を展開し…

(元記事 2019/02/25)Rを使うことになったので新しくカテゴリを作った ヘビーユーザーにはならないと信じてるので、凝ったことはしない - インストール Mac用のRは https://cran.r-project.org の「Download R for (Mac) OS X」からインストールできる追加の…