R/abdl.r

In edseab/LeadershipModel: Evolutionary simulation of Leadership Strategies

abdl <- function(pvol,pvolHiL=0,grpsz=5, L=3,LL=0.8, Lavg = 3, P=1, Pavg=1, inv = 0.6,volcost=0.1, Lcost=0.2, E=0, Eavg=0, Ecoef = 0.9, HiLprop = 0,HiL=4,PHiL= 0.5, Baseline=5){

  Lavg <- L*(1-HiLprop)+HiL*HiLprop
  Pavg <- P*(1-HiLprop)+PHiL*HiLprop
  volavg <- pvol*(1-HiLprop)+pvolHiL*HiLprop
  
  
  EHiL <- min((LL + (HiL-LL)*inv)*grpsz,E)
  E <- min((LL + (L-LL)*inv)*grpsz,E)
  Eavg <- min((LL + (Lavg-LL)*inv)*grpsz,E)
  
  HiL_return <- ((LL + (HiL-LL)*inv)*grpsz - EHiL)/grpsz + Baseline# Leadership of ego
  L_return <- ((LL + (L-LL)*inv)*grpsz - E)/grpsz + Baseline# Leadership of ego
  L_return_avg <-((LL + (Lavg-LL)*inv)*grpsz - Eavg)/grpsz + Baseline
  

HiL_ldr_stay <- HiL_return  + PHiL - Lcost - inv*PHiL + Ecoef*E
LoL_ldr_stay <- L_return  + P - Lcost - inv*P + Ecoef*E
HiL_ldr_gone <- (1-pavg)^(grpsz-1)*(LL+P + Baseline) + (1-(1-pavg)^(grpsz-1))*(L_return_avg + PHiL)
LoL_ldr_gone <- (1-pavg)^(grpsz-1)*(LL+P + Baseline) + (1-(1-pavg)^(grpsz-1))*(L_return_avg + P)


HiL_ldr_adv <- HiL_ldr_stay-HiL_ldr_gone
LoL_ldr_adv <- LoL_ldr_stay-LoL_ldr_gone
return(c(HiL=HiL_ldr_adv,LoL=LoL_ldr_adv))
}