inst/data-raw/process/PMID23139348_Korpos-2012/process.R

In avucoh/nPOD: What the Package Does (One Line, Title Case)

# The input data is from Supp Table 3

# Calculate weighted percentages for 6052 
islet.n <- c(32, 41, 58, 45)
InsPos <- c(0, 4, 2, 3)
InsPos.prct <- InsPos/islet.n * 100

InsPos.prct <- round(weighted.mean(InsPos.prct, islet.n), 1)
InsNeg.prct <- 100 - InsPos.prct

# Percentage for 6070, which is more uncertain given that only one sample was examined
InsPos.prct.6070 <- 2/21 * 100
InsNeg.prct.6070 <- 19/21 * 100