Conditional Logit with Capacity Constraints

conditional_logit_with_capacity_constraints.R

set.seed(1)

N <- 800
J <- 80

id_firm = 1:J
delta = c(1:J)/J

capacity <- c(rep(10, 7*J/8), rep(100, J/8))
enrollment <- rep(0, J)

for(i in 1:N){
	placed <- FALSE
	temp_list <- delta + EnvStats::revd(J)
	while(!placed){
		pref <- which(temp_list == max(temp_list))
		if(capacity[pref] <= enrollment[pref]){
			temp_list[pref] <- -Inf
		} else {
			enrollment[pref] <- enrollment[pref] + 1
			placed <- TRUE
		}
	}
}

shares_with_constraint <- enrollment / N
shares_with_constraint