Go to the RMD, R, PDF, or HTML version of this file. Go back to fan’s REconTools research support package, R4Econ examples page, PkgTestR packaging guide, or Stat4Econ course page.
Build up a list of strings, where the strings share common components. Iteratre over lists to generate variations in elements of the string list.
# common string components
st_base_name <- 'snwx_v_planner_docdense'
st_base_middle <- 'b1_xi0_manna_88'
# numeric values to loop over
ar_st_beta_val <- c('bt60', 'bt70', 'bt80', 'bt90')
ar_st_edu_type <- c('e1lm2', 'e2hm2')
# initialize string list
ls_snm <- vector(mode = "list", length = length(ar_st_beta_val)*length(ar_st_edu_type))
# generate list
it_ctr = 0
for (st_beta_val in ar_st_beta_val) {
for (st_edu_type in ar_st_edu_type) {
it_ctr = it_ctr + 1
# snm_file_name <- 'snwx_v_planner_docdense_e2hm2_b1_xi0_manna_88_bt90'
snm_file_name <- paste(st_base_name, st_edu_type, st_base_middle, st_beta_val, sep ='_')
ls_snm[it_ctr] <- snm_file_name
}
}
# print
for (snm in ls_snm) {
print(snm)
}## [1] "snwx_v_planner_docdense_e1lm2_b1_xi0_manna_88_bt60"
## [1] "snwx_v_planner_docdense_e2hm2_b1_xi0_manna_88_bt60"
## [1] "snwx_v_planner_docdense_e1lm2_b1_xi0_manna_88_bt70"
## [1] "snwx_v_planner_docdense_e2hm2_b1_xi0_manna_88_bt70"
## [1] "snwx_v_planner_docdense_e1lm2_b1_xi0_manna_88_bt80"
## [1] "snwx_v_planner_docdense_e2hm2_b1_xi0_manna_88_bt80"
## [1] "snwx_v_planner_docdense_e1lm2_b1_xi0_manna_88_bt90"
## [1] "snwx_v_planner_docdense_e2hm2_b1_xi0_manna_88_bt90"# if string in string
grepl('snwx_v_planner', snm)## [1] TRUESave a list of matrixes. Retrieve Element of that list via loop.
# Define an array to loop over
ar_fl_mean <- c(10, 20, 30)
# store restuls in named list
ls_mt_res = vector(mode = "list", length = length(ar_fl_mean))
ar_st_names <- paste0('mean', ar_fl_mean)
names(ls_mt_res) <- ar_st_names
# Loop and generat a list of dataframes
for (it_fl_mean in seq(1, length(ar_fl_mean))) {
fl_mean = ar_fl_mean[it_fl_mean]
# dataframe
set.seed(it_fl_mean)
tb_combine <- as_tibble(
matrix(rnorm(4,mean=fl_mean,sd=1), nrow=2, ncol=3)
) %>%
rowid_to_column(var = "id") %>%
rename_all(~c(c('id','var1','varb','vartheta')))
ls_mt_res[[it_fl_mean]] = tb_combine
}
# Retrieve elements
print(ls_mt_res[[1]])
print(ls_mt_res$mean10)
print(ls_mt_res[['mean10']])
# Print via Loop
for (it_fl_mean in seq(1, length(ar_fl_mean))) {
tb_combine = ls_mt_res[[it_fl_mean]]
print(tb_combine)
}# Define Lists
ls_num <- list(1,2,3)
ls_str <- list('1','2','3')
ls_num_str <- list(1,2,'3')
# Named Lists
ar_st_names <- c('e1','e2','e3')
ls_num_str_named <- ls_num_str
names(ls_num_str_named) <- ar_st_names
# Add Element to Named List
ls_num_str_named$e4 <- 'this is added'Initiate an empty list and add to it
# Initiate List
ls_abc <- vector(mode = "list", length = 0)
# Add Named Elements to List Sequentially
ls_abc$a = 1
ls_abc$b = 2
ls_abc$c = 'abc\'s third element'
# Get all Names Added to List
ar_st_list_names <- names(ls_abc)
# Print list in a loop
print(ls_abc)## $a
## [1] 1
##
## $b
## [1] 2
##
## $c
## [1] "abc's third element"for (it_list_ele_ctr in seq(1,length(ar_st_list_names))) {
st_list_ele_name <- ar_st_list_names[it_list_ele_ctr]
st_list_ele_val <- ls_abc[it_list_ele_ctr]
print(paste0(st_list_ele_name,'=',st_list_ele_val))
}## [1] "a=1"
## [1] "b=2"
## [1] "c=abc's third element"The function below ffi_lst2str is also a function in REconTools: ff_sup_lst2str.
# list to String printing function
ffi_lst2str <- function(ls_list, st_desc, bl_print=TRUE) {
# string desc
if(missing(st_desc)){
st_desc <- deparse(substitute(ls_list))
}
# create string
st_string_from_list = paste0(paste0(st_desc, ':'),
paste(names(ls_list), ls_list, sep="=", collapse=";" ))
if (bl_print){
print(st_string_from_list)
}
}
# print full
ffi_lst2str(ls_num)## [1] "ls_num:=1;=2;=3"ffi_lst2str(ls_str)## [1] "ls_str:=1;=2;=3"ffi_lst2str(ls_num_str)## [1] "ls_num_str:=1;=2;=3"ffi_lst2str(ls_num_str_named)## [1] "ls_num_str_named:e1=1;e2=2;e3=3;e4=this is added"# print subset
ffi_lst2str(ls_num[2:3])## [1] "ls_num[2:3]:=2;=3"ffi_lst2str(ls_str[2:3])## [1] "ls_str[2:3]:=2;=3"ffi_lst2str(ls_num_str[2:4])## [1] "ls_num_str[2:4]:=2;=3;=NULL"ffi_lst2str(ls_num_str_named[c('e2','e3','e4')])## [1] "ls_num_str_named[c(\"e2\", \"e3\", \"e4\")]:e2=2;e3=3;e4=this is added"Generate a multiple dimensional list:
List allows for different data types to be stored together.
Note that element specific names in named list are not preserved when the list is reshaped to be two dimensional. Two dimensional list, however, could have row and column names.
# Dimensions
it_M <- 2
it_Q <- 3
it_N <- it_M*it_Q
# Initiate an Empty MxQ=N element list
ls_2d_flat <- vector(mode = "list", length = it_N)
ls_2d <- ls_2d_flat
# Named flat
ls_2d_flat_named <- ls_2d_flat
names(ls_2d_flat_named) <- paste0('e',seq(1,it_N))
ls_2d_named <- ls_2d_flat_named
# Reshape
dim(ls_2d) <- c(it_M, it_Q)
# named 2d list can not carry 1d name after reshape
dim(ls_2d_named) <- c(it_M, it_Q)Print Various objects generated above, print list flattened.
# display
ffi_lst2str(ls_2d_flat_named)## [1] "ls_2d_flat_named:e1=NULL;e2=NULL;e3=NULL;e4=NULL;e5=NULL;e6=NULL"# print(ls_2d_flat_named)
ffi_lst2str(ls_2d_named)## [1] "ls_2d_named:=NULL;=NULL;=NULL;=NULL;=NULL;=NULL"print(ls_2d_named)## [,1] [,2] [,3]
## [1,] NULL NULL NULL
## [2,] NULL NULL NULLSelect element from list:
# Select Values, double bracket to select from 2dim list
print('ls_2d[[1,2]]')## [1] "ls_2d[[1,2]]"print(ls_2d[[1,2]])## NULLFor naming two dimensional lists, rowname and colname does not work. Rather, we need to use dimnames. Note that in addition to dimnames, we can continue to have element specific names. Both can co-exist. But note that the element specific names are not preserved after dimension transform, so need to be redefined afterwards.
How to select an element of a two dimensional list:
Neither dimnames nor names are required, but both can be used to select elements.
# Dimensions
it_M <- 3
it_Q <- 4
it_N <- it_M*it_Q
# Initiate an Empty MxQ=N element list
ls_2d_flat_named <- vector(mode = "list", length = it_N)
dim(ls_2d_flat_named) <- c(it_M, it_Q)
# Fill with values
for (it_Q_ctr in seq(1,it_Q)) {
for (it_M_ctr in seq(1,it_M)) {
# linear index
ls_2d_flat_named[[it_M_ctr, it_Q_ctr]] <- (it_Q_ctr-1)*it_M+it_M_ctr
}
}
# Replace row names, note rownames does not work
dimnames(ls_2d_flat_named)[[1]] <- paste0('row',seq(1,it_M))
dimnames(ls_2d_flat_named)[[2]] <- paste0('col',seq(1,it_Q))
# Element Specific Names
names(ls_2d_flat_named) <- paste0('e',seq(1,it_N))
# Convert to Matrix
tb_2d_flat_named <- as_tibble(ls_2d_flat_named) %>% unnest()
mt_2d_flat_named <- as.matrix(tb_2d_flat_named)Print various objects generated above:
# These are not element names, can still name each element
# display
print('ls_2d_flat_named')## [1] "ls_2d_flat_named"print(ls_2d_flat_named)## col1 col2 col3 col4
## row1 1 4 7 10
## row2 2 5 8 11
## row3 3 6 9 12
## attr(,"names")
## [1] "e1" "e2" "e3" "e4" "e5" "e6" "e7" "e8" "e9" "e10" "e11" "e12"print('tb_2d_flat_named')## [1] "tb_2d_flat_named"print(tb_2d_flat_named)
print('mt_2d_flat_named')## [1] "mt_2d_flat_named"print(mt_2d_flat_named)## col1 col2 col3 col4
## [1,] 1 4 7 10
## [2,] 2 5 8 11
## [3,] 3 6 9 12Select elements from list:
# Select elements with with dimnames
ffi_lst2str(ls_2d_flat_named[['row2','col2']])## [1] "ls_2d_flat_named[[\"row2\", \"col2\"]]:=5"# Select elements with element names
ffi_lst2str(ls_2d_flat_named[['e5']])## [1] "ls_2d_flat_named[[\"e5\"]]:=5"# Select elements with index
ffi_lst2str(ls_2d_flat_named[[5]])## [1] "ls_2d_flat_named[[5]]:=5"There are two discrete random variables, generate some random discrete probability mass, name the columns and rows, and then convert to matrix.
set.seed(123)
# Generate prob list
it_Q <- 2
it_M <- 2
ls_2d <- vector(mode = "list", length = it_Q*it_M)
dim(ls_2d) <- c(it_Q, it_M)
# Random joint mass
ar_rand <- runif(it_Q*it_M)
ar_rand <- ar_rand/sum(ar_rand)
# Fill with values
it_ctr <- 0
for (it_Q_ctr in seq(1,it_Q)) {
for (it_M_ctr in seq(1,it_M)) {
# linear index
ls_2d[[it_M_ctr, it_Q_ctr]] <- ar_rand[(it_Q_ctr-1)*it_M+it_M_ctr]
}
}
# Replace row names, note rownames does not work
dimnames(ls_2d)[[1]] <- paste0('E',seq(1,it_M))
dimnames(ls_2d)[[2]] <- paste0('A',seq(1,it_Q))
# rename
ls_prob_joint_E_A <- ls_2d
mt_prob_joint_E_A <- matrix(unlist(ls_prob_joint_E_A), ncol=it_M, byrow=F)
print('ls_prob_joint_E_A')## [1] "ls_prob_joint_E_A"print(ls_prob_joint_E_A)## A1 A2
## E1 0.1214495 0.1727188
## E2 0.3329164 0.3729152print(mt_prob_joint_E_A)## [,1] [,2]
## [1,] 0.1214495 0.1727188
## [2,] 0.3329164 0.3729152Create conditional probabilities: \(F=P(A_1|E_1)\), \(B=P(A_1|E_2)\), \(C=P(E_1|A_1)\), \(D=P(E_1|A_2)\)
fl_F <- mt_prob_joint_E_A[1,1]/sum(mt_prob_joint_E_A[1,])
fl_B <- mt_prob_joint_E_A[2,1]/sum(mt_prob_joint_E_A[2,])
fl_C <- mt_prob_joint_E_A[1,1]/sum(mt_prob_joint_E_A[,1])
fl_D <- mt_prob_joint_E_A[1,2]/sum(mt_prob_joint_E_A[,2])
print(paste0('fl_F=', fl_F, ',fl_B=',fl_B,',fl_C=',fl_C,',fl_D=',fl_D))## [1] "fl_F=0.412857205138471,fl_B=0.471665472604598,fl_C=0.267294503388642,fl_D=0.316546995323062"