算額あれこれ

Julia の RCall が動く環境

JuliaとRのバージョンの組み合わせで動かないことがある。現在のところ，以下の組み合わせでは問題ない。Julia:Version1.7.3R:Version4.2.1x86_64M1用のRだと動かない。JuliaのRCallが動く環境

2022/06/30 11:03

Julia: コマンドラインの特殊記号の入力法

sedを使って単語の置換を行い，結果をファイルに書き出そうと，run(`sed-e's/\t//g'foo.txt>out.txt`)としたら，LoadError:parsingcommand`sed-e's/\t//g'foo.txt>out.txt`:specialcharacters"#{}()[]<> &*?~;"mustbequotedincommandsとエラーを吐いた。一体どうやってquoteしたらよいか，ぐぐったら，run(`sed-e's/\t//g'foo.txt$(>)out.txt`)とすればよいのだった。Julia:コマンドラインの特殊記号の入力法

2022/06/27 15:47

Python のインストール

やれHomebrewだのpyenveだのcondaだのAnacondaだのMinicondaだの挙句の果てに，>Pythonを直接インストールする場合、ターミナルを利用せずにインストールできる点はメリットといえるでしょう。しかし、それ以外のメリットはありません。Homebrew+pyenvや、MinicondaのようにPythonのインストールだけでなく、今後のパッケージ管理の容易さなどを考えると、Pythonの直接インストールはあまりおすすめできません。しかし、「とにかくPython3系をインストールして使ってみたい」という方には、簡単にインストールするための手段として有効です。「それ以外のメリット」ってなんなん。「パッケージ管理」にあくせくする必要なんかないですよね。「Pythonの公式サイトにアクセ...Pythonのインストール

2022/06/26 21:00

数値から文字列への変換

数値から文字列への変換https://qiita.com/dzonesasaki/items/cc33ede78ca679f17f6dJuliaの例が記載されていなかったので...Pythonの場合>>>fVal=3.141592653589793>>>str(fVal)'3.141592653589793'>>>'{:01.8f}'.format(fVal)'3.14159265'>>>f'{fVal:01.8f}''3.14159265'Juliaではjulia>fVal=piπ=3.1415926535897...julia>string(fVal)"π"ulia>string(1pi)"3.141592653589793"julia>usingPrintfjulia>@printf("%.8F",...数値から文字列への変換

2022/06/17 21:35

R-4.2-branch 4.2.1 RC 久々に Status OK になった

https://mac.r-project.org/もう，2ヶ月近くか，StatusOKでなかったのだけど，今日見たら，（4.3.0も含め）x86_64用もarm64用もOKになっていた。R-4.2-branch4.2.1RC久々にStatusOKになった

2022/06/16 22:00

じゃあ!四捨五入問題，どうすりゃいいのよ！！（結論：学校で習った四捨五入は忘れましょう！！）

まあ，説明するのも面倒なのですが，学校で習った四捨五入は「不適切な方法」なので，忘れましょう。JISZ8401（数値の丸め方）（PDF、75KB）https://www.erca.go.jp/fukakin/before/pdf/suuchi.pdfあなたが思ったのと違った結果になっても，それは，「あなた（学校教育あるいは文部科学省）が間違っている（不適切だ）から」なのです（残念！！！）。じゃあ!四捨五入問題，どうすりゃいいのよ！！（結論：学校で習った四捨五入は忘れましょう！！）

2022/06/13 22:35

R, Python, Julia での round の挙動（結論：同じだ）

Rのroundのexample>example(round)round>round(.5+-2:4)#IEEE/IECrounding:-2002244[1]-2002244round>##(thisis*good*behaviour--do*NOT*reportitasbug!)round>round>(x1<-seq(-2,4,by=.5))[1]-2.0-1.5-1.0-0.50.00.51.01.52.02.53.03.54.0round>round(x1)#--IEEE/IECrounding![1]-2-2-10001222344round>x1[trunc(x1)!=floor(x1)][1]-1.5-0.5round>x1[round(x1)!=floor(x1+.5)][1]-1.50....R,Python,Juliaでのroundの挙動（結論：同じだ）

2022/06/13 15:34

round() は四捨五入ではないよという話（何回目?）

Python公式ドキュメントに>注釈浮動小数点数に対するround()の振る舞いは意外なものかもしれません:例えば、round(2.675,2)は予想通りの2.68ではなく2.67を与えます。これはバグではありません:これはほとんどの小数が浮動小数点数で正確に表せないことの結果です。詳しくは浮動小数点演算、その問題と制限を参照してください。と書いてあるそうだ。この説明は正確ではない。2.5も3.5も「浮動小数点数で正確に表せる」が，roundの結果は異なる。>>>round(2.5)2>>>round(3.5)42.5は16進表現で(2.8)16で，10進数に変換するのは2*16^0+8*16^(-1)=2+0.5である。3.5は16進表現で(3.8)16で，10進数に変換するのは3*16^0+8*16^...round()は四捨五入ではないよという話（何回目?）

2022/06/13 13:37

AS 177（正規所数値の期待値） を Julia に翻訳

「正規所数値の期待値」とは聞き慣れない訳語かもしれないが，Expectedvaluesofnormalorderstatisticsのこと。H.LeonHarter（1961）ExpectedValuesofNormalOrderStatistics,Biometrika,48,151-165.https://www.jstor.org/stable/2333139大本に近いのは以下AlgorithmAS177:ExpectedNormalOrderStatistics(ExactandApproximate)http://imaging.mrc-cbu.cam.ac.uk/statswiki/FAQ/Corder?action=AttachFile&do=get&target=royston.pdfFO...AS177（正規所数値の期待値）をJuliaに翻訳

2022/06/12 20:30

ダメ出し：日付の等差数列 (*^^*)

「2020年1月から2022年5月まで、連続の月次データを生成する」ということで，持って回ったプログラムが紹介されている。library(tidyverse)library(lubridate)year_month<-seq(as.Date("2020-01-01"),as.Date("2022-05-01"),by="1day")%>%floor_date(unit="months")%>%unique()seq(as.Date("2020-01-01"),as.Date("2022-05-01"),by="1month")で十分だし，seq(as.Date("2020-01-01"),as.Date("2022-05-01"),by="1week")seq(as.Date("2020-01-01"),...ダメ出し：日付の等差数列(*^^*)

2022/06/05 21:36