axjack's blog

axjack is said to be an abbreviation for An eXistent JApanese Cool Klutz.

2018-12-23から1日間の記事一覧

カイ二乗分布をRで図示する

ソース # df個のN(0,1)な確率変数の二乗和 f <- function(df) sum( rnorm(n=df,m=0,sd =1)^2 ) # 繰り返し回数 N <- 10000; # カイ二乗分布をN回計算してヒストグラムを書く showHistOfChiSqrd <- function (df){ cs <- sapply(rep(df,N),f) hist(cs, freq =…

axjack is said to be an abbreviation for An eXistent JApanese Cool Klutz.