powered by
Synthetic dataset of 8 Rb-Sr analysis that form a 1Ga isochron.
data(rbsr,package='geostats') plot(rbsr[,'RbSr'],rbsr[,'SrSr']) fit <- lm(SrSr ~ RbSr,data=rbsr) abline(fit)
Run the code above in your browser using DataLab