Twitter APIの仕様変更のため、「いいね」の新規取得を終了いたしました

黒木玄 Gen Kuroki

@genkuroki

  • いいね数 389,756/311,170
  • フォロー 995 フォロワー 14,556 ツイート 293,980
  • 現在地 (^-^)/
  • Web https://genkuroki.github.io/documents/
  • 自己紹介 私については https://twilog.org/genkuroki と https://genkuroki.github.io と https://github.com/genkuroki と https://github.com/genkuroki/public を見て下さい。
並び順 : 新→古 | 古→新

2020年04月25日(土)

数学女子 @phasetrbot

20年4月25日

もうJuliaは本を買うか。置く場所ないからPDFがいいのだが、本当に社会は厳しい。

タグ:

posted at 23:42:28

出羽もがみ @DewaMogami

20年4月25日

山形県の偏見地図を作りました pic.twitter.com/CkRpoe3y67

タグ:

posted at 23:37:09

数学女子 @phasetrbot

20年4月25日

Python の数値計算、numpy のベクトル計算(?)がちょっと大変ではあるものの原始的で馬鹿でも書けるという感じがするのだが、Julia は何か謎の @ がすぐ出てきて市民にはとても厳しい代物になっている。そして古めのコードだと動かないようなのでその辺も厳しく、ここでも愚鈍な市民に社会は厳しい。

タグ:

posted at 23:30:04

数学女子 @phasetrbot

20年4月25日

Julia、@ animate とかそういうのが頻発してきて厳しい。

タグ:

posted at 23:21:22

数学女子 @phasetrbot

20年4月25日

optie.hatenablog.com/entry/2018/03/... このページのコードを動かしたいのだが何もわからない。linspaceで怒られて range がいいのではと言われるので変えたがまだ怒られるし何にしたらいいのかわからないのでしらべている。

タグ:

posted at 23:12:09

数学女子 @phasetrbot

20年4月25日

Julia全く分からない。

タグ:

posted at 23:04:02

Matt Henderson @matthen2

20年4月25日

combining two cellular automata- rule 30 inside the central rectangle feeds the game of life outside it pic.twitter.com/vE9HTVKU5O

タグ:

posted at 23:03:26

徳 @trouble_mr

20年4月25日

暴走族はコロナに関係なく走ってるな。
#仙台市

タグ: 仙台市

posted at 22:56:37

KokyuHatuden @breathingpower

20年4月25日

後手後手の感はあるけれど、やらないよりは遥かにまし。他のソーシャルネットワークサービスもニセ科学排除を徹底してほしい。 / “Facebook、ターゲット広告の分類から「疑似科学」を削除 | TechCrunch Japan” htn.to/3YnKVopHbr

タグ:

posted at 22:29:32

黒木玄 Gen Kuroki @genkuroki

20年4月25日

#Julia言語 行列やベクトルを組み合わせて行列を作ることもできます。例えば、

A = [
1 2
3 4
]

b = [
5
6
]

c = [7, 8]'

d = 9

B = [
A b
c d
]

の結果は

B = [
1 2 5
3 4 6
7 8 9
]

と同じ。

タグ: Julia言語

posted at 21:48:58

ねとらぼ @itm_nlab

20年4月25日

200RT 香川ゲーム条例、パブコメ原本(約4000ページ)を公開します nlab.itmedia.co.jp/nl/articles/20... pic.twitter.com/3WsN8kn2hP

タグ:

posted at 21:46:00

黒木玄 Gen Kuroki @genkuroki

20年4月25日

#Julia言語 Float64型成分の行列を

A = Float64[
1 2 3
4 5 6
]

とか

A = Float64[1 2 3; 4 5 6]

で作れます。セミコロンと改行はほぼ同じ意味。

以上でFloat64を略すと整数成分の行列ができる。行列と2次元配列は同じ。縦ベクトルと1次元配列は同じ。

タグ: Julia言語

posted at 21:45:07

TBDD(モデルナ3回接種済) @TBDD_yahoo

20年4月25日

@SciCafeShizuoka しかも、労組が提案しているところもあるとか…

twitter.com/tbdd_yahoo/sta... twitter.com/tbdd_yahoo/sta...

タグ:

posted at 21:44:47

黒木玄 Gen Kuroki @genkuroki

20年4月25日

#Julia言語 縦ベクトルに双対の横ベクトルは

w = [1, 2, 3]'

とか

w = transpose([1, 2, 3])

で作れます。

MATLABやそのクローンに慣れた人は、Juliaでの横ベクトルは縦ベクトルの双対であり、1×n行列ではないことに注意する必要があります。これ重要!

タグ: Julia言語

posted at 21:41:04

ねとらぼ @itm_nlab

20年4月25日

無効分を除いた全ページをスキャンし、PDF化しました

【香川ゲーム条例】パブコメ原本(約4000ページ)を公開します nlab.itmedia.co.jp/nl/articles/20... @itm_nlab pic.twitter.com/ciwmq6TU4T

タグ:

posted at 21:38:34

黒木玄 Gen Kuroki @genkuroki

20年4月25日

#Julia言語 のベクトルや行列の取り扱いのスタイルはMATLABに強く影響を受けています。

v = [1, 2, 3] は整数の縦ベクトル。同じ縦ベクトルを

v = [
1
2
3
]



v = [1; 2; 3]

でも作れます。

タグ: Julia言語

posted at 21:37:42

黒木玄 Gen Kuroki @genkuroki

20年4月25日

#Julia言語

A = [3, 1, 0, 7] はArray{Int, 1}型(整数成分の1次元配列)になります。

A = Float64[3, 1, 0, 7] とすると成分がFloat64型に変換されて、Array{Float64, 1}型になります。

A = Float64[] とすると、Float64型成分の1次元空配列が作成されます。

タグ: Julia言語

posted at 21:33:40

黒木玄 Gen Kuroki @genkuroki

20年4月25日

#Julia言語 Float64型のゼロで埋め尽くされた配列Aの作成は

A = zeros(3, 4) # 3×4の2次元配列

