裏 RjpWiki

Julia ときどき R, Python によるコンピュータプログラム,コンピュータ・サイエンス,統計学

ナイチンゲールの鶏頭図--再び

2022年02月23日 | ブログラミング

2019/12/10 にもナイチンゲールの鶏頭図について書いたが,それは R であった。

そのときには,HistData パッケージの Nightingale は知らなかった。

中澤先生が example(Nightingale) で描けると言及していたので,やってみた。

melt がないとか,reshape がないとか,ggplot2 がないとかうるさかった割に,ちょっと残念な出来栄えだった。

ので,Julia で描いたみた。ggplot の呪いがないので,スッキリ書けたと思う。プログラムは下の方に。

2 枚の図は描画範囲を一致させている。ちょっと,色がどぎつい。

using RDatasets, DataFrames
Nightingale = dataset("HistData", "Nightingale");

using Colors, Plots, Dates

function spidernet(df, begin_)
    θ = range(0, 360, length=1200)
    for deg in 250:250:1250
        r = sqrt(deg)
        plot!(r * cosd.(θ), r * sind.(θ), color=:gray70)
        annotate!(r * cosd(150), r * sind(150),
            text(string(deg), 6, :gray70, "times"))
    end
    for i in 1:12
        pos = sqrt(1200)
        x = pos*cosd(begin_[i] - 15)
        y = pos*sind(begin_[i] -15)
        plot!([0, x], [0, y], color=:gray70)
        annotate!(x, y,
            text(Dates.format.(df[i, 1], "yyyy/mm"), 8,
                :gray70, "times", 4<=i<=9 ? :left : :right))
    end
end

function legend_(df, color)
    for (y, col, leg) in zip(-14:-3:-20, color,
                             ["■ Disease", "■ Other", "■ Wounds"])
        annotate!(32, y, text(leg, 8, :left, col, "times"))
    end
end

function nightingale(df)
    begin_ = 180:-30:-210
    # color = [colorant"#b1bbc1", colorant"#e6b4b7", colorant"#6f625e"];
    color = [colorant"#f6716d", colorant"#6298ff", colorant"#1fbf38", ];
    gr(grid=false, showaxis=false, ticks=false, aspect_ratio=1, lw=0.2,
        titlefont = (10, "times"), label="")
    p = plot(title="Diagram of the Causes of Mortality in the Army in the East")
    spidernet(df, begin_)
    legend_(df, color)
    for i in 1:12
        deg = range(begin_[i], begin_[i+1], length = 100)
        for j in sortperm(Vector(df[i, 8:10]), rev = true)
            r = sqrt(df[i, j+7])
            x = vcat(0, r * cosd.(deg))
            y = vcat(0, r * sind.(deg))
            plot!(x, y, seriestype=:shape, fillcolor = color[j], label="")
        end
    end
    p
end

nightingale(Nightingale[1:12, :])
nightingale(Nightingale[13:24, :])

 

コメント (1)
  • X
  • Facebookでシェアする
  • はてなブックマークに追加する
  • LINEでシェアする

PVアクセスランキング にほんブログ村

PVアクセスランキング にほんブログ村