黒木玄 Gen Kuroki(@genkuroki)/2019年10月07日

Re: RT

#Julia言語で音入りもmp4ファイルを楽に作るにはどうすれば良いのだろうか？

ffmpegが定跡？

twinklesmile.blog42.fc2.com/blog-entry-187...
ffmpeg で連続画像ファイルと音声ファイルを結合して動画ファイルにする

* pngファイル群作成
* wavファイル作成
* ffmpegを呼ぶ

誰かJuliaで実装してないかな？

タグ： Julia言語

posted at 23:49:50

#統計以下のリンク先の応用例2については

genkuroki.github.io/documents/Stat...
Ridge正則化とStein推定量

の3.8節を参照。

twitter.com/genkuroki/stat...

タグ：統計

posted at 23:38:54

ぱおとんず @paotongzi

野党の皆さん！今回の入試改革に関して、政界における出発点であるこの方もそろそろ追求して下さい！
文科省が改革強行に拘る背景に、この人物とB社(他にもT中H蔵など)が関与しているはずです！ m.newspicks.com/news/1295889/

タグ：

posted at 23:23:34

ラプラス・ルンゲ・レンツベクトルについての資料要望があったので

僕の演習ノート：sites.google.com/site/naotoshir...　のp46から
国場さんの解説記事：webpark1739.sakura.ne.jp/atsuo/Lenz2007...
（立川さんの講義ノートはなぜか消えてた…）

僕のノートは国場さんの解説の最初の方だけを薄くしたような感じです。

タグ：

posted at 23:23:34

Teiji Kunihiro @KunihiroTeiji

佐藤光著「群と物理」(丸善)、拙著「量子力学」(東京図書)でも説明されている。また、B.L. van der Waerden ‘Sources of Quantum Mechanics’(Dover Books on Physics)では当のパウリの原論文が収録されている。古典力学レベルの話は、少なくとも第2版ごろまでの、ゴールドスタイン「古典力学」にある。 twitter.com/perfect_inside...

タグ：

posted at 23:20:24

#統計自由度に関する間違った説明関連

自由度の説明は、帰無仮説に対応する確率分布にパラメーターが入らなければ易しくなるが(レベル1)、2×2の分割表の独立性検定のようにパラメーターが入る場合には難しくなる(レベル2)。

twitter.com/genkuroki/stat...

タグ：統計

posted at 23:13:34

#統計ところが、

a b
c d

でa+b, c+d, a+c, b+dをすべて固定すると、aの値を決めるだけで残りのb,c,dの値が決まることになり、aの1個の分だけの自由度しか残らないので、自由度は1になる

というような

間違った説明

が教科書によく書いてある。

レベル1の回答にミスリードされている感じ。

タグ：統計

posted at 23:13:33

#統計例：2×2の分割表

a b
c d

で a+b と c+d を固定したもの全体上の2項分布×2はd₁=2個のパラメーターを持つ。独立性の帰無仮説を満たす確率分布にはd₀=1個のパラメーターが入る。だから、対応するカイ二乗検定で使うカイ二乗分布の自由度はd₁-d₀=1になる。

タグ：統計

posted at 23:13:33

#統計例：2×2の分割表

a b
c d

で a+b+c+d=n を固定したもの全体上の4項分布はd₁=3個のパラメーターを持つ。独立性の帰無仮説を満たす確率分布にはd₀=2個のパラメーターが入る。だから、対応するカイ二乗検定で使うカイ二乗分布の自由度はd₁-d₀=1になる。

タグ：統計

posted at 23:13:32

#統計

「カイ二乗分布の自由度とは何か？」への回答レベル2

カイ二乗分布の自由度＝パラメーター空間の部分空間の余次元

d₁次元空間のd₀次元部分空間の余次元の定義はd₁-d₂になる。

タグ：統計

posted at 23:13:32

#統計自由度の話 Part 2

「カイ二乗分布の自由度とは何か？」への回答レベル1

カイ二乗分布の自由度＝対応する多変量正規分布の台の次元

例えば、r項分布は中心極限定理によってr-1次元空間上の多変量正規分布で近似される。r項分布のカイ二乗検定で使うカイ二乗分布の自由度はr-1になる。

