今週号で連載三回目の「異法人」がちょっと気になる。>モーニング
Ritch Hicky。スライドをちょっと見ただけだけど面白そう
Design, Composition and Performance
ロジスティック曲線を描いてみる。
式はまあ色々あるんだけど、簡単な?ところで。
> logi <- function(x, b, c) 1/(1+(b*exp(1)**(-c*x)))
> logi(0,1,1)
[1] 0.5
> seq(-6,6,0.5)
[1] -6.0 -5.5 -5.0 -4.5 -4.0 -3.5 -3.0 -2.5 -2.0 -1.5 -1.0 -0.5 0.0 0.5
[15] 1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 5.0 5.5 6.0
> logi(seq(-6,6,0.5),1,1)
[1] 0.002472623 0.004070138 0.006692851 0.010986943 0.017986210
[6] 0.029312231 0.047425873 0.075858180 0.119202922 0.182425524
[11] 0.268941421 0.377540669 0.500000000 0.622459331 0.731058579
[16] 0.817574476 0.880797078 0.924141820 0.952574127 0.970687769
[21] 0.982013790 0.989013057 0.993307149 0.995929862 0.997527377
> z <- seq(-6,6,0.1)
> plot(z, logi(z,1,1), ann=F, type="l")
> par(new=T)
> plot(z, logi(z,0.3,1), ann=F, type="l")
> par(new=T)
> plot(z, logi(z,3,1), ann=F, type="l")
>
しまった。同じ色、線種だからどの曲線がどのパラメータで描画されたものかわからんw
上の例ではb の値だけを変化させましたが、c の値を0に近づけるとだんだん直線に近くなります。
たとえば
plot(z, logi(z,1,0.2), ann=F, type="l")
で描画させた結果はこう。
さて話をどう広げよう。