0

我的 y 轴上有一个对数刻度,值为负。我正在使用“细分”功能来绘制我的误差线。误差线出现在我的 x 轴上,但不是 y。

这是我的脚本:

massaxis<-read.table("mass_axis.dat", header=TRUE)

x<-log((massaxis$mass))
y<-log((massaxis$semimajoraxis))
plot(x, y, type="p",
     main="Log(semi-major axis) vs Log(mass) for all known exoplanets",
     ylab="log(Semi Major Axis) (AU)",
     xlab="log(mass) (mass in jupiter mass (Mj))")

fit<-lsfit(x, y, intercept=TRUE)
print(summary(fit))

a<-log(0.431)
b<-log(0.0556)
sigma<-0.095
delta<-0.0126
alo<-a-sigma
ahi<-a+sigma
blo<-b-delta
bhi<-b+delta

segments(alo, b, ahi, b, col="red")
segments(bhi, a, blo, a, col="red")
points(a, b, col="red")

我有一个绘制了很多点的图表,然后我过度绘制了一个特定点的错误。一位朋友建议我对我的轴设置限制,但这并没有解决我的问题。

请帮忙!!

4

0 回答 0