のように作れます。Int8型成分の配列は

B = zeros(Int8, 9, 9)

のようにして作れる。ゼロで埋め尽くす時間がもったいないなら

C = Array{Int8, 2}(undef, 10^3, 10^3)

のように面倒な書き方をする。

タグ: Julia言語

posted at 21:29:25

黒木玄 Gen Kuroki @genkuroki

20年4月25日

#Julia言語 例えば

using Distributions

bin = Binomial(100, 0.4)

とすると二項分布オブジェクト bin ができます。

その乱数や確率函数や累積分布函数は

rand(bin)
rand(bin, n) # サイズnのサンプルの作成
pdf(bin, x)
cdf(bin, x)

などの形式で使用できます。

タグ: Julia言語

posted at 21:23:02

黒木玄 Gen Kuroki @genkuroki

20年4月25日

#Julia言語

型Aのオブジェクトaのメソッドfを

a.f(x)

のスタイルで使うことの、Julia版は

function f(a::A, x)

end

f(a, x)

です。このスタイルは特にDistributions.jlで明瞭なので、そのドキュメントを参考にすると良いと思います。

github.com/JuliaStats/Dis...

タグ: Julia言語

posted at 21:19:29

Noboru Hagino @Noboru_Hagino

20年4月25日

医者は割と「新型コロナ感染症を早期診断できても現時点では救命できる保障なんて全く無い(いまの人類には無理)」とか、色々言ってきたつもりなんですけど、いわゆる「言論人」にはほぼ届いてなかった実感はありますね。
むしろ「政治的な問題」に回収されていました。日本以外のどこかなら可能と。 twitter.com/hazuma/status/...

タグ:

posted at 21:18:18

黒木玄 Gen Kuroki @genkuroki

20年4月25日

#Julia言語 割り算は要注意で、Int分のIntの 6/3 はFloat64型の 2.0 になります。inv(-1)はFloat64型の-1.0になる。

タグ: Julia言語

posted at 21:11:44

ゴルゴ・サーディーン @golgo_sardine

20年4月25日

なるほど。
そういうパターンもありうるわけか。#掛算 twitter.com/if_itself/stat...

タグ: 掛算

posted at 21:09:10

黒木玄 Gen Kuroki @genkuroki

20年4月25日

#Julia言語 数式中では、型が自動的にpromoteして、値がその型に自動変換されるのですが、その辺のことは

docs.julialang.org/en/v1/manual/c...

に書いてあります。例えば、3はInt型で4.0はFloat64型なのですが、3 + 4.0 はFloat64型の7.0になる。

この辺も理解しておかないと、型の伝搬を把握できなくなる。

タグ: Julia言語

posted at 21:07:33

黒木玄 Gen Kuroki @genkuroki

20年4月25日

#Julia言語

大域変数を

g::Float64 = 9.8

のように定義する誤解をたまに見かけます。これは無意味です。

const g = 9.8

と書かないと速度劣化を引き起こします。const で「型を変更できない大域変数」を定義できます(型は変更不可能だが、値は警告付きで変更可能)。

タグ: Julia言語

posted at 21:00:49

舞田敏彦 @tmaita77

20年4月25日

自殺と同様,強盗も失業と相関するのよね。

「コロナで失職して,腹が減った」と,スーパーに押し入った事件が報じられたが。 pic.twitter.com/t8ifgaSh8C

タグ:

posted at 20:56:47

黒木玄 Gen Kuroki @genkuroki

20年4月25日

#Julia言語 函数の定義で引数の型指定は基本的に必要ないことを強調しましたが、structで型を定義するときには、

docs.julialang.org/en/v1/manual/p...

に書いてあることに注意が必要です。

私はお気楽モードで使う場合には

struct Foo{A, B, C}
a::A
b::B
c::C
end

と書いて手抜きをよくする。

タグ: Julia言語

posted at 20:56:25

黒木玄 Gen Kuroki @genkuroki

20年4月25日

#Julia言語

* 函数の定義において、引数の型指定はmultiple dispatchを利用したい場合を除けば必要ない。下手に引数の型指定をすると、函数の適用可能範囲は無用に狭めまる。

* 型の伝搬に気を使う。zero(x)のような書き方をする。型の自動promoteを上手に使うと楽をできる。

タグ: Julia言語

posted at 20:51:15

黒木玄 Gen Kuroki @genkuroki

20年4月25日

#Julia言語 続き。しかし、すでにAがArray{Float64,1}型のとき、

A[1] = 3

とすると、右辺の3は自動的にFloat64の3に変換されて、A[1]はFloat64の3になってくれます。Array{Float64,1}型の配列の成分に整数を代入できたりしません。これと

a = 3

の違いを認識しておかなければいけない。

タグ: Julia言語

posted at 20:47:52

月世界への招待 @MoonWorld_jp

20年4月25日

夕方の西空で 細い月と金星がとても綺麗でした。(1枚目)
明日は更に近づいて見えます。
2020/4/25 月齢2.3 pic.twitter.com/HNEcevfJXM

タグ:

posted at 20:10:53

motoko gram @momogramo

20年4月25日

もう10年程前か、今では親友となった生粋ロンドナーの彼が日本出張し一回りでかくなり帰って来た「飯がうますぎる、うどん&ご飯なんて炭水化物の塊じゃないか、だが逆らえない」それを聞いた夫「お前人の事いえるか」3人で向かったフィッシュ&チップス店でみた衝撃的なチップスバティを未だに忘れない pic.twitter.com/A53tBBfnRa

タグ:

posted at 19:53:34

Erik Engheim @erikengheim

20年4月25日

@cgarciae88 @JuliaComputing Why would that be hard? Put package install instructions into a Julia script and run that. Look at how nextjournal setup their Julia docker containers: nextjournal.com/nextjournal/ju...

Not sure why the choice was made but I greatly prefer it. Much nicer to work with real Julia functions

タグ:

posted at 19:42:40

tea_hashi @tea_hashi

20年4月25日

