Homec4science
Diffusion hcv_modeling (HCV_v2) Browse

R3744/FP_Scenarios.R
424c6b36b4adHCV_v2

Last Change
Raw File

FP_Scenarios.R
View Options

FP_scenario <- 1
source("index_block_fs_stage.r") # it identify that in which block and stage we are in
######################################################################################
##########################################################################################
if (FP_scenario == 1)
{
### FIBROSIS PROGRESSION VOGEL
fp.fun = function(t, bl, history, r)
{
index.fun = stage_identifier(r, bl, history)
index_block = index.fun[1]
index_fs = index.fun[2]
state = index.fun[3]
#######################################
if (is.element(state, c(31:36)) )
index_fs = 6 # DC
######################################
if (is.element(state, c(37:42)) )
index_fs = 7 # HCC
####################################
if (is.element(state, c(43:48)) )
index_fs = 8 #LT
####################################
#####################################
res<- 0 + 0 * t
rr<- 0
{
ifelse(index_block == 1 || index_block == 2 ||index_block == 3 ||index_block == 4 || index_block == 5 , rr<- 1, rr<- 0.1 ) # detectable VL/undetectable VL
{
###################
r = switch (index_fs,
c(1,1.16,1.33),
c(1,1.3,2.22),
c(1,1.3,2.22),
c(1,1.16,4)
)
#r = c(1,1,1.33)
###################
if (index_fs == 1)
{
if (bl["Gender"] == 1) #male
{
res[(t + sum(history) + bl["Age"])< 20] <- { (rr * 0.045 * r[bl["Alcohol"] + 1]) + 0 * t }
res[(t + sum(history)+bl["Age"])>= 20] <- { (rr*0.037 * r[bl["Alcohol"] + 1]) + 0 * t }
res[(t + sum(history)+bl["Age"])>= 30] <-{ (rr*0.027 * r[bl["Alcohol"] + 1]) + 0 * t }
res[(t + sum(history)+bl["Age"])>= 40] <- { (rr * 0.099 * r[bl["Alcohol"]+1]) + 0 * t }
res[(t + sum(history)+bl["Age"])>= 50] <- { (rr * 0.121 * r[bl["Alcohol"] + 1]) + 0 * t }
res[(t + sum(history)+bl["Age"])>= 60] <- { (rr * 0.138 * r[bl["Alcohol"] + 1]) + 0 * t }
res[(t + sum(history)+bl["Age"])>= 70] <- { (rr * 0.155 * r[bl["Alcohol"] + 1]) + 0 * t }
res[(t + sum(history)+bl["Age"])>= 80] <- { (rr * 0.127 * r[bl["Alcohol"] + 1]) + 0 * t }
}
if (bl["Gender"] == 0) # female
{
res[(t + sum(history) + bl["Age"])< 20] <- { (rr * 0.038 * r[bl["Alcohol"] + 1]) + 0 * t }
res[(t + sum(history)+bl["Age"])>= 20] <- { (rr*0.031 * r[bl["Alcohol"] + 1]) + 0 * t }
res[(t + sum(history)+bl["Age"])>= 30] <-{ (rr*0.022 * r[bl["Alcohol"] + 1]) + 0 * t }
res[(t + sum(history)+bl["Age"])>= 40] <- { (rr * 0.082 * r[bl["Alcohol"]+1]) + 0 * t }
res[(t + sum(history)+bl["Age"])>= 50] <- { (rr * 0.102 * r[bl["Alcohol"] + 1]) + 0 * t }
res[(t + sum(history)+bl["Age"])>= 60] <- { (rr * 0.115 * r[bl["Alcohol"] + 1]) + 0 * t }
res[(t + sum(history)+bl["Age"])>= 70] <- { (rr * 0.129 * r[bl["Alcohol"] + 1]) + 0 * t }
res[(t + sum(history)+bl["Age"])>= 80] <- { (rr * 0.106 * r[bl["Alcohol"] + 1]) + 0 * t }
}
}#(index_fs == 1)
######################
if (index_fs == 2)
{
if (bl["Gender"] == 1)
{
res[(t + sum(history) + bl["Age"])< 20] <- { (rr * 0.033 * r[bl["Alcohol"] + 1]) + 0 * t }
res[(t + sum(history)+bl["Age"])>= 20] <- { (rr*0.027 * r[bl["Alcohol"] + 1]) + 0 * t }
res[(t + sum(history)+bl["Age"])>= 30] <-{ (rr*0.019 * r[bl["Alcohol"] + 1]) + 0 * t }
res[(t + sum(history)+bl["Age"])>= 40] <- { (rr * 0.072 * r[bl["Alcohol"]+1]) + 0 * t }
res[(t + sum(history)+bl["Age"])>= 50] <- { (rr * 0.088 * r[bl["Alcohol"] + 1]) + 0 * t }
res[(t + sum(history)+bl["Age"])>= 60] <- { (rr * 0.100 * r[bl["Alcohol"] + 1]) + 0 * t }
res[(t + sum(history)+bl["Age"])>= 70] <- { (rr * 0.112 * r[bl["Alcohol"] + 1]) + 0 * t }
res[(t + sum(history)+bl["Age"])>= 80] <- { (rr * 0.130 * r[bl["Alcohol"] + 1]) + 0 * t }
}
if (bl["Gender"] == 0)
{
res[(t + sum(history) + bl["Age"])< 20] <- { (rr * 0.028 * r[bl["Alcohol"] + 1]) + 0 * t }
res[(t + sum(history)+bl["Age"])>= 20] <- { (rr*0.022 * r[bl["Alcohol"] + 1]) + 0 * t }
res[(t + sum(history)+bl["Age"])>= 30] <-{ (rr*0.016 * r[bl["Alcohol"] + 1]) + 0 * t }
res[(t + sum(history)+bl["Age"])>= 40] <- { (rr * 0.060 * r[bl["Alcohol"]+1]) + 0 * t }
res[(t + sum(history)+bl["Age"])>= 50] <- { (rr * 0.074 * r[bl["Alcohol"] + 1]) + 0 * t }
res[(t + sum(history)+bl["Age"])>= 60] <- { (rr * 0.083 * r[bl["Alcohol"] + 1]) + 0 * t }
res[(t + sum(history)+bl["Age"])>= 70] <- { (rr * 0.094 * r[bl["Alcohol"] + 1]) + 0 * t }
res[(t + sum(history)+bl["Age"])>= 80] <- { (rr * 0.077 * r[bl["Alcohol"] + 1]) + 0 * t }
}
}#(index_fs == 2)
######################
if (index_fs == 3)
{
if (bl["Gender"] == 1)
{
res[(t + sum(history) + bl["Age"])< 20] <- { (rr * 0.047 * r[bl["Alcohol"] + 1]) + 0 * t }
res[(t + sum(history)+bl["Age"])>= 20] <- { (rr*0.038 * r[bl["Alcohol"] + 1]) + 0 * t }
res[(t + sum(history)+bl["Age"])>= 30] <-{ (rr*0.028 * r[bl["Alcohol"] + 1]) + 0 * t }
res[(t + sum(history)+bl["Age"])>= 40] <- { (rr * 0.102 * r[bl["Alcohol"]+1]) + 0 * t }
res[(t + sum(history)+bl["Age"])>= 50] <- { (rr * 0.124 * r[bl["Alcohol"] + 1]) + 0 * t }
res[(t + sum(history)+bl["Age"])>= 60] <- { (rr * 0.141 * r[bl["Alcohol"] + 1]) + 0 * t }
res[(t + sum(history)+bl["Age"])>= 70] <- { (rr * 0.159 * r[bl["Alcohol"] + 1]) + 0 * t }
res[(t + sum(history)+bl["Age"])>= 80] <- { (rr * 0.130 * r[bl["Alcohol"] + 1]) + 0 * t }
}
if (bl["Gender"] == 0)
{
res[(t + sum(history) + bl["Age"])< 20] <- { (rr * 0.039 * r[bl["Alcohol"] + 1]) + 0 * t }
res[(t + sum(history)+bl["Age"])>= 20] <- { (rr*0.031 * r[bl["Alcohol"] + 1]) + 0 * t }
res[(t + sum(history)+bl["Age"])>= 30] <-{ (rr*0.023 * r[bl["Alcohol"] + 1]) + 0 * t }
res[(t + sum(history)+bl["Age"])>= 40] <- { (rr * 0.085 * r[bl["Alcohol"]+1]) + 0 * t }
res[(t + sum(history)+bl["Age"])>= 50] <- { (rr * 0.104 * r[bl["Alcohol"] + 1]) + 0 * t }
res[(t + sum(history)+bl["Age"])>= 60] <- { (rr * 0.118 * r[bl["Alcohol"] + 1]) + 0 * t }
res[(t + sum(history)+bl["Age"])>= 70] <- { (rr * 0.132 * r[bl["Alcohol"] + 1]) + 0 * t }
res[(t + sum(history)+bl["Age"])>= 80] <- { (rr * 0.109 * r[bl["Alcohol"] + 1]) + 0 * t }
}
}#(index_fs == 3)
######################
if (index_fs == 4)
{
if (bl["Gender"] == 1)
{
res[(t + sum(history) + bl["Age"])< 20] <- { (rr * 0.006 * r[bl["Alcohol"] + 1]) + 0 * t }
res[(t + sum(history)+bl["Age"])>= 20] <- { (rr*0.018 * r[bl["Alcohol"] + 1]) + 0 * t }
res[(t + sum(history)+bl["Age"])>= 30] <-{ (rr*0.040 * r[bl["Alcohol"] + 1]) + 0 * t }
res[(t + sum(history)+bl["Age"])>= 40] <- { (rr * 0.063 * r[bl["Alcohol"]+1]) + 0 * t }
res[(t + sum(history)+bl["Age"])>= 50] <- { (rr * 0.034 * r[bl["Alcohol"] + 1]) + 0 * t }
res[(t + sum(history)+bl["Age"])>= 60] <- { (rr * 0.070 * r[bl["Alcohol"] + 1]) + 0 * t }
res[(t + sum(history)+bl["Age"])>= 70] <- { (rr * 0.136 * r[bl["Alcohol"] + 1]) + 0 * t }
res[(t + sum(history)+bl["Age"])>= 80] <- { (rr * 0.136 * r[bl["Alcohol"] + 1]) + 0 * t }
}
if (bl["Gender"] == 0)
{
res[(t + sum(history) + bl["Age"])< 20] <- { (rr * 0.004 * r[bl["Alcohol"] + 1]) + 0 * t }
res[(t + sum(history)+bl["Age"])>= 20] <- { (rr*0.015 * r[bl["Alcohol"] + 1]) + 0 * t }
res[(t + sum(history)+bl["Age"])>= 30] <-{ (rr*0.033 * r[bl["Alcohol"] + 1]) + 0 * t }
res[(t + sum(history)+bl["Age"])>= 40] <- { (rr * 0.053 * r[bl["Alcohol"]+1]) + 0 * t }
res[(t + sum(history)+bl["Age"])>= 50] <- { (rr * 0.028 * r[bl["Alcohol"] + 1]) + 0 * t }
res[(t + sum(history)+bl["Age"])>= 60] <- { (rr * 0.059 * r[bl["Alcohol"] + 1]) + 0 * t }
res[(t + sum(history)+bl["Age"])>= 70] <- { (rr * 0.113 * r[bl["Alcohol"] + 1]) + 0 * t }
res[(t + sum(history)+bl["Age"])>= 80] <- { (rr * 0.113 * r[bl["Alcohol"] + 1]) + 0 * t }
}
}#(index_fs == 4)
#######################
}
}
return(res)
} # FIBROSIS PROGRESSION VOGEL
}# end if (FP_scenario == 1)
#######################################################################################
##########################################################################################
c4science · Help