0

過半数の一員であることは、動物グループのリーダーとして浮上するのにどのように役立つかを評価したいと思います.

リーダーが多数派か少数派かを評価した 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 の二項分布を含めるにはどうすればよいかということです。

4

0 に答える 0