A ordered dataset created from pairData of the outbreak of 100 individuals
including the relative transmission probabilities for each pair estimated using the function
nbProbabilities. The code to recreate this dataset from pairData
is shown below.
nbResultsA data frame with 9900 rows and 24 variables:
A pair-level ID variable (the individual IDs separated by an '_').
The label for the run, here "SNPs".
The mean transmission probability for the pair over all runs.
The standard deviation of the transmission probability for the pair over all runs.
The mean relative transmission probability for the pair over all runs: pAvg scaled so that the probabilities for all infectors per infectee add to 1.
The rank of the probability of the the pair out of all pairs for that infectee (in case of ties all values have the minimum rank of the group).
The number of probability estimates that contributed to pAvg. This represents the number of prediction datasets this pair was included in over the 10x1 cross prediction repeated 50 times.
The ID of the potential "infector".
The ID of the potential "infectee".
Did individual.1 truly infect individual.2?
The number of SNPs between the individuals.
The date and time of infection of individualID.1.
The date and time of infection of individualID.2.
The date and time of sampling of individualID.1.
The date and time of sampling of individualID.2.
The number of days between sampleDate.1 and sampleDate.2.
The number of days between infectionDate.1 and infectionDate.2.
The number of years between infectionDate.1 and infectionDate.2.
A categorical representation of infectionDiff: <1y, 1-2y, 2-3y, 3-4y, 4-5y, >5y.
Pair-level covariate derived from X1: 1 if match, 0 if not match.
Pair-level covariate derived from X2: 1 if match, 0 if not match.
Pair-level covariate derived from X3: 1 if a-a, 2 if b-b, 3 if a-b, 4 if b-a.
Pair-level covariate derived from X4: 1 if match, 2 if adjacent, 2 otherwise.
Logical value indicating if a pair is a probable link. TRUE if the pair has fewer than 3 SNPs, FALSE if the pair has more than 12 SNPs, NA otherwise
# ## NOT RUN ##
# ## This is the code used to create this dataset ##
# orderedPair <- pairData[pairData$infectionDiff > 0, ]
# orderedPair$snpClose <- ifelse(orderedPair$snpDist < 3, TRUE,
# ifelse(orderedPair$snpDist > 12, FALSE, NA))
# set.seed(0)
# covariates = c("Z1", "Z2", "Z3", "Z4", "timeCat")
# resGen <- nbProbabilities(orderedPair = orderedPair,
# indIDVar = "individualID",
# pairIDVar = "pairID",
# goldStdVar = "snpClose",
# covariates = covariates,
# label = "SNPs", l = 1,
# n = 10, m = 1, nReps = 50)
# nbResults <- merge(resGen[[1]], orderedPair, by = "pairID", all = TRUE)