読み込んだファイルをお手軽ランダムシャッフル． x0 <- read.table("Forward_Number_Sherlock.txt") x <- x0[,1] set.seed(10) y <- sample(x) write.table(y,"Shuffle10_Number_Sherlock.txt",append=FALSE,row.names=FALSE,col.names=FALSE,quote=FALSE)…

レビューにまわっていた論文がrejectされた． ひとりのレビュアーからは厳しい言葉ばかりの回答だった．「an elephant in the room」とまで皮肉を言われた． これは著者が7人もいるのに，重要な問題を見逃しているという手痛い指摘で，こちらが恥ずかしくな…

軸名をテキスト形式にした，バーを作りたい場合，以下のようなwaveを作って， User Ticks from WavesでLabelにテキスト，Locationsにwave1を入れる． textWave1 wave1 Mon. 0.5 Tue. 1.5 Wed. 2.5 Thu. 3.5 Fri. 4.5 Sat. 5.5 Sun. 6.5

単位根検定 library(tseries) a <- PP.test(x) b <- adf.test(x) ave <- round(mean(x),digits=5) std <- round(sd(x),digits=5) a2 <- PP.test(z) b2 <- adf.test(z) ave2 <- round(mean(z),digits=5) std2 <- round(sd(z),digits=5) out <- paste(commandA…

雑記．大阪市長が従軍慰安婦(comfort woman)は申し訳ないけど必要だったとか，風俗もっと利用しろ言っている件で． 彼がが風俗好きなのは分かった．弁護士時代からの仕事で風俗業の人と付き合いがあって，彼らが「ちゃんと」仕事しているのを認めたいのも分…

空飛ぶ広報室／有川浩テレビドラマが面白かったので，一足先に原作を読んだ．ドラマより恋愛要素は少ないものの，そういえばこういう若い恋愛も絡んだ物語も久しぶり．空幕広報室に一人のりこむ，力の入った稲ぴょんが2006年頃の自分とダブって見えて仕方が…

匿名ダイアリーに投稿したネタがちょっとだけバズっていた． 正確に言うと，私がバズっていたネタにコメントして，その私のコメントを引用する形で，他の人ネタがバズっているという．不思議な経験だった．とりあえず，人目に触れるのは大事なのは分かった．…

対数の底がnの計算をするには，log_n(X)=log(X)/log(n)を使って計算すればよい．IGOR Pro リファレンス