我用盖革（Geiger）计数器测量了X射线管的光谱。这些数据不是正态分布的，而是遵循Poisson distribution。我想对测量数据拟合一条曲线：
a = 1
b = 1
f(x) = a*(x/b - 1)/x**2
fit f(x) 'data.txt' via a,b
但是，这将执行最小二乘回归，在这种情况下，理论上是错误的。相反，我需要执行Poisson regression。我该如何实现？
泊松分布的link function为ln（μ）。有没有一种方法可以使用此信息，例如。用gnuplot的link command？