• 検索結果がありません。

chapter 5 フリーソフト「R」を使った自然史情報の統計処理入門 Shun (R) TAKAGI

N/A
N/A
Protected

Academic year: 2018

シェア "chapter 5 フリーソフト「R」を使った自然史情報の統計処理入門 Shun (R) TAKAGI"

Copied!
8
0
0

読み込み中.... (全文を見る)

全文

(1)

♧菙 ⻉简䕎幉さ

ٌرٕ

ؿٔ٦اؿزչ3պ׾⢪׏׋荈搫〷䞔㜠ך窟鎘Ⳣ椚ⰅꟌ

痥鿇

Ⱔ䏧源甧➂ה荈搫ך⽆暟긫 灇瑔㆞

넝加 ⤥

UBLBHJ!IJUPIBLVKQ

幉さٌرٕד䪔ֲ׮ך

暴חꅾ銲זأٓ؎سחכ

ӹو٦ؙ׃ג֮׶תׅ

Ⳣ椚ך⸬卓זוך䕦갟׾鋅׋ְ

㔿㹀⸬卓

ח⸇ִג

⦐⡤䊴זו铎䊴ה׃ג罋䣁׃׋ְ

ٓٝتي⸬卓

׾鐰⣣

Č

㔿㹀⸬卓Ĕ

J

–9

J

铎䊴ė

Č

㔿㹀⸬卓Ĕ

J

–9

J

Č

ٓٝتي⸬卓C

K

铎䊴ė

♧菙简䕎ٌرٕ

♧菙简䕎

幉さ

ٌرٕ

幉さٌرٕ

הٓٝتي⸬卓

٥➙תדךٌرٕכر٦ةך杝甧䚍׾⟎㹀

٥⩝ך粸׶鵤׃瘝杝甧䚍ָ弫׋ׁ׸זְ朐屣׮

ر٦ةך杝甧䚍

0.

0

0.

5

1.

0

1.

5

2.

0

귤" 귤#

A B C D

0.

0

0.

5

1.

0

1.

5

2.

0

귤" 귤#

⢽귤珏ָ傸赇ך䧭ꞿꆀח♷ִ׷䕦갟׾鋅׋ְ

⼏ך䎒赇׾ꅿ㢩ַ׵꧊׭ג㹋꿀

̔ر٦ةכ杝甧

⼏ך鋵ַ׵⼏׆א䎒赇׾꧊׭ג㹋꿀

̔ず鋵ر٦ةכ杝甧דזְ

⩝ך粸׶鵤׃

(2)

鋵꟦ Ⳣ椚ⰻךⴓ侔

Ⳣ椚꟦ךⴓ侔

Ⳣ椚Jך⚥ךK殢湡ך鋵ך⸬卓

Ⳣ椚Jך⸬卓

䧭ꞿⳢ椚

鋵꟦ך䊴

铎䊴

5PUBMל׵אֹ 

Ⳣ椚꟦ל׵אֹ

Ⳣ椚ⰻ鋵꟦ל׵אֹ

⦐⡤꟦铎䊴

Ⳣ椚꟦ל׵אָֹ鋵꟦ל׵אֹח嫰ץג⼧ⴓ㣐ֹֽ׸לⳢ椚ך⸬卓ָ֮׷

幉さٌرٕך罋ִ倯

BPWח״׷鍑匿

"HEXE VIEHGWZHEXEGWZ8

"HEXETEVIRX EWJEGXSVHEXETEVIRX

"WYQQEV]ESZKVS[XLbXVX)VVSVTEVIRXHEXE

)VVSVTEVIRX

(J7YQ7U1IER7U*ZEPYI4V"*

XVX

6IWMHYEPW

)VVSV;MXLMR

(J7YQ7U1IER7U*ZEPYI4V"*

6IWMHYEPW

EWJEGXSV़ॾঙইঽ৏ၜᙋऽबऻ⟿⡤परॠ

XVXTEVIRXKVS[XL G G G G G G G G G G X X X X X X X X X X

QSHIP ESZ=b<)VVSV়৆ও঳൮ច

WYQQEV]QSHIP

