黒木玄 Gen Kuroki
- いいね数 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 を見て下さい。
2017年11月21日(火)
研究不正があったら大学の評価を下げると言われたって、不正する人は大学の評価なんか気にしてないわけで、なんの対策にもならないわな。「不正を隠蔽したら評価を下げる」のほうが、まだしも効果あるかもねえ
タグ:
posted at 23:35:08
九九の唱え方が標準どおりでないと「ひどい」のだろうか?
( 標準どおりでない読み方を教えて一字一字までソックリにすることを指導していたら「ひどい」が、その元同僚(後輩)氏の場合、それはあるまい。)#掛算
タグ:
posted at 22:50:17
非公開
タグ:
posted at xx:xx:xx
#統計 WAICの大小によるモデル選択はWTの大小によるモデル選択と同じになると思って構いません(私が試した場合にはそうなるが一般には誤差が出る)。WTは予測誤差KLと正反対方向に揺らぎます。予測誤差が大きくなると、WTは小さくなり、モデル選択に失敗するかもしれない。
pic.twitter.com/ShNBFs2LfV
タグ: 統計
posted at 22:11:26
#統計
WT = WAIC - (真の分布を含む確率モデル共通の揺らぎ成分)
KL = (汎化損失GL) - (真の分布で決まる定数) = (予測分布の予測誤差)
pic.twitter.com/ShNBFs2LfV
タグ: 統計
posted at 22:07:16
@Essence_sci @UFOprofessor @kamo_hiroyasu 引用文献に出てくる、金田茂裕氏や浪川幸彦氏は超算数推奨者。こうやって超算数が継承されるのでしょうね。
タグ:
posted at 22:06:58
@Essence_sci @UFOprofessor @kamo_hiroyasu (1)4 メートルのひもと 8 メートルのひもをあわせると、何メートルの長さになりますか。
(2)20 cm の棒があります。その棒に、さらに 20 cm の棒 4 本をまっすぐにくっつけると、何 cm
になりますか。
これらだって、つなぎ目の長さが書いていないから回答不可能、ということだって可能。
タグ:
posted at 22:04:58
@Essence_sci @UFOprofessor @kamo_hiroyasu 「不合理問題」とされるものも、出題されたら非現実的と思っても計算すると思う。
タグ:
posted at 22:04:21
#統計 WAICの精度を正確に評価するために作ったグラフが既出の添付画像。渡辺澄夫著『ベイズ統計の理論と方法』p.119の定理15のβ=1の場合が再現されています。WAIC - T_true の揺らぎは予測分布の予測誤差を意味するKL情報量KLの揺らぎとちょうど反対向きになっています。
pic.twitter.com/ShNBFs2LfV
タグ: 統計
posted at 22:04:20
@Essence_sci @UFOprofessor @kamo_hiroyasu 問題も酷くて「4 個のボールと 8 個のボールがあります。かけると何個でしょう。」と聞かれたら、32個と答えるのは普通だと思う。意味だって、4個と8個から1個ずつ取り出す組み合わせの個数と解釈可能。
タグ:
posted at 22:02:31
@OpenGate_AL はじめましてです。夏頃に訪れた某県庁所在地のバスターミナルでこんな光景に絶句いたしました。これを設計した建築デザイナー氏や施設管理者氏に、弱視体験メガネ+白杖でココを歩いていただきたいものです。 pic.twitter.com/D2U3D8C8Gz
タグ:
posted at 21:58:06
#JuliaLang その辺のことについて私が最も参考になったのが、次のリンク先のJulia言語のJupyter notebookです。計算してその結果を見ながら、特殊函数の数値計算をするコード(連分数展開とテイラー展開を使う)を最適化して行く様子を見ることができます。
nbviewer.jupyter.org/github/steveng...
タグ: JuliaLang
posted at 21:34:28
#JuliaLang 具体的には、ぱぱっと計算してくれること(計算が速いこと)だけではなく、計算結果をすぐにプロットできることはとても大事。計算するためのコードを書くときには、計算結果を見ながらコードを積み重ねて行くので、計算結果をすぐに見れないと苦しい。
タグ: JuliaLang
posted at 21:30:02
#JuliaLang ユーザーインターフェースが気軽に使えるように設計されていること、ソースコードど計算結果を気軽にワンボタンで公開できるようになっていること、各種ライブラリを気軽に使用できること、十分な計算速度を確保できること(速いは正義!)、はどれも必須の条件だと思う。
タグ: JuliaLang
posted at 21:24:03
#統計 #JuliaLang 数値計算をぱぱっと行ってプロットした結果を見せてくれる人が社会の中に稠密に存在すれば、数学が関わる事柄について色々スムーズに行くようになると思う。
そして、数値計算をぱぱっとできるだけではなく、厳密な証明も説明できる人は貴重なので社会的に大事にされる必要がある。
posted at 21:20:27
@sekibunnteisuu @UFOprofessor @kamo_hiroyasu 示していただいた直前の一文「作文学力を伸ばすことが文章題学力を高めることに繋がると考える方が妥当であり、逆の方向での因果関係は考えにくい。」も酷いですよ。これと、「測定時期の前後」の話から「作文学力を伸ばすことによって算数文章題学力も向上する」と結論しているとかひどすぎる。
タグ:
posted at 21:19:12
#統計 #JuliaLang Julia言語を初めて使う人は
next.juliabox.com
を無料で使わせてもらうのがよいと思います。自分のパソコンに入れたい人は
nbviewer.jupyter.org/gist/genkuroki...
で公開している私によるインストールの記録を参考にして下さい。こういうことができる環境が普及するとよいと思う。
posted at 21:17:41
@sekibunnteisuu @UFOprofessor @kamo_hiroyasu 読んでみました。いやこれは酷い。誤字脱字の多さには目をつぶったとしても、酷い。要旨だと因果関係の認定まで踏み込んでいないことから推測するに、だめだと分かってて考察に入れているように見えました。
タグ:
posted at 21:16:36
#統計 #JuliaLang あと、私のような統計もど素人でかつ数値計算もど素人であっても、以上のような計算を気軽に実行して、グラフをプロットして、ソースコードも公開し、SNSで得られた情報を気軽に拡散できるわけです。そのための道具として、Jupyter notebook + Julia言語は極めて優秀だと思います。
posted at 21:16:08
#統計 添付画像は
WAIC-真の分布の経験損失
KL情報量KL=汎化損失GL-真の分布のShannon情報量
を比較するための散布図です。誤差は小さいですが、逆相関がはっきり見えています。これらの揺らぎはちょうど逆向きになります。(渡辺澄夫『ベイズ統計の理論と方法』p.119の定理15のβ=1の場合) pic.twitter.com/ShNBFs2LfV
タグ: 統計
posted at 21:11:02
#統計 添付画像はWAICと汎化損失GLの散布図です。汎化損失GLは予測誤差の正確な指標です。WAICはGLの推定値なのですが、大きくずれていることがわかり、WAICを使いたい人は心配になると思います。しかし、ずれの成分の大部分はモデル選択に影響を与えないずれです。 pic.twitter.com/lcseheLC7d
タグ: 統計
posted at 21:06:51
#統計 添付画像はLOOCVとWAICの散布図。ほぼ45度線に乗っています。すなわちそれらはほとんど等しい。LOOCVは工夫していない計算の仕方をしました。私が扱ったケースではWAICとLOOCVのモデル選択の結果はほぼ同じになりました。 pic.twitter.com/XGJSjjWWuN
タグ: 統計
posted at 21:04:58
#統計 #JuliaLang まだ確認していないのは、私の方法(MCMCはβ=1の1回だけしかせずに、数値積分は既存パッケージにまかせる方法)で数値計算した、2F(1)の近似値の誤差がどの程度であるか。F(1)の値を手計算で求められるケースで実験すればいいのですが、まだやっていない。誰かやってくれると嬉しい。
posted at 20:45:33
#統計 #JuliaLang 添付画像は上で紹介したJupyter notebookからの引用。
WBIC - FreeEnergy のプロット。サイズn=128のサンプル千個の各々について計算した結果です。WBICはFreeEnergy=2∫_0^1 E[β,H]dβより値が大きくなる傾向があるようです。
Jupyter notebookには他にもたくさんの情報がある。 pic.twitter.com/ieIWIRDgn5
posted at 20:43:11
#統計 #JuliaLang 計算には7時間程度かかる。
添付画像は私の方法で計算した FreeEnergy=2F(1) を横軸にWBICを縦軸に取った散布図です。サンプルサイズn=8で千個のサンプルについて2F(1)とWBICを計算しました。MCMCによる誤差と数値積分による誤差が入っています。しかし、ほぼ45度線に乗っている。 pic.twitter.com/WTnaeIJrJV
posted at 20:37:44
#統計 #JuliaLang そのシミュレーション結果を数表にしたり、グラフにしたりすることは以下のJupyter notebookでやっています。n=インプットの標準正規分布のサンプルサイズ
nbviewer.jupyter.org/gist/genkuroki...
n=8
nbviewer.jupyter.org/gist/genkuroki...
n=32
nbviewer.jupyter.org/gist/genkuroki...
n=128
posted at 20:32:43
#統計 #JuliaLang 私がやったことは、標準正規分布で生成したサンプルを次の3つのモデルで推定して、WAIC、LOOCV、汎化損失GL、WBIC、FreeEnergy=2F(1)などを求めることです。
normal=パラメーター数2の正規分布モデル
normal1=分散1の正規分布モデル
mixnormal=分散1の山が2個の混合正規分布モデル
posted at 20:29:56
#統計 #JuliaLang 続き。渡辺澄夫著『ベイズ統計の理論と方法』の記号で H(w) は nL_n(w) と書かれています。だから
nbviewer.jupyter.org/gist/genkuroki...
のIn [7]で E[β,H] を E2nLn と書いてあります。2倍を意味する2が入っているのは伝統的なAICとBICのスケールに合わせるためです。
posted at 20:27:20
#統計 続き。逆温度βごとにMCMCを走らせる必要があるかのような説明がググると普通に見付かるので私はびっくりしたわけです。私は何か勘違いしているのでしょうか?
私は以上の方針でF(1)=FreeEnergyとWBICを計算するコードを
nbviewer.jupyter.org/gist/genkuroki...
で公開しています。その In [7].
タグ: 統計
posted at 20:24:25
「水をかき混ぜ続けると熱湯へ変わる」と聞くと疑う人もいるだろうが、実際に成功させる動画。
そもそも科学者ジュールはこれで熱の仕事当量を導いた歴史的実験であり、理科の教科書にも載っているが忘れがち
youtu.be/GjcOobt9Ef8
タグ:
posted at 20:10:00
斉藤ひでみ・現職教師(西村祐二) @kimamanigo0815
ホワイト部活で有名な岐阜県多治見市の不都合な真実。
「部活動は全生徒加入制にしている」
校長よ、新聞で胸張って語るな。これは退職ものだ。
さらに調べたら、多治見=東濃の強制入部の割合が…なんと100%であることが!
多治見市よ、教員の勤務がホワイトだったらそれでいいのか。
拡散希望です。 pic.twitter.com/Adj3XlBVfz
タグ:
posted at 20:02:26
@togetter_jp あ、ちなみに医師免許持ってても、田中佳(脳外科医;ameblo.jp/emkanayoshi/en...)や
母里啓子(感染症専門だが臨床医ではない;www.hmv.co.jp/artist_%E6%AF%...)みたいに、『劇薬』の定義を知らなかった人もいるので、ご安心を。薬理学の授業受けてなかったのかな?この人たち
タグ:
posted at 17:58:06
@togetter_jp 薬剤に詳しい医療関係者なら、『劇薬=危険』という意味でないことは自明(www.fizz-di.jp/archives/10075...)。劇薬とは『過量に使用するときわめて危険性が高い医薬品』なので解熱剤でも劇薬となる。当然Seki_yoさん、知っててワクチンの不安を煽ってるんだよね?
twitter.com/seki_yo/status...
タグ:
posted at 17:53:10
@kikyokikyo519 #超算数 も実際はごく一部のぶっ飛んだ人がやっていると思ったら、市町村全体でそろえてやっているのを見て、戦慄しました
タグ: 超算数
posted at 17:15:46
非公開
タグ:
posted at xx:xx:xx
非公開
タグ:
posted at xx:xx:xx
あるいは、アメリカのとあるウラン鉱周辺で“ただ一度”だけ起きたケロシン漏れによる鳥類への被害を、世界中のウラン鉱でその産出量に応じて同様の被害が毎年起こるかのように計上しているとか、その他。 この指摘が正しいなら、あれはとても信用出来ない数字。
atomicinsights.com/nukes-kill-mor...
タグ:
posted at 10:46:45
あの数字には生物学者などからも異論が(風力の影響が低すぎと)出ている模様です。
一例→ “Sovacool’s (2009) estimate of the average number of birds killed per GWh of wind power is incorrect and omits a large body of easily accessible, published data.”
www.uregina.ca/science/biolog... twitter.com/scicafeshizuok...
タグ:
posted at 10:26:02
息子(6)が「猫ちゃんとは言うけど犬ちゃんとは言わんよねえ。ニャンちゃんとは言わんけどワンちゃんとは言うよねえ。逆だねぇ」と言っていた。わしはずっと気づかんかったし面白い。言語的現象に名前なんか説明があるんかもなー。
タグ:
posted at 06:58:26
@musicisthebest_ 数研出版に問い合わせして回答をもらいました。
8254.teacup.com/kakezannojunjo...
8254.teacup.com/kakezannojunjo...
タグ:
posted at 05:18:32
@Essence_sci @UFOprofessor @kamo_hiroyasu この論文も因果関係と相関関係の理解があやしい
glim-re.glim.gakushuin.ac.jp/bitstream/1095...
p90
【また、本研究の手続きにおいては、2 学期に作文学力を測定し、3 学期に文章題学力を測定した。時系列的にみても、先に測定した作文学力が後に測定した文章題学力を予測したとする、因果関係を想定することができる。】
タグ:
posted at 04:48:30
粘菌コンピュータというものがあり、日本地図の上に住んでいる人工の分だけ餌を盛り、あるバクテリアを放ちます。バクテリアは効率よく餌を得るため餌同士を繋ぐように繁殖します。その線が現在の日本の路線図とほぼ一致し、ずれているところは人類が間違えて建設している非合理的な輸送経路なのです。 pic.twitter.com/Gwjzv2SHA1
タグ:
posted at 00:37:00
非公開
タグ:
posted at xx:xx:xx
困ったな。ぼくは量子相転移、イジングモデルについてはガチの専門家なのだが、この注 1 の説明はまったく理解できない。ものすごく文章が下手な人が書いたのかもしれないが、キーワードを羅列したニセ科学系文書みたいに見えてしまう。
どこかにちゃんとした解説はないかな?
www.ntt.co.jp/news2016/1610/...
タグ:
posted at 00:15:23