Hi hope this helps
#Just enter only these 4 figures and these codes will get you all statistics relating two by two table of diagnostic test data
TP <- 20
FP <- 180
FN <- 10
TN <- 210
TOTP <- TP + FP
TOTN <- FN + TN
TOTDI <- TP + FN
TOTND <- FP + TN
TOTSAM <- TP + FP + FN + TN
mConfusionMatrix <- matrix(c(TP, FN, FP, TN), nrow = 2)
mfisher <- fisher.test(mConfusionMatrix)
mPrevalence <- (TP+FN)/TOTSAM #Prevalence
mAccuracy <- (TP+TN)/TOTSAM # Accuracy
mPPV <- TP/(TP+FP) #Precision positive predictive value
mNPV <- TN/(FN+TN) #negative predictive value
mFDR <- FP/(TP+FP) # False Discovery Rate
mFOR <- FN/(FN+TN) # False omission rate
mTPR <- TP/(TP+FN) #sensitivity or recall
mFPR <- FP/(FP+TN) #fall-out; or false alarm
mFNR <- FN/(TP+FN) #False Negative rate or Miss rate
mTNR <- TN/(FP+TN) #Specificity
mPosLR <- (TP/(TP+FN))/(FP/(FP+TN)) # Positive likelihood ratio
mNegLR <- (FN/(TP+FN)) /(TN/(FP+TN)) # negative likelihood ratio
mDOR <- ((TP/(TP+FN))/(FP/(FP+TN)))/((FN/(TP+FN)) /(TN/(FP+TN)))
mPrevalenceTH <- (sqrt(mTPR*(1-mTNR)) +mTNR -1) / (mTPR + mTNR-1) #Prevalence Threshold
mF1Score <- 2 * (mPPVmTPR) / (mPPV+mTPR) # F1 score
mCHISQ <- (TOTSAM * ((FPFN) - (TPTN))^2) /(TP+FP) / (FN+TN) / (TP+FN) / (FP+TN)
mYULESQ <- ((FPFN) - (TPTN)) / ((FPFN) + (TP*TN))
#find Chi-Square critical value
mChisq_Critical <- qchisq(p=.05, df=1, lower.tail=FALSE)
mPvalChiSq <-pchisq(mCHISQ, df=1, lower.tail=FALSE)
library(epitools)
mepiOR <- epitab(c(TP,FP,FN,TN),method = 'oddsratio')
#odds ratio and confidence interval
mOR <- mepiOR$tab[10] #Odds Ratio
mORCIL <- mepiOR$tab[12]
mORCIU <- mepiOR$tab[14]
mORPValue <-round(mepiOR$tab[16],4)
#Relative risk and confidence interval
mRR <- (TP / (TP+FP)) / (FN / (FN+TN)) #Relative Risk
mRRSE <- sqrt((1/TP - 1/(TP+FP))+(1/FN - 1/(FN+TN)) ) #Standard Error for Odds Ratio
mRRCIL <- exp(log(mRR)- (1.95mRRSE))
mRRCIU <- exp(log(mRR)+ (1.95mRRSE))
mRRZSTAT <- log(mRR)/mRRSE
mRRPValue <- round(2*pnorm(-abs(mRRZSTAT)),4)
mEER <- TP / (TP+FP) #Experimental Event Rate
mCER <- FN / (FN+TN) #Control Event Rate
mRRR <- (mCER-mEER)/mCER #Relative Risk Reduction
mRRI <- (mEER-mCER)/mCER #Relative Risk Increase
mARR <- mCER-mEER #Absolute Risk Reduction
mARI <- mEER-mCER #Absolute Risk Increase
mNNT <- 1/mARR #Number Needed to Treat
mNNH <- 1/mARI #Number Needed to Harm