黒木玄 Gen Kuroki(@genkuroki)/2021年10月24日

@n_hroko @genkuroki 英語圏には教育が分権化しているところが多いとは聞いています。カリキュラムだけではなく、おそらく個々の教員の裁量の重みも違うのではなかと想像しました。

タグ：

posted at 23:34:25

解答略 @kaitou_ryaku

天むす名古屋 Temmus @temmusu_n

インストールで詰まった人はこれを読むんや〜 twitter.com/kaitou_ryaku/s...

タグ：

posted at 23:29:55

#超算数現地校と補習校の方針の違いで困っている保護者の方がいるんだね。割って入っていって邪魔扱いされている人が観測されたtwitter.com/flute23432/sta...。

タグ：超算数

posted at 23:05:49

積分定数 @sekibunnteisuu

#超算数　「ハジキ、クモワじゃなく図で教える」

二重数直線のことなら、ろくでもない。 twitter.com/TeacherhaGreat...

タグ：超算数

posted at 22:44:37

©︎waku @wiimuuwooooo

これを超える写真に出会えない。 pic.twitter.com/xthm7Lu3W5

タグ：

posted at 22:31:37

解答略 @kaitou_ryaku

みんな笑うだろうけど、僕はこの未来が本当に来るんじゃないかと思ってるぜ
twitter.com/kaitou_ryaku/s...

タグ：

posted at 22:24:04

解答略 @kaitou_ryaku

Windows11をインストールして、PowerShellでwsl --install -d Ubuntuすると、全自動でLinuxが立ち上がった。標準ファイラーのNautilusやxeyesも動いた

ここまでくると完全なLinux DesktopとしてWindowsを使えちゃうな。今後、Windowsのネイティブアプリはどんどん消滅し、全てがLinuxになる気がするぜ pic.twitter.com/wW3v2ZJccs

タグ：

posted at 22:22:27

非公開

タグ：

posted at xx:xx:xx

Hiroyasu Kamo @kamo_hiroyasu

そういえば、書名が『幾何学大辞典』であって『幾何学大事典』でないのはなぜなんでしょう？

タグ：

posted at 21:04:00

Hiroyasu Kamo @kamo_hiroyasu

比較できる類書がありません。これは初等幾何について事典です。この分野で他に事典と呼べる網羅性のある本は他に見たことがありません。
#マシュマロを投げ合おう
 marshmallow-qa.com/messages/4107a... pic.twitter.com/pgGYOCe6f2

タグ：マシュマロを投げ合おう

posted at 21:03:00

非公開

タグ：

posted at xx:xx:xx

@hasundue @KB_satou #Julia言語あと eachcol を使うこともできます。 pic.twitter.com/CzuIsAO73G

タグ： Julia言語

posted at 20:03:18

清　史弘 @f_sei

【即興曲】
この前のスペースのときに弾いた即興演奏です。
どうぞお聴き下さい。 pic.twitter.com/9K4tPsw3V1

タグ：

posted at 19:55:20

@KB_satou @hasundue #Julia言語

成分ごとの作用はdot syntaxによるブロードキャストまたはmapで、行列の列または行ごとの作用にはmapslicesを使えます。(もっと一般的なことをできる。) pic.twitter.com/dZwgmo0BtD

タグ： Julia言語

posted at 19:35:40

dubhe @dubhejp

元はこんなこです。これは一年前の写真ですが。。。孤高のねこで、滅多にニャーとも言いません。媚びません。 pic.twitter.com/F79cmLgjva

タグ：

posted at 19:24:58

避雷 @lucknknock

GPUでグラフを扱えるようになったので物流ネットワークみたいな表現が出来るようになりました！GPUでグラフ上の近傍探索を並列化しています…！ pic.twitter.com/Ni6T0eKOTN

タグ：

posted at 19:06:24

兎です。(FAKE) @Soviet_Usako

この円グラフ好き pic.twitter.com/vcUCi7oIUs

タグ：

posted at 19:04:45

dubhe @dubhejp

ねこを可愛いく撮影するってとてもむずかしいな。。。
なんで皆さんのねこは可愛く撮れるのだろう。 pic.twitter.com/6p7ijXM83u

