お久しぶりです。
この記事はTDN覚書なので、スルーしてどうぞ。
最近は大学からの脱出のために論文を書いているのですが、これまで(盲目的に)扱っていた分散分析ANOVAから一般化線形モデルGLMに移行しなさいという命令が出たので勉強しています。
・・・
とある虫4種(A、B、C、D)がいて、各種の♀成虫の持っている卵の数を数えてデータにしました。グラフにすると、なんとなく種Bが一番卵を持ってそうだ。ということで、今までは分散分析をしてくれる有料ソフトに突っ込んで、P値をだし、P<0.05なら種Bが他より有意に卵を多くもっているといえる(多重比較とかやって、おわり)。というストーリーでした。
GLMとはいったい・・うごごご。
1晩いろんなウェブサイトを廻ってみましたが、よくわからない。どうも今までやってきたものとは毛色が違う様子?
図書館に走って本を適当に借りてきました。「データ解析のための統計モデリング入門 著:久保拓弥」の第3章まで読んでなんとなく理解。
今回の例で言うと、卵の数を計算するためのモデルづくり(観測データにうまいこと当てはまるようなモデルづくり)のための武器なのかな?
今までの線形モデルLM(分散分析ふくむ)では、前提条件で等分散性だったり正規分布だけだったり?なるほど聞いた覚えはある。といった感じ。GLMにレベルアップするとそういうのは考えなくてよい(自分で設定できる)?らしい。聞こえはシンプルでイイネ・
はてなマークばかりだけども、GLMを使うために統計ソフトRくんに計算してもらう。
Speciesは4種(A、B、C、D)あって、卵の数Eggのモデルづくりは、
>glm(Egg~Species, family="poisson",data=元データ) と入力
すると、
切片とB C D の係数がでた。Aはどこいったの?(←ここで1日悩んだ)
どうも基準として種Aの係数を0(ゼロ)にしてから、他の種の係数を計算しているそうです。あたしきいてない!
とにかくこれで、各種に対する卵の数のモデルが組み立てられた感じ。
で、モデルを作ってどうするんだろう・・・
がイマココです。
他にも虫の齢期などの要因があるので、他のモデルを作ってAICの比較をするつもりだけども・・・
さらに、このglm()の結果をaov()に放り込んだりするとそれっぽいのが出るけど、変なことやってそう。
勉強したいけどそろそろ提出期限ががが
おしまい
この記事はTDN覚書なので、スルーしてどうぞ。
最近は大学からの脱出のために論文を書いているのですが、これまで(盲目的に)扱っていた分散分析ANOVAから一般化線形モデルGLMに移行しなさいという命令が出たので勉強しています。
・・・
とある虫4種(A、B、C、D)がいて、各種の♀成虫の持っている卵の数を数えてデータにしました。グラフにすると、なんとなく種Bが一番卵を持ってそうだ。ということで、今までは分散分析をしてくれる有料ソフトに突っ込んで、P値をだし、P<0.05なら種Bが他より有意に卵を多くもっているといえる(多重比較とかやって、おわり)。というストーリーでした。
GLMとはいったい・・うごごご。
1晩いろんなウェブサイトを廻ってみましたが、よくわからない。どうも今までやってきたものとは毛色が違う様子?
図書館に走って本を適当に借りてきました。「データ解析のための統計モデリング入門 著:久保拓弥」の第3章まで読んでなんとなく理解。
今回の例で言うと、卵の数を計算するためのモデルづくり(観測データにうまいこと当てはまるようなモデルづくり)のための武器なのかな?
今までの線形モデルLM(分散分析ふくむ)では、前提条件で等分散性だったり正規分布だけだったり?なるほど聞いた覚えはある。といった感じ。GLMにレベルアップするとそういうのは考えなくてよい(自分で設定できる)?らしい。聞こえはシンプルでイイネ・
はてなマークばかりだけども、GLMを使うために統計ソフトRくんに計算してもらう。
Speciesは4種(A、B、C、D)あって、卵の数Eggのモデルづくりは、
>glm(Egg~Species, family="poisson",data=元データ) と入力
すると、
切片とB C D の係数がでた。Aはどこいったの?(←ここで1日悩んだ)
どうも基準として種Aの係数を0(ゼロ)にしてから、他の種の係数を計算しているそうです。あたしきいてない!
とにかくこれで、各種に対する卵の数のモデルが組み立てられた感じ。
で、モデルを作ってどうするんだろう・・・
がイマココです。
他にも虫の齢期などの要因があるので、他のモデルを作ってAICの比較をするつもりだけども・・・
さらに、このglm()の結果をaov()に放り込んだりするとそれっぽいのが出るけど、変なことやってそう。
勉強したいけどそろそろ提出期限ががが
おしまい