ffp_opt_sobin_target_row: Theorem 1 and 2, solves the binary targeting queue with Ai_i,...

Description Usage Arguments Value Author(s) References Examples

View source: R/ffp_opt_sobin.R

Description

This does not solve for the full solution, just the targeting queue, individual position on queue.

Usage

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
ffp_opt_sobin_target_row(
  ls_row,
  fl_rho,
  ar_A,
  ar_alpha,
  ar_beta,
  svr_A_i = "A",
  svr_alpha_i = "alpha",
  svr_beta_i = "beta"
)

Arguments

ls_row

list a row from dataframe tibble row where there are variables for A, alpha and beta

fl_rho

float preference for equality for the planner, negative infinity to 1

ar_A

float array of expected outcome without provision for all N

ar_alpha

float array of expected effect of provision for each i of N

ar_beta

float array of planner bias

svr_A_i

string name of the A_i variable, any nonlinear or linear evaluated expected outcome

svr_alpha_i

string name of the alpha_i variable, individual specific effects of treatment

svr_beta_i

string name of the beta_i variable, relative preference weight for each child

Value

a list with an integer and an array

Author(s)

Fan Wang, http://fanwangecon.github.io

References

https://fanwangecon.github.io/PrjOptiAlloc/reference/ffp_opt_sobin_target_row.html https://fanwangecon.github.io/PrjOptiAlloc/articles/ffv_opt_sobin_rkone_allrw_car.html

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
library(tibble)
library(dplyr)
library(tidyr)
ar_alpha <- c(0.1,1.5,2.5,4)
ar_A <-     c(0.5,1.5,3.5,6.5)
ar_beta <- c(0.25,0.25,0.25,0.25)
mt_alpha_A <- cbind(ar_alpha, ar_A, ar_beta)
ar_st_varnames <- c('alpha', 'A', 'beta')
tb_alpha_A <- as_tibble(mt_alpha_A) %>% rename_all(~c(ar_st_varnames))
tb_alpha_A

ar_rho <- c(-100, -1.1, 0.01, 0.10, 0.9)
for (it_rho_ctr in seq(1,length(ar_rho))) {
  fl_rho <- ar_rho[it_rho_ctr]
  ls_ranks <- ffp_opt_sobin_target_row(tb_alpha_A[1,], 0.1, ar_A, ar_alpha, ar_beta)
  it_rank <- ls_ranks$it_rank
  ar_it_rank <- ls_ranks$ar_it_rank
  ar_fl_rank_val <- ls_ranks$ar_fl_rank_val
  cat('fl_rho:', fl_rho, 'it_rank:', it_rank, '\n')
  cat('ar_it_rank:', ar_it_rank, '\n')
  cat('ar_fl_rank_val:', ar_fl_rank_val, '\n')
}

FanWangEcon/PrjOptiAlloc documentation built on Jan. 25, 2022, 6:55 a.m.