I have a model where the volume ($V$) of a finger is normally distributed, with mean $\mu = \beta_0 L^{\beta_1}D^{\beta_2}$ (where $L=$ length, $D=$ diameter and $\beta_i \in \Bbb R$ for $i=0,1,2$) and some variance $\sigma^2$.
I was thinking that if this is a GLM (although I am not sure if we have an exponential family here), a log link function would be appropriate since $\ln (\mu)$ is a linear function of the natural log of the predictors: $$\ln (\mu) = \ln(\beta_0 L^{\beta_1}D^{\beta_2}) = \ln(\beta_0) + \beta_1\ln(L) + \beta_2\ln(D) $$
However, a Poisson distribution is not ideal since $V, L, D \in \Bbb R$ are not counts.
So, I am not sure if I should proceed with this, or if another link function would be more appropriate. Any suggestions would be appreciated.
glm(V ~ log(L) + log(D), family=gaussian(link="log"))
. But I would expectlm(log(V) ~ log(L) + log(D))
(assuming that $V$ conditional on $L$ and $D$ follows a lognormal distribution) to give a better fit to the data. This is also (for good reasons) how most allometric data such as this is analysed, see en.wikipedia.org/wiki/Allometry $\endgroup$glm(V ~ log(L) + log(D), family=gaussian(link="log"))
, but you think it would be better fit bylm(log(V) ~ log(L) + log(D))
? $\endgroup$