This function generates an amplitude plot of two fluorescence channels as found in droplet digital PCR.

plot_vic_fam(vic, fam, col_vic = "green", col_fam = "blue", circle = TRUE)

Arguments

vic
Amplitudes of the VIC channel - object of class dpcr.
fam
Amplitudes of the FAM channel - object of class dpcr.
col_vic
Color of the VIC channel.
col_fam
Color of the FAM channel.
circle
If TRUE circles are drawn, if FALSE not. If "numeric", specifies the radius of circles.

Details

Droplet digital PCR experiments consist of three steps (droplet generation, clonal amplification, droplet amplitude analysis). Typically 20000 nano-sized droplets are analyzed and separated into amplification-positive and amplification-negative droplets. An example of such system is the Bio-Rad QX100 and QX200 (Pinheiro et al. 2012). Such systems have applications in the detection of rare DNA target copies, the determination of copy number variations (CNV), detection of mutation, or expression analysis of genes or miRNA. Each droplet is analyzed individually using a virtual two-color detection system. The channels are treated separately but finally aligned (e.g., FAM and VIC or FAM and HEX).

References

Pinheiro, L.B., Coleman, V.A., Hindson, C.M., Herrmann, J., Hindson, B.J., Bhat, S., and Emslie, K.R. (2012). Evaluation of a droplet digital polymerase chain reaction format for DNA copy number quantification. Anal. Chem. 84, 1003 - 1011.

Examples

# Generate an amplitude plot for the first fluorescence channel (e.g., FAM) fluos1 <- sim_dpcr(m = 16, n = 30, times = 100, pos_sums = FALSE, n_exp = 1, fluo = list(0.1, 0)) # Generate an amplitude plot for the second fluorescence channel (e.g., VIC) fluos2 <- sim_dpcr(m = 16, n = 30, times = 100, pos_sums = FALSE, n_exp = 1, fluo = list(0.1, 0)) # Plot the amplitudes of both fluorescence channel in an aligned fashion plot_vic_fam(fam = fluos1, vic = fluos2)
# Same as above but different colors plot_vic_fam(fam = fluos1, vic = fluos2, col_vic = "red", col_fam = "yellow")
# Same as above without circles plot_vic_fam(fam = fluos1, vic = fluos2, col_vic = "red", col_fam = "yellow", circle = FALSE)
# Generate two channels in one object and plot them fluos_both <- sim_dpcr(m = 16, n = 30, times = 100, pos_sums = FALSE, n_exp = 2, fluo = list(0.1, 0))
#> The assumed volume of partitions in each run is equal to 1.
#> The assumed volume uncertainty in each run is equal to 0.
plot_vic_fam(extract_run(fluos_both, 1), extract_run(fluos_both, 2))