タグ：

posted at 18:52:19

せきゅーん @integers_blog

@tsujimotter そして、1.5によればそれは(Z/pZ)^*の勝手な原始根 or それ+pのいずれかは満たすことが示されていると思います。なので、non-Wieferichであれば原始根を共有できるということですかね。

タグ：

posted at 18:24:38

せきゅーん @integers_blog

@tsujimotter integers.hatenablog.com/entry/2016/07/...
の1.6によればg\not\equiv 1\pmod{p^2}なmodpの原始根gは(Z/p^kZ)^*の原始根にもなることが示されているように思います。

タグ：

posted at 18:22:16

せきゅーん @integers_blog

@tsujimotter つまり、(Z/pZ)^*の原始根と(Z/p^kZ)^*の原始根を共有できるかという問題でしょうか・

タグ：

posted at 18:19:32

#Julia言語よりすすんだことについては

youtu.be/wXRMwJdEjX4
Package development: improving engineering quality & latency | Workshop | JuliaCon 2021

タグ： Julia言語

posted at 18:16:04

#Julia言語 Meta.@ dump の結果は見難いので、Meta. show_sexpr(:(コード)) でS式っぽい表示の方を見た方が分かりやすいことが多い。

もしくは私の MetaUtils.jl を使う。
github.com/genkuroki/Meta...

Juliaの式への理解はメタプログラミングで必須。

タグ： Julia言語

posted at 18:12:35

#Julia言語が裏で何をやっているかを積極的に見るようにすると理解が進み易い。

:( ) で囲む → パーサーの挙動の確認
Meta.@ dump → パースの結果どのようなASTが得られるか
@ code_warntype → 型伝搬の様子
@ time → 計算時間とメモリ割当
@ code_llvm, @ code_native → コンパイル結果

など

タグ： Julia言語

posted at 18:09:05

tsujimotter 日曜数学者 @tsujimotter

#Julia言語型安定性の確保とメモリ割当の節約がJuliaでの最適化のイロハのイなのですが、その辺に関する具体的な処方箋については

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

@ code_warntype による型伝搬の様子の観察と、@ time によるメモリ割当のモニターが基本中の基本。

タグ： Julia言語

posted at 18:03:51

昨日の #日曜数学会の発表にインスパイアされて、循環小数を調べていたら面白い定理を発見・証明できました！
（多分知られていない性質だと思うので）ぜひご覧ください！

（自由研究）1/p^k型循環小数のフルサイクル性について - tsujimotterのノートブック
tsujimotter.hatenablog.com/entry/full-rep...

タグ：日曜数学会

posted at 17:56:28

実際の使用時の多くの場面では、#Julia言語の最適化では重要だとされる「型安定性」を無視したコードを気楽に書いても動くことが非常に便利です。

型安定性が重要な場面での基本は、函数の @ code_warntype や @ code_typed の結果を何はともあれ見てみること。これは基本中の基本。

タグ： Julia言語

posted at 17:56:07

#Julia言語【型推論の仕組みがJulia本体にあるはずだし、それを使ったLinter的なものがあったら最高】と言われると、Julia界の新ヒーローになっているaviateskさんのJET.jlおよび関連のツールの試みが思い浮かびます。我々日本人は応援できたら良いと思います。

github.com/aviatesk/JET.jl twitter.com/meows_whiskers...

タグ： Julia言語

posted at 17:48:57

