I identify the following claims in the comments of Savage, Enquist and West(hereafter referred to as SEW) (Savage et al., 2007) about my manuscript(Chaui-Berlinck, 2006):

  1. that the first term in parentheses in equation 5a is not 0/0 and,therefore, the correct limit corroborates their results;

  2. that I have performed an analysis `in a regime where the equation does not hold';

  3. that a geometric constant is absorbed into the Lagrange multiplierλ k, and, therefore, the distinction between a cube or a sphere is irrelevant;

  4. that I have rehashed mistaken ideas from other authors and ignored the answers given to these authors.

In relation to (1), I agree with SEW that I have mistranscribed the first term in the equation, and I apologize for this. Therefore, that term is not 0/0, only the other two are, and the limit that SEW present is correct.

However, these were not the problems I addressed in my manuscript. In their original paper (West et al.,1997), the authors state that: `Although the resultβ k=n–1/3 is independent of k, it is not area-preserving and therefore does not give a=3/4 when used in Eq. 5; instead,it gives a=1.' Thus, the authors are discussing a general case ofβ=n–1/3, not the specific one related to two distinct rules for this ratio (i.e. equation 5a in Chaui-Berlinck, 2006). If one analyses the general case and, as I state in my manuscript, `take into account the possibility that the product nγβ2 could be equal to 1', then applying the sum of a power series to obtain Vb from equation 4 of West et al.(West et al., 1997), the net result is Vb=VcnN(N+1),as shown below:
\[\ r_{\mathrm{k}}=\frac{r_{\mathrm{k}+1}}{{\beta}}{\iff}r_{\mathrm{N-j}}=\frac{r_{\mathrm{N}}}{{\beta}^{\mathrm{j}}},\]
(1a)
\[\ l_{\mathrm{k}}=\frac{l_{\mathrm{k}+1}}{{\gamma}}{\iff}l_{\mathrm{N-j}}=\frac{l_{\mathrm{N}}}{{\gamma}^{\mathrm{j}}},\]
(1b)
where rN and lN are the radius and the length of a capillary, respectively, and j is a counter from the capillaries (j=0) towards the great vessels. Inserting r and l into equation 4 from West et al.(West et al., 1997) gives:
graphic
View largeDownload slide

As given by West et al. (West et al.,1997), the volume Vc of a capillary isπ rN2lN and, by the assumption of the case nγβ2=1, the result is the one presented above: Vb=VcnN(N+1). This prevents West et al. from obtaining a linear relationship of Vb with body mass in a general case. It has nothing to do with applying or not applying L'hospital's rule. Moreover, in contrast to what SEW state in their comments, this computation of Vbis not `the most critical claim made by Chaui-Berlinck'. This is,simply, another source of inconsistency in the West et al.(West et al., 1997) model.

In relation to (2), all I can say is that in my manuscript(Chaui-Berlinck, 2006) I present a step-by-step reasoning based on what West et al. stated in their paper (West et al., 1997). Every step in the analysis can be easily tracked by means of figure 1 in the manuscript (Chaui-Berlinck,2006). The impedance for each type of flow was analyzed by the equations put forward by West et al. themselves(West et al., 1997). For example, in my section `Impedances and resistances to flow'(Chaui-Berlinck, 2006), it is clearly demonstrated that West et al. cannot obtain the necessary ratio between radii, β=n–1/2, a conclusion that other authors, in a much deeper analysis, have drawn as well (e.g. Painter et al., 2006).

Item (3) is related to the energy minimization procedure taken by West et al. (West et al., 1997). On the one hand, SEW are correct in that a constant term (3π/2) could or could not be incorporated in a given Lagrange multiplier, just because it is a constant. On the other hand, West et al. directly assert that: `For a fixed mass M,the auxiliary Lagrange function F, which incorporates the constraints, must be minimized with respect to all variables for the entire system (rk, lk, and n).'(underlines have been added by me); then, they put their augmented function F:
\[\ F(r_{\mathrm{k}},l_{\mathrm{k}},n)=W(r_{\mathrm{k}},l_{\mathrm{k}},n_{\mathrm{k}},M)+{\lambda}V_{\mathrm{b}}(r_{\mathrm{k}},l_{\mathrm{k}},n_{\mathrm{k}},M)+{{\sum}_{k=0}^{N}}{\lambda}_{k}N_{\mathrm{k}}l_{\mathrm{k}}^{3}+{\lambda}_{\mathrm{M}}M.\]
(3)
Therefore, up to this point, one should understand that body mass Mis to be taken as a constraint in the model, which has, as variables, the radius r, the length l and the branching ratio n. So, the question is not whether we are dealing with spheres or cubes. The question concerns how West et al. constructed a way to obtain their desired results: `Now varying M and minimizing F in Eq. 7(∂F/∂M = 0) leads to Vb ∝ M, which is just the relation needed to derive Eq. 5.'(West et al., 1997). Therefore,they took the derivative of the augmented function F in relation to body mass, i.e. the derivative in relation to a constraint of their model, a completely heterodox/unexplained step in optimization via Lagrange multipliers.

In relation to (4), mea culpa. I agree with the criticism that I did not cite their response to other authors. The following explanation to the`service volume' issue can be found in Savage et al.(Savage et al., 2004): `As pointed out in the quotations below, WBE clearly state that only the characteristics of the capillaries themselves are assumed to be invariant. Nevertheless, K & K incorrectly interpreted this size-invariance to mean that each capillary must supply a constant volume of tissue' (underlines have been added by me). However, West et al. (West et al.,1997) state: `The network must branch so that a group of cells, referred to here as a “service volume,” is supplied by each capillary. Because rklk and the total number of branchings N is large, the volume supplied by the total network can be approximated by the sum of spheres whose diameters are that of a typical kth-level vessel, namely 4/3π(lk/2)3Nk.'(underlines have been added by me). Thus, what the authors didn't realize is that when they concluded that `the volume serviced by each capillary must scale as M1/4...' (West et al., 1997), they made themselves free from their own assumptions. Simply, they have two non-coincident explanations for `service volume'.

As I show here, the problem is not with my mathematical literacy. And it is not with my geometric skills. It is also not with those who reviewed my manuscript. The problem lies in a model that was, and still is, presented to the audience as a complete and general structure for dealing with almost all topics related to biological scaling. Such a model is fully discussed in my manuscript, which presents a number of fundamental issues that have been avoided more than answered by SEW.

In conclusion, the request for the retraction of `A critical understanding of the fractal model of metabolic scaling'(Chaui-Berlinck, 2006) is, at least, delusional.

Chaui-Berlinck, J. G. (
2006
). A critical understanding of the fractal model of metabolic scaling.
J. Exp. Biol.
209
,
3045
-3054.
Painter, P. R., Edén, P. and Bengtsson, H.-U.(
2006
). Pulsatile blood flow, shear force, energy dissipation and Murray's Law.
Theor. Biol. Med. Mod.
3
,
31
.
Savage, V. M., Gillooly, J. F., Woodruff, W. H., West, G. B.,Allen, A. P., Enquist, B. J. and Brown, J. H. (
2004
). The predominance of quarter-power scaling in biology.
Funct. Ecol.
18
,
257
-282.
Savage, V. M., Enquist, B. J. and West, G. B.(
2007
). Comment on `A critical understanding of the fractal model of metabolic scaling'.
J. Exp. Biol.
210
,
3873
-3874.
West, G. B., Brown, J. H. and Enquist, B. J.(
1997
). A general model for the origin of allometric scaling laws in biology.
Science
276
,
122
-126.
© The Company of Biologists Limited 2007
2007