@aoi_mokei 国勢調査員が漏らす場合もあります。仕事の都合でクッソど田舎に一年だけ住んでいたとき、田舎大嫌いなので、無視してたのですが、隣の家が大家でうるさくて、大ウソ書いて提出しました。
※職業にはAV女優と書きました。

タグ:

posted at 19:35:42

黒木玄 Gen Kuroki @genkuroki

20年4月25日

#Julia言語 そういう注意を要するのは新たに変数を作る場合であり、数式中の計算ではJuliaが適当に型をpromoteしてくれるので結構テキトーに書いても大丈夫。

例えば、xがFloat64型のとき x + 3 の1はFloat64型の3に自動的に変換される。

xがFloat64型であってほしいのに

x = 3

と書くにはアウト。

タグ: Julia言語

posted at 18:52:46

黒木玄 Gen Kuroki @genkuroki

20年4月25日

#Julia言語 例えば、引数の型と同じ型のゼロを返す函数 zero(x) を使うことが結構あります。引数Xが配列でXの成分と同じ型のゼロが欲しければ zero(eltype(X)) を使ったりします。

タグ: Julia言語

posted at 18:49:16

黒木玄 Gen Kuroki @genkuroki

20年4月25日

#Julia言語 こういう仕組みなので、以下のようなコードの書き方が基本になります。

* 函数の引数の型を指定しない。(指定する場合にはその函数の適用可能範囲を無用に狭めないように気を使う。)

* 函数内にコードは引数の型が他の変数の型に適切に伝搬するように書く。

タグ: Julia言語

posted at 18:44:18

黒木玄 Gen Kuroki @genkuroki

20年4月25日

#Julia言語 大事なことなので繰り返し説明しますが、

* 函数は実行時にネイティブコードにコンパイルされる。
* 函数実行時には函数の引数の型は決まっている。

そして

* 函数のコンパイルは引数の型を種にした型推論に基いて実行される。

続く

タグ: Julia言語

posted at 18:40:33

黒木玄 Gen Kuroki @genkuroki

20年4月25日

#Julia言語 では、函数は実行時にコンパイルされます。例えば、函数

f(x) = x^2



f(1)



f(2.4)

の形式で実行するときには、引数xの型はそれぞれIntとFloat64に確定しています。実行時に確定する引数の型を指定することは、その函数の適用可能範囲を狭めることになります。

タグ: Julia言語

posted at 18:36:54

黒木玄 Gen Kuroki @genkuroki

20年4月25日

#Julia言語 では、函数の定義時に引数の型を指定しても計算速度は速くなりません。

multiple dispatchを利用するために函数の引数の型指定を適切にやっているとき以外には、函数の引数の型指定はその函数の適用範囲を狭くすることになります。

部品の汎用性を無駄に下げるコードを書くべきではない。

タグ: Julia言語

posted at 18:33:19

TaKu @takusansu

20年4月25日

@sekibunnteisuu @metameta007 @temmusu_n @KabuTaro @LimgTW 学研(1964年)の「言葉の式」のはしりのようなもの
twitter.com/temmusu_n/stat...
>[1人ぶんの代金]×[人数]がぜんたいの代金だから

1972年には、明確に「言葉の式」を使っているものがあります。
8254.teacup.com/kakezannojunjo...
>(ひとりぶんのあめの数)×(人数)=(全体のあめの数)
>この言葉の式で、

タグ:

posted at 18:32:53

黒木玄 Gen Kuroki @genkuroki

20年4月25日

#Julia言語 の型システムやmultiple dispatchを十分に理解するまでは

* 函数の定義に限っては型指定を行わない

と覚えておけば良いと思います。例えば

function f(x::Float64)
cos(x)*sin(x)
end

function g(X::Array{Float64,1})
sum(exp(-x^2) for x in X)
end

のように書かない。 twitter.com/stattan/status...

タグ: Julia言語

posted at 18:28:44

TaKu @takusansu

20年4月25日

@sekibunnteisuu @metameta007 @temmusu_n @KabuTaro @LimgTW #超算数 「言葉の式」の影響と考えられませんかね。
1977年度
>もとの かず  かける かず  こたえ

1980年度
>1はこぶんの かず  はこの かず  ぜんたいの かず

「かける かず」とかより具体的な状況を「言葉の式」で表したのが、区別を強化する原因の一端と思うのですが。

タグ: 超算数

posted at 18:20:18

非公開

タグ:

posted at xx:xx:xx

非公開

タグ:

posted at xx:xx:xx

甲乙丙丁 @villageofsound

20年4月25日

Julia, めっちゃ気が効くやんけ・・・
Jupyter上でmagic command入れたら「Juliaには似たのがあるから、こっち使ってね」とご紹介を受けた・・・  #Julia #julialang pic.twitter.com/w2uopkiAKH

タグ: Julia julialang

posted at 16:53:28

Massimo @Rainmaker1973

20年4月25日

The ambiguous illusion objects invented by mathematician Kokichi Sugihara now include a Batman logo, printed and sold by 3DNY Design [more: buff.ly/2yHS4yB] pic.twitter.com/s8eifWKCpg

タグ:

posted at 16:38:21

黒木玄 Gen Kuroki @genkuroki

20年4月25日

以上のようなことを、ど素人の私であっても、気楽に書いて実行して、動画を作るところまで行ける点が #Julia言語 の優れた点だと思います。

気楽さと計算速度が両立しているので、微分方程式の数値解を求めるコードを書く気になれる。

ソルバの部分はC++で書けとか言われてもやる気になれない。

タグ: Julia言語

posted at 16:35:55

黒木玄 Gen Kuroki @genkuroki

20年4月25日

#Julia言語 添付動画は自前コードで計算した「自由境界条件の正五角形上の波動方程式」の数値解です。正五角形に限らず、領域の形は自由に設定できるように設計されています。

nbviewer.jupyter.org/gist/genkuroki... pic.twitter.com/Aw76lHnuG0

タグ: Julia言語

posted at 16:31:48

黒木玄 Gen Kuroki @genkuroki

20年4月25日

#Julia言語 「実質編微分方程式の数値解」→「実質的に偏微分方程式の数値解」

