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")

多くの点がプロットされたグラフがあり、特定の 1 つの点の誤差を過剰にプロットしています。友人は、軸に制限を設定することを提案しましたが、これは私の問題を解決しませんでした.

助けてください!!

4

0 に答える 0