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)