twitter.com/genkuroki/stat...

タグ: Julia言語

posted at 16:13:13

黒木玄 Gen Kuroki @genkuroki

20年4月25日

#Julia言語 DifferentialEquations.jlの使用例

以下のリンク先のリンク先

nbviewer.jupyter.org/gist/genkuroki...

ではオープン戸田格子のHamiltonianを与えるだけで、対応する正準方程式の初期値問題の数値解を求めてくれるDifferentialEquations.jlの機能を利用しています。

twitter.com/genkuroki/stat...

タグ: Julia言語

posted at 16:09:37

黒木玄 Gen Kuroki @genkuroki

20年4月25日

#Julia言語 DifferentialEquations.jlの使用例

「跳ね返り」の類もシンプルに記述できます。

これに限らず、DifferentialEquations.jl はものすごく多機能で、私はその全貌を全く把握できていません。めちゃくちゃ、すごいです。

twitter.com/genkuroki/stat...

タグ: Julia言語

posted at 16:03:02

黒木玄 Gen Kuroki @genkuroki

20年4月25日

#Julia言語 DifferentialEquations.jl のミニマルな使用例

twitter.com/genkuroki/stat...

タグ: Julia言語

posted at 16:00:23

黒木玄 Gen Kuroki @genkuroki

20年4月25日

#Julia言語 DifferentialEquationsの使用例

Lorenz方程式の初期値を微小に変えたN個の解の同時プロットのアニメーション。

ぶわぁーってなる(笑)

twitter.com/genkuroki/stat...

タグ: Julia言語

posted at 15:58:08

竹田陽介@病院マーケティングサミットJA @yosuke_takeda

20年4月25日

どうか拡散を。

乳がん術後の放射線治療で免疫力が下がる(≒骨髄抑制)ことはありません。
照射は局所です。
安心して受けてください。
乳癌を治すために大切な治療です。

正しく理解して正しく怖れることが重要です。

放射線治療医の先生からの情報(学会のアナウンス)をシェアします。 pic.twitter.com/YMzojFY49O

タグ:

posted at 15:56:48

黒木玄 Gen Kuroki @genkuroki

20年4月25日

#Julia言語 DifferentialEquations.jl の使用例

N個の質点がバネで繋がっている系の常微分方程式の数値解は実質編微分方程式の差分化の数値解になる。

twitter.com/genkuroki/stat...

タグ: Julia言語

posted at 15:55:41

ぷにチコ @punipunichikori

20年4月25日

高校数学なんて大人になっても役に立たないと思っている学生諸君、キュアマーメイドのスパイラルヒラヒラの型紙を作るためには三角関数、微分積分の知識が必要になるからよく勉強しとけよ pic.twitter.com/PSvu6tEn2e

タグ:

posted at 15:53:10

黒木玄 Gen Kuroki @genkuroki

20年4月25日

#Julia言語 における優れものの巨大パッケージに

DifferentialEquations.jl
github.com/SciML/Differen...

があるのですが、PythonとRのラッパーもあります。

github.com/SciML/diffeqpy

github.com/SciML/diffeqr

例えば、微分方程式をRで扱う必要が生じたら、使用を検討してみるといいかも。

タグ: Julia言語

posted at 15:48:45

のらんぶる @nolimbre

20年4月25日

coinvariantsがcoronavirusに見えたがcontravariantがcoronavirusに見えたに見えた……

タグ:

posted at 15:33:58

Akinori Ito @akinori_ito

20年4月25日

Rから速度が必要な時はRcppをつかってるけどJuliaのほうが書くのが楽だな

タグ:

posted at 15:33:40

Katsushi Kagaya @katzkagaya

20年4月25日

さらに例えば、黒木さんが紹介されていたZerotierを利用してVPNをつくればリモートでrstudio-serverを使うことができます。twitter.com/genkuroki/stat...

タグ:

posted at 15:23:22

TBDD(モデルナ3回接種済) @TBDD_yahoo

20年4月25日

労組が提案、だ、と…

ちょっと何言ってるかわからない

「給付の10万円、県職員は寄付を」労組が提案 神奈川 [新型コロナウイルス]:朝日新聞デジタル www.asahi.com/articles/ASN4S...

タグ:

posted at 15:20:47

Katsushi Kagaya @katzkagaya

20年4月25日

{tidyverse}も{rstan}もいちいちコンパイルが走ってバイナリで入れる方法が分かってません。なので、けっこうインストールに時間がかかります。Windowsだとバイナリインストールされて瞬殺なので少し忍耐がいります。

タグ:

posted at 15:20:05

if @if_itself

20年4月25日

2年生で「かけられる数」「かける数」の概念が身についてしまうと、3年生で割り算する時に混乱しませんか?
うちの子、「かけられる数」「かける数」の両方が「わる数」になり得るということに納得してないんだけど、これは「交換法則も習ったでしょ」でいいのかな。

タグ:

posted at 15:08:55

香川愛生 Manao Kagawa ☗ @MNO_shogi

20年4月25日

もう番長って呼ばれてかれこれ6年半くらい経つので、新しいあだ名募集してみます!

タグ:

posted at 14:56:46

アオイ模型 土-東エ38b @aoi_mokei

20年4月25日

親の話を聞いていると、老人の間でどういう感じでデマや噂が広がっていくか見て取れるのですが、
ダスキンやヤクルトといった昔ながらの宅配販売が「✕✕さんの家にコロナが出たらしい」と噂話を広めていくというのは、中々に盲点だった

地元密着ゆえの厄介さかも知れない

タグ:

posted at 14:51:52

モト@PPMMPP @29silicon

20年4月25日

パチンコ店を強制的に閉めさせはしないが、公務員個人の給付金は半強制的に召し上げるってどうよ?

タグ:

posted at 14:49:51

Katsushi Kagaya @katzkagaya

20年4月25日

Stanを使う場合、WSLの方が速いらしいので。
twitter.com/hankagosa/stat... twitter.com/hankagosa/stat...

タグ:

posted at 14:32:29