3؝٦س

➖ⰱ⠊኉१࿜ीबऴ%23:%✈

ᦗ኉タ➖ೃৎᄽⰱ⠊኉ダ

鸐䌢ךⴓ侔ⴓ匿הך麩ְ

"WYQQEV]ESZKVS[XLbXVXHEXE

(J7YQ7U1IER7U*ZEPYI4V"*

XVXI

6IWMHYEPW

"WYQQEV]ESZKVS[XLbXVX)VVSVTEVIRXHEXE

)VVSVTEVIRX

(J7YQ7U1IER7U*ZEPYI4V"*

XVX

6IWMHYEPW

)VVSV;MXLMR

(J7YQ7U1IER7U*ZEPYI4V"*

6IWMHYEPW

婍䊴˘΍

Ⳣ椚ד铡僇דֹ׋ל׵אֹ

婍䊴鋵꟦˘Ύ

Ⳣ椚ד铡僇דֹ׋ל׵אֹ

婍䊴鋵ⰻ˘Ώ

Ⳣ椚ד铡僇דֹזַ׏׋婍䊴΍鋵꟦婍䊴Ύ鋵ⰻ婍䊴Ώחⴓ鍑ׁ׸׋

'嫰׾鎘皾ׅ׷ⴓ嫡ך铎䊴ָ㢌׻׏׋΍̔Ύךד1⦼׮㢌⻉

ⴓ嫡ך荈歋䏝׮̔ח

/FTUFE

Ⰵ׸㶨

3BOEPNJ[FEDPNQMFUFCMPDL

VOSFQMJDBUFE ✉㝄

$PNQMFUFMZSBOEPNJ[FE

杝甧رؠ؎ٝ

$PNQMFUFMZSBOEPNJ[FE

GBDUPSJBM 湫❛

4QMJUQMPU

ⴓⶴ

杝甧رؠ؎ٝ⟃㢩דכ

ٓٝتي⸬卓罋䣁ך䗳銲

罋䣁ך➬倯׮رؠ؎ٝד

殯ז׷

(3)

" ءؕךְ׷㖑㚖ךءؕ叟ⰻה叟㢩ד⦐⡤ך埠넝庠㹀

# ءؕךְ׷㖑㚖ךءؕ叟ⰻה叟㢩ד⦐⡤ך埠넝庠㹀

$ ءؕךְ׷㖑㚖ծְזְ㖑㚖ד⦐⡤׆א埠넝庠㹀

% ءؕךְ׷㖑㚖ծְזְ㖑㚖ד⦐⡤׆א埠넝庠㹀

–㖑㚖

–㖑㚖

–㖑㚖

–㖑㚖

"湫❛

$Ⰵ׸㶨

#✉㝄

–㖑㚖

–㖑㚖

%杝甧

ءؕ٥؝ٝزٗ٦ٕךرؠ؎ٝד׮Ⰻ鿇麩ֲ

⢽ءؕך䕦갟׾׫׷אך㹋꿀رؠ؎ٝ

EFFSTJUF

–

"

EFFSTJUF

–

#

EFFSTJUF

–

–

$

EFFSTJUF

–

–

%

3ח׻׋ׅر٦ةך鎸鶢岀

HEXEEGWZ

HEXEFGWZ

HEXEGGWZ

HEXEHGWZ

ءؕך㖈♶㖈ծ㖑挿JEכ侧㶵דⰅ׏גְ׷ךדծֿ׸׾ؕذ؞ٔ

ٕؕ㢌侧ה׃ג钠陎ׇׁ׷䗳銲ָ֮׷

⠐॔⧥॔

HEXEE VIEHGWZHEXEEGWZ8

HEXEF VIEHGWZHEXEFGWZ8

HEXEG VIEHGWZHEXEGGWZ8

HEXEH VIEHGWZHEXEHGWZ8

HIIV१ॾঙইঽ৏ၜᙋऽबऻ⟿⡤परॠ

HEXEEHIIV EWJEGXSVHEXEEHIIV

HEXEFHIIV EWJEGXSVHEXEFHIIV

HEXEGHIIV EWJEGXSVHEXEGHIIV

HEXEHHIIV EWJEGXSVHEXEHHIIV

WMXI१ॾঙইঽ৏ၜᙋऽबऻ⟿⡤परॠ

HEXEEWMXI EWJEGXSVHEXEEWMXI

HEXEFWMXI EWJEGXSVHEXEFWMXI

HEXEGWMXI EWJEGXSVHEXEGWMXI

HEXEHWMXI EWJEGXSVHEXEHWMXI

鍑匿⵸ח㢌侧ך㘗㢌䳔

–

WYQQEV]ESZLIMKLXbHIIV)VVSVWMXIHIIVHEXEE

㔿㹀⸬卓

㢌ꆀ⸬卓

婍䊴

TJUFד铡僇ׁ׸׷ל׵אֹ

EFFSד铡僇ׁ׸׷ל׵אֹ EFFSTJUFד铡僇ׁ׸׷ל׵אֹ

婍׶⦐⡤䊴

)VVSVWMXI (J7YQ7U1IER7U*ZEPYI4V"* 6IWMHYEPW )VVSVWMXIHIIV (J7YQ7U1IER7U*ZEPYI4V"* HIIV 6IWMHYEPW )VVSV;MXLMR (J7YQ7U1IER7U*ZEPYI4V"* 6IWMHYEPW

㢌ꆀ⸬卓

❛✼⡲欽

$PNQMFUFMZSBOEPNJ[FEGBDUPSJBMEFTJHO

(4)

㔿㹀⸬卓 㢌ꆀ⸬卓 婍䊴

TJUFד铡僇ׁ׸׷ל׵אֹ

EFFSד铡僇ׁ׸׷ל׵אֹ 婍׶⦐⡤䊴ַא⸬卓ך㖑㚖䊴

–

)VVSVWMXI

(J7YQ7U1IER7U*ZEPYI4V"* 6IWMHYEPW

)VVSV;MXLMR

(J7YQ7U1IER7U*ZEPYI4V"* HIIV 6IWMHYEPW

WYQQEV]ESZLIMKLXbHIIV)VVSVWMXIHEXEF

#✉㝄رؠ؎ٝ

3BOEPNJ[FEDPNQMFUFCMPDLEFTJHO

㔿㹀⸬卓 㢌ꆀ⸬卓

婍䊴

TJUFד铡僇ׁ׸׷ל׵אֹ ַא⸬卓ך㖑㚖䊴 EFFSד铡僇ׁ׸׷ל׵אֹ

婍׶⦐⡤䊴

)VVSVWMXI

(J7YQ7U1IER7U*ZEPYI4V"* HIIV 6IWMHYEPW

)VVSV;MXLMR

(J7YQ7U1IER7U*ZEPYI4V"* 6IWMHYEPW

WYQQEV]ESZLIMKLXbHIIV)VVSVWMXIHEXEG

