はじめに この記事はStan Advent Calendar2018,12月04日のエントリー記事です。 12月01日のアドカレで,スーパーマラドーナが優勝す
Home | about | Works | Notes | Support | YUEP | BLOG |
はじめに この記事はStan Advent Calendar2018,12月04日のエントリー記事です。 12月01日のアドカレで,スーパーマラドーナが優勝す
はじめに この記事はStan Advent Calendar2018,12月01日のエントリー記事です。 またこの記事は,昨年のStan Advent Calendar20
授業履修前に次のコードを実行して準備しておいてください。 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
左右対称で単峰な分布であれば,平均値,中央値,最頻値は一致しますが,歪んでいると三つはズレます。 このことを正規分布とカイ二乗分布で図示しまし
不偏分散(\(N-1\)で割る方)が標本分散よりも分散の推定値として良いことを可視化します。 HUGOMORE42 set.seed(20180704) library(tidyverse) ## ─ Attaching packages ─────────────────
標準正規分布を描くコード。領域を塗りつぶします。 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),
「標本平均の平均が母平均に一致する」,これを説明するために,次のような例を示して教えてます。 100個のデータを作って「これを母集団とする」と
標準正規分布のテイルを塗りつぶす。 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 =