左巻健男(サマキタケオ) @samakikaku

20年4月25日

こういう事態はEM菌の衰退を示しており、EM菌側が批判者に非常に攻撃的になっている理由だと思う。
もっとも活動的と思ったNPO法人北海道EM普及協会は、情報誌が92号 (2019年3月発行)以来サイトに出ていない。会員数の減少で年4回から3回に減らしたがそれも出ていないのか?

タグ:

posted at 14:31:49

左巻健男(サマキタケオ) @samakikaku

20年4月25日

EM菌は大丈夫か?
・会社案内 | EM研究機構|※以上 平成28年6月27日 現在と数年間更新なし。
・EM菌推進のEM研究機構顧問だったDND出口俊一氏(ヒノキヤグループ社外取締役)サイトで連載コラムはEM菌比嘉照夫氏だけ。比嘉氏肩書きは何年も前のもの。

タグ:

posted at 14:27:55

黒木玄 Gen Kuroki @genkuroki

20年4月25日

#統計 話題がどんどん脱線するが、「分割表のFisher検定は正確だが、χ²はその近似に過ぎない。サンプルサイズが小さい場合には(補正無しの)χ²検定を使わずに、Fisher検定を使うべきである」というような__不合理__な考え方の普及は、つまらないプレッシャーを現場にかけて、研究活動を阻害したと思う。

タグ: 統計

posted at 14:27:41

黒木玄 Gen Kuroki @genkuroki

20年4月25日

#統計 Rのbinom.test, fisher.test, chisq.testの仕様が十分に合理的でないことは、不合理な伝統の維持の現れの1つだと思う。

「所詮は丼勘定」であることを常に意識していれば実害はないかもしれませんが。

タグ: 統計

posted at 14:22:38

黒木玄 Gen Kuroki @genkuroki

20年4月25日

#統計 今回は結論だけを述べますが、2×2の分割表の独立性検定で「Fisher検定は正確な検定であり、χ²検定はその近似に過ぎない」という考え方は誤り。

多くの解説がおかしい。どうして訂正されないのか?

「統計学入門」という分野は不合理な伝統を維持している不合理な分野になっていると思う。

タグ: 統計

posted at 14:18:44

黒木玄 Gen Kuroki @genkuroki

20年4月25日

#統計 「所詮は丼勘定」なので実践的には細かいことにこだわるのはあまり合理的ではないのですが、Rのbinom.test, fisher.test, chisq.testの仕様が十分に合理的でない点はちょっと気になります。

chisq.testがデフォルトであのろくでもないYatesの連続性補正を行う仕様は特に納得できない。

タグ: 統計

posted at 14:14:09

黒木玄 Gen Kuroki @genkuroki

20年4月25日

#統計 Rのfisher.testでもP値と信頼区間が整合的でなくなることがある件については例えば以下のリンク先を参照。

exact2x2::fisher.exactならば整合的なP値と信頼区間を表示してくれます。

twitter.com/genkuroki/stat...

タグ: 統計

posted at 14:10:59

黒木玄 Gen Kuroki @genkuroki

20年4月25日

#統計 2つ目の問題の答え

Rのbinom.testが出力するp-valueと95%信頼区間は異なるP値函数で計算されているので、それらの間は整合的でなくなる。

P値函数の定義はただ一つではなく、複数あります。

同じP値函数でp-valueと信頼区間を計算する仕様にしていない理由が不明。私は合理的でないと思う。 pic.twitter.com/vuVShBo5bG

タグ: 統計

posted at 14:03:46

黒木玄 Gen Kuroki @genkuroki

20年4月25日

#統計 添付画像は私が書いた #Julia言語 での二項検定のP値函数およびP値函数を与えると信頼区間を計算してくれる函数の実装例です。このようにシンプルに書けます。

P値函数が与えられると95%信頼区間はP値が5%以上になるパラメータpの範囲として計算されます。

nbviewer.jupyter.org/gist/genkuroki... pic.twitter.com/Im7ZWVyq3N

タグ: Julia言語 統計

posted at 13:51:39

黒木玄 Gen Kuroki @genkuroki

20年4月25日

#統計 表示されるP値と信頼区間の間に整合性が無くても、「所詮はどんぶり勘定」なので大して気にする必要はないと思いますが、Rのbinom.testとfisher.testのコードを書いている人達の理解度については若干の不安が生じてしまいます。

次の問題:Rのbinom.testでP値と信頼区間が不整合になる理由は?

タグ: 統計

posted at 13:47:11

Katsushi Kagaya @katzkagaya

20年4月25日

ちょっと前のですが、WSLでRStudio server の手順、僕も書いたのがあるので御参考になれば。Tidyverse入れるとき、ちょっとトリッキーなところがありました。note.com/katzkagaya/n/n...

タグ:

posted at 13:45:33

黒木玄 Gen Kuroki @genkuroki

20年4月25日

#統計 問題の答え

95%信頼区間はP値が5%以上のパラメータpの範囲に一致して欲しいが、P値が5.5%になるパラメータの値p=0.1469が表示された95%信頼区間[0.0165044, 0.1458632]に含まれていない。

実はRのfisher.testでも同様のことが起こります。 pic.twitter.com/Oqad9ZkMiu

タグ: 統計

posted at 13:42:49

黒木玄 Gen Kuroki @genkuroki

20年4月25日

#統計 以下のリンク先ツイート関連の余談

Rのbinom.test(4, 67, p=0.1469)が出力するP値と95%信頼区間の組み合わせに整合性がないという話。

添付画像は #Julia言語 RCall.jl経由で使ったRの出力のスクショ。

問題:どのようにP値と95%信頼区間に整合性がないか? twitter.com/philomyu/statu... pic.twitter.com/H9jdLFxqUd

タグ: Julia言語 統計

posted at 13:42:47

森 勇一 @ymori117

20年4月25日

NATROMのブログ:「論座」の新型コロナ感染症の記事から学べること - natrom.hatenablog.com/entry/2020/04/... 「佐藤章氏が間違うのも仕方がない」「上昌広・医療ガバナンス研究所理事長も、陽性的中割合の計算を間違ったぐらいだ」ちょww

