Nothing
## ---- include = FALSE---------------------------------------------------------
knitr::opts_chunk$set(
collapse = TRUE,
comment = "#>"
)
## ----setup--------------------------------------------------------------------
library(CRTConjoint)
## -----------------------------------------------------------------------------
data("immigrationdata")
## -----------------------------------------------------------------------------
form = formula("Y ~ FeatEd + FeatGender + FeatCountry + FeatReason + FeatJob +
FeatExp + FeatPlans + FeatTrips + FeatLang + ppage + ppeducat + ppethm + ppgender")
left = colnames(immigrationdata)[1:9]
right = colnames(immigrationdata)[10:18]
left; right
## ---- eval = FALSE------------------------------------------------------------
# education_test = CRT_pval(formula = form, data = immigrationdata, X = "FeatEd",
# left = left, right = right, non_factor = "ppage", B = 100, analysis = 2)
# education_test$p_val
## ---- eval = FALSE------------------------------------------------------------
# constraint_randomization = list() # (Job has dependent randomization scheme)
# constraint_randomization[["FeatJob"]] = c("Financial analyst","Computer programmer",
# "Research scientist","Doctor")
# constraint_randomization[["FeatEd"]] = c("Equivalent to completing two years of
# college in the US", "Equivalent to completing a graduate degree in the US",
# "Equivalent to completing a college degree in the US")
## ---- eval = FALSE------------------------------------------------------------
# job_test = CRT_pval(formula = form, data = immigrationdata, X = "FeatJob",
# left = left, right = right, design = "Constrained Uniform",
# constraint_randomization = constraint_randomization, non_factor = "ppage", B = 100)
# job_test$p_val
## ---- eval = FALSE------------------------------------------------------------
# profileorder_test = CRT_profileordereffect(formula = form, data = immigrationdata,
# left = left, right = right, B = 100)
# profileorder_test$p_val
## ---- eval = FALSE------------------------------------------------------------
# resample_func_immigration = function(x, seed = sample(c(0, 1000), size = 1), left_idx, right_idx) {
# set.seed(seed)
# df = x[, c(left_idx, right_idx)]
# variable = colnames(x)[c(left_idx, right_idx)]
# len = length(variable)
# resampled = list()
# n = nrow(df)
# for (i in 1:len) {
# var = df[, variable[i]]
# lev = levels(var)
# resampled[[i]] = factor(sample(lev, size = n, replace = TRUE))
# }
#
# resampled_df = data.frame(resampled[[1]])
# for (i in 2:len) {
# resampled_df = cbind(resampled_df, resampled[[i]])
# }
# colnames(resampled_df) = colnames(df)
#
# #escape persecution was dependently randomized
# country_1 = resampled_df[, "FeatCountry"]
# country_2 = resampled_df[, "FeatCountry_2"]
# i_1 = which((country_1 == "Iraq" | country_1 == "Sudan" | country_1 == "Somalia"))
# i_2 = which((country_2 == "Iraq" | country_2 == "Sudan" | country_2 == "Somalia"))
#
# reason_1 = resampled_df[, "FeatReason"]
# reason_2 = resampled_df[, "FeatReason_2"]
# levs = levels(reason_1)
# r_levs = levs[c(2,3)]
#
# reason_1 = sample(r_levs, size = n, replace = TRUE)
#
# reason_1[i_1] = sample(levs, size = length(i_1), replace = TRUE)
#
# reason_2 = sample(r_levs, size = n, replace = TRUE)
#
# reason_2[i_2] = sample(levs, size = length(i_2), replace = TRUE)
#
# resampled_df[, "FeatReason"] = reason_1
# resampled_df[, "FeatReason_2"] = reason_2
#
# #profession high skill fix
# educ_1 = resampled_df[, "FeatEd"]
# educ_2 = resampled_df[, "FeatEd_2"]
# i_1 = which((educ_1 == "Equivalent to completing two years of college in the US" |
# educ_1 == "Equivalent to completing a college degree in the US" |
# educ_1 == "Equivalent to completing a graduate degree in the US"))
# i_2 = which((educ_2 == "Equivalent to completing two years of college in the US" |
# educ_2 == "Equivalent to completing a college degree in the US" |
# educ_2 == "Equivalent to completing a graduate degree in the US"))
#
#
# job_1 = resampled_df[, "FeatJob"]
# job_2 = resampled_df[, "FeatJob_2"]
# levs = levels(job_1)
# # take out computer programmer, doctor, financial analyst, and research scientist
# r_levs = levs[-c(2,4,5, 9)]
#
# job_1 = sample(r_levs, size = n, replace = TRUE)
#
# job_1[i_1] = sample(levs, size = length(i_1), replace = TRUE)
#
# job_2 = sample(r_levs, size = n, replace = TRUE)
#
# job_2[i_2] = sample(levs, size = length(i_2), replace = TRUE)
#
# resampled_df[, "FeatJob"] = job_1
# resampled_df[, "FeatJob_2"] = job_2
#
# resampled_df[colnames(resampled_df)] = lapply(resampled_df[colnames(resampled_df)], factor )
#
# return(resampled_df)
# }
## ---- eval = FALSE------------------------------------------------------------
# carryover_df = immigrationdata
# own_resamples = list()
# B = 100
# for (i in 1:B) {
# newdf = resample_func_immigration(carryover_df, left_idx = 1:9, right_idx = 10:18, seed = i)
# own_resamples[[i]] = newdf
# }
## ---- eval = FALSE------------------------------------------------------------
# J = 5
# carryover_df$task = rep(1:J, nrow(carryover_df)/J)
#
# carryover_test = CRT_carryovereffect(formula = form, data = carryover_df, left = left,
# right = right, task = "task", supplyown_resamples = own_resamples, B = B)
# carryover_test$p_val
## ---- eval = FALSE------------------------------------------------------------
# fatigue_df = immigrationdata
# fatigue_df$task = rep(1:J, nrow(fatigue_df)/J)
# fatigue_df$respondent = rep(1:(nrow(fatigue_df)/J), each = J)
#
# fatigue_test = CRT_fatigueeffect(formula = form, data = fatigue_df, left = left,
# right = right, task = "task", respondent = "respondent", B = 100)
# fatigue_test$p_val
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.