l n ( M c , p , d , r ) = C c , p , d + β 1 ∗ l n ( Y γ + B e t a 2 ∗ l n ( H C γ ) + B e t a 3 ∗ ( l n ( Y γ ) ) 2 + B e t a 4 ∗ T + B 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+\mathrm {B} _{5}*ln(SI_{c,k,p,r})}