タグ:

posted at 12:56:43

numachi @numachinomajo

20年4月25日

「論座」の新型コロナ感染症の記事から学べること - NATROMのブログ natrom.hatenablog.com/entry/2020/04/...
あいかわらず良い解説だなあ。ほんといつも問題点を判りやすく書いてくださる。

タグ:

posted at 12:39:48

はっぴぃ @shiawaseahiru

20年4月25日

10ccをすぐ飲んでしまったので追加で少しあげたのにまたすぐ飲んじゃって未練たらしいうぃるそん君。 pic.twitter.com/iqIva5eDAM

タグ:

posted at 12:28:05

飯田泰之 @iida_yasuyuki

20年4月25日

変わる……以前に

各制度のNPO法人への適用のように「元々排除する意図はないが,明示しなかったので対象じゃないように現場がうけとってしまった」というケースもある

気づいた人が声を上げる必要があるよね

タグ:

posted at 11:46:04

bechori @bechori777

20年4月25日

静止画もお納めください(定例)これは爽やかすぎてしまっている…!!

#文具の杜
#杜の四季インク
#青葉晴嫋
#インク沼 pic.twitter.com/U7VfJiZB9D

タグ: インク沼 文具の杜 杜の四季インク 青葉晴嫋

posted at 11:45:41

TEN/他@情報整理屋さん @tenfirefighter

20年4月25日

経産省の業種別支援リーフレットもめっちゃわかりやすい(医療業界向けもあるよ)

www.meti.go.jp/covid-19/leafl... twitter.com/tenfirefighter... pic.twitter.com/aHMQOhlFR2

タグ:

posted at 11:43:50

bechori @bechori777

20年4月25日

文具の杜さんの待望の新色がやってきた!ので早速書いてみました✍🏼明るいターコイズのような、ミントグリーンのような、これからの季節にぴったりな爽やかな色には間違いない…!☘️☘️文具の杜さんありがとうございます!静止画に続きます→

#文具の杜
#青葉晴嫋
#インク沼

@office_vender pic.twitter.com/h5hgyRRy1g

タグ: インク沼 文具の杜 青葉晴嫋

posted at 11:43:03

飯田泰之 @iida_yasuyuki

20年4月25日

給付金といい,事情ある別居者への個別支給といい,

そして雇用調整助成金の要件緩和と拡充……

どれも素早く大変よい対応.みんなの声で政策かわるぞ!

(はじめからやっときゃ全然印象違ったのに)

→休業手当100%国が補助、小規模企業向け雇用調整助成金: www.nikkei.com/article/DGXMZO...

タグ:

posted at 11:35:35

数学女子 @phasetrbot

20年4月25日

JuliaでPDE解くのはともかく、それをどう動画化すればいいかを調べに行く必要がある。

タグ:

posted at 11:17:18

白樺香澄 @kasumishirakaba

20年4月25日

嘘みたいなハサミの使い方、宙を舞う卵の殻、安定の高笑い、安定の目分量、切れてないニラ、乱れ飛ぶニラ、びしょびしょで直置きのビニール袋、なのに出来上がりはちゃんと美味しそう。「私たちの見たい平野レミ」の全てがここにある。 twitter.com/remi_hirano/st...

タグ:

posted at 11:11:13

数学女子 @phasetrbot

20年4月25日

Julia の GRUtils あたりで plot するといちいちビューアー的なものが立ち上がるっぽいのだが、単に画像だけ生成したい場合どうすればいいのだろうか。

タグ:

posted at 11:07:10

統計たん @stattan

20年4月25日

ゆるぼ:いい感じのJulia教材 twitter.com/mgg03300/statu...

タグ:

posted at 10:38:30

しらんがな @mgg03300

20年4月25日

: がいろんな意味に使われるのが未だに慣れない… twitter.com/stattan/status...

タグ:

posted at 10:26:25

数学女子 @phasetrbot

20年4月25日

Julia 単純な動画サンプル GRUtils youtu.be/U34q8iOqjQo @YouTubeさんから
Julia での初投稿。単にサンプルコードをコピペしただけ。

タグ:

posted at 10:25:13

統計たん @stattan

20年4月25日

Julia,型の指定の感じが慣れない。。。

タグ:

posted at 10:23:47

平野レミ @Remi_Hirano

20年4月25日

皆さん、おうちゴハン頑張ってますか?皆さんの力になりたくて、私もノーメイクで頑張りました(笑)。自宅で撮った「ふわふわのニラ玉」です。文字のレシピはこちらをどうぞ→remy.jp/recipe/061.html pic.twitter.com/IBEeWr6rTK

タグ:

posted at 09:57:22

非公開

タグ:

posted at xx:xx:xx

ふろ仙人 @doctorhirosan

20年4月25日

@kumadesukedooo ご無事で何よりです…😭😭
胸骨圧迫本当に大事ですよね!躊躇わず!

タグ:

posted at 09:36:37

非公開

タグ:

posted at xx:xx:xx

黒木玄 Gen Kuroki @genkuroki

20年4月25日

#Julia言語 nbviewerで閲覧しながら、ブラウザ操作でリロードすると、プロット部分が表示されたり(一部分だけ表示されることも多い)、プロットされていた部分が消えたりと不安定な挙動を示す。

上手く表示される場合もあるので、適切に訂正すればnbviewerでマウスぐりぐりプロットが可能になるはず。

タグ: Julia言語

posted at 09:18:27

黒木玄 Gen Kuroki @genkuroki

20年4月25日

#Julia言語

Plots.jl + plotlyjs() + nbviewer

の組み合わせで上手く行く条件が実はよく分からない。

添付画像は上手く行っているばあい。

nbviewer.jupyter.org/gist/genkuroki...

ブラウザは iPad Chromeで「PC版サイトを見る」にした直後。再度読み込むとプロット部分が消える。

twitter.com/nextstep2osx/s... pic.twitter.com/1BvOlrxPqK

タグ: Julia言語

