「ぎょーむ日誌」目次に戻る
|
KuboWeb top に戻る
|
twilog
|
atom
ぎょーむ日誌 2017-06-22
苦情・お叱りは, たいへんお手数かけて恐縮ですが, 久保 (
kubo@ees.hokudai.ac.jp
) までお知らせください.
本日
(
kubolog20170622
) |
次の日
|
1 日前
|
7 日前
|
31 日前
|
365 日前
|
top
2017 年 06 月 22 日 (木)
0610 起床. 朝飯の準備その他ばたばた. 0815 保育園着. 0845 研究室着.
院生が「なんで3個体のときは個体差のランダム効果を推定できないんですか?」とつめよってくる…えーと,個体差が平均ゼロの正規分布にしたがう場合,その標準偏差の推定がむちゃくちゃになるから,とご説明する作図.院生と教員の役割が逆転…?いやいや,教員もこうやって鍛えられるべきかも…
pic.twitter.com/9zqQB5WvSh
10:13:21
RT
@berobero11
: StanのWikiにGelman先生による事前分布のオススメの選び方(Prior Choice Recommendations)のページができました。重要です。
github.com/stan-dev/stan/…
10:21:05
さてさて…階層モデルにおける「個体数が少ないときには,個体差のばらつきの推定が難しいかも」問題は「個体差のばらつき」の事前分布の設定とかにも関係するわけですが…このあたり Gelman 先生の意見
github.com/stan-dev/stan/…
を参考にしますか…というかんじで…
10:26:49
RT
@berobero11
: 階層事前分布の標準偏差の事前分布には半コーシー分布がよさげ(2006)、とのことですが、最近ではそれでは弱すぎる、自由度4ぐらいの半t分布がよい、とGelman先生が以下のページで書いてます。ふぁぼるだけじゃなくて読んで下さい>_<
twitter.com/berobero11/sta…
10:28:36
RT
@KuboBook
: 「グループ」の個数が 3, 4, 5 のときは半コーシー分布がよい…知らなかった…うむ,勉強しよう.さすがに,いまどき逆ガンマな事前分布を使ってるヒトは…じつは,いないわけでもなかったり…
www.slideshare.net/hoxo_m/ss-5941…
10:29:10
個体数=3 の「個体差の標準偏差」,その事前分布は半コーシーとかもっと保守的なやつとか,いろいろ議論あるけど…「真の困難」は環境科学な院生の多くが JAGS,Stan 利用を嫌がっていて,階層事前分布とか自由に指定できない…という点にあるかも.みなさん,lme4 とかもうヤメて〜
10:37:22
GLMM 推定関数では「古参」の glmmML,Broström さんは今もきちんと保守していて…しかも階層事前分布が gaussian, logistic, cauchy から選べます…これは SD の分布ではなく「個体差」(とか地域差)の分布の指定ですが…
pic.twitter.com/vG1OeUhyUA
10:46:41
R の library(lme4) は好きではないのだが,近年 Bolker さんたちががんばって(Bates さんは Julia に移動?)保守・改訂したりドキュメントを整理してくれたりしているわけですが…
cran.r-project.org/web/packages/l…
10:54:46
RT
@KuboBook
: R の library(lme4) の論文が出た.Bolker 先生など著者が増えたおかげで,かつてなくわかりやすいような…
twitter.com/bolkerb/status…
10:54:59
【音声配信】「毒性のヒアリ、その生態と対策。そして、特定外来生物との付き合い方」五箇公一×勝田吉彰…
www.tbsradio.jp/158118
昨晩の放送,いまごろ聴いていますが…五箇さんいわく,大田区・横浜市で侵入してきたアルゼンチンアリを根絶,世界初の事例とのこと.
11:13:13
五箇さん,10年以内に国内のアルゼンチンアリを根絶させる…とのことなので,根絶される前にさっさとアルゼンチンアリ vs 在来アリ群集のデータを解析して論文にしなくては…
pic.twitter.com/rH5EIWRcRa
11:19:59
五箇さんによると,アリの専門家でもヒアリと在来アリの区別は難しいとのこと. "ヒアリ"はほかのアリとどう違う?東京都が見分け方を公開 - LINE NEWS
news.line.me/issue/science-…
#linenews
@news_line_me
さんから
11:28:04
RT
@Anthunter1212
: これから神戸で「ヒアリ等対策会議」に。等がつきましたが、これはこれで、外来種対策としては正しいかと思います
11:35:18
RT
@konotarogomame
: 文科省の窓口にぜひ、ご連絡を入れてください。手ぐすねひいて待ってますから。
twitter.com/kikumaco/statu…
12:16:33
平均ゼロ・標準偏差 (SD) 1 の正規乱数を使って,サンプルサイズ n が小さいときに SD がばらつきますよ図示,n = 3 (左) では左右非対称になってるな…と今さら気づいたのですが,これは「不偏ではない」推定値だから (n-1 ではなくn で割ってるから) ですかね…
pic.twitter.com/6OE23JdNJo
12:29:19
@ibaibabaibai
はい,そうでした…
15:52:29
標準正規乱数の標準偏差 (SD) の分布ってカイ二乗分布
ja.wikipedia.org/wiki/%E3%82%AB…
でしょ…とご指摘いただき「あ,そうか」状態…
16:09:11
いちおう可視化というか作図…サンプルサイズ 3 のときに 3 で割る(上の図)と 2 で割る(下の図)…まあ,後者のモードは 1 に少し近づきますが…まあ非対称.乱数生成時の SD / 平均値が小さければ,サンプルサイズが小さいときでも左右対照になるんでしょうね.
pic.twitter.com/pfZdmSuBMT
16:13:54
Methods in Ecology and Evolution ではいまだに Akaike Weight 議論が続いてるもよう… A farewell to the sum of Akaike weights: …
onlinelibrary.wiley.com/doi/10.1111/20…
16:26:57
RT
@willyoes
: 米国大学院のRA制度とほぼ同じ。学費の一部または全部免除があればなお良しなのだが。
www.nikkan.co.jp/articles/view/…
16:36:06
1710 保育園着. とらを散髪につれていく. そのあと, かねてより「行きたい」とくりかえし要望していたパスタ屋へ. 帰宅して風呂. 2140 とら寝かしつけ ……
本日
(
kubolog20170622
) |
次の日
|
1 日前
|
7 日前
|
31 日前
|
365 日前
|
top
KuboLog
|
KuboWeb