l n ( M c , p , d , r ) = C c , p , d + β 1 ∗ l n ( Y γ + β 2 ∗ l n ( H C γ ) + β 3 ∗ ( l n ( Y γ ) ) 2 + β 4 ∗ T + β 5 ∗ l n ( S I c , k , p , r ) {\displaystyle ln(M_{c,p,d,r})=C_{c,p,d}+\beta _{1}*ln(Y_{\gamma }+\beta _{2}*ln(HC_{\gamma })+\beta _{3}*(ln(Y_{\gamma }))^{2}+\beta _{4}*T+\beta _{5}*ln(SI_{c,k,p,r})}