posted at 09:15:21

数学女子 @phasetrbot

20年4月25日

JuliaでのPNG生成、Plots/GRと GRUtils はどちらがいいのだろう。最終的にはPDEの時間発展のような動画が作りたくて、 とにかく速度がほしい。

タグ:

posted at 09:00:20

Nyoho @NeXTSTEP2OSX

20年4月25日

やってみました。しかし plotly() のときと同じように plotlyjs() でも Plotly が見つからないエラーが出て動きませんでした。requirejs で CDN から plotly を読むようになっているのになぜ動いてないんだろう…… twitter.com/genkuroki/stat...

タグ:

posted at 08:29:46

Cristian Garcia @cgarciae88

20年4月25日

Can somebody in the Julia community kindly explain why Pkg isn't a script like in all other package managers (npm, pip)? Having a hard time imagining how to take Julia code into production...

@JuliaComputing #JuliaLang

タグ: JuliaLang

posted at 07:32:19

カタバミ @oroti56562

20年4月25日

そのリプ欄は見る気がしないがこれ思い出した twitter.com/samoedon/statu... pic.twitter.com/MctkjwLwVA

タグ:

posted at 07:21:10

非公開

タグ:

posted at xx:xx:xx

黒木玄 Gen Kuroki @genkuroki

20年4月25日

大昔にやっていたとあるゲームでは、最初にダガーを手に入れて、ダガーを装備した状態で祭壇で祈りまくって、ダガーに麻痺+1がつけば勝ちというパターンがあったな。

麻痺で動きが止まるモンスターで最強のものと戦うと、楽勝できてかつレベルが一挙に上がる。

タグ:

posted at 06:43:21

黒木玄 Gen Kuroki @genkuroki

20年4月25日

#数楽 Cholesky分解以外の行列分解も、行列の対角化を記述しているとみなせる可積分系の離散化だとみなされるものがあります。

こんな感じで、行列分解、可積分系、離散可積分系、行列の対角化などなどが全部繋がる。

タグ: 数楽

posted at 06:33:21

黒木玄 Gen Kuroki @genkuroki

20年4月25日

#数楽 オープン戸田格子は質点が反発して離れて行く状況を記述する常微分方程式の一種なのですが、時間がたつと質点の間の距離が大きくなり、運動量は一定に近付きます。これをLax表示で見直すと行列の数値的対角化になっている。

その離散版は行列の対角化の数値計算のアルゴリズム!

タグ: 数楽

posted at 06:31:36

黒木玄 Gen Kuroki @genkuroki

20年4月25日

#数楽 オープン戸田更新は連続時間における常微分方程式による時間発展なのですが、その離散化は本質的に「Cholesky分解の積の順序交換」の繰り返しです。「Cholesky分解の積の順序交換」の繰り返しでオープン戸田格子のLax表示のL行列が数値的に対角化されます。

タグ: 数楽

posted at 06:28:43

黒木玄 Gen Kuroki @genkuroki

20年4月25日

#数楽 オープン戸田格子のLax表示において、時間発展するL行列の非対角成分はq_iとq_{i+1}の差が大きくなると0に収束する形になっており、対角成分はp_iです。

添付画像からオープン戸田格子の時間発展によって、Lax表示のL行列が数値的に対角化されることがわかる。

nbviewer.jupyter.org/gist/genkuroki... pic.twitter.com/m6AfVDVSDj

タグ: 数楽

posted at 06:26:34

黒木玄 Gen Kuroki @genkuroki

20年4月25日

#数楽

添付画像は何も知らないと面白くも何ともない画像。

実線はオープン戸田格子の常微分方程式を数値的に解いた結果のグラフで、整数時刻におけるドットは本質的に「Cholesky分解の積の順序交換」の繰り返しのプロットです。ぴったり一致しています!!!

nbviewer.jupyter.org/gist/genkuroki... pic.twitter.com/A8RdQ6CSaT

タグ: 数楽

posted at 06:19:53

黒木玄 Gen Kuroki @genkuroki

20年4月25日

#数楽 Cholesky分解1つを取っても、数学的に非常に楽しめるネタになっています。「Cholesky分解の積の順序の交換」によって定義される行列の離散時間発展は、オープン戸田格子の誤差ゼロの離散化になっています!!!

戸田格子は色々な意味で良い数学ネタになっています。

nbviewer.jupyter.org/gist/genkuroki... pic.twitter.com/qsW86drBy5

タグ: 数楽

posted at 06:16:42

黒木玄 Gen Kuroki @genkuroki

20年4月25日

#数楽 数学畑でこれから線形代数の解説を書こうと思っている人は、 #Julia言語 のマニュアルの行列分解の節

docs.julialang.org/en/v1/stdlib/L...

も見ておいた方がいいかも。そこに書いてある行列分解達が便利な道具として普通に使われています。(実際にはそこのリストにない分解もよく使われている。)

タグ: Julia言語 数楽

posted at 06:11:20

黒木玄 Gen Kuroki @genkuroki

20年4月25日

#数楽 数学に「線形代数」という分野はなくて、「線形代数」は何となく線形代数と呼んでよさそうな便利な代数的道具達を集積しておく場所になっています。

だから、「線形代数」の名のもとで解説される道具の難易度の幅は大きい。滅茶苦茶易しいものから、非常に難しい道具まで色々。

タグ: 数楽

posted at 05:40:17

黒木玄 Gen Kuroki @genkuroki

20年4月25日

#数楽 線形代数に限らず、「易しいことから順番に」が基本。

線形代数では「以下同様」が数学的に通用し易いので、最初の方で考え方のコツがつかめれば、それ以後は「以下同様」の感覚で楽に理解できます。

2次元や3次元で図を描ける場合にたくさんの絵を描いておいて、「以下同様」で高次元に拡張。

タグ: 数楽

posted at 05:36:11

黒木玄 Gen Kuroki @genkuroki

20年4月25日

#数楽 「標準メニュー」の段階で四苦八苦している場合には梶原健さんによる解説がきっと役に立つと思います。

twitter.com/genkuroki/stat...

