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

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

Read more

Share

t分布の描画

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

Read more

Share

三つの代表値がずれる

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

Read more

Share

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

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

Read more

Share

標準正規分布を描く

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

Read more

Share

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

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

Read more

Share

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

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

Read more

Share

任意の相関係数の散布図を描く

任意の相関係数をもつ散布図を描くコード。 library(tidyverse) ## ─ Attaching packages ──────────────────────────────────────── tidyverse 1.2.1 ─ ## ✔

Read more

Share

共分散(相関)の説明に使えそうな図

共分散や相関係数を説明するのに使えそうな図のセットです。 準備 # 準備 library(tidyverse) ## ─ Attaching packages ────────────────────────────────

Read more

Share

正規分布の図を描く

標準偏差の違う正規分布を色々描いてみる。 library(ggplot2) library(gridExtra) old = theme_set(theme_gray(base_family = "HiraKakuProN-W3")) p1 <- ggplot(data = data.frame(X = c(0, 100)), aes(x = X)) + stat_function(fun = dnorm, args = list(mean = 90, sd = 8), color = "red") + xlab("") + ylab("英語&

Read more

Share