S (ﾂｲｰﾄはｽﾚｯﾄﾞ全体をご確認く @esumii

::Float64とかcollectはなくても大丈夫。 #Julia言語

-3.0:0.1:3.0はそのままで縦ベクトル扱いされるAbstractVectorです。

引数の型名を書いていない函数の方が適用範囲が広く、結果的に思わぬバグの発生を防ぎ易くなる。::Float64と書いた途端に自動微分などの便利なパッケージを適用不可能に。 twitter.com/sophia8wendel/...

タグ： Julia言語

posted at 17:33:44

あれ、GmailのiPhoneアプリがアップデートで、To:やCc:先をドラッグ&ドロップで移動できなくなった？　むー…

タグ：

posted at 17:32:41

#統計

 github.com/genkuroki/publ...
2×2の分割表のオッズ比の信頼区間のシンプルな #Julia言語での実装例

各種のP値函数と信頼区間函数を数行で書けます‼️

どれも数行なのでコードを読む気になれる分量。

χ²検定に付随する信頼区間の背景にはLagrangeの未定乗数法を使った最尤法の計算が隠れている。 pic.twitter.com/qqHdvtfDPv

タグ： Julia言語統計

posted at 17:17:11

#統計だから、データのオッズ比の対数を見てパラメータのオッズ比の信頼区間を求める方法では、正規分布を仮定していません。正規分布による近似は結論です。

しかし、データのオッズ比の対数が近似的に正規分布に従うことを使う方法は誤差が大きいので、使わない方が良いと思います。

タグ：統計

posted at 17:17:08

#統計 2×2の分割表のモデル化としてよく使われているのは

* 4つのPoisson分布
* 4項分布
* 2つの2項分布
* 超幾何分布(独立性(パラメータのオッズ比=1)を仮定しない場合はFisherの非心超幾何分布)

の4種類です。データのオッズ比の対数が近似的に正規分布に従うことはこれらのモデルの帰結。続く

タグ：統計

posted at 17:17:07

#統計

信頼区間=P値がα以上になるパラメータの範囲

なのでオッズ比(これもモデルのパラメータ)の信頼区間の定義はP値函数の定義の種類だけ沢山あります。

対数オッズ比がサンプルサイズ→∞で漸近的に正規分布に従うことを利用する信頼区間の求め方もありますが、正規分布性はモデルの結論。続く twitter.com/kota222studyin...

タグ：統計

posted at 17:17:06

kashi @mkashi

いざ作ってみると微妙だな。もうちょっと見やすく調整すればよくなるのか？

タグ：

posted at 16:26:06

kashi @mkashi

ソースコードを書いてボタンを押すと、syntax highlightされたコードのページに飛ばされるが、そのページをtwitterに貼るとtwitter cardの機能で画像が見える、はず。
kashi.zenno.info/exw/twicode/d/...

タグ：

posted at 16:24:29

kashi @mkashi

ごまふあざらし(GomahuAzaras @MathSorcerer

twitter.com/mkashi/status/...
これを試しに作ってみた。ただしたぶんすぐ消す。
kashi.zenno.info/exw/twicode/tw...

タグ：

posted at 16:21:47

壁ドンするとな円周率わかるんやで
緑玉が質量１００で
紫が質量１としている
#Julia言語 pic.twitter.com/ouJdD1EUZY

タグ： Julia言語

posted at 16:19:45

非公開

タグ：

posted at xx:xx:xx

竜生九子博物館 @kyouinoheya

これはこれはッ！！！実物が観れるとは！！ pic.twitter.com/vjZ5SXvyN4

タグ：

posted at 13:31:48

#統計 2×2の分割表のχ²検定を「超幾何分布の正規分布近似」によって理解したつもりになるのは危ないです。

それだと、2×2の分割表のχ²検定が、4つのPoisson分布、多項分布、2つの二項分布において相当に正確になることを理解できなくなり、Yatesさんと同様の誤解に陥る可能性が出て来ます。

タグ：統計

posted at 13:28:02

Finbergさんの指摘

* 条件付き確率分布として得られる超幾何分布ではなく、マージンを固定しない多項分布や2つの二項分布においては、無補正のPearsonのχ²統計量の分布は自由度1のχ²分布にかなり近くなる。

* その場合にはYatesさんの主張とは違ってχ²検定はconditional testにはなっていない。

タグ：

posted at 13:23:19

elastic @sophia8ds

勾配降下法の更新プロセスのグラフ
これブログになかったので自作した💪
#Julia言語 pic.twitter.com/QndOYTZNPU

タグ： Julia言語

posted at 13:23:09

#統計 Yatesさんの「分割表のχ²検定はFisher検定と同様に条件付き確率分布を基礎とする検定である」という基本前提が単純な誤解に過ぎないことは、私が指摘するまでもなく、同論文に付属の

rss.onlinelibrary.wiley.com/doi/pdfdirect/...
Discussion of Dr Yates's Paper

でもしっかり指摘されています。
↓ pic.twitter.com/BXVeDz2QOW

タグ：統計

posted at 13:18:07

アヒル93 @SomeKindOfDuck2

Turingで複数のチェインを並列で走らせることができた。自分の環境ではconfig.fishに`set -x JULIA_NUM_THREADS = 4`を追加したら問題なく動いた。ちゃんと`plot`でそれぞれのチェインの結果を重ねて出してくれるしもう文句ないわな。
discourse.julialang.org/t/multi-thread... pic.twitter.com/ootaEesKSH

タグ：

posted at 13:15:54

#統計分割表のχ²検定がFisher検定と同様に条件付き確率分布を基礎とする検定(conditional test)だと誤解している偉い人がいるので要注意。Yates補正のYatesさんはもろにそういう誤解をしていました。

rss.onlinelibrary.wiley.com/doi/pdfdirect/...
F. Yates, Tests of significance for 2× 2 contingency tables, 1984
↓ pic.twitter.com/lLCu4fyoB3

タグ：統計

posted at 12:47:24

uep @uepuep

昭和基地の夜。星がほぼ真上で回っています。天の川、マゼラン星雲とか見えます。
#南極 pic.twitter.com/9Lncd8BMGo

タグ：南極

posted at 12:46:35

少し雑な言い方になるけど，こういうツイート群に見られる「安直さ」が「論理国語」の方向性の危うさの一つだと当初から思っている．それは「文学」と「論理」を分けることと同じくらい大きい問題なのではないか．

タグ：

posted at 12:45:27

新井氏は記事中の但し書きを引用していない．（海外の臨床試験の結果や腎臓や肝臓機能が低下している人への使用注意．）当初の方針や今の方針の正しさや適切さについて論じるには様々なデータを取り寄せたり，方針決定の経緯についての情報をみて判断しなければならない．この記事だけではそれは無理． pic.twitter.com/64Z8xthrs7

タグ：

posted at 12:28:58

#統計縦横の合計が全固定されているという非現実的な制限による条件付き確率分布(独立性の帰無仮説の下で超幾何分布になる)で、Pearsonのχ²統計量を使った近似をするとひどく不正確になっていても、それを足し上げて元の分布に戻れば十分正確になる場合が多い。

タグ：統計

posted at 12:28:09

#統計 2×2の分割表の超幾何分布でのPearsonのχ²統計量の分布を足し上げて、2つの二項分布でのPearsonのχ²統計量の分布を作る計算の様子のアニメーションが以下のリンク先にある。

足し上げる前のt=1はひどく不正確だが、足し上げると正確になって行く。 twitter.com/genkuroki/stat...

タグ：統計

posted at 12:24:35

しかし新井氏のツイートは「答え」がないので，新井氏が「論理的に自分の意見を持てるようになって欲しい」とか「身につけたい国語力」として何を想定しているのかがわからない．（この記事だけから，正しさや適切さについての判断を行うことは，それ自体”論理的”ではないと私は思う．） pic.twitter.com/SCLGYeo7kU

タグ：

posted at 12:23:26

朝日の記事「（新型コロナ）治療薬、増える選択肢　他の病気向けの薬、転用」の記事を素材に新井紀子氏が「こういう解説文を正確に読んだり書いたりする力は，ぜひとも身につけたい国語力」「学校という場所で指導すべき」と述べ，”モデル問題”的なツイートを2つ付けている．
digital.asahi.com/articles/DA3S1... pic.twitter.com/02zA5Bhf1L

タグ：

posted at 12:19:16

こーた丼 @rental_math

手計算のアルゴリズムを、そのまま実装しても速度が落ちないからJulia言語は強いですね🥰 twitter.com/genkuroki/stat...

タグ：

posted at 12:17:15

個人的な印象では、一般に、手でする計算にも強い人の方が、コンピュータによる計算も上手にしている感じ。 twitter.com/genkuroki/stat...

タグ：

posted at 12:16:09

#統計より最近のFisher検定と独立性のχ²検定に関する解説
↓ twitter.com/genkuroki/stat...

タグ：統計

posted at 12:12:54

#統計より最近のFisher検定と独立性のχ²検定に関する解説
↓ twitter.com/genkuroki/stat...

タグ：統計

posted at 12:12:53

#統計真に科学的であろうとするために最も重要なことは、権威を根拠にしないことです。権威を根拠にしないことは科学の出発点です。

教科書や講義や論文でみんなそう言っているから〇〇は正しいと主張することは最も非科学的な行為だと思います。

タグ：統計

posted at 11:57:02

#統計第一の過誤が起こる確率はコンピュータを使ったモンテカルロシミュレーションで容易に計算できます。

教科書や大学の講義では「小サンプルでは不正確なカイ二乗検定ではなく、正確なFisher検定を使おう」となっていても、自分でコンピュータをまわせば、誰でもおかしいと気付けます。

タグ：統計

posted at 11:54:35

#統計本当はFisher's exact testは小サンプルでびっくりするほど正確ではなく、カイ二乗検定の方が第一種の過誤の確率と有意水準がより近いという意味で正確な場合が多いので、小サンプルの場合にトレードオフの問題にせずに一方的にFisher検定を勧める伝統的流儀は間違っています。

タグ：統計

posted at 11:50:55

#統計 α以下の値達の荷重平均はαよりずっと小さくなりがちという当たり前のことがFisher検定の検出力を下げている。特に小サンプルでその問題はひどくなります。

しかし、これはトレードオフの問題で、Fisher検定では第一の過誤の確率を確実に名目有意水準α以下にできます。

タグ：統計

posted at 11:47:36

#統計しかし、すべての条件付き確率分布でP値がα以下になる確率をα以下にしたとき、元の制限されていない確率分布でのP値がα以下になる確率はα以下の値の荷重平均として計算され、αよりずっと小さな値になりがちになります。

これがFisher検定が大きく不正確になりがちな真の理由です。

タグ：統計

posted at 11:44:03

#統計続き。これによって、Fisher検定では、第一種の過誤が起こる確率が確実に名目有意水準α以下になるようにしています。続く

タグ：統計

posted at 11:44:02

#統計続き。③のP値の構成では、任意の有意水準αについて、条件付き確率分布内でP値がα以下になる確率がα以下になるようにします。すべての条件付き確率分布でP値がα以下になる確率がα以下になるようにしているので、元の確率分布(そこでは縦横の合計は全固定されていない)でも同様のことが成立。続く

タグ：統計

posted at 11:44:02

#統計このスレッドの上の方での説明の仕方は以下のリンク先の意味で良くないです。

Fisher検定：

①縦横の合計をデータと同じに固定す流ことによって条件付き確率分布に移る。

②独立性の帰無仮説の元でその条件付き確率分布は超幾何分布になる。

③その条件付き確率分布でP値を構成。 twitter.com/genkuroki/stat...

タグ：統計

posted at 11:32:12

@kota222studying #統計その辺の私の書き方を今読むと、Fisher検定が条件付き確率分布を使ってP値を構成する十分に面白い方法であることに十分な敬意を払っていない点が非常に良くないと自分を非難したくなりました。カイ2乗検定の側の説明はフェアだと思います。自己引用RTに続けるつもり。

タグ：統計

posted at 11:28:49

K.B.砂糖 @KB_satou

dot syntax も map もよく理解してないことがわかった。(特に dot syntax)

タグ：

posted at 11:03:18

K.B.砂糖 @KB_satou

A を計算したい縦ベクトルを横に並べた行列とすると
A' * A で内積の計算自体はできるけど余分な計算しちゃう
ちょうど (A' * A) の対角要素だけをベクトルとして出力してくれるとちょうどいい twitter.com/KB_satou/statu...

タグ：

posted at 11:01:15

K.B.砂糖 @KB_satou

Julia で行列の行(列) ごとにノルムをとったりしたいんだけど for 文使わずにかける?

タグ：

posted at 10:51:08

elastic @sophia8ds

勾配のグラフ
quiverはなぜかsurfの時と同じようにやると変なグラフになる🤔
X軸とY軸はMatrixにしておいた方が無難? pic.twitter.com/Y0UMnoLkry

タグ：

posted at 10:02:43

大須賀覚 @SatoruO