タグ: 数楽

posted at 05:33:30

黒木玄 Gen Kuroki @genkuroki

20年4月25日

#数楽 沢山ある行列の標準形と沢山ある行列の積への分解について色々知っていると、「真の線形代数ユーザー」(笑)になれた感じがして来るはずです。

しかし、もっと基本的な基礎概念について理解していないと、「ケース・バイ・ケースで複雑なことをやるだけ」のような感じになって苦しくなる。

タグ: 数楽

posted at 05:29:36

黒木玄 Gen Kuroki @genkuroki

20年4月25日

#数楽 沢山ある「行列分解」についても詳しく書かれた線形代数の入門書があれば、線形代数ユーザーにとってありがたいかもしれませんが、そういう本は知らない。

ある程度以上の数学を知っていると、各種分解についての視界を広げることができる。

docs.julialang.org/en/v1/stdlib/L...

タグ: 数楽

posted at 05:27:22

黒木玄 Gen Kuroki @genkuroki

20年4月25日

#数楽

梶原健さんは線形代数の本を何冊も書いていますが、私の手元にあるのは

『線形代数のコツ』
『基礎からわかる!しっかりわかる!!線形代数ゼミ』

の2冊です。

タグ: 数楽

posted at 05:21:41

黒木玄 Gen Kuroki @genkuroki

20年4月25日

#数楽 あと、まだ紹介したことがなかったと思いますが、梶原健さんによる線形代数の解説も魂がこもっている感じで非常にいいです。たぶん、長谷川本よりも梶原本の方が初心者向けかもしれない。

www.amazon.co.jp/s?k=%E6%A2%B6%...

タグ: 数楽

posted at 05:18:25

黒木玄 Gen Kuroki @genkuroki

20年4月25日

#数楽 長谷川浩司『線型代数』は非常によい本です。

線形代数の本の多くは、教える側が了解している「標準メニュー」を単に分かり易く見えるように解説するだけの本になっています。そういう本を他人には勧めにくい。

長谷川本は全然違っていて魂がこもっています。

twitter.com/MathSorcerer/s...

タグ: 数楽

posted at 05:13:13

黒木玄 Gen Kuroki @genkuroki

20年4月25日

#Julia言語 パッケージ Foo.jl について

using Foo: Foo

として、Fooで定義されているfを

Foo.f(x)

のように利用できることを知っていると便利な場合があります。

お気楽に使う場合には

using Foo

を使い、名前空間が重なるのがどうしても嫌な場合には

using Foo: Foo

を使うという感じ。

タグ: Julia言語

posted at 05:02:43

黒木玄 Gen Kuroki @genkuroki

20年4月25日

#Julia言語

using Plots
pyplot()
using PyPlot: PyPlot, plt

のようにして、Plots.jl経由でのPyPlotの利用と直接的なPyPlotの利用を両立できます。Plots.jlではPyPlotの全ての機能が使えるわけではないので、こうした方が便利な場合があります。

twitter.com/genkuroki/stat...

タグ: Julia言語

posted at 05:00:06

黒木玄 Gen Kuroki @genkuroki

20年4月25日

#Julia言語 他の例(@ の後の空白は除く)

using SymPy: SymPy, sympy, Sym, @ vars, @ syms, oo, PI

として、

@ vars x
(x^120 - 1).factor()

とか

@ vars a positive=true
@ vars x
sympy.integrate(exp(-a*x^2), (x, -oo, oo))

のようにしてSymPyを使う。

Sym(π) は PI になります。 pic.twitter.com/2Xue6NDCOT

タグ: Julia言語

posted at 04:57:59

黒木玄 Gen Kuroki @genkuroki

20年4月25日

#Julia言語 には「お気楽に使ってよい」という雰囲気があると思う。だから、

using パッケージ名

がよく使われる。名前がぶつかる場合には例えば、

using PyPlot: PyPlot, plt

のようにパッケージを読み込んで、

PyPlot.plot



plt.plot

を使うようにすればよい。続く

タグ: Julia言語

posted at 04:50:11

黒木玄 Gen Kuroki @genkuroki

20年4月25日

#Julia言語 試してみようと思いましたが、データのcsvファイル達の入手法が説明だけだとよくわからなかったので断念。

plotly() ではなく plotlyjs() ならば

nbviewer.jupyter.org

経由でマウスでぐりぐりできるグラフを閲覧できる場合があります(かなり不安定)。

twitter.com/NeXTSTEP2OSX/s...

タグ: Julia言語

posted at 04:39:07

黒木玄 Gen Kuroki @genkuroki

20年4月25日

メモリ16GBは人権

速いは正義

タグ:

posted at 04:00:33

Anmol Rana @ref_sys

20年4月25日

I knew it👀! ParallelKmeans.jl in #julialang is orders of magnitude faster than Py's scikit-learn and R implentations, and even Clustering.jl. And it all started as a heated discussion on #julialang ❤️ discourse abt
optimizng code!
pydatablog.github.io/ParallelKMeans... pic.twitter.com/jvozI9Equj

タグ: julialang

posted at 03:24:34

マナマの生マナマコ @chor0sK

20年4月25日

"Integers and Floating-Point Numbers · The Julia Language" docs.julialang.org/en/v1/manual/i...

コンピュータの浮動小数点の話はここにええ感じの資料がまとまっているなーということしか分からない
#julialang

タグ: julialang

posted at 01:46:39

ごまふあざらし(GomahuAzaras @MathSorcerer

20年4月25日

知らなかった.Julia ラッパーがあるのね

github.com/plotly/Dash.jl

タグ:

posted at 01:35:06

非公開

タグ:

posted at xx:xx:xx

@genkurokiホーム
スポンサーリンク
▲ページの先頭に戻る
ツイート  タグ  ユーザー

User

» More...

Tag

» More...

Recent

Archive

» More...

タグの編集

掛算 統計 超算数 Julia言語 数楽 JuliaLang 十分 と教 モルグリコ 掛け算

※タグはスペースで区切ってください

送信中

送信に失敗しました

タグを編集しました