過半数の一員であることは、動物グループのリーダーとして浮上するのにどのように役立つかを評価したいと思います.
リーダーが多数派か少数派かを評価した 10 のケースがあるとします。
Leader <- c(1,1,1,1,0,1,1,1,0,1,0,0,0,0,1,0,0,0,1,0)
Case <- as.factor(c(1,2,3,4,5,6,7,8,9,10,1,2,3,4,5,6,7,8,9,10))
Majority <- as.factor(c("Maj","Maj","Maj","Maj","Maj","Maj","Maj","Maj","Maj","Maj",
"Min","Min","Min","Min","Min","Min","Min","Min","Min","Min"))
leadMaj <- data.frame(Leader,Case,Majority)
binomial.glmer <- glmer(Leader ~ Majority + (1|Case),
family = binomial, data = leadMaj)
summary(binomial.glmer)
結果は、マイノリティ出身であると、グループをリードする可能性が大幅に低下することを示しています
Generalized linear mixed model fit by maximum likelihood (Laplace Approximation) [glmerMod]
Family: binomial ( logit )
Formula: Leader ~ Majority + (1 | Case)
Data: leadMaj
AIC BIC logLik deviance df.resid
26 29 -10 20 17
Scaled residuals:
Min 1Q Median 3Q Max
-2.0 -0.5 0.0 0.5 2.0
Random effects:
Groups Name Variance Std.Dev.
Case (Intercept) 0 0
Number of obs: 20, groups: Case, 10
Fixed effects:
Estimate Std. Error z value Pr(>|z|)
(Intercept) 1.3863 0.7906 1.754 0.0795 .
MajorityMin -2.7726 1.1180 -2.480 0.0131 *
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
Correlation of Fixed Effects:
(Intr)
MajorityMin -0.707
ただし、グループは大多数が 8 人、少数が 2 人で構成されていました。ケースの 80% で過半数が主導したことがわかります。これは予想どおりです。
問題は、p=0.5 ではなく p=0.8 の二項分布を含めるにはどうすればよいかということです。