はじめに この記事はStan Advent Calendar2018，12月04日のエントリー記事です。 12月01日のアドカレで，スーパーマラドーナが優勝す

Read more

はじめに この記事はStan Advent Calendar2018，12月01日のエントリー記事です。 またこの記事は，昨年のStan Advent Calendar20

Read more

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

Read more

授業履修前に次のコードを実行して準備しておいてください。 HUGOMORE42 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.xdom

Read more

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

Read more

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

Read more

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

Read more

標準正規分布を描くコード。領域を塗りつぶします。 HUGOMORE42 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),

Read more

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

Read more

標準正規分布のテイルを塗りつぶす。 HUGOMORE42 標準偏差+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 =

Read more