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)