R による離散選択モデルの推定方法メモ 平成 21 年1月 26 日 東京海洋大学 兵藤 哲朗 1 目 目 目 目 次 次 次 次 1. 1. 1. 1. 『 『 『 『LOGIT.FOR LOGIT.FOR LOGIT.FOR LOGIT.FOR』 』 』 』と と と と四半世紀 四半世紀 四半世紀 四半世紀................................ ................................ ................................ ................................................................ ................................ ................................ .............
RでMarkov chain Monte Carlo(MCMC、マルコフ連鎖モンテカルロ)を行うには大きく2つの方法があります。 1つ目はパッケージR2WinBUGSを使って、Rの後ろでWinBugsを走らせる方法。 もう1つはパッケージMCMCpackを使って直接R上で計算する方法。 ちなみにMCMCの利点は、混合効果モデルのように複雑なモデリングになった場合の複雑なパラメータ推定を安定にさせるというもの。 ロジスティック回帰で交互作用がいっぱいのときも安定するらしい。 それではパッケージの説明です。 MCMCpack #MCMCpack library(MCMCpack) y <- c(1,3,3,3,5) x1 <- c(1,2,3,5,5) x2 <- c(2,3,1,5,2) x <- matrix(c(x1,x2),ncol=2,nrow=5) result3.sim <- M
今回は多変量解析についてです。その前にそもそもですが、「多変量解析」という言葉は様々な意味で使えるので、なるべく使うのを止めましょう。私が経験してきた中で、このような意味で使われていました。重回帰、一般線形モデル一般化線形モデル変数選択(ステップワイズ法)変数縮小(主成分分析) どの手法も目的がまったく違っています。「多変量解析をやりたいのですが、、、」と相談されると、こちらとしては「多変量解析」が何を意味するのかを探るところから始めます。 具体的には、解析手法はこのように使い分けます。何かの結果変数を説明するモデルを作る→重回帰同じ目的で説明変数が連続値以外→一般化線形モデル(GLIM、ぐりむと発音) 実は「重回帰」も「GLIM」もほとんど同じ意味ですが、ニュアンスとして重回帰は一般線形モデル(GLM、じーえるえむ)を指す事が多いです。正確には「重回帰」は「単回帰」と対になる言葉で、説明
13. 基本コード オプション 効果量 eta・・・「eta = T」とすることで,分散分析表にイータ二乗(eta squared)を追加 peta・・・「peta = T」とすることで,分散分析表に偏イータ二乗(partial eta squared)を追加 geta・・・「geta = T」とすることで,分散分析表に一般化イータ二乗(generalized eta squared)を追加eps・・・ 「eps = T」とすることで,分散分析表にイプシロン二乗(epsilon squared)を追加 peps・・・「peps = T」とすることで,分散分析表に偏イプシロン二乗(partial epsilon squared)を追加 gpes・・・「gpes = T」とすることで,分散分析表に一般化イプシロン二乗(generalized epsilon squared)を追加 omega・
2007-03-20 2007-05-14 ちょっと修正 生態学会松山大会 (ESJ 54) 自由集会「データ解析で出会う統計的問題 – ベイズ統計学の考えかた・使いかた」 簡単な例題で理解できる 階層ベイズ + MCMC 計算 http://hosho.ees.hokudai.ac.jp/∼kubo/ce/EcoSj2007.html 講釈: 久保拓弥 kubo@ees.hokudai.ac.jp 2007-03-20 1/40 今日のハナシ: 階層ベイズモデル使いになろう! 1. 登場篇: 階層ベイズモデルが必要な状況 そもそも階層ベイズモデルって何なの? 2. 実例篇: 階層ベイズモデルの作りかた 尤度 ∼ 事前分布 ∼ 超事前分布 ∼ ♪ 3. MCMC 篇: 階層ベイズモデルの推定計算 尤度が高そうなあたりをさまよう? p( | )Data prior posterior ~
技術力の低い作品ばかりを集めた、ダメなMakerFaireみたいなやつやりたい(アンケート募集中) - nomolkのブログ そこそこ話題にはなったものの肝心の出演者が集まらず今のところ開催は難しそうな状況なんだけど、こういうときはもったいぶらずにどんどん情報公開していくべきと心得ておりますので、集まった作品の中からこれはと思ったものを紹介します。あと応募きたわけじゃないけど追加で思い出したやつもいくつか。 自転する寿司(マンスーンさん) 【自宅で楽しめる】卓上回転寿司マシンを作った。 | ハイエナズクラブ まず、前回から揚げ君ロボを紹介したmansooonさんから回答をいただきました。ありがとうございます。 卓上回転寿司マシン2 【Tabletop Conveyor belt sushi machine】 - YouTube から揚げ君よりも工作のクオリティは高いけど、それでも寿司が転げ
この記事は検証可能な参考文献や出典が全く示されていないか、不十分です。出典を追加して記事の信頼性向上にご協力ください。(このテンプレートの使い方) 出典検索?: "京大天皇事件" – ニュース · 書籍 · スカラー · CiNii · J-STAGE · NDL · dlib.jp · ジャパンサーチ · TWL(2015年4月) 1951年の昭和天皇来学時、京都大学吉田分校(現在の吉田南構内キャンパス)正門前に立てられた看板 / 同学会の「公開質問状」と比べさほど問題にならなかった。 京大天皇事件(きょうだいてんのうじけん)は、1951年11月12日の昭和天皇の京都大学来学に際して発生した混乱、およびこれを発端とする学生処分事件で、「京大事件」あるいは「天皇事件」とも称される。 概要[編集] 1951年11月12日、関西巡幸途上の昭和天皇が京都大学に来学したとき、多数の学生が正門付近に
このたび、東京大学総合研究博物館と京都大学総合博物館は、東京駅前丸の内側JPタワー内の「インターメディアテク」の企画展示スペース「モデュール」において、共催事業「『驚異の部屋』京都大学ヴァージョン」東京展を開催することになりました。 本展は、京都大学総合博物館の所蔵する学術標本コレクションのなかから、旧制第三高等学校の遺産を中心とする歴史的な標本約五十点を選び出し、一般に公開しようとするものです。 ご存じの通り、京都大学と東京大学の両校は、いずれも官制大学を前身とする国立大学法人ですが、これまで学術標本を用いた博物館レベルでの協働企画展をおこなう機会に恵まれませんでした。また、大学博物館等協議会に加盟している国内諸館のあいだでも本展のような、二館協働の本格的な特別企画が実現したためしはありません。現代社会における大学博物館の存在意義と学術的使命を広く社会へ伝えるために、そしてまた大学博物館
こちら葛飾区亀有公園前派出所、つまりあれだ、ダメ警官の両さんの漫画だ。 こういう書き方は誤解を招くからみんな1巻から読むといいよ んで、ラーメン屋で古い単行本を読んだ時に、ひどく記憶に残ったシーンがある。 確か、NHKともう一チャンネルしか映らないような離島で、小学校の教員を両さんが代行する。 さんざダメな大人について語った後に 『でも、そんなやつでも何とかなっちまうから人生面白い』 と続けて、教室みんなで笑う。 こういう「ダメな大人だな~」って言いながらも、それを笑えるのは、教育で出来ると思う。 両さんは紛れも無くダメな大人で、大原部長が怒り心頭で探しまわったり、島流しされたりするオチが多いけど、 次週にはしれっと現れる。 駄目だ駄目だと怒鳴られながらも、頼りにされたり、馬鹿にされたり、信頼されたりする。 そういうのを、もっとリアルな授業でもやれば良いと思う。 全国を渡り歩きながら絵本を
なんて表現したらいいかわかんなくて、開発支援系サービスって謎表現したけど…。なんつーか、開発支援向けのサービス?クラウドってやつ?ってかいわゆる外部がやってくれる系のサービス(モニタリング/ホスティング/etc)が充実してますよね。んで、一介のWebエンジニアのおれがこの先生きのこるにはどうするかを真剣に考えていたところだった。きのこ。何割かはネタ。 思いついたものを挙げてみる。AWSやGitHubは割愛。言うまでもねーだろ…。 New Relic http://newrelic.com/ 有名なNew Relic。これも説明するまでもないかな。今のチームでコレのお金払う版を使ってるんだけど、「外部APIとの通信個所とDBとの通信個所が遅いように思えるので調査しますわ」→「それNew Relicで見れるよ」とか「各テーブルへのアクセス頻度集計しますわ」→「それNew Relicで見れるよ」
リリース、障害情報などのサービスのお知らせ
最新の人気エントリーの配信
j次のブックマーク
k前のブックマーク
lあとで読む
eコメント一覧を開く
oページを開く