SEQanswers

SEQanswers (http://seqanswers.com/forums/index.php)
-   General (http://seqanswers.com/forums/forumdisplay.php?f=16)
-   -   Using R package to implement Bayesian phase I/II dose-finding design for three outcom (http://seqanswers.com/forums/showthread.php?t=100829)

Ran20 06-07-2021 10:17 AM

Using R package to implement Bayesian phase I/II dose-finding design for three outcom
 
Dear all,

I am trying to implement the work done by Suyu Liu, "A Bayesian Phase I/II Trial Design for Immunotherapy", using R, since the code attached with that work takes a lot of time (more than 20 hours, and the code not complete). So that I tried to use trialr package since it used rstan, but this package allowed me to use two outcomes ( toxicity, efficacy ) and the work of Liu used three outcomes (immune response, toxicity, and efficacy).

I tried to use the R package trialr for two outcomes using the utility for sensitivity analysis written in R code below (table 1 in the article), I want to see if I used the correct utility and to see how to add a third outcome ( immune response ) to the model

and thanks in advance.

Here is the work of Liu:

https://www.tandfonline.com/doi/supp...260?scroll=top

````

### My code###
rm(list = ls())
library(trialr)

#Utility
Uti <- array(0,c(2,3,2)) # order: tox, eff, immuno
Uti[,,1] <- matrix(c(0,0,50,10,80,35),nrow=2)
Uti[,,2] <- matrix(c(5,0,70,20,100,45),nrow=2)
N.max= 60 # patients
outcomes <- '1NNN 2NNT 3NNT 4NNN 5NTN'
doses = c(.1,.3,.5,.7,.9)


fit <- stan_efftox(outcomes,
real_doses =doses,
efficacy_hurdle = 0.5, toxicity_hurdle = 0.3,
p_e = 0.1, p_t = 0.1,
eff0 = 0.5, tox1 = 0.65,
eff_star = 0.7, tox_star = 0.25,
alpha_mean = -7.9593, alpha_sd = 3.5487,
beta_mean = 1.5482, beta_sd = 3.5018,
gamma_mean = 0.7367, gamma_sd = 2.5423,
zeta_mean = 3.4181, zeta_sd = 2.4406,
eta_mean = 0, eta_sd = 0.2,
psi_mean = 0, psi_sd = 1,
seed = 123)

ndoses <- length(fit$prob_tox)
plot(1:ndoses, fit$prob_tox, type="b", pch=19, xlab="Dose level", ylab="Probability of toxicity", ylim=c(0,max(fit$prob_tox) + 0.15), col="green")
points(1:ndoses,fit$prob_eff, type="b", pch=18, col="blue")
abline(h=0.3, lwd=2, lty=4, col = "red")
legend(1, 0.4, legend=c("Toxicity", "Effecacy"),
col=c("green", "blue"), lty=1:2, cex=0.8)

````

monicall89 06-23-2021 03:44 AM

Thank you so much for this information; itís really a very good information and very authentic and useful information for us. I really appreciate you have share this information.
If you know more about this information just click here ptcl speed test and check more information with more reasonable and authentic way.


All times are GMT -8. The time now is 01:22 AM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2021, vBulletin Solutions, Inc.