–

–

$Ⰵ׸㶨رؠ؎ٝ

/FTUFEEFTJHO

㔿㹀⸬卓 婍䊴

EFFSד铡僇ׁ׸׷ל׵אֹ 婍׶⦐⡤䊴ַא㖑㚖 䊴ַא⸬卓ך㖑㚖䊴

(J7YQ7U1IER7U*ZEPYI4V"* HIIV 6IWMHYEPW

WYQQEV]ESZLIMKLXbHIIVHEXEH

–

–

%杝甧رؠ؎ٝ

$PNQMFUFMZSBOEPNJ[FEEFTJHO

–

ءؕך⸬卓

㖑㚖ך⸬卓

ءؕך⸬卓ך㖑㚖䊴

⦐⡤䊴

–

ءؕך⸬卓

㖑㚖ך⸬卓

⦐⡤䊴ءؕך⸬卓ך㖑㚖䊴

–

–

ءؕך⸬卓

㖑㚖ך⸬卓ءؕך⸬卓ך㖑㚖䊴

⦐⡤䊴

–

–

ءؕך⸬卓

㖑㚖ך⸬卓ءؕך⸬卓ך㖑㚖䊴⦐⡤䊴

3BOEPNJ[F

EDPNQMFUF

CMPDL

✉㝄

$PNQMFUFMZ

SBOEPNJ[FE

GBDUPSJBM

湫❛

/FTUFE

Ⰵ׸㶨

$PNQMFUFMZ

SBOEPNJ[FE

杝甧

(5)

٥姻鋉ⴓ䋒⟃㢩׾⟎㹀ׅ׷

♧菙⻉简䕎幉さٌرٕ

(-..

׮

♧菙⻉简䕎幉さٌرٕ

ৎKPQIVMRPMFVEV]PQI

HMNד⢪ֲ״ֲזⴓ䋒ח㼎䘔ծ醱侧ךٓٝتي⸬卓׮䪔ִ׷կ

姻鋉ⴓ䋒דװ׷הֹכ

PQIV

ৎKPQQEHQFMRPMFVEV]KPQQ%(1&

刿ח圫ղזⴓ䋒ד(-.٥(-..ך鍑匿〳腉կ

تؐٝٗ٦س؟؎زך鼅䫛3'PSHF ַ׵؎ٝأز٦ٕ〳腉

♧菙简䕎幉さٌرٕ

♧菙⻉简䕎幉さٌرٕ

Ⳣ椚ך⸬卓 㜥䨽꟦ל׵אֹ

ESZ]bXVX)VVSVWMXI

PQIV]bXVX`WMXI

KPQIV]bXVX`WMXIJEQMP]!TSMWWSR

KPQQEHQF]bXVX`WMXIJEQMP]!׎TSMWWSR׏

♧菙٥♧菙⻉简䕎幉さٌرٕ׾䪔ֲꟼ侧

ⴓ侔ⴓ匿邌אֻ׷ז׵ ٌرٕ鼅䫛הַ׮׃׋ְז׵

⦐ղך⸬卓ך䓼ׁח莆㄂ָ剣׷̔

㔿㹀⸬卓

ל׵אֹ׾罋䣁׃׋ְ׌ֽ̔

ٓٝتي⸬卓

㖑㚖זו瑞꟦涸ז׮ך׌ֽדזֻծ⦐⡤䊴ծ儗꟦זו

׮ず圫חٌرؚٔٝ〳ず♧⦐⡤ך⿾䗁庠㹀זו

䠐㔳׃׋鍑匿חׅ׷׋׭חכծر٦ةךⰅ׸倯٥ٓٝ

تي⸬卓ךⰅ׸倯ח岣䠐

תה׭♧菙⻉简䕎幉さٌرٕ

׮׏ה莆㄂֮׷➂כֿך劤֮׋׶ד˘

&YUFOEJOHUIFMJOFBSNPEFMXJUI3

(FOFSBMJ[FE-JOFBS.JYFE&GGFDUTBOE

/POQBSBNFUSJD3FHSFTTJPO.PEFMT

⿫罋项俱

(6)

/FTUFEEFTJHO

銲㔓"ך⚥ח銲㔓#ָⰅ׸㶨ծ粸׶鵤׃֮׶

˘

˘

"GJYFE #SBOEPN

㖑挿꟦ךⴓ侔 ؼءWTꟚ佝꟦ךⴓ侔 ؼءך⸬卓 㖑挿꟦铎䊴 㖑挿ⰻ㜥䨽꟦铎䊴

⢽

ؼء䌒㖑挿٥Ꟛ佝宏꬗㖑挿ָ֮

׶ծ׉׸׊׸ך㖑挿ד皘䨽׆א庠

錁庠侧 4JUF 4JUF 4JUF 4JUF

0QFO

5SBQB

/FTUFEEFTJHO

㖑挿꟦ךⴓ侔

ؼءWTꟚ佝宏꬗꟦ךⴓ侔

⢽

ؼء䌒㖑挿٥Ꟛ佝宏꬗㖑挿ָ֮

׶ծ׉׸׊׸ך㖑挿ד皘䨽׆א庠

0QFO

5SBQB

4JUF 4JUF 4JUF

4JUF

3BOEPNJ[FECMPDLEFTJHO

✉㝄

銲㔓"ה銲㔓#כ湫❛ծ粸׶鵤׃ז׃

"GJYFE #SBOEPN

ؼءך⸬卓 ـٗحؙ꟦铎䊴ְ׆׸ד׮铡僇ׁ׸זְ婍䊴

⢽

ـٗحָؙ֮׶ծ׉׸׊׸ךـٗح

ؙדؼءⴘ׶《׶٥ⴘ׶婍׃㖑挿׆

א庠㹀

Ⳣ椚٥ـٗحؙד铡僇 ׁ׸זְⴓ侔婍䊴 ⴘ׶《׶WTⴘ׶婍׃꟦ךⴓ 侔

錁庠侧 #MPDL

#MPDL #MPDL #MPDL

ⴘ《

0

ⴘ婍5

3BOEPNJ[FECMPDLEFTJHO

铡僇דֹזְ铎䊴

ⴘ׶婍׃WTⴘ׶《׶꟦ךⴓ侔

⢽

ـٗحָؙ֮׶ծ׉׸׊׸ךـٗح

ؙדⴘ׶《׶٥ⴘ׶婍׃㖑挿׆א庠

#ML #ML #ML

#ML

5

0

0

5

0

5

0

5

(7)

'BDUPSJBMEFTJHO

銲㔓"ה銲㔓#כ湫❛ծぐ穈׫さ׻ׇח粸׶鵤׃֮׶

˘

˘

"GJYFE #SBOEPN

⢽

ـٗحָؙ֮׶ծ׉׸׊׸ךـٗح

ؙדⴘ׶《׶٥ⴘ׶婍׃㖑挿׆א庠

ؼءך⸬卓 ـٗحؙ꟦铎䊴

ְ׆׸ד׮铡僇ׁ׸זְ婍䊴

ؼءך⸬卓ךـٗحؙ꟦铎䊴

ⴘ׶《׶ך⸬卓ךـٗح ؙ꟦דךל׵אֹ ⴘ׶《׶WTⴘ׶婍׃꟦ךⴓ 侔

'BDUPSJBMEFTJHO

⢽

ـٗحָؙ֮׶ծ׉׸׊׸ךـٗح

ؙדⴘ׶《׶٥ⴘ׶婍׃㖑挿׆א庠

錁庠侧 #MPDL

#MPDL

ⴘ《

0

ⴘ婍5

#MPDL

#MPDL

5

0

0

5

ⴘ׶婍׃WTⴘ׶《׶꟦ךⴓ侔

ⴘ׶《׶⸬卓ךـٗحؙ꟦דךⴓ 侔

㹋꿀رؠ؎ٝח״׷"/07"ך麩ְ

$PNQMFUFMZ

SBOEPNJ[FE

杝甧

/FTUFE

Ⰵ׸

㶨

3BOEPNJ[F

ECMPDL

✉㝄

'BDUPSJBM

湫❛

Ⳣ椚ך⸬卓Ⳣ椚ד铡僇דֹזְ婍䊴

Ⳣ椚ך⸬卓Ⳣ椚הـٗحؙך⸬卓׾ꤐְ׋婍䊴

ـٗحؙⰟ鸐ךⳢ椚ך⸬卓ـٗحؙ׀הךⳢ椚ך⸬卓ךל׵א

Ⳣ椚 婍䊴

Ⳣ椚ך⸬卓Ⳣ椚ⰻךـٗحؙ꟦ךל׵אֹ

Ⳣ椚ⰻـٗحؙ

ـٗحؙ

Ⳣ椚ـٗحؙ❛✼⡲ 欽

⢽겗ءְָؕ׷ה،ؔؗכ幾׷ַ

ءؕךⴓ䋒㖑㚖ծꬊⴓ䋒㖑㚖דぐ㖑㚖吳ח

ꟼ׃גワ㔲ךԅ吳ך侧׾庠㹀/FTUFEEFTJHO

ꬊ頾侭侧⦼

family=poisson

atago fudago godai hiratuka kimikame

0

5

10

15

site

(8)

glmer.model5.3<- glmer(female~deer+(1|site),data5.3,family=poisson) summary(glmer.model5.3)

Generalized linear mixed model fit by maximum likelihood ['glmerMod'] Family: poisson ( log )

Formula: female ~ deer + (1 | site) Data: data5.3

AIC BIC logLik deviance 125.0005 128.6572 -59.5003 119.0005

Random effects:

Groups Name Variance Std.Dev. site (Intercept) 0.1544 0.393 Number of obs: 25, groups: site, 5

Fixed effects:

Estimate Std. Error z value Pr(>|z|) (Intercept) 1.0085 0.2760 3.654 0.000258 *** deernd 1.2864 0.4043 3.182 0.001463 ** ---

Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Correlation of Fixed Effects: (Intr)

deernd -0.683

⤘侧ך䱿㹀⦼ח㼎ׅ׷8BME嗚㹀ח״׷1⦼ 㢌侧ך⸬卓׾׫׷ז׵㽍䏝嫰嗚㹀٥

1BSBNFUSJD#PPUTUSBQךקֲָ葺ְ

䱿㹀

ءؕ㖈劍䖉㺘䏝 d = exp(1.0085)

ءؕ♶㖈劍䖉㺘䏝 nd = exp(1.00851.2864)

嗚㹀

ؕ؎✳⛦鵚⡂׾欽ְ׋㽍䏝嫰嗚㹀

1BSBNFUSJDCPPUTUSBQ岀ח״׷鵚⡂ח״׵זְ嗚㹀

glmer.model5.3.0<- glmer(female~1+(1|site),data5.3,family=poisson) anova(glmer.model5.3.0,glmer.model5.3)

Data: data5.3 Models:

glmer.model5.3.0: female ~ 1 + (1 | site) glmer.model5.3: female ~ deer + (1 | site)

Df AIC BIC logLik deviance Chisq Chi Df Pr(>Chisq) glmer.model5.3.0 2 128.35 130.79 -62.177 124.36 glmer.model5.3 3 125.00 128.66 -59.500 119.00 5.3543 1 0.02067 * ---

Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ̒鋅׋ְ⸬卓Ⰵ׶EFFS̒鋅׋ְ⸬卓䫙ֹOVMM

–㼎侧㽍䏝ך䊴㽍䏝嫰

չ鋅׋ְ⸬卓׾䫙ְ׋ٌرٕך䱿㹀⦼ך⯋דر٦ة׾涪欰ׁ ׇծ㹋ر٦ةהず圫ך鍑匿׾遤ְծ㽍䏝嫰׾鎘皾ׇׁ׷պ׾ 粸׶鵤׃遤ְ㔐הַծ鵚⡂ח״׵זְ㽍䏝嫰ךⴓ䋒 ׾䖤׷倯岀կ鎘皾ח儗꟦ַַָ׷կ鑫׃ְ倯岀כأؙٔف ز⿫撑

ءؕך⸬卓֮׶׉ֲ

参照

関連したドキュメント

There is a bijection between left cosets of S n in the affine group and certain types of partitions (see Bjorner and Brenti (1996) and Eriksson and Eriksson (1998)).. In B-B,

(The Elliott-Halberstam conjecture does allow one to take B = 2 in (1.39), and therefore leads to small improve- ments in Huxley’s results, which for r ≥ 2 are weaker than the result

[r]

“Breuil-M´ezard conjecture and modularity lifting for potentially semistable deformations after

lines. Notice that Theorem 4 can be reformulated so as to give the mean harmonic stability of the configuration rather than that of the separate foliations. To this end it is

S., Oxford Advanced Learner's Dictionary of Current English, Oxford University Press, Oxford

At the end of the section, we will be in the position to present the main result of this work: a representation of the inverse of T under certain conditions on the H¨older

支払方法 支払日 ※② 緊急時連絡先等 ※③.