    install.packages('IsoplotR')

    # load the IsoplotR package:
    library(IsoplotR)

    # navigate to the system directory that stores the built-in data files:
    setwd(system.file(package='IsoplotR'))

    # Rb-Sr isochron:
    RbSr <- read.data('RbSr1.csv',method='Rb-Sr',format=1)
    isochron(RbSr)

    # Weighted mean plot:
    meandat <- read.data('LudwigMean.csv',method='other')
    weightedmean(meandat)

    # Cumulative Age Distribution:
    densdat <- read.data('LudwigKDE.csv',method='other')
    cad(densdat)

    # Kernel Density Estimate:
    mixture <- read.data('LudwigMixture.csv',method='other')
    kde(densdat,pch='|')

    # Radial plot with minimum age estimate:
    radialplot(mixture,k='min',bg='yellow')

    # 40Ar/39Ar release spectrum:
    ArAr <- read.data('ArAr3.csv',method='Ar-Ar',format=3)
    agespectrum(ArAr)

    # Concordia diagram of common-Pb corrected data with concordia age:
    UPb <- read.data('UPb6.csv',method='U-Pb',format=6)
    concordia(UPb,common.Pb=2,show.age=1,exterr=TRUE)

    # U-series evolution diagram:
    ThU <- read.data('ThU1.csv',method='Th-U',format=1)
    evolution(ThU,levels=ThU$x[,'U238Th232'],
              clabel=expression(paste(""^"238","U/"^"232","Th")))

    # U-series age vs. initial 234U/238U diagram:
    evolution(ThU,transform=TRUE,detrital=TRUE,
              ellipse.col=rgb(1,0,0,0.2),
              show.numbers=TRUE,isochron=TRUE)

    # Logratio plot of U-Th-Sm-He data:
    UThSmHe <- read.data('UThSmHe.csv',method='U-Th-He')
    helioplot(UThSmHe,model=3,
              levels=log10(UThSmHe[,'Sm']),
              clabel=expression("log[Sm]"))

    # Ternary diagram of U-Th-Sm-He data:
    helioplot(UThSmHe,model=1,logratio=FALSE,ellipse.col='lightblue')

    # Multidimensional Scaling (MDS) configuration of detrital ages:
    DZ <- read.data('DZ.csv',method='detritals')
    mds(DZ)