Jecon-jpa

2018-10-21 · 130 words · 1 minute read

心理学用bibtexスタイルファイルをアップデートしました。中の人の助けを借りて，より正解に近くなりました。スタイルファイルだけでは対応でき

山口大学教育学部・集中講義準備のコード

2018-08-07 · 162 words · 1 minute read

授業履修前に次のコードを実行して準備しておいてください。 install.packages('tidyverse') install.packages('rstan') install.packages('bayesplot') install.packages('summarytools') install.packages('gridExtra') install.packages('GGally') install.packages('loo') source('http://riseki.php.xdoma

t分布の描画

2018-07-06 · 359 words · 1 minute read

Notes

R

ある自由度のt分布とその棄却域，p値の塗りつぶし，片側と両側の違いを比較した図です。自由度18のt分布における棄却域と，t(18)=2.79

三つの代表値がずれる

2018-07-06 · 300 words · 1 minute read

Notes

R

左右対称で単峰な分布であれば，平均値，中央値，最頻値は一致しますが，歪んでいると三つはズレます。このことを正規分布とカイ二乗分布で図示しまし

不偏分散と標本分散の比較

2018-07-04 · 451 words · 1 minute read

Notes

R

不偏分散(\(N-1\)で割る方)が標本分散よりも分散の推定値として良いことを可視化します。 set.seed(20180704) library(tidyverse) ## ─ Attaching packages ──────────────────

標準正規分布を描く

2018-06-30 · 80 words · 1 minute read

Notes

R

標準正規分布を描くコード。領域を塗りつぶします。 library(ggplot2) limitRange <- function(fun, min, max) { function(x) { y <- fun(x) y[x < min | x > max] <- NA return(y) } } ggplot(data.frame(x=c(-4, 4)), aes(x=x)) + stat_function(fun = dnorm) + stat_function(fun = limitRange(dnorm, -1, 1), geom="area", fill="blue", alpha=0.2)+ stat_function(fun = limitRange(dnorm, -2, 2), geom="area",

標本平均の平均が母平均に一致することを表す例

2018-06-29 · 775 words · 2 minute read

Notes

R

「標本平均の平均が母平均に一致する」，これを説明するために，次のような例を示して教えてます。 100個のデータを作って「これを母集団とする」と

標準正規分布のテイルを塗りつぶす

2018-06-29 · 103 words · 1 minute read

標準正規分布のテイルを塗りつぶす。標準偏差+2以上のテイル library(ggplot2) library(gridExtra) limitRange <- function(fun, min, max) { function(x) { y <- fun(x) y[x < min | x > max] <- NA return(y) } } # p(x>2) p <- ggplot(data.frame(x=c(-4, 4)), aes(x=x)) p + stat_function(fun = dnorm) + stat_function(fun = limitRange(dnorm,