我用盖革(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?