小编cjm*_*s85的帖子

我试图绘制一个U-Pb等时线,它基本上是一个xy图,在每个数据点周围有一个误差椭圆,表示该数据点的精度.这是我想要实现的一个例子,唯一的区别是我有更多的数据点:

我对x和y空间中的每个数据点都有对称的+和 - 错误,我想用它来指示错误椭圆的形状.

下面是我到目前为止编写的代码.这会绘制我的数据,但不是错误省略号,我目前有错误栏.

# Select data file.
fname<-file.choose()
# Rename imported data file.
upbiso <- read.table(file=fname, header= TRUE)
# Attaches data file as default data source.
attach(upbiso)
# Reads and display column headers in console.
names(upbiso)
# Sets margins around plot.
par(mar=c(5,5,4,2))
# Plots 238U/206Pb vs 207Pb/206Pb with superscript notation for isotopic masses, xlim and ylim sets min and max limit for axes, las rotates y axis labels, cex.lab increases the size of the axis labels.
plot(UPb, PbPb, …