shortASAQConstructsScore_combinedData,
shortASAQscores_combinedData,
probs10th,
constructASAQdesc$CABR)
probs20th <- c(.20, .40, .60, .80)
minSampleConfIntervLongASAQConstructScore_combinedData_p20th <-
calcMinSampleConfIntervASAQConstructScore(
name29ASAs,
longASAQConstructsScore_combinedData,
longASAQscores_combinedData,
probs20th,
constructASAQdesc$CABR)
minSampleConfIntervShortASAQConstructScore_combinedData_p20th <-
calcMinSampleConfIntervASAQConstructScore(
name29ASAs,
shortASAQConstructsScore_combinedData,
shortASAQscores_combinedData,
probs20th,
constructASAQdesc$CABR)
probs25th <- c(.25, .50, .75)
minSampleConfIntervLongASAQConstructScore_combinedData_p25th <-
calcMinSampleConfIntervASAQConstructScore(
name29ASAs,
longASAQConstructsScore_combinedData,
longASAQscores_combinedData,
probs25th,
constructASAQdesc$CABR)
minSampleConfIntervShortASAQConstructScore_combinedData_p25th <-
calcMinSampleConfIntervASAQConstructScore(
name29ASAs,
shortASAQConstructsScore_combinedData,
shortASAQscores_combinedData,
probs25th,
constructASAQdesc$CABR)
# Calculate maximum sample size for evaluating an ASA based on the confident interval
maxSampleSizeConfidenceIntervalASAQConstructScore_combinedData <-
getMaxSampleSizeConfidenceIntervalASAQConstructScore(
minSampleConfIntervLongASAQConstructScore_combinedData_p5th,
minSampleConfIntervShortASAQConstructScore_combinedData_p5th,
minSampleConfIntervLongASAQConstructScore_combinedData_p10th,
minSampleConfIntervShortASAQConstructScore_combinedData_p10th,
minSampleConfIntervLongASAQConstructScore_combinedData_p20th,
minSampleConfIntervShortASAQConstructScore_combinedData_p20th,
minSampleConfIntervLongASAQConstructScore_combinedData_p25th,
minSampleConfIntervShortASAQConstructScore_combinedData_p25th)
maxSampleSizeConfidenceIntervalASAQConstructScore_combinedData
# registered participants active in the past 90 days (retrieved on 17-07-2024)
activeTotal = 216627
man = 65603 # Man including Trans Male/Trans Man
woman = 94475 # Woman including Trans Female/Trans Woman
nonBinary = 3340 #Non-binary (would like to give more detailed)
notSay = 1138 #Rather not say
genderUnknown = activeTotal - man - woman - nonBinary - notSay
age18to25 = 59734
age26to35 = 79101
age36to45 = 40597
age46to55 = 22222
age56to65 = 10780
ageOlder65 = 4063
ageUnknown = activeTotal-age18to25-age26to35-age36to45-age46to55- age56to65-ageOlder65
calcTotalLongASAQscore(longASAQscores_combinedData)
showTable(relFreqLongASAQConstructsScore_combinedData, constructASAQdesc$CABR,
constructASAQdesc$CNAME, 0, 1, 1, 1, 1, 1, 1)
showTable(relFreqShortASAQConstructsScore_combinedData, constructASAQdesc$CABR,
constructASAQdesc$CNAME, 0, 1, 1, 1, 1, 1, 1)
# Get the number of sample of each agent
getSampleASA(name29ASAs, combinedData)
showTable(longASAQscores_combinedData, constructASAQdesc$CABR,
constructASAQdesc$CNAME, 1, 1, 1, 1, 1, 1, 1)
showTable(shortASAQscores_combinedData, constructASAQdesc$CABR,
constructASAQdesc$CNAME, 1, 1, 1, 1, 1, 1, 1)
showTable(percentileLongASAQScore_combinedData, constructASAQdesc$CABR,
constructASAQdesc$CNAME, 0, 1, 1, 1, 1, 1, 1)
showTable(percentileShortASAQScore_combinedData, constructASAQdesc$CABR,
constructASAQdesc$CNAME, 0, 1, 1, 1, 1, 1, 1)
showTable(percentileDiffLongTwoASAQscore_combinedData, constructASAQdesc$CABR,
constructASAQdesc$CNAME, 0, 1, 1, 1, 1, 1, 1)
showTable(percentileDiffShortTwoASAQscore_combinedData, constructASAQdesc$CABR,
constructASAQdesc$CNAME, 0, 1, 1, 1, 1, 1, 1)
showTable(percentileEffectSizeLongASAQConstructScore_combinedData,
constructASAQdesc$CABR, constructASAQdesc$CNAME, 0, 1, 1, 1, 1, 1, 1)
showTable(percentileEffectSizeShortASAQConstructScore_combinedData,
constructASAQdesc$CABR, constructASAQdesc$CNAME, 0, 1, 1, 1, 1, 1, 1)
minSamplePowerAnalysisASAQConstructScore_combinedData <-
cbind(minSamplePowerAnalysisLongASAQConstructScore_combinedData,
minSamplePowerAnalysisShortASAQConstructScore_combinedData[, 2:ncol(
minSamplePowerAnalysisShortASAQConstructScore_combinedData)])
names(minSamplePowerAnalysisASAQConstructScore_combinedData) <-
c("CID", "LONG_Small_q25th", "LONG_Medium_q50th", "LONG_Large_q75th",
"SHORT_Small_q25th", "SHORT_Medium_q50th", "SHORT_Large_q75th")
showTable(minSamplePowerAnalysisASAQConstructScore_combinedData,
constructASAQdesc$CABR, constructASAQdesc$CNAME, 0, 1, 0, 0, 0, 0, 1)
minSampleConfIntervASAQConstructScore_combinedData_p5th <-
cbind(minSampleConfIntervLongASAQConstructScore_combinedData_p5th,
minSampleConfIntervShortASAQConstructScore_combinedData_p5th[, 2:ncol(
minSampleConfIntervShortASAQConstructScore_combinedData_p5th)])
names(minSampleConfIntervASAQConstructScore_combinedData_p5th) <-
c("CID", "LONG_90CI", "LONG_95CI", "LONG_97.5CI", "LONG_99CI",
"SHORT_90CI","SHORT_95CI", "SHORT_97.5CI", "SHORT_99CI")
showTable(minSampleConfIntervASAQConstructScore_combinedData_p5th,
constructASAQdesc$CABR, constructASAQdesc$CNAME, 0, 1, 0, 0, 0, 0, 1)
minSampleConfIntervASAQConstructScore_combinedData_p10th <-
cbind(minSampleConfIntervLongASAQConstructScore_combinedData_p10th,
minSampleConfIntervShortASAQConstructScore_combinedData_p10th[, 2:ncol(
minSampleConfIntervShortASAQConstructScore_combinedData_p10th)])
names(minSampleConfIntervASAQConstructScore_combinedData_p10th) <-
c("CID", "LONG_90CI", "LONG_95CI", "LONG_97.5CI", "LONG_99CI",
"SHORT_90CI","SHORT_95CI", "SHORT_97.5CI", "SHORT_99CI")
showTable(minSampleConfIntervASAQConstructScore_combinedData_p10th,
constructASAQdesc$CABR, constructASAQdesc$CNAME, 0, 1, 0, 0, 0, 0, 1)
minSampleConfIntervASAQConstructScore_combinedData_p20th <-
cbind(minSampleConfIntervLongASAQConstructScore_combinedData_p20th,
minSampleConfIntervShortASAQConstructScore_combinedData_p20th[, 2:ncol(
minSampleConfIntervShortASAQConstructScore_combinedData_p20th)])
names(minSampleConfIntervASAQConstructScore_combinedData_p20th) <-
c("CID", "LONG_90CI", "LONG_95CI", "LONG_97.5CI", "LONG_99CI",
"SHORT_90CI","SHORT_95CI", "SHORT_97.5CI", "SHORT_99CI")
showTable(minSampleConfIntervASAQConstructScore_combinedData_p20th,
constructASAQdesc$CABR, constructASAQdesc$CNAME, 0, 1, 0, 0, 0, 0, 1)
minSampleConfIntervASAQConstructScore_combinedData_p25th <-
cbind(minSampleConfIntervLongASAQConstructScore_combinedData_p25th,
minSampleConfIntervShortASAQConstructScore_combinedData_p25th[, 2:ncol(
minSampleConfIntervShortASAQConstructScore_combinedData_p25th)])
names(minSampleConfIntervASAQConstructScore_combinedData_p25th) <-
c("CID", "LONG_90CI", "LONG_95CI", "LONG_97.5CI", "LONG_99CI",
"SHORT_90CI","SHORT_95CI", "SHORT_97.5CI", "SHORT_99CI")
showTable(minSampleConfIntervASAQConstructScore_combinedData_p25th,
constructASAQdesc$CABR, constructASAQdesc$CNAME, 0, 1, 0, 0, 0, 0, 1)
data_mean <-shortASAQscores_combinedData
d_sum<- calcTotalLongASAQscore(shortASAQscores_combinedData)
abrConstructs = constructASAQdesc$CABR
nameConstructs =constructASAQdesc$CNAME
d_samples <- getSampleASA(name29ASAs, combinedData)
d_mean <- convertToNumeric(shortASAQscores_combinedData, 2)
d_sd <- convertToNumeric(shortASAQscoreSD_combinedData, 2)
relFreqLongASAQConstructsScore_combinedData
View(relFreqLongASAQConstructsScore_combinedData)
showTable(relFreqLongASAQConstructsScore_combinedData, constructASAQdesc$CABR,
constructASAQdesc$CNAME, 0, 1, 1, 1, 1, 1, 1)
t=showTable(relFreqLongASAQConstructsScore_combinedData, constructASAQdesc$CABR,
constructASAQdesc$CNAME, 0, 1, 1, 1, 1, 1, 1)
showTable <- function(data_ori, abrConstructs, nameConstructs, type,
showDesc, showMean, showSD, showMedian, showMin, showMax){
dat <- data.frame(matrix(nrow = length(abrConstructs), ncol = 0))
data <- data.frame()
if (type == 0)
data <- convertToNumeric(data_ori, 2)
else if (type == 1){
headers <- data_ori[, 1]
data_ori <- convertToNumeric(data_ori, 2)
for(i in 1:length(headers)){
for (j in 1:length(abrConstructs)){
data[j, headers[i]] <- data_ori[i, j]
}
}
colnames(data) <- headers
}
if (type == 2){
data <- convertToNumeric(data_ori, 2)
n = ncol(data)/2
h1 <- paste("LONG", colnames(data)[1:n], sep = "_")
n = n + 1
h2 <- paste("SHORT", colnames(data)[n:ncol(data)], sep = "_")
colnames(data) <- c(h1,h2)
}
data <- round(data,2)
dat[, "CABR"] <- abrConstructs
dat[, "CNAME"] <- nameConstructs
dat <- cbind(dat, format(data, nsmall=2))
print(dat)
if (showDesc == 1){
# Do data description
i = 1
desc <- data.frame(matrix(nrow = 0, ncol = ncol(data) + 2))
if (showMean == 1){
desc [i, ] <- c("","Mean:", format(round(apply(data, 2, mean, na.rm=TRUE), 2)),
nsmall = 2)
i = i + 1
}
if (showSD == 1){
desc [i, ] <- c("","SD:", format(round(apply(data, 2, sd, na.rm=TRUE), 2)),
nsmall = 2)
i = i + 1
}
if (showMedian == 1){
desc [i, ] <- c("","Median:", format(round(apply(data, 2, median, na.rm=TRUE), 2)),
nsmall = 2)
i = i + 1
}
if (showMin == 1){
desc [i, ] <- c("","Min:", format(apply(data, 2, min, na.rm=TRUE), nsmall = 2))
i = i + 1
}
if (showMax == 1){
desc [i, ] <- c("","Max:", format(apply(data, 2, max, na.rm=TRUE), nsmall = 2))
i = i + 1
}
# Display data description
d <- desc[,2:ncol(desc)]
colnames(d) <- c("Description", colnames(data))
print(d)
}
}
t=showTable(relFreqLongASAQConstructsScore_combinedData, constructASAQdesc$CABR,
constructASAQdesc$CNAME, 0, 1, 1, 1, 1, 1, 1)
showTable <- function(data_ori, abrConstructs, nameConstructs, type,
showDesc, showMean, showSD, showMedian, showMin, showMax){
dat <- data.frame(matrix(nrow = length(abrConstructs), ncol = 0))
data <- data.frame()
if (type == 0)
data <- convertToNumeric(data_ori, 2)
else if (type == 1){
headers <- data_ori[, 1]
data_ori <- convertToNumeric(data_ori, 2)
for(i in 1:length(headers)){
for (j in 1:length(abrConstructs)){
data[j, headers[i]] <- data_ori[i, j]
}
}
colnames(data) <- headers
}
if (type == 2){
data <- convertToNumeric(data_ori, 2)
n = ncol(data)/2
h1 <- paste("LONG", colnames(data)[1:n], sep = "_")
n = n + 1
h2 <- paste("SHORT", colnames(data)[n:ncol(data)], sep = "_")
colnames(data) <- c(h1,h2)
}
data <- round(data,2)
dat[, "CABR"] <- abrConstructs
dat[, "CNAME"] <- nameConstructs
dat <- cbind(dat, format(data, nsmall=2))
print(dat)
return(dat)
}
t=showTable(relFreqLongASAQConstructsScore_combinedData, constructASAQdesc$CABR,
constructASAQdesc$CNAME, 0, 1, 1, 1, 1, 1, 1)
showTable <- function(data_ori, abrConstructs, nameConstructs, type,
showDesc, showMean, showSD, showMedian, showMin, showMax){
dat <- data.frame(matrix(nrow = length(abrConstructs), ncol = 0))
data <- data.frame()
if (type == 0)
data <- convertToNumeric(data_ori, 2)
else if (type == 1){
headers <- data_ori[, 1]
data_ori <- convertToNumeric(data_ori, 2)
for(i in 1:length(headers)){
for (j in 1:length(abrConstructs)){
data[j, headers[i]] <- data_ori[i, j]
}
}
colnames(data) <- headers
}
if (type == 2){
data <- convertToNumeric(data_ori, 2)
n = ncol(data)/2
h1 <- paste("LONG", colnames(data)[1:n], sep = "_")
n = n + 1
h2 <- paste("SHORT", colnames(data)[n:ncol(data)], sep = "_")
colnames(data) <- c(h1,h2)
}
data <- round(data,2)
dat[, "CABR"] <- abrConstructs
dat[, "CNAME"] <- nameConstructs
dat <- cbind(dat, format(data, nsmall=2))
print(dat)
if (showDesc == 1){
# Do data description
i = 1
desc <- data.frame(matrix(nrow = 0, ncol = ncol(dat)))
if (showMean == 1){
desc [i, ] <- c("","Mean:", format(round(apply(data, 2, mean, na.rm=TRUE), 2)),
nsmall = 2)
i = i + 1
}
if (showSD == 1){
desc [i, ] <- c("","SD:", format(round(apply(data, 2, sd, na.rm=TRUE), 2)),
nsmall = 2)
i = i + 1
}
if (showMedian == 1){
desc [i, ] <- c("","Median:", format(round(apply(data, 2, median, na.rm=TRUE), 2)),
nsmall = 2)
i = i + 1
}
if (showMin == 1){
desc [i, ] <- c("","Min:", format(apply(data, 2, min, na.rm=TRUE), nsmall = 2))
i = i + 1
}
if (showMax == 1){
desc [i, ] <- c("","Max:", format(apply(data, 2, max, na.rm=TRUE), nsmall = 2))
i = i + 1
}
# for return purpose
colnames(desc) <- names(dat)
dat <- rbind(dat, desc)
# for printing purpose
colnames(desc) <- c("", "Description", headers)
print(desc)
}
return(dat)
}
t=showTable(relFreqLongASAQConstructsScore_combinedData, constructASAQdesc$CABR,
constructASAQdesc$CNAME, 0, 1, 1, 1, 1, 1, 1)
showTable <- function(data_ori, abrConstructs, nameConstructs, type,
showDesc, showMean, showSD, showMedian, showMin, showMax){
dat <- data.frame(matrix(nrow = length(abrConstructs), ncol = 0))
data <- data.frame()
if (type == 0)
data <- convertToNumeric(data_ori, 2)
else if (type == 1){
headers <- data_ori[, 1]
data_ori <- convertToNumeric(data_ori, 2)
for(i in 1:length(headers)){
for (j in 1:length(abrConstructs)){
data[j, headers[i]] <- data_ori[i, j]
}
}
colnames(data) <- headers
}
if (type == 2){
data <- convertToNumeric(data_ori, 2)
n = ncol(data)/2
h1 <- paste("LONG", colnames(data)[1:n], sep = "_")
n = n + 1
h2 <- paste("SHORT", colnames(data)[n:ncol(data)], sep = "_")
colnames(data) <- c(h1,h2)
}
data <- round(data,2)
dat[, "CABR"] <- abrConstructs
dat[, "CNAME"] <- nameConstructs
dat <- cbind(dat, format(data, nsmall=2))
print(dat)
if (showDesc == 1){
# Do data description
i = 1
desc <- data.frame(matrix(nrow = 0, ncol = ncol(dat)))
if (showMean == 1){
desc [i, ] <- c("","Mean:", format(round(apply(data, 2, mean, na.rm=TRUE), 2)),
nsmall = 2)
i = i + 1
}
if (showSD == 1){
desc [i, ] <- c("","SD:", format(round(apply(data, 2, sd, na.rm=TRUE), 2)),
nsmall = 2)
i = i + 1
}
if (showMedian == 1){
desc [i, ] <- c("","Median:", format(round(apply(data, 2, median, na.rm=TRUE), 2)),
nsmall = 2)
i = i + 1
}
if (showMin == 1){
desc [i, ] <- c("","Min:", format(apply(data, 2, min, na.rm=TRUE), nsmall = 2))
i = i + 1
}
if (showMax == 1){
desc [i, ] <- c("","Max:", format(apply(data, 2, max, na.rm=TRUE), nsmall = 2))
i = i + 1
}
# for return purpose
#colnames(desc) <- names(dat)
#dat <- rbind(dat, desc)
# for printing purpose
colnames(desc) <- c("", "Description", names(data))
print(desc)
}
#return(dat)
}
t=showTable(relFreqLongASAQConstructsScore_combinedData, constructASAQdesc$CABR,
constructASAQdesc$CNAME, 0, 1, 1, 1, 1, 1, 1)
showTable <- function(data_ori, abrConstructs, nameConstructs, type,
showDesc, showMean, showSD, showMedian, showMin, showMax){
dat <- data.frame(matrix(nrow = length(abrConstructs), ncol = 0))
data <- data.frame()
if (type == 0)
data <- convertToNumeric(data_ori, 2)
else if (type == 1){
headers <- data_ori[, 1]
data_ori <- convertToNumeric(data_ori, 2)
for(i in 1:length(headers)){
for (j in 1:length(abrConstructs)){
data[j, headers[i]] <- data_ori[i, j]
}
}
colnames(data) <- headers
}
if (type == 2){
data <- convertToNumeric(data_ori, 2)
n = ncol(data)/2
h1 <- paste("LONG", colnames(data)[1:n], sep = "_")
n = n + 1
h2 <- paste("SHORT", colnames(data)[n:ncol(data)], sep = "_")
colnames(data) <- c(h1,h2)
}
data <- round(data,2)
dat[, "CABR"] <- abrConstructs
dat[, "CNAME"] <- nameConstructs
dat <- cbind(dat, format(data, nsmall=2))
print(dat)
if (showDesc == 1){
# Do data description
rows = showMean + showSD + showMedian + showMin + showMax
i = 1
desc <- data.frame(matrix(nrow = rows, ncol = ncol(dat)))
if (showMean == 1){
desc [i, ] <- c("","Mean:", format(round(apply(data, 2, mean, na.rm=TRUE), 2)),
nsmall = 2)
i = i + 1
}
if (showSD == 1){
desc [i, ] <- c("","SD:", format(round(apply(data, 2, sd, na.rm=TRUE), 2)),
nsmall = 2)
i = i + 1
}
if (showMedian == 1){
desc [i, ] <- c("","Median:", format(round(apply(data, 2, median, na.rm=TRUE), 2)),
nsmall = 2)
i = i + 1
}
if (showMin == 1){
desc [i, ] <- c("","Min:", format(apply(data, 2, min, na.rm=TRUE), nsmall = 2))
i = i + 1
}
if (showMax == 1){
desc [i, ] <- c("","Max:", format(apply(data, 2, max, na.rm=TRUE), nsmall = 2))
i = i + 1
}
# for return purpose
colnames(desc) <- names(dat)
dat <- rbind(dat, desc)
# for printing purpose
colnames(desc) <- c("", "Description", names(data))
print(desc)
}
return(dat)
}
t=showTable(relFreqLongASAQConstructsScore_combinedData, constructASAQdesc$CABR,
constructASAQdesc$CNAME, 0, 1, 1, 1, 1, 1, 1)
showTable <- function(data_ori, abrConstructs, nameConstructs, type,
showDesc, showMean, showSD, showMedian, showMin, showMax){
dat <- data.frame(matrix(nrow = length(abrConstructs), ncol = 0))
data <- data.frame()
if (type == 0)
data <- convertToNumeric(data_ori, 2)
else if (type == 1){
headers <- data_ori[, 1]
data_ori <- convertToNumeric(data_ori, 2)
for(i in 1:length(headers)){
for (j in 1:length(abrConstructs)){
data[j, headers[i]] <- data_ori[i, j]
}
}
colnames(data) <- headers
}
if (type == 2){
data <- convertToNumeric(data_ori, 2)
n = ncol(data)/2
h1 <- paste("LONG", colnames(data)[1:n], sep = "_")
n = n + 1
h2 <- paste("SHORT", colnames(data)[n:ncol(data)], sep = "_")
colnames(data) <- c(h1,h2)
}
data <- round(data,2)
dat[, "CABR"] <- abrConstructs
dat[, "CNAME"] <- nameConstructs
dat <- cbind(dat, format(data, nsmall=2))
print(dat)
return(dat)
}
t=showTable(relFreqLongASAQConstructsScore_combinedData, constructASAQdesc$CABR,
constructASAQdesc$CNAME, 0, 1, 1, 1, 1, 1, 1)
rows = 5
i = 1
desc <- data.frame(matrix(nrow = rows, ncol = ncol(dat)))
dat<-showTable(relFreqLongASAQConstructsScore_combinedData, constructASAQdesc$CABR,
constructASAQdesc$CNAME, 0, 1, 1, 1, 1, 1, 1)
rows = 5
i = 1
desc <- data.frame(matrix(nrow = rows, ncol = ncol(dat)))
if (showMean == 1){
desc [i, ] <- c("","Mean:", format(round(apply(data, 2, mean, na.rm=TRUE), 2)),
nsmall = 2)
i = i + 1
}
data_ori = relFreqLongASAQConstructsScore_combinedData
abrConstructs = constructASAQdesc$CABR
nameConstructs = constructASAQdesc$CNAME
type = 0
showDesc = 1
showMean = 1
showSD= 1
showMedian= 1
showMin= 1
showMax= 1
dat <- data.frame(matrix(nrow = length(abrConstructs), ncol = 0))
data <- data.frame()
if (type == 0)
data <- convertToNumeric(data_ori, 2)
else if (type == 1){
dat <- data.frame(matrix(nrow = length(abrConstructs), ncol = 0))
data <- data.frame()
if (type == 0)
data <- convertToNumeric(data_ori, 2)
else if (type == 1){
if (type == 0) data <- convertToNumeric(data_ori, 2)
else if (type == 1){
if (type == 0) {
data <- convertToNumeric(data_ori, 2)
}
else if (type == 1){