タグ：統計

posted at 23:13:27

Astellon @astellon_music

Cutoff/Qやってみました pic.twitter.com/rvU7sVUpkq

タグ：

posted at 23:13:04

特殊関数を使って解くと、エネルギー縮退がn^2になるのが「たまたま解いたらそうなった」程度にしか説明されない。ラプラス・ルンゲ・レンツベクトルを用いた代数的方法を使うと、その点も綺麗に説明される。

タグ：

posted at 23:02:16

そもそも水素原子が初めて解かれたのは、シュレディンガー方程式が提案されるより前（！）なので、特殊関数で解かれたはずはない。だが、後者に言及した量子力学のテキスト・講義ノートは、（自分の演習を除くと）知ってる範囲では立川さんの講義ノートだけだと思う。

タグ：

posted at 22:56:24

量子力学の講義ではほぼ必ず水素原子が出てきて「水素原子が解けたのは歴史的にも重要」と説明されることが多い。だが水素原子が初めて解かれた際は、大体の講義でやる特殊関数を用いた方法ではなく、ラプラス・ルンゲ・レンツベクトルを用いた代数的方法で解かれたのに、後者への言及はほとんどない。

タグ：

posted at 22:54:11

しぎょういつみ @s_itsme

高3の秋になって急に大学受験を決心しても、今の制度なら、ちゃんとした指導者のもとでそれなりの努力をすれば、現役で志望大学に合格することも夢ではありません。英語民間試験を活用する制度になると、高3の秋では民間試験を受ける手続きが間に合わず、浪人してもう1年勉強しなければなりません。

タグ：

posted at 22:49:05

勝又清和 @katsumata

「消えた戦法の謎ができたー」
「よう、私は24年後から来たお前だ。矢倉は5手目▲６六歩から▲７七銀に戻り、５二飛戦法はゴキゲン中飛車の名前で大流行、石田流も角換わり早繰り銀も指されている。
そして29才名人対17才七段の将棋は17才のひねり飛車に対して名人がたこ金戦法で勝った。」
「うそ」 pic.twitter.com/CHZrzs4uIn

タグ：

posted at 22:41:48

左巻健男（サマキタケオ） @samakikaku

#数楽

例3: 分配函数

Z(β) = ∫ exp(-βH(w))φ(w)dw

の漸近挙動が対応するゼータ函数

ζ(-s) = ∫ H(w)^{-s} φ(w) dw

の極から得られる。(ベイズ統計における分配函数の漸近挙動はこの方法で得られる。)

タグ：数楽

posted at 22:33:07

暮らしのなかのニセ科学 (平凡社新書) 左巻健男 www.amazon.co.jp/dp/4582858473/... @amazonJPさんから
＊うわっ！平凡社新書で1位になって「一時的に在庫切れ; 入荷時期は未定です。」
すぐに入荷されると思うけど、昨日の産経記事のおかげかな。

タグ：

posted at 22:29:39

