set.seed(1234)
## Parameters of demand function
alpha <- 4
beta <- c(-5, -4)
gamma <- c(1, 0.5)
## Make a list of factor level
char_levels <- list(char1 = c(0,1), char2 = c(0,1))
## Make list of all potential products
XP <- expand.grid(char_levels)
J <- dim(XP)[1]
XP <- data.frame(XP, prodname=1:J)
a1 <- sample(c(0,1), 4, replace = TRUE)
a2 <- sample(c(0,1), 4, replace = TRUE)
Xi <- rnorm(1, 0, 1)
Omega <- rnorm(1, 0, 1)
my_prod_table <- data.frame(j = XP$prodname, delta = (cbind(XP$char1, XP$char2) %*% beta), mc = (cbind(XP$char1, XP$char2) %*% gamma))
my_mkt_prods <- rbind(my_prod_table[which(a1==1),], my_prod_table[which(a2==1),])
my_mkt_prods <- data.frame(my_mkt_prods, firm = c(1,1,1,0,0))
my_ldm_objR6 <- Logit_Demand_Market$new(my_mkt_prods$firm, my_mkt_prods$j, my_mkt_prods$delta, Mc_fixed = my_mkt_prods$mc, Struct_error = Xi, U_out_opt = 1, Deriv_price = alpha)
my_ldm_objR6$getMarket()
my_ldm_objR6$share()
my_ldm_objR6$getMarket()
my_ldm_objR6$getOwnership()
my_ldm_objR6$getDs()
my_ldm_objR6$Ds_fun()
my_ldm_objR6$getDs()
#my_ldm_objR6_2 <- Logit_Demand_Market$new(my_mkt_prods$firm, my_mkt_prods$j, my_mkt_prods$delta, Mc_fixed = my_mkt_prods$mc, Struct_error = Xi, U_out_opt = 1, Deriv_price = alpha)
#my_ldm_objR6_2$getMarket()
#my_ldm_objsR6 <- purrr::map(1:500, .f=~my_ldm_objR6_2$clone)
my_ldm_objR6$zeta_fixed_point(1e-6)
st_err_fun <- function(J){
rnorm(J, 0, 1)
}
mc_err_fun <- function(J){
runif(J, 0, 1)
}
test1 <- my_ldm_objR6$exp_profits(mc_error_fun = mc_err_fun, struct_error_fun = st_err_fun, draws = 500, tol=1e-6)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.