Julia で統計解析 2022/8/14 にバージョンアップ https://r-de-r.github.io/stats/Julia-stats1.htmlhttps://r-de-r.github.io/stats/Julia-stats2.htmlhttps://r-de-r.github.io/stats/Julia-stats3.htmlhttps://r-de-r.github.io/stats/Julia-stats4.htmlh…

ナイチンゲールの鶏頭図よりは折れ線グラフ...という R のプログラムも，あまり良くない気がして，Julia で書いて（描いて）みた。 R の example(Nightingale) の図 折れ線に重ねるデータポイントが，文字通り「間が抜けている」 Julia で描いた図。自画自賛…

2019/12/10 にもナイチンゲールの鶏頭図について書いたが，それは R であった。 そのときには，HistData パッケージの Nightingale は知らなかった。 中澤先生が example(Nightingale) で描けると言及していたので，やってみた。 melt がないとか，reshape …

"""DataFrames の describe() は，sum(), var() を計算しない。 sum(), var() に対応する関数 describe2() を定義しておく。 指定できるのは :sum, :mean, :var, :std, :min, :q25, :median, :q75, :max, :nmissing である。 すべての統計量を計算するときは…

""" Julia の Clustering.hclust() は StatsPlots パッケージで普通の（垂直方向の）デンドログラムを描く。機能を持たない。 以下のプログラムは，R の plot() と同じであるが水平方向のデンドログラムを描く。""" using Plotsfunction plot_hclust_horizon…

using Plotsz = randn(1000, 100);x = sum(z[:, 1:70], dims=2);y = sum(z[:, 30:100], dims=2);scatter(x, y, grid=false, # なくてもいい tick_direction=:out, # なくてもいい label="", # なくてもいい )

using Plotsx = randn(10000)histogram(x, grid=false, # なくてもいい tick_direction=:out, # なくてもいい label="", # なくてもいい )

これらの文書群は github で管理することとした 最新バージョン 2022-02-12 22:38 以下を参照のこと https://r-de-r.github.io/stats/Julia-stats3.html 1. データを使用するための準備 1.1. 既存のデータを使用する 1.2. 自前のデータ2. データフレームの概…

これらの文書群は github で管理することとした 最新バージョン 2022-02-12 22:23 以下を参照のこと https://r-de-r.github.io/stats/Julia-stats7.html 1. HypothesisTests に含まれる検定関数の使用法 1.1. 検定関数関数の呼び出し方 1.2. 検定関数関数に…

うれしいなぁ https://julialang.org/downloads/#current_stable_release

しばらく前からというのも，Julia 1.7.1 になったり，mac OS Monterey になったり，マシンが M1 チップ搭載の Mac mini になったりで，わたしの環境で RCall が動かなくなっていた。 今日，たまたまなんの気無しにやってみたら，動いた。ばんざい！！ julia>…

SymPy で n^10 までと，n^50 をやってみました julia> summation(k^50, (k, 1, n)) |> factor |> print n*(n + 1)*(2*n + 1)*(429*n^48 + 10296*n^47 + 75504*n^46 - 118404*n^45 - 3433716*n^44 + 5209776*n^43 + 177855964*n^42 - 269388834*n^41 - 879055…

思ったより，簡単で，びっくりした。

HypothesisTests の二項検定の p 値が駄目だ。正確に言うと，母比率 ≒ 0.5 のときの両側検定の場合に限るのだが。 母比率 ≒ 0.5のときは，確率分布が歪んでいる。 例えば，x = 18, n = 24, p = 0.68 だと，以下のようになる。 julia> using Plots julia> bar…

Julia の HypothesisTests.PowerDivergenceTest() は問題がある。 対象とする分割表に，度数が 0 のセルがあると不適切な（役に立たない）結果を返す。 julia> x = [4 5 2 0; 0 7 6 1; 1 0 3 1]3×4 Matrix{Int64}: 4 5 2 0 0 7 6 1 1 0 3 1 julia> using Hyp…

SymPy でやった。簡単だった。

Julia の HypothesisTests.ChisqTest() は色々問題がある。 今までもいくつか指摘したところであるが，今回は独立性の検定（χ二乗検定）に引数として 2 つのベクトルを与える場合について書く。 ChisqTest() の定義で，引数に 2 変数のベクトルを与える場合…

それが起こったのは，たぶん 2022 年 1 月末のある日。 以前にはなんのエラーもなく実行できていたプログラムでエラーが生じる。FreqTables.freqtable() で作成したクロス集計表の列や行を抽出できなくなっていた。 再現すると以下のようである。 まず以下の…

これらの文書群は github で管理することとした 最新バージョン 2022-02-02 17:05 以下を参照のこと https://r-de-r.github.io/stats/Julia-stats6.html 1. 数値データの可視化 1.1. ヒストグラム 1.1.1. 一標本の場合 1.3.2. 複数標本の場合 1.4. ボックス…

これらの文書群は github で管理することとした 最新バージョン 2022-02-02 16:54 以下を参照のこと https://r-de-r.github.io/stats/Julia-stats5.html 1. 離散データの可視化 1.1. 例示に使用するデータセット 1.1.1. カテゴリーデータ 1.2. 棒グラフ 1.2.…

これらの文書群は github で管理することとした 最新バージョン 2022-02-02 16:51 以下を参照のこと https://r-de-r.github.io/stats/Julia-stats4.html 1. 一変量統計 1.1. 一変数の場合 1.1.1. 基礎統計量 1.1.2. パーセンタイル値 1.1.3. 度数分布 1.2. …