http://www.sr.bham.ac.uk/~ajrs/R/gallery/plot_cluster_z_vs_kT.txt
中的source("http://www.sr.bham.ac.uk/~ajrs/R/scripts/errorbar.R")
#--Add errorbars (using function in sourced errorbar.R script):
with(A, errorbar(z, 0, 0, Tx, -Tx.le, Tx.ue, bnd=F, col=col))        # bnd=F means input data are errors not confidence bounds

自己研究下

Output: