估计R中形式为z = k（x ^ a）（y ^ b）的模型

Question

给定数据x，y，z的数组我需要估计

中的常数k，a和b

z = k x^a y^b 

某些z数据包含零，这使得双方的日志变得困难。在对cross validated进行讨论后，我认为我需要拟合具有泊松链接函数的广义线性模型。但是如何在R中表达这个？我试过了

glm(z ~ x*y,data,family=poisson)

但是这会分别输出x，y和x * y的估计值，因此我怀疑它们与我想要估算的幂项无关。同时

glm(log z ~ log x + log y, data, family=gaussian)

适合NaN / inf的抱怨，大概是因为z中的零。

如何正确编码此模型？

Answer 1