knitr::opts_chunk$set(echo = TRUE) library(FiDEL)
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()
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.