ごまふあざらし(GomahuAzaras @MathSorcerer

#数楽

例1: スターリングの公式

log N! = Σ_{n=1]^N log n

の漸近挙動の発散部分は、Riemannのゼータ函数の導函数の-1倍 -ζ'(s)のs=1での2位の極から出て来る。

例2: x以下の素数の個数の漸近挙動は

-ζ'(s)/ζ(s)

の極(=ζ(s)の極と零点)から得られる。

nbviewer.jupyter.org/github/genkuro...

タグ：数楽

posted at 22:28:23

これ不思議よねーって教科書読んで思った。 twitter.com/genkuroki/stat...

タグ：

posted at 22:08:46

田中健一／講師業、著述業 @TNK_KNCH

#数楽逆Mellin変換の応用関連

仮にゼータ函数の定義を「何かのMellin変換」と定義する。そのとき、もとの何かを「ゼータ函数の逆Mellin変換」で表示することによって、ゼータ函数の極の情報からもとの何かの漸近挙動が分かる。

これがゼータ函数の最も基本的な使い方。

twitter.com/genkuroki/stat...

タグ：数楽

posted at 22:02:04

田中健一／講師業、著述業 @TNK_KNCH

文科省は各大学に試験問題を公表するよう求めていて、各大学はもちろん公開しています。しかし英語民間試験業者は「我々は問題を公表しません。なぜなら我々は問題を公表しないからです」みたいな態度で、これに対しても文科省は「いやいや、公開しなさいよ」と重ねて要請することはありません。

タグ：

posted at 21:43:04

非公開

タグ：

posted at xx:xx:xx

左巻健男（サマキタケオ） @samakikaku

本日印象的だったのは

民間試験に不備があります、
文科省は改善要請を出しました、
試験業者の対応策が不十分です、
文科省は再度の改善要請は出しません、

みたいなケースが沢山あって、試験業者がやりましたと言えばやったことになり、試験業者が大丈夫だと言えば大丈夫なことになる構図です。

タグ：

posted at 21:35:59

非公開

タグ：

posted at xx:xx:xx

非公開

タグ：

posted at xx:xx:xx

ニセ科学批判をすればそれに命をかけているニセ科学側から強い攻撃をされる。
これまで研究室に来られて「月夜の晩ばかりではない」などと言われたこともある。
裁判もされる。
でももう残り少ない人生だし、ガハハと笑って受けて立ちたい。
何事も楽しめるといいな。

タグ：

posted at 21:25:32

非公開

タグ：

posted at xx:xx:xx

とおん舎 @toonsha_m

息子(小5)が素数を調べる宿題があって、自分で調べるのがめんどくさいからと言って、スクラッチで素数を調べるプログラミングを作っていた！

天才か！！😆✨ pic.twitter.com/w8G70Sp4is

タグ：

posted at 21:14:28

Karolina Fortin (Sho @oneye

山本健人(外科医けいゆう, Takehi @keiyou30

Newest Shogi Project for Fans Abroad! ShogiHub Reactivation youtu.be/2LFa4kzCr-Q via @YouTube

タグ：

posted at 21:14:20

ほんの２年ほど前までは、こうしたニュースがあるたび、「急性胆管炎ってどんな病気？調べてみました」的な、誤った内容のブログ記事が大量に現れ、検索結果の上位を占めていました。身をよじるほど悔しい思いでした。
今は上位をガイドラインや学会関連の記事が占め、私のこの記事も１ページ目です。

タグ：

posted at 21:06:59

教育新聞 @kyoiku_shimbun

【英語民間試験】東京私大教連「即時中止を求める」
www.kyobun.co.jp/news/20191007_...

タグ：

posted at 21:04:00

FUJI Goro @__gfx__

まじか…qmailってかdjbは改めてすごいな。

"他人のコードや設計を見て1番これはあり得ないだろと思う実装はありますか？への回答" jp.quora.com/%E4%BB%96%E4%B...

タグ：

posted at 20:56:12

非公開

タグ：

posted at xx:xx:xx

abap34 @abap34

山本健人(外科医けいゆう, Takehi @keiyou30

JuliaBox無料化廃止で悲しい人はGoogle ColabでJuliaを使いましょう

タグ：

posted at 20:23:31

Yahoo!ニュースで記事を公開。
元プロ野球選手の金田正一氏が逝去されました。
原因とされている急性胆管炎は、私たちが扱う消化器疾患の中でも一二を争うほど緊急性が高く、重症化するとあっという間に命の危機に瀕することがあります。
その理由について解説しました。
news.yahoo.co.jp/byline/yamamot...

タグ：

posted at 19:56:56

羽藤由美 @KITspeakee

もう一度言いますね。野党が英語民間試験の導入延期法案提出を検討してくださるとのこと。「やっとっこまで来たか」と本当に有り難いですが，一方では多数決で否決されればそこでオシマイ。今，どうしても必要なのは，それを許さない世論の盛り上がりと「できない」という証拠です。 twitter.com/KITspeakee/sta...

タグ：

posted at 19:45:53

岡崎渉 @z_karaagen

残念ながら…

「逃げていくのは、教職の現状や社会的な位置を自ら調べて考える、教員になってほしい学生だ」

（記者解説）教員のなり手が減少　「ブラック職場」敬遠、対策も後手　編集委員・氏岡真弓：朝日新聞デジタル www.asahi.com/articles/DA3S1...

タグ：

posted at 19:14:14

数学女子 @phasetrbot

日本で女子の理系進出が少ないの、理系が奴隷のような扱いしか受けないから、親も行かせたがらないし、当人も嫌がると言うだけの話なのではないかという気もする。

タグ：

posted at 18:49:06

マナマの生マナマコ @chor0sK

今日はむっちゃ多重ディスパッチとmap do文の恩恵を受けた。#julialang

タグ： julialang

posted at 18:36:15

漏洩がないことの徹底を確認しようもない）。
ほかにも無数の問題が露呈。公平性や実態を踏まえない参加要件が策定、そこに記載されていない以上「試験団体に要請します」を繰り返すほかないかの文科省。
受験の未来を考えて設計したとは到底思えない制度。延期法案を提出しなければならない事態です。

タグ：

posted at 18:15:32

試験実施主体が講座や問題集でお金を儲けることに利益相反があるのではないかという問いに、「各業者に試験問題を作成する部署とそうでないところは部署を分けていて、問題の漏洩がない対策をしていることを確認している」と噛み合わない答弁（そもそも、文科省が張り付いているわけにもいかないから

タグ：

posted at 18:13:20

高３も合否が決まる前に申し込まなければいけない。合格後、返金してもらえないのも不合理。
障がいのある受験生につき「合理的な配慮をしていることを公表する」しかし、九州の障がいのある高校生が、業者から「配慮があるのは東京の会場だから東京で受験を」といわれたら？文科省はよしとするのか。

タグ：

posted at 18:11:41

大学入試問題の公表は義務づけられている。しかし、民間試験は公表されない。そうすると、出題ミスがあっても発見できない、どうするのか、という問いには、「事業者に判断していただく」（？）。受験生が出題ミスに気づいたらどうすればいいのか。事業者が公表しないなら、出題ミスを検証できません。

タグ：

posted at 18:07:47

２分の１の補助を概算要求。予算組みしない自治体の受験生は補助されない。また、交通費や宿泊費がかかるのは、離島の受験生だけではない。
　地域格差、所得格差についての認識を問うても、「会場設置の協力の要請を重ねている」「国立大学には会場費の減免を要請している」との要領を得ない答弁。

タグ：

posted at 18:06:04

受験料の軽減措置を受けられる経済的に厳しい受験生についての質問に資料はなし、口頭のみ。大学入試センターに申請する。大学入試センターが認定する。生活保護者や奨学金受給者などを想定。軽減分は業者が負担。
試験会場から遠い生徒の交通費や宿泊費については？離島の受験生についてのみ想定。

タグ：

posted at 18:03:56

受験の早期化を文科省はよしとしているのかとの問いには、「高校大学の関係者との協議も経て実施要項を出した」と繰り返すのみ。
　経済的に厳しい生徒に対する軽減策として、たとえば、英検準１級９，８００円→９，３１０円など、数百円レベル。２万円を超す検定料のところも。

タグ：

posted at 18:02:59

予備校の先生方からも、必ずしも出願資格としない大学でも３７％程度が合否判定の資料として加点するとしており、家庭の経済力が大学入試の合否を決定することになる、検定料だけでなく民間業者のテキスト代も高いこと、など切実なご指摘。
様々な質問に対する文科省担当者の説明は残念なもの…。

タグ：

posted at 17:59:29

大学入学テストで英語民間試験の活用への懸念が高まっています。本日の第3回英語民間試験導入問題野党合同ヒアリングで、現時点でも深刻な問題が多いことが一層明らかに。複数の高校生から、将来が企業の利潤のためにふりまわされる事態や、高校が予備校化する懸念等切実なご意見をいただきました。

タグ：

posted at 17:54:29

荒木健太郎 @arakencloud

台風19号の開眼の様子．
大型で非常に強い勢力となった台風19号は今後北上して，今週末の12～13日に本州にかなり接近するおそれがあります．今後の動向に本気でお気をつけ下さい． pic.twitter.com/BJf38WITZf

タグ：

posted at 17:48:23

理解の限界に近い領域で実際にさまよった経験がないと、そういう経験がどのように役に立つかがピンと来ないということはあるかも。

タグ：

posted at 14:55:22

数学科に限らず、数学的に複雑なことを学んだ経験がある人達は、理解することに関しては人類の限界に近い領域の経験者。

今現在持っている知識よりも、自分が理解できる限界領域に近いところでさまよった経験を評価して欲しいと思います。

関連↓
twitter.com/mathsorcerer/s...

タグ：

posted at 14:55:21

名取宏（なとろむ） @NATROM

ポイントは、よしんば盗用で不正だとしても、論文の内容自体に異議が唱えられたわけではなく、STAP細胞のように研究内容が疑わしいわけではない点。過剰診断の議論を理解していない/理解したくない人は、論文の内容の是非には触れず、「Welchは研究不正で大学を辞職した」とだけ言い立てるので注意。

タグ：

posted at 13:36:43

名取宏（なとろむ） @NATROM

大学は告発者である準教授を論文著者に入れるように提案したが、Welchはこの提案を蹴って大学を辞職した。Welchの言い分は「研究に参加していない人物の名前を論文に記載するほうが不正である。不正してまで大学に留まりたくない」。NEJM誌は「盗用ではない。論文撤回に該当しない」と判断した。

タグ：

posted at 13:35:55

名取宏（なとろむ） @NATROM

余談だが、Gilbert Welchは大学は辞職したが論文は撤回されていない。盗用が疑われたのだが、Welchの見解は「私たちが以前から言ってきたことから自然に導かれるアイディアであり盗用ではない」。 twitter.com/miakiza2010090...

タグ：

posted at 13:35:04

RochejacMonmo @RochejacMonmo

英語民間試験の利益相反事案。高校会場を確保するために優先的に当該高の学生の座席を確保したり、自団体の講座等の受講で受験料を割り引くなどの手法について現時点では文科省として規制できないとのこと。このような方法は、実質的に特定の受験者の差別的取り扱いにすら該当するのではないか。

タグ：

posted at 12:57:41

RochejacMonmo @RochejacMonmo

英語民間試験の問題と解答の公表に関して。文科省担当者は各団体の判断に任せるとのこと。しかし、異なる実施時期の試験を比較するためには、問題非公開＋過去問の再利用等でスコアの標準化が必須。「問題の公開」と「実施時期の異なる試験のスコアを相対評価で比較」とは本質的に両立しえない。 pic.twitter.com/ouODX65b7w

タグ：

posted at 12:50:09

カミケ @KKami1115

JuliaBox無料版使えなくなるのか

タグ：

posted at 11:16:08

#Julia言語

b = zero(Int8) # Int8の0
b = 4 # Int64の4

とすると、Int8の0が格納された変数bがInt64の4が格納された変数に上書きされます。

1つ前のツイートでaは配列であったことに注意！

タグ： Julia言語

posted at 11:14:19

#Julia言語

Juliaは初心者がやりがちなことをやっても適切に動くようにするために相当な手間をかけている節が見られます。

例えば、

a = zeros(Int8, 3) # Int8の長さ3の1次元配列
a[1] = 4 # Int8の配列の第1成分にInt64の4を代入

とすると、a[1]はInt8の4になってくれます。

タグ： Julia言語

posted at 11:11:49

#Julia言語あ、もしかしたら、以上を読んで「Juliaは難しいからパス！」と思う人がいたら、誤解を招く解説をしてしまったことになり、申し訳ない感じです。

ほとんどの数値的な計算はFloat64とInt64の範囲内で済むので、細かい仕様を知らなくてもJuliaは便利に使えます。

タグ： Julia言語

posted at 11:03:33

須山敦志 Suyama Atsushi @sammy_suyama

#Julia言語数学用語を使えば、たぶん、

(半)環 R 上の加群 M について、x∈M について one(x) は R の要素としての 1 になり、M がR上ランク1の自由加群で基底 u が決められているならば、x∈M=Ru に関する oneunit(x) は u になる

ということなのだと思います。実際にそうなら色々わかりやすい。

タグ： Julia言語

posted at 10:59:09

Juliaboxの無料ver.が10月末までらしいので有料版に切り替えた

タグ：

posted at 10:58:58

#Julia言語さらに

addone(x) = x + one(x)

よりも

addoneunit(x) = x + oneunit(x)

の方が適用可能範囲が広がります。なぜならば、

Day(3)

のように(日数の)次元を持つ量もJuliaで扱えるからです。 one(Day(3)) は Int64 の 1 になってしまいます。oneunit(Day(3)) ならば Day(1) になる。

タグ： Julia言語

posted at 10:55:14

#Julia言語

add1(x) = x + 1

よりも、

addone(x) = x + one(x)

の方がより広い範囲で使い回しが効く函数になっています。

add1(x)だと1がInt64なので、x=Int8(3)のとき、add1(x)はInt64の4になってしまいます。add1(x)だとInt8を使ってメモリ効率を上げることができない。

addone(x)ならOK.

タグ： Julia言語

posted at 10:55:13

#Julia言語問題はそういう自動的なpromoteで対処できない場合にどうするかです。

そういう場合には函数の引数の型が適切に伝搬するようにコードを書く必要が生じます。

型の伝搬を記述するための仕組みをJuliaは豊富に提供してくれています。これ、Juliaでは本質的に重要。

タグ： Julia言語

posted at 10:55:13

#Julia言語自動的に型をpromoteしてくれる機能のおかげで、「Float64の計算なのだから、1ではなく1.0と書かなければいけない」というように神経質になる必要がなくなっています。

数式内でFloat64とInt64が混ざっていれば、Int64がFloat64の側に自動的に変換された後にコンパイルされる。

タグ： Julia言語

posted at 10:55:12

#Julia言語以上でも以下でも64bit環境を仮定。

Juliaで1はInt64の1を意味し、3.2はFloat64の3.2を意味します。このとき、

3.2 + 1

の1は自動的にFloat64の1.0にpromoteされ、実際には、

3.2 + 1.0

が計算されます。

タグ： Julia言語

posted at 10:55:12

しぶちん⏳ @pmagshib

舞台裏が見える瞬間。コンピュータって詳細を知らなくていいように設計しても、なんかの拍子に、やっぱり知ってないとわからないことが出てくるのよね。テクノロジー全般にそうだけど。 twitter.com/genkuroki/stat...

タグ：

posted at 10:50:23

#Julia言語

one(x) = xの型におけるかけ算の意味での1
oneunit(x) = xと同じ次元を持つ1

などが(他にもこの手の函数が沢山ある)、型の伝搬を適切に記述するために使われます。

この辺は結構ややこしくて、promoteの仕組みを理解していないとうまく使いこなせないことが多い。

タグ： Julia言語

posted at 10:42:05

#Julia言語しかし、函数に引数を与えて引数の型が確定したときのJITコンパイル時に函数内の変数の型が適切に確定しないと効率の悪いコードが生成されてしまいます。

だから、Juliaでは函数の引数の型が適切に伝搬するように書きます。

タグ： Julia言語

posted at 10:42:05

ごまふあざらし(GomahuAzaras @MathSorcerer

#Julia言語解説：Juliaでは函数の引数の型を一切指定せずに書くことが多く、引数の型を指定する場合にも、その函数の適用可能範囲を狭めないようにゆるく指定することが普通。

同一の函数をより多くの場面で使い回しできた方が便利。

タグ： Julia言語

posted at 10:42:05

ほー twitter.com/genkuroki/stat...

タグ：

posted at 10:34:22

#Julia言語ちょっと面白いことを発見。

nbviewer.jupyter.org/gist/genkuroki...
Julia v1.3.0-RC1での1,one,oneunit の比較

add1(x)=x+1
addone(x)=x+one(x)
addoneunit(x)=x+oneunit(x)

x=Float16(3.2)のとき、v1.3.0-RC1なら、これらのコンパイル結果llvmレベルで同じ。

しかし、v1.2.0では同じにならない！

タグ： Julia言語

posted at 10:32:51