knitr::opts_chunk$set(echo = TRUE)
library(FiDEL)

R Markdown

auclist <- (2:48)*0.01 + 0.5
rholist <- (2:18)*0.05
N <- 1000
res <- create_beta_mu(auclist, rholist, N=1)
rhoN <- length(unique(res$rho))
AUCN <- length(unique(res$AUC))

rho <- unique(res$rho)
AUC <- unique(res$AUC)

betaN <- matrix(res[['beta*N']], nrow=rhoN, ncol=AUCN)
muN <- matrix(res[['mu/N']], nrow=rhoN, ncol=AUCN)
#pdf('results/3dplot.pdf', width=12, height=6)
par(mfrow=c(1, 2))
persp(rho, AUC, betaN, theta = 30, phi = 15, shade=.3, ticktype='detailed', expand=.8, scale=T)
persp(rho, AUC, muN, theta = 30, phi = 15, shade=.3, ticktype='detailed', expand=.8, scale=T)
#dev.off()


sungcheolkim78/FiDEL documentation built on Nov. 13, 2024, 7:58 a.m.