# Patent

## Title
Bioimpedance methods and apparatus

## Abstract

Methods and apparatus for providing bioimpedance analysis are provided. In certain aspects, equivalent circuit frequency response models are provided which lead to improved correlations with MRI data. The frequency response models take account of body composition, including the fat component of a body segment. Data obtained by performing bioimpedance spectroscopy (BIS) and MRI on the calves of subjects illustrates the improved correlations achieved compared to single frequency analyses at 50 kilohertz and analyses performed using the conventional Cole-Cole model.

## Background
CROSS REFERENCE TO RELATED APPLICATIONS
This application is a division of U.S. application Ser. No. 10/570,368 filed Mar. 2, 2006 now U.S. Pat. No. 7,945,317 which is the U.S. national phase under 35 USC §371 of International Application No. PCT/US2004/29711, filed Sep. 10, 2004, which was published in English under PCT Article 21(2) on Mar. 31, 2005 as International Publication No. WO 2005/027717. The contents of U.S. application Ser. No. 10/570,368 and International Application No. PCT/US2004/29711 are incorporated herein in their entireties
This application claims the benefit under 35 USC §119(e) of U.S. Provisional Application No. 60/502,483, filed Sep. 12, 2003, U.S. Provisional Application No. 60/580,166, filed Jun. 16, 2004, and U.S. Provisional Application No. 60/587,652, filed Jul. 14, 2004, the contents of all of which are incorporated herein in their entirety.


FIELD OF THE INVENTION
This invention relates to bioimpedance methods and apparatus and, more particularly, to the use of such methods and apparatus to determine the muscle content, fat content, and/or extracellular fluid content of a segment of the body of a subject, where a segment can include essentially the subject's entire body.
BACKGROUND OF THE INVENTION
The electrical properties of biologic tissue has been of scientific interest for a substantial period time. Many developments and new devices based on a knowledge of bioelectricity have been used in biology and the biomedical area in the last century.
Bioelectrical impedance analysis is one of the interesting and challenging subjects in this area. Bioimpedance has been studied in many areas of medicine because of its potential ability to measure body composition with noninvasive, simple, inexpensive, and portable methods. In particular, bioimpedance has been employed in clinical research for many years. For example, clinical applications using bioimpedance were reported at an early stage by Nyboer [3] and Patterson [4].
Many authors have investigated the nature of the electrical properties of living tissue [2, 4, 5]. Schwan et al. described the relationship between the dielectric properties of the cell membrane using multi-frequency currents [2]. A basic theory to explain electrical properties of tissue in the body has been well established by Cole [1]. In particular, Cole successfully developed an equivalent circuit model (hereinafter the “Cole-Cole model”) to explain the electrical response of cells and their membranes to AC current.
A method of using bioimpedance spectroscopy (BIS) has been suggested to measure extracellular (ECV) and intracellular (ICV) fluid volumes based on the Cole-Cole model and the Hanai method [5, 6, 10]. The methodology of multi-frequency bioimpedance analysis can now provide some information about extracellular and intracellular fluid volume in the total or segmental body compartment [6].
However, the accuracy of bioimpedance analysis, including BIS, is a major point of concern by the clinical user [7, 8]. Even though many studies have reported the use of bioimpedance analysis to estimate body fluids, the current techniques have not been accepted widely in clinical practice because of questions regarding reliability, validity, and accuracy.
There are many factors which adversely affect the accuracy of the measurement and analysis using currently available bioimpedance techniques. In accordance with certain of its aspects, the present invention is concerned with one of those factors, namely, the model used to analyze bioimpedance data. The bioimpedance model commonly used to date to calculate electrical properties of different tissues has a basic assumption that fat has a high resistivity compared to fat free mass, and that therefore, fat mass can be ignored. However, a recent study found that bioimpedance measurements at 50 kilohertz are affected when subjects have large amounts of adipose tissue [7]. As discussed in detail below, in accordance with the present invention it has been found that the amount of fat mass is one of the major factors affecting the accurate measurement of body composition by BIS for subjects having a variety of body mass index (BMI) values.
SUMMARY OF THE INVENTION
In accordance with a first aspect, the invention provides a method for analyzing bioimpedance data for a body segment of a subject, said body segment having an external skin surface, said method comprising:
(a) applying alternating current at a plurality of frequencies to at least two points on the external skin surface so as to cause current to pass through the segment;
(b) for each frequency, recording the voltage difference between at least two other points on the external skin surface, said recorded voltage differences comprising both magnitude and phase information (i.e., magnitude and phase values or, equivalently, resistance and reactance values); and
(c) using the recorded voltage differences at the plurality of frequencies to determine at least one numerical value indicative of the muscle, fat, and/or extracellular fluid content of the segment, said numerical value being determined using an impedance model for the segment which at least comprises three parallel paths, one of which consists of a capacitor CM and a resistor RI in series which represent primarily the muscle component of the segment, one of which consists of a capacitor CF and a resistor RF in series which represent primarily the fat component of the segment, and one of which consists of a resistor.
In accordance with a second aspect, the invention provides a method for analyzing bioimpedance data for a body segment of a subject, said body segment having an external skin surface, said method comprising:
(a) applying alternating current at a plurality of frequencies to at least two points on the external skin surface so as to cause current to pass through the segment;
(b) for each frequency, recording the voltage difference between at least two other points on the external skin surface, said recorded voltage differences comprising both magnitude and phase information (i.e., magnitude and phase values or, equivalently, resistance and reactance values); and
(c) using the recorded voltage differences at the plurality of frequencies to determine at least one numerical value indicative of the fat and/or extracellular fluid content of the segment, said numerical value being determined using an impedance model for the segment which at least comprises two parallel paths, one of which consists of a capacitor CF and a resistor RF in series which represent primarily the fat component of the segment and the other of which consists of a resistor which primarily represents the extracellular fluid component of the segment;
wherein:
(i) the two parallel paths are the only parallel paths of the impedance model which represent the composition of the segment internal to the skin; and
(ii) each of the frequencies applied in step (a) is less than or equal to 10 kilohertz.
In accordance with certain embodiments of the foregoing aspects of the invention, a correlation equation is used which transforms a model parameter (e.g., RF) to a physiological value (e.g., a fat content value) for the segment. The correlation equation can be obtained by:
(i) performing steps (a), (b), and (c) on a plurality of calibration subjects to obtain a model parameter value for each of said subjects;
(ii) performing a measurement on the segment for the plurality of calibration subjects (e.g., a measurement of fat content using magnetic resonance imaging) to obtain a physiological value for the segment for each of said subjects; and
(iii) performing a regression analysis on the values obtained in steps (i) and (ii) to obtain the correlation equation.
Preferably, the plurality of calibration subjects includes at least one subject having a body mass index less than 20 and at least one subject having a body mass index greater than 35. More preferably, the plurality of calibration subjects includes at least one subject having a body mass index less than 20 and at least one subject having a body mass index greater than 40.
In accordance with a third aspect, the invention provides a method for determining the circumference of a portion of a body segment covered by skin comprising:
(a) applying a series of electrodes around said portion, said series of electrodes having a first electrode and a last electrode, the circumferential distances between all electrodes in the series being known, except for the distance between the first and last electrodes;
(b) determining the resistance between at least two electrodes of the series, other than the first and last electrode, by applying a low frequency current which does not substantially penetrate the skin;
(c) determining a resistivity value per unit length for the skin from the resistance determined in step (b) and the known circumferential distance between the two electrodes;
(d) determining the resistance between the first and last electrodes of the series by applying a low frequency current which does not substantially penetrate the skin; and
(e) calculating the distance between the first and last electrodes of the series from the resistance measured in step (d) and the resistivity value per unit length determined in step (c).
It is to be understood that both the foregoing general description and the following detailed description are merely exemplary of the invention, and are intended to provide an overview or framework for understanding the nature and character of the invention.
Additional features and advantages of the invention are set forth in the detailed description which follows, and in part will be readily apparent to those skilled in the art from that description or recognized by practicing the invention as described herein. The accompanying drawings are included to provide a further understanding of the invention, and are incorporated in and constitute a part of this specification.



BRIEF DESCRIPTION OF THE DRAWINGS
FIG. 1 shows a body composition model in accordance with the invention, both in cross-sectional form (FIG. 1A) and three-dimensional form (FIG. 1B).
FIGS. 2A and 2B are equivalent circuit frequency response models for the body composition model of FIG. 1.
FIG. 3 is an equivalent circuit frequency response model for the body composition model of FIG. 1 under low frequency conditions, e.g., frequencies less than 10 kilohertz.
FIG. 4 shows an electrode configuration that can be used in the practice of the invention.
FIG. 5 shows a resistor model for the electrode configuration of FIG. 4.
FIG. 6 shows the electrode placement used to generate the bioimpedance data of the examples presented below.
FIG. 7 shows correlations between measured and calculated p values obtained using Equations 12 and 13. The steeper slope of the correlation obtained using Equation 12 illustrates the value of including information regarding fat content in equivalent circuit frequency response models for bioimpedance measurements.
FIG. 8 shows p values of a t-test between groups of subjects having different fat contents.
FIG. 9 shows curve fitting of impedance for three subjects (patients) having different fat contents.
FIG. 10 shows impedance magnitude (FIG. 10A) and impedance phase (FIG. 10B) for a subject (patient) with a 0.588 kg fat mass.
FIG. 11 shows impedance magnitude (FIG. 11A) and impedance phase (FIG. 11B) for a subject (patient) with a 0.134 kg fat mass.
FIG. 12 shows a correlation between muscle (MS) estimation by MRI and by a BIS technique of the invention.
FIG. 13 shows a correlation between muscle mass measured by MRI and resistance measured at 5 kHz.
FIG. 14 compares (i) a correlation between muscle mass measured by MRI and intracellular resistance (RI) determined using the multi-frequency, three parallel path model of the invention (FIG. 14A) with (ii) a correlation between muscle mass measured by MRI and impedance at 50 KHz (FIG. 14B).
FIG. 15 compares (i) a correlation between fat mass measured by MRI and fat resistance (RF) determined using the multi-frequency, three parallel path model of the invention (FIG. 15A) with (ii) a correlation between fat mass measured by MRI and impedance at 50 KHz (FIG. 15B).
FIG. 16 compares (i) a correlation between muscle mass measured by MRI and intracellular resistance (RI) determined using the multi-frequency, three parallel path model of the invention (FIG. 16A) with (ii) a correlation between muscle mass measured by MRI and intracellular resistance determined using the Cole-Cole model (FIG. 16B).
FIG. 17 compares (i) a correlation between fat mass measured by MRI and fat resistance (RF) determined using the multi-frequency, three parallel path model of the invention (FIG. 17A) with (ii) a correlation between fat mass measured by MRI and intracellular resistance determined using the Cole-Cole model (FIG. 17B).



DETAILED DESCRIPTION OF THE INVENTION AND ITS PREFERRED EMBODIMENTS
As discussed above, the present invention relates to bioimpedance methods and apparatus which can provide information regarding the fat, muscle, and/or extracellular fluid content of a segment of a body of a subject, i.e., a human or an animal.
The segment will typically be a part or all of a subject's limb, e.g., a subject's calf or forearm, but can also be all or part of the subject's torso, e.g., part or all of the subject's abdomen. Similarly, the segment can be as small as a specific muscle or a part of a muscle and as large as essentially the subject's entire body.
The location and size of the segment will depend on the placement of current applying electrodes on the subject's skin. In particular, the segment will constitute the portion of the subject's body through which substantial current passes when the current applying electrodes are activated. By suitable choices of the number, location, and polarity of the current applying electrodes, a variety of current patterns within the subject's body can be achieved. Indeed, by varying the polarity of selected current applying electrodes, more than one segment can be analyzed without the need to move the current applying electrodes.
As known in the art, the potential recording electrodes will typically be located inboard of the current applying electrodes, i.e., the potential recording electrodes will typically be located on skin which surrounds the portion of the patient's body through which substantial current passes when the current applying electrodes are activated.
The application of current and the recording of potentials can be performed with bioimpedance equipment, including current applying and recording electrodes, now known or subsequently developed, e.g., commercially available equipment such as the 4000B Bio-Impedance Spectrum Analyzer System (Xitron Technologies, Inc., San Diego, Calif.) used in the examples discussed below. Alternatively, customized equipment can be used in the practice of the invention.
Processing of the data obtained by the bioimpedance equipment can be performed entirely within the equipment or can be performed on-line using a separate computer. Alternatively, the data can be stored and processed subsequent to the time of measurement.
Preferably, the bioimpedance equipment will include a microprocessor programmed to perform at least a portion of the analysis procedures of the present invention. For example, the bioimpedance equipment can apply a regression equation obtained from a study on calibration subjects to measured impedance data for a subject and thus directly report the subject's fat, muscle, and/or extracellular fluid content(s) to the user of the equipment, e.g., the subject himself and/or other personnel, e.g., a health care provider, interested in the information. The fat, muscle, and/or extracellular fluid content can be reported graphically, numerically (e.g., as a fat and/or muscle percentage value), by color (e.g., red for high fat content), or the like.
In certain embodiments of the invention, alternating current at a plurality of frequencies is applied to the subject's skin. Preferably, at least 10 frequencies are used and more preferably, approximately 50 frequencies are used. Preferably, the plurality of frequencies comprises frequencies between 5 and 1000 kilohertz, although frequencies over larger or smaller ranges can be used if desired. Most preferably, the frequencies are logarithmically distributed over the frequency range.
As discussed above, the bioimpedance model commonly used to date to calculate electrical properties of different tissues has suffered from a number of problems. One of those problems is the basic assumption that fat has a high resistivity compared to fat free mass, and that therefore, fat mass can be ignored. Another is related to the model's calculation of extracellular and intracellular resistances. Although not accounted for in the existing model, that model, which is based on the tissue under the skin, is influenced by the skin and amount of adipose tissue.
The present invention is directed to reducing various of the problems with the present analysis approach by examining the response of different constituent body components, such as skin, fat, muscle, to electrical input. More particularly, the invention, in accordance with certain of its aspects, provides models to describe the effects of different components of body tissue using equivalent circuits.
Specifically, to improve current bioimpedance techniques, in accordance with certain aspects of the present invention, an improved electrical model is provided which is able to explain the electrical properties with different proportions of body tissues, and the effect of a wide range of current frequencies. The data presented below in Example 1 specifically evaluates the relationship of the resistivity at 5 kHz measured at the skin to the resistivity by calculation with the model.
FIG. 1 shows a segmental body composition model in accordance with the invention which can be used to describe the components of conductivity in the limbs. In particular, the figure shows conductive components in a limb segment which form the basis for the electrical models of the invention. For reference, Table 1 shows the electrical properties (permittivity and resistivity) of different tissues using a 10 kHz current frequency as reported in a previous study [8].
An equivalent circuit model to FIG. 1 is shown in FIG. 2A, where RG represents resistance of measurement by two electrodes on the skin and RS, RF, RE, RI, and RB represent resistance of the skin, fat mass, extracellular volume, intracellular volume, and bone, respectively. In FIG. 2A, CIN, CF, and CM represent capacitance between skin and electrode, capacitance of fat mass, and capacitance of membrane of cells, respectively.
From the model in FIG. 2A, the total electrical current (IG) can be given by:

IG=IS+IF+IE+II+IB  Eq. 1

Potential across RG can be obtained by:




  
    
      
        
          
            R
            G
          
          ⁢
          
            I
            G
          
        
        =
        
          
            
              R
              E
            
            ⁢
            
              I
              E
            
          
          +
          
            
              2
              ⁢
              
                I
                G
                ′
              
            
            
              ω
              ⁢
              
                
              
              ⁢
              
                C
                IN
              
            
          
        
      
    
    
      
        Eq
        .
        
          
        
        ⁢
        2
      
    
  




where I′G is the current pass through capacitor CIN, and is given by:

I′G=IF+IE+II+IB  Eq. 3

Because skin resistance (RS) is much higher than in other tissue IS is very low, so that we have:

I′G≈IG  Eq. 4

Since current at low frequency (e.g., 5 kHz) will not pass through the intracellular space, the equivalent circuit can be modified as shown in FIG. 3.
In FIG. 3, E represents the potential between two measuring electrodes (across RG) and EE represents the potential across RE. Therefore, the following potential equation can be written:




  
    
      
        E
        =
        
          
            E
            E
          
          +
          
            
              2
              ⁢
              
                I
                G
              
            
            
              jω
              ⁢
              
                
              
              ⁢
              
                C
                IN
              
            
          
        
      
    
    
      
        Eq
        .
        
          
        
        ⁢
        5
      
    
  




where the parameter





  
    2
    ⁢
    
      I
      G
    
  
  
    jω
    ⁢
    
      
    
    ⁢
    
      C
      IN
    
  




represents the potential across the capacitor CIN. This potential can be ignored at high current frequency or when CIN is large.

According to the parallel circuit of FIG. 3, the resistance RG of measurement from skin electrodes can be calculated by:




  
    
      
        
          R
          G
        
        =
        
          
            
              
                R
                E
              
              ⁡
              
                (
                
                  
                    1
                    
                      jω
                      ⁢
                      
                        
                      
                      ⁢
                      
                        C
                        F
                      
                    
                  
                  +
                  
                    R
                    F
                  
                
                )
              
            
            
              
                R
                E
              
              +
              
                1
                
                  jω
                  ⁢
                  
                    
                  
                  ⁢
                  
                    C
                    F
                  
                
              
              +
              
                R
                F
              
            
          
          +
          
            2
            
              jω
              ⁢
              
                
              
              ⁢
              
                C
                IN
              
            
          
        
      
    
    
      
        Eq
        .
        
          
        
        ⁢
        6
      
    
  




Reactance of fat mass (XF) and reactance between skin and electrode (XIN) are as follows:




  
    
      X
      F
    
    =
    
      -
      
        j
        
          ω
          ⁢
          
            
          
          ⁢
          
            C
            F
          
        
      
    
  
  ,
  
    
      X
      IN
    
    =
    
      -
      
        j
        
          ω
          ⁢
          
            
          
          ⁢
          
            C
            IN
          
        
      
    
  




Thus, RG can be simplified to:




  
    
      
        
          
            R
            G
          
          =
          
            
              
                R
                E
              
              (
              
                1
                +
                
                  
                    
                      R
                      E
                    
                    /
                    
                      X
                      F
                    
                  
                  
                    1
                    +
                    
                      
                        (
                        
                          
                            R
                            E
                          
                          +
                          
                            R
                            F
                          
                        
                        )
                      
                      /
                      
                        X
                        F
                      
                    
                  
                
              
              )
            
            -
            
              2
              ⁢
              
                X
                IN
              
            
          
        
        ⁢
        
          

        
        ⁢
        and
      
    
    
      
        Eq
        .
        
          
        
        ⁢
        7
      
    
  
  
    
      
        
          R
          G
        
        =
        
          
            
              R
              E
            
            (
            
              1
              +
              
                
                  A
                  F
                
                
                  
                    
                      
                        ρ
                        F
                      
                      
                        ρ
                        E
                      
                    
                    ⁢
                    
                      A
                      E
                    
                  
                  +
                  
                    A
                    F
                  
                
              
            
            )
          
          -
          
            2
            ⁢
            
              X
              IN
            
          
        
      
    
    
      
        Eq
        .
        
          
        
        ⁢
        8
      
    
  




where AF and AE represent the cross sectional area of fat mass from MRI and ECV in this segment, respectively. According to a previous study, the ratio of resistivity in fat mass to resistivity in ECV is approximately as follows [6].





  
    
      
        
          
            ρ
            F
          
          
            ρ
            K
          
        
        ≈
        
          3
          -
          5
        
      
    
    
      
        Eq
        .
        
          
        
        ⁢
        9
      
    
  




Equation 8 can be further simplified to:




  
    
      
        
          ρ
          G
        
        =
        
          
            
              ρ
              E
            
            ⁢
            
              
                A
                G
              
              
                A
                E
              
            
            ⁢
            
              (
              
                1
                +
                
                  
                    A
                    F
                  
                  
                    A
                    G
                  
                
              
              )
            
          
          -
          
            2
            ⁢
            
              X
              IN
            
            ⁢
            
              
                A
                G
              
              L
            
          
        
      
    
    
      
        Eq
        .
        
          
        
        ⁢
        10
      
    
  




where ρG is the resistivity measured on the skin and ρE is the resistivity in ECV space with a constant value. Letting





  
    
      
        
          k
          E
        
        =
        
          (
          
            1
            +
            
              
                A
                F
              
              
                A
                G
              
            
          
          )
        
      
    
    
      
        Eq
        .
        
          
        
        ⁢
        11
      
    
  




and assuming that





  2
  ⁢
  
    X
    IN
  
  ⁢
  
    
      A
      G
    
    L
  




is small, we have:





  
    
      
        
          ρ
          
            G
            ,
            Cal
          
        
        =
        
          
            k
            E
          
          ⁢
          
            ρ
            E
          
          ⁢
          
            
              A
              G
            
            
              A
              E
            
          
        
      
    
    
      
        Eq
        .
        
          
        
        ⁢
        12
      
    
  




If kE=1, resistivity can be calculated by:




  
    
      
        
          ρ
          
            G
            ,
            Cal
          
          *
        
        =
        
          
            ρ
            E
          
          ⁢
          
            
              A
              G
            
            
              A
              E
            
          
        
      
    
    
      
        Eq
        .
        
          
        
        ⁢
        13
      
    
  




The application of these equations and, in particular, Equations 12 and 13 to experimental data obtained from a set of calibration subjects is set forth below in Example 1. The results described therein indicate that the volume of fat mass is an important factor influencing the estimation of body composition using the standard four-electrode bioimpedance technique.
The experimental data of Example 1 was obtained using one frequency for the applied current, namely, 5 kilohertz. The invention's equivalent circuit models of segmental body components are preferably used to describe responses to multi-frequency current. Example 2 shows the results of such multi-frequency testing. In particular, this example shows different decreases in impedance with increases in current frequency for subjects with different fat contents.
FIG. 2B is a modified version of the equivalent circuit model of FIG. 2A. In FIG. 2B, ZTotal represents total impedance of measurement by two electrodes at A and B. As in FIG. 2A, RS, RF, RE, RI, and RB represent resistance of the skin, fat mass, extracellular and intracellular volume, and bone, respectively, and CIN, CF, and CM represent capacitance between skin and electrode, capacitance of fat mass and the capacitance of the membranes of the cells, respectively. Ri1 and Ri2 represent resistance of segmental body between injecting and measuring electrodes.
To obtain a relationship based on frequency for the model of FIG. 2, equations were used to simplify the model so that it can be standardized. First, we let RP represent the parallel resistance of RE and RB:




  
    
      
        
          R
          P
        
        =
        
          
            
              R
              E
            
            ×
            
              R
              B
            
          
          
            
              R
              E
            
            +
            
              R
              B
            
          
        
      
    
    
      
        Eq
        .
        
          
        
        ⁢
        14
      
    
  




Secondly, we used the following parameters a, b, c, and d to simplify the calculation:

a=RIRFCICF  Eq. 15

b=CICF(RI+RF)  Eq. 16

c=CFCI  Eq. 17

d+CF+CI  Eq. 18

Total impedance (Ztotal) can be obtained from the circuit of FIG. 2 and is given by the following Eq. 19:




  
    
      
        
          Z
          Total
        
        =
        
          
            
              
                
                  
                    
                      aC
                      IN
                    
                    ⁢
                    
                      R
                      P
                    
                    ⁢
                    
                      
                        
                          R
                          S
                        
                        ⁡
                        
                          (
                          jπ
                          )
                        
                      
                      3
                    
                  
                  +
                
              
            
            
              
                
                  
                    
                      
                        R
                        S
                      
                      ⁡
                      
                        [
                        
                          
                            
                              (
                              
                                2
                                +
                                
                                  C
                                  IN
                                
                              
                              )
                            
                            ⁢
                            
                              bR
                              P
                            
                          
                          +
                          
                            2
                            ⁢
                            a
                          
                        
                        ]
                      
                    
                    ⁢
                    
                      
                        (
                        jπ
                        )
                      
                      2
                    
                  
                  +
                
              
            
            
              
                
                  
                    R
                    S
                  
                  [
                  
                    
                      
                        (
                        
                          
                            
                              C
                              IN
                            
                            ⁢
                            
                              cR
                              P
                            
                          
                          +
                          
                            2
                            ⁢
                            b
                          
                          +
                          
                            2
                            ⁢
                            dRp
                          
                        
                        ]
                      
                      ⁢
                      
                        (
                        jπ
                        )
                      
                    
                    +
                    
                      2
                      ⁢
                      
                        cR
                        S
                      
                    
                  
                
              
            
          
          
            
              
                
                  
                    
                      
                        C
                        IN
                      
                      ⁡
                      
                        [
                        
                          
                            
                              R
                              S
                            
                            ⁡
                            
                              (
                              
                                a
                                +
                                
                                  bR
                                  P
                                
                              
                              )
                            
                          
                          +
                          
                            aR
                            P
                          
                        
                        ]
                      
                    
                    ⁢
                    
                      
                        (
                        jπ
                        )
                      
                      3
                    
                  
                  +
                
              
            
            
              
                
                  
                    
                      [
                      
                        
                          
                            R
                            S
                          
                          ⁢
                          
                            
                              C
                              IN
                            
                            ⁡
                            
                              (
                              
                                
                                  dR
                                  P
                                
                                +
                                b
                              
                              )
                            
                          
                        
                        +
                        
                          
                            bR
                            P
                          
                          ⁡
                          
                            (
                            
                              2
                              +
                              
                                C
                                IN
                              
                            
                            )
                          
                        
                        +
                        
                          2
                          ⁢
                          a
                        
                      
                      ]
                    
                    ⁢
                    
                      
                        (
                        jπ
                        )
                      
                      2
                    
                  
                  +
                
              
            
            
              
                
                  
                    
                      [
                      
                        
                          
                            cR
                            S
                          
                          ⁢
                          
                            C
                            IN
                          
                        
                        +
                        
                          
                            (
                            
                              
                                cC
                                IN
                              
                              +
                              
                                2
                                ⁢
                                d
                              
                            
                            )
                          
                          ⁢
                          
                            R
                            P
                          
                        
                        +
                        
                          2
                          ⁢
                          b
                        
                      
                      )
                    
                    ⁢
                    
                      (
                      jπ
                      )
                    
                  
                  +
                  
                    2
                    ⁢
                    c
                  
                
              
            
          
        
      
    
    
      
        Eq
        .
        
          
        
        ⁢
        19
      
    
  




To standardize Eq. 19, we let

μ=CIN[RS(a+bRP)+aRP]  Eq. 20

ZTotal is then given by Eq. 21 as follows:




  
    
      
        
          Z
          Total
        
        =
        
          
            
              
                
                  
                    
                      (
                      
                        1
                        μ
                      
                      )
                    
                    ⁢
                    
                      aC
                      IN
                    
                    ⁢
                    
                      R
                      P
                    
                    ⁢
                    
                      
                        
                          R
                          S
                        
                        ⁡
                        
                          (
                          jπ
                          )
                        
                      
                      3
                    
                  
                  +
                
              
            
            
              
                
                  
                    
                      1
                      μ
                    
                    ⁢
                    
                      
                        R
                        S
                      
                      ⁡
                      
                        [
                        
                          
                            
                              (
                              
                                2
                                +
                                
                                  C
                                  IN
                                
                              
                              )
                            
                            ⁢
                            
                              bR
                              P
                            
                          
                          +
                          
                            2
                            ⁢
                            a
                          
                        
                        ]
                      
                    
                    ⁢
                    
                      
                        (
                        jπ
                        )
                      
                      2
                    
                  
                  +
                
              
            
            
              
                
                  
                    
                      R
                      S
                    
                    μ
                  
                  [
                  
                    
                      
                        (
                        
                          
                            
                              C
                              IN
                            
                            ⁢
                            
                              cR
                              P
                            
                          
                          +
                          
                            2
                            ⁢
                            b
                          
                          +
                          
                            2
                            ⁢
                            dRp
                          
                        
                        ]
                      
                      ⁢
                      
                        (
                        jπ
                        )
                      
                    
                    +
                    
                      
                        2
                        ⁢
                        
                          cR
                          S
                        
                      
                      μ
                    
                  
                
              
            
          
          
            
              
                
                  
                    
                      (
                      jπ
                      )
                    
                    3
                  
                  +
                  
                    
                      
                        1
                        μ
                      
                      ⁡
                      
                        [
                        
                          
                            
                              
                                
                                  
                                    R
                                    S
                                  
                                  ⁢
                                  
                                    
                                      C
                                      IN
                                    
                                    ⁡
                                    
                                      (
                                      
                                        
                                          dR
                                          P
                                        
                                        +
                                        b
                                      
                                      )
                                    
                                  
                                
                                +
                              
                            
                          
                          
                            
                              
                                
                                  
                                    bR
                                    P
                                  
                                  ⁡
                                  
                                    (
                                    
                                      2
                                      +
                                      
                                        C
                                        IN
                                      
                                    
                                    )
                                  
                                
                                +
                                
                                  2
                                  ⁢
                                  a
                                
                              
                            
                          
                        
                        ]
                      
                    
                    ⁢
                    
                      
                        (
                        jπ
                        )
                      
                      2
                    
                  
                  +
                
              
            
            
              
                
                  
                    
                      
                        
                          1
                          μ
                        
                        [
                        
                          
                            
                              cR
                              S
                            
                            ⁢
                            
                              C
                              IN
                            
                          
                          +
                          
                            
                              (
                              
                                
                                  cC
                                  IN
                                
                                +
                                
                                  2
                                  ⁢
                                  d
                                
                              
                              )
                            
                            ⁢
                            
                              R
                              P
                            
                          
                          +
                          
                            2
                            ⁢
                            b
                          
                        
                        )
                      
                      ]
                    
                    ⁢
                    
                      (
                      jπ
                      )
                    
                  
                  +
                  
                    
                      2
                      ⁢
                      c
                    
                    μ
                  
                
              
            
          
        
      
    
    
      
        Eq
        .
        
          
        
        ⁢
        21
      
    
  




Quantitative analysis of frequency response data for individual subjects is preferably performed using a group of parameters to represent the function of components in the equivalent circuit. Such a group of parameters helps in the identification of differences in the body composition of individual subjects by facilitating the use of a standard transfer function for the analysis.
Equation 21 can, for example, be normalized using the following relationship, where P1, P2, P3, P4, Q1, Q2, Q3 are the group of parameters:




  
    
      
        
          Z
          Total
        
        =
        
          
            
              
                
                  P
                  1
                
                ⁡
                
                  (
                  jπ
                  )
                
              
              3
            
            +
            
              
                
                  P
                  2
                
                ⁡
                
                  (
                  jπ
                  )
                
              
              2
            
            +
            
              
                P
                3
              
              ⁡
              
                (
                jπ
                )
              
            
            +
            
              P
              4
            
          
          
            
              
                (
                jπ
                )
              
              3
            
            +
            
              
                
                  Q
                  1
                
                ⁡
                
                  (
                  jπ
                  )
                
              
              2
            
            +
            
              
                Q
                2
              
              ⁡
              
                (
                jπ
                )
              
            
            +
            
              Q
              3
            
          
        
      
    
    
      
        Eq
        .
        
          
        
        ⁢
        22
      
    
  




Values for RI, RF, CI, CF, RP, RS, CIN for individual patients (individual subjects) are calculated using the following equations:




  
    
      
        
          
            
              
                P
                1
              
              =
              
              ⁢
              
                
                  (
                  
                    1
                    μ
                  
                  )
                
                ⁢
                
                  aC
                  IN
                
                ⁢
                
                  R
                  P
                
                ⁢
                
                  R
                  S
                
              
            
          
        
        
          
            
              =
              
              ⁢
              
                
                  
                    aC
                    IN
                  
                  ⁢
                  
                    R
                    P
                  
                  ⁢
                  
                    R
                    S
                  
                
                
                  
                    C
                    IN
                  
                  ⁡
                  
                    [
                    
                      
                        
                          R
                          S
                        
                        ⁡
                        
                          (
                          
                            a
                            +
                            
                              bR
                              P
                            
                          
                          )
                        
                      
                      +
                      
                        aR
                        P
                      
                    
                    ]
                  
                
              
            
          
        
        
          
            
              =
              
              ⁢
              
                
                  
                    aR
                    P
                  
                  ⁢
                  
                    R
                    S
                  
                
                
                  
                    
                      R
                      S
                    
                    ⁡
                    
                      (
                      
                        a
                        +
                        
                          bR
                          P
                        
                      
                      )
                    
                  
                  +
                  
                    aR
                    P
                  
                
              
            
          
        
      
    
    
      
        Eq
        .
        
          
        
        ⁢
        23
      
    
  
  
    
      
        
          
            
              
                P
                2
              
              =
              
              ⁢
              
                
                  1
                  μ
                
                ⁢
                
                  
                    R
                    S
                  
                  ⁡
                  
                    [
                    
                      
                        
                          (
                          
                            2
                            +
                            
                              C
                              IN
                            
                          
                          )
                        
                        ⁢
                        
                          bR
                          P
                        
                      
                      +
                      
                        2
                        ⁢
                        a
                      
                    
                    ]
                  
                
              
            
          
        
        
          
            
              =
              
              ⁢
              
                
                  
                    R
                    S
                  
                  ⁡
                  
                    [
                    
                      
                        
                          (
                          
                            2
                            +
                            
                              C
                              IN
                            
                          
                          )
                        
                        ⁢
                        
                          bR
                          P
                        
                      
                      +
                      
                        2
                        ⁢
                        a
                      
                    
                    ]
                  
                
                
                  
                    C
                    IN
                  
                  ⁡
                  
                    [
                    
                      
                        
                          R
                          S
                        
                        ⁡
                        
                          (
                          
                            a
                            +
                            
                              bR
                              P
                            
                          
                          )
                        
                      
                      +
                      
                        aR
                        P
                      
                    
                    ]
                  
                
              
            
          
        
      
    
    
      
        Eq
        .
        
          
        
        ⁢
        24
      
    
  
  
    
      
        
          
            
              
                P
                3
              
              =
              
              ⁢
              
                
                  
                    R
                    S
                  
                  μ
                
                ⁡
                
                  [
                  
                    
                      
                        C
                        IN
                      
                      ⁢
                      
                        cR
                        P
                      
                    
                    +
                    
                      2
                      ⁢
                      b
                    
                    +
                    
                      2
                      ⁢
                      dRp
                    
                  
                  ]
                
              
            
          
        
        
          
            
              =
              
              ⁢
              
                
                  
                    R
                    S
                  
                  ⁡
                  
                    [
                    
                      
                        
                          C
                          IN
                        
                        ⁢
                        
                          cR
                          P
                        
                      
                      +
                      
                        2
                        ⁢
                        b
                      
                      +
                      
                        2
                        ⁢
                        dRp
                      
                    
                    ]
                  
                
                
                  
                    C
                    IN
                  
                  ⁡
                  
                    [
                    
                      
                        
                          R
                          S
                        
                        ⁡
                        
                          (
                          
                            a
                            +
                            
                              bR
                              P
                            
                          
                          )
                        
                      
                      +
                      
                        aR
                        P
                      
                    
                    ]
                  
                
              
            
          
        
      
    
    
      
        Eq
        .
        
          
        
        ⁢
        25
      
    
  
  
    
      
        
          P
          4
        
        =
        
          
            
              2
              ⁢
              
                cR
                S
              
            
            μ
          
          =
          
            
              2
              ⁢
              
                cR
                S
              
            
            
              
                C
                IN
              
              ⁡
              
                [
                
                  
                    
                      R
                      S
                    
                    ⁡
                    
                      (
                      
                        a
                        +
                        
                          bR
                          P
                        
                      
                      )
                    
                  
                  +
                  
                    aR
                    P
                  
                
                ]
              
            
          
        
      
    
    
      
        Eq
        .
        
          
        
        ⁢
        26
      
    
  
  
    
      
        
          
            
              
                Q
                1
              
              =
              
              ⁢
              
                
                  1
                  μ
                
                ⁡
                
                  [
                  
                    
                      
                        R
                        S
                      
                      ⁢
                      
                        
                          C
                          IN
                        
                        ⁡
                        
                          (
                          
                            
                              dR
                              P
                            
                            +
                            b
                          
                          )
                        
                      
                    
                    +
                    
                      
                        bR
                        P
                      
                      ⁡
                      
                        (
                        
                          2
                          +
                          
                            C
                            IN
                          
                        
                        )
                      
                    
                    +
                    
                      2
                      ⁢
                      a
                    
                  
                  ]
                
              
            
          
        
        
          
            
              =
              
              ⁢
              
                
                  [
                  
                    
                      
                        R
                        S
                      
                      ⁢
                      
                        
                          C
                          IN
                        
                        ⁡
                        
                          (
                          
                            
                              dR
                              P
                            
                            +
                            b
                          
                          )
                        
                      
                    
                    +
                    
                      
                        bR
                        P
                      
                      ⁡
                      
                        (
                        
                          2
                          +
                          
                            C
                            IN
                          
                        
                        )
                      
                    
                    +
                    
                      2
                      ⁢
                      a
                    
                  
                  ]
                
                
                  
                    C
                    IN
                  
                  ⁡
                  
                    [
                    
                      
                        
                          R
                          S
                        
                        ⁡
                        
                          (
                          
                            a
                            +
                            
                              bR
                              P
                            
                          
                          )
                        
                      
                      +
                      
                        aR
                        P
                      
                    
                    ]
                  
                
              
            
          
        
      
    
    
      
        Eq
        .
        
          
        
        ⁢
        27
      
    
  
  
    
      
        
          
            
              
                
                  Q
                  2
                
                =
                
                ⁢
                
                  
                    1
                    μ
                  
                  [
                  
                    
                      
                        cR
                        S
                      
                      ⁢
                      
                        C
                        IN
                      
                    
                    +
                    
                      
                        (
                        
                          
                            cC
                            IN
                          
                          +
                          
                            2
                            ⁢
                            d
                          
                        
                        )
                      
                      ⁢
                      
                        R
                        P
                      
                    
                    +
                    
                      2
                      ⁢
                      b
                    
                  
                  )
                
              
              ]
            
          
        
        
          
            
              =
              
              ⁢
              
                
                  
                    
                      [
                      
                        
                          
                            cR
                            S
                          
                          ⁢
                          
                            C
                            IN
                          
                        
                        +
                        
                          
                            (
                            
                              
                                cC
                                IN
                              
                              +
                              
                                2
                                ⁢
                                d
                              
                            
                            )
                          
                          ⁢
                          
                            R
                            P
                          
                        
                        +
                        
                          2
                          ⁢
                          b
                        
                      
                    
                    )
                  
                  ]
                
                
                  
                    C
                    IN
                  
                  ⁡
                  
                    [
                    
                      
                        
                          R
                          S
                        
                        ⁡
                        
                          (
                          
                            a
                            +
                            
                              bR
                              P
                            
                          
                          )
                        
                      
                      +
                      
                        aR
                        P
                      
                    
                    ]
                  
                
              
            
          
        
      
    
    
      
        Eq
        .
        
          
        
        ⁢
        28
      
    
  
  
    
      
        
          Q
          3
        
        =
        
          
            
              2
              ⁢
              c
            
            μ
          
          =
          
            
              2
              ⁢
              c
            
            
              
                
                  C
                  IN
                
                ⁡
                
                  [
                  
                    
                      
                        R
                        S
                      
                      ⁡
                      
                        (
                        
                          a
                          +
                          
                            bR
                            P
                          
                        
                        )
                      
                    
                    +
                    
                      aR
                      P
                    
                  
                  ]
                
              
              ]
            
          
        
      
    
    
      
        Eq
        .
        
          
        
        ⁢
        29
      
    
  




where a, b, c, and d are as defined above in Equations 15-18.

In these equations, RI, RF, CI, CF, RP, RS, CIN are the variables to be determined and P1, P2, P3, P4 Q1, Q2, Q3 are data values for an individual patient obtained by curve fitting of the frequency response data for the patient. In particular, the values of the seven variables RI, RF, CI, CF, RP, RS, CIN for an individual patient are obtained by solving the above seven equations (Eq. 23 to Eq. 29) using the P1, P2, P3, P4 Q1, Q2, Q3 values for that patient.
In Example 2 below, the curve fitting program of the MATLAB toolbox (The MathWorks, Inc., Natick, Mass.) was used to obtain the values of the parameters P1, P2, P3, P4, Q1, Q2, Q3 for individual patients. Simulation of frequency response curves was obtained using the standard signal processing program of the MATLAB toolbox.
Although Eq. 22 is a preferred curve fitting equation for use with the type of frequency response data obtained during a bioimpedance experiment, as will be evident to persons skilled in the art from the present disclosure, other curve fitting equations can be used if desired. Similarly, other math processing programs besides the MATLAB programs can be used in the practice of the invention.
Segmental length and circumference can be measured manually. Alternatively, in accordance with certain embodiments of the invention, the circumference and the length of a segment can be measured using groups of electrodes instead of single electrodes (see FIG. 4). The resistance at low frequency (e.g., <10 Hz, preferably <1 Hz) can reflect the electrical properties of the skin. Since the distance between any adjacent two electrodes in each pair is known, the circumference of the measurement area can be calculated as,




  
    
      
        C
        =
        
          
            D
            ·
            
              (
              
                n
                -
                1
              
              )
            
          
          +
          
            
              R
              n
            
            λ
          
        
      
    
    
      
        Eq
        .
        
          
        
        ⁢
        30
      
    
  




where C represents the circumference, n is the number of electrodes, D is a known distance between two adjacent electrodes, Rn (n=8 in FIG. 5) is the resistance between the first and last electrode when the electrodes form a circle on any segment, and λ is the ratio of resistivity to the area of cross sectional area which can be calculated by





  
    
      
        λ
        =
        
          
            1
            
              n
              -
              1
            
          
          ⁢
          
            
              ∑
              
                i
                =
                1
              
              
                n
                -
                1
              
            
            ⁢
            
              
                R
                i
              
              
                D
                i
              
            
          
        
      
    
    
      
        Eq
        .
        
          
        
        ⁢
        31
      
    
  




where Ri is the resistance between any adjacent two electrodes except the resistance between the first and end, Di is the distance between any adjacent two electrodes except the distance between the first and the last electrodes.

Thus, specific resistivity in skin of this segment is

ρ′=λ×A  Eq. 32

where A is the cross sectional area in the segment. Since A=C2/(4*π), ρ′ can be calculated by

ρ′=λ*C2/(4*π)  Eq. 33

Once the skin resistivity is determined, the length can be calculated by the equation

L=A*RL/ρ′  Eq. 34

where L is the length and A is the cross sectional area of the segment and RL is the resistance between two cross sectional areas along the vertical axis of the segment.

In FIG. 4, E1 to E8 are electrodes, D1 to D7 represent the distance between two adjacent electrodes respectively, and D1=D2= . . . =D7. In FIG. 5, R1 to R7 represents the skin resistance between adjacent two electrodes. R8 is the resistance between, as discussed above, the first (E1) and the last electrode (E8). The distance between E1 and E8 is unknown but it, as discussed above, can be calculated by R8/λ and therefore, the circumference (C) can be calculated by Eq. 30.
The mathematical operations described herein can be performed using a variety of computers and software. For example, those operations can be performed using the commercially available MATLAB program discussed above and a personal computer configured to run that program in accordance with the program manufacturer's specifications. The program as customized to perform the mathematical operations of the invention can be embodied as an article of manufacture comprising a computer usable medium, such as a magnetic disc, an optical disc, or the like, upon which the program is encoded. The data generated by the programs can similarly be stored on various types of storage media.
Without intending to limit it in any manner, the present invention will be more fully described by the following examples.
Example 1
Effect of Fat Mass on Bioimpedance Analysis
This example shows that fat mass is a major factor affecting the accurate measurement of body composition, including the accurate measurement of body fluid in clinical practice in dialysis patients, using bioimpedance analysis.
28 chronic hemodialysis patients were studied before treatment. Table 2 sets forth the relevant parameters for this group of subjects (patients).
Calf bioimpedance measurements were performed using BIS device Xitron 4200 (Xitron Technologies, San Diego, Calif.). The procedure of placement of electrodes was published in a previous study [9]. As described therein, to place the electrodes in the normal position, the distance (H) between the central part of the patella and the center of the lateral malleolus was measured while the subject was in the sitting position (see FIG. 6). With the patient recumbent, two measuring electrodes were placed, one 5 cm above the mid point M and another 5 cm below. Two current injecting electrodes were placed, one 3 cm above and the other 3 cm below the measuring electrodes. The circumference of the calf at the two measuring electrodes was measured using a tape with 0.1 cm precision.
Calf resistance and reactance were measured using the BIS device with 50 logarithmically distributed frequencies from 5 kHz to 1 MHz. Extracellular resistance (Re) and intracellular resistance (Ri) were calculated and ECV and ICV were obtained using a program provided by the BIS device [9, 10]. Fat mass, muscle mass, and bone were separately measured in the same area by MRI.
Geometric volume (VG) was defined as the cross sectional area (AG) times the length (L) of the segment. AG can be calculated by AG=C2/(4π), where C is average value of circumferences in the calf. Resistivity (ρG,Mea) was calculated by the value of measurement shown as follows:




  
    
      
        
          ρ
          
            G
            ,
            Mea
          
        
        =
        
          
            
              R
              G
            
            ×
            
              A
              G
            
          
          L
        
      
    
    
      
        Eq
        .
        
          
        
        ⁢
        35
      
    
  




where RG is resistance measured at 5 kHz, AG is cross sectional area and length was fixed at 10 cm for all subjects.

In order to compare the effect of volume of fat mass on the value of resistivity from both measurement and calculation, the 28 patients were divided into two groups according to the ratio of fat mass (VF) to geometric volume (VG): the high fat group was defined as VF/VG>0.2 and the normal fat group was defined as VF/VG≦0.2. Data is presented as mean value±SD, and the student t-test was used to compare the data from different groups. A difference between groups was considered significant if the p value<0.05.
Table 3 shows results of comparison of the resistivity between the two groups. There were significant differences in ρG,Mea (p<0.05), ρG,Cal (p<0.005) and Δρ (p<0.05), however, no significant difference was seen in ρ*G,Cal between the two groups (Table 3). The difference (p<0.05) between ρG,Mea and ρ*G,cal was significant but no significant difference was seen between ρG,Mea and ρG,Cal.
FIG. 7 shows correlations between calculated and measured resistivity values using Equations 12 and 13. The figure shows a high degree of correlation between the resistivity (ρG,Mea) determined from Eq. 35 which is based on the measured value and the resistivity (ρG,Cal) calculated using Eq. 12. In FIG. 7, the diamond symbols represent the correlation of ρG,Mea with ρ*G,Cal by Eq. 13 (ρ*G,Cal=0.66ρG,Mea+76.7, R2=0.96) and the symbols of solid circles represent the correlation between ρG,Mea and ρG,Cal with individual various kE calculated by Eq. 12 (ρG,Cal=0.93ρG,Mea+41.1, R2=0.9). Even though the correlation is high between ρ*G,Cal and ρG,Mea, the slope (diamond symbol) is only 0.66. After calibration with individual kE using Eq. 12, the slope is 0.93 (solid circles), i.e., significantly better.
Knowledge of the relationship of resistivity measured at the skin to the various subcutaneous tissues is essential to understand the relevant bioelectric phenomena. Many early studies have shown the resistivity of various biologic tissues and organs by direct measurement [6]. However, since clinical measurements are performed on the surface of skin and the resistivity of different tissues is calculated by a number of theoretical models, errors may occur due to: (1) individual different impedance interfaces of electrode and skin; (2) different volumes of adipose tissue between individuals and (3) theoretical models that do not include the variables of differences in skin impedance and fat mass between individuals.
The present invention provides a model which includes the variable of impedance of fat mass and may include the variables of resistance of skin and interface between capacitance. Table 3 shows that there is significant difference in resistivity (Δp) at 5 kHz between calculation by the model and by measurement as the amount of fat mass increases. This suggests that the value of resistivity at low frequency (5 kHz) obtained from measurement on the skin should include the contribution not only from extracellular fluid compartment but also from the fat mass compartment. FIG. 7 shows that after calibration with individual kE the slope of the curve as shown by the solid circles is almost equal to one, which indicates that the calculated value is close to the measurement values. It is important to understand that the variability of resistivity in healthy subjects depends mainly on the individual volume of fat mass. Therefore using this model, correct information concerning ECV can be obtained by measurement of, for example, calf resistivity. This information is important to provide, among other things, an accurate parameter of hydration for adjusting the variably abnormal hydration in dialysis patients; this will permit targeting appropriate body weight by removal of excess body fluid during dialysis.
In summary, this example illustrates some of the aspects of the invention that provide a model which can be used as a basis to calculate the resistivity in different tissues from the measurement of resistivity at the skin surface. That is, this example illustrates some of the aspects of the invention that provide an electrical resistivity model of segmented body composition for use in bioimpedance analysis, where the model describes the relationship of the resistivity with measurement on the skin surface and by calculation of tissues composing a limb.
In particular, in this example, a specific analysis was performed in a group of subjects with 5 kHz current injected at the surface of the calf. The study showed that calculated resistivity using the model was highly correlated with the values of resistivity by actual measurement. In the next example, the response of the system to current of different frequencies is investigated. The results of that example show that the resistivity of cell membranes can be tested using the three parallel path model of the invention.
Example 2
Multi-Frequency Analysis
To improve existing BIS techniques, a correct electrical model should be able to explain the electrical properties with different proportions of body tissues over a wide range of current frequencies. The aim of this study was to evaluate frequency response of the equivalent circuit model from 5 kHz to 1 MHz measured at the skin of the calf. The equivalent circuit model is shown in FIG. 2B.
The same 28 chronic hemodialysis patients as in Example 1 were studied before treatment (see Table 2). Again, calf bioimpedance measurements were performed using BIS device Xitron 4200 (Xitron Technologies, San Diego, Calif.) and the placement of electrodes was as described above in Example 1 (see FIG. 6).
Calf resistance and reactance were measured using BIS device with 50 logarithmically distributed frequencies from 5 kHz to 1 MHz. Fat mass, muscle mass and bone were separately measured in the same area of BIS measurement by MRI. Patients were divided into three groups according to the amount of fat mass in the calf with G1: fat>0.4 kg, G2: 0.25<fat<0.4 kg and G3: fat<0.25 kg.
Significant differences in resistance, reactance and impedance between G1 and G3 were found; however, the significant difference between G1 and G2 show only with frequencies higher than 40 kHz (FIG. 8). Results of curve fitting to the raw data (circle, triangle and square) of impedance with different frequencies are shown in FIG. 9. The data shown by circles is for a patient having a fat mass of 0.59 kg as estimated by MRI, while the triangles and squares show data for two other patients having fat masses of 0.29 kg and 0.23 kg, respectively, again estimated by MRI.
The frequency response curve (A) and phase response (B) are shown in FIGS. 10 and 11 from two patients: one had 0.55 kg fat mass, 1.046 kg muscle mass and 0.021 kg bone (FIG. 10) as estimated by MRI and the other had 0.134 kg fat mass, 1.022 kg muscle mass and 0.022 kg bone (FIG. 11), again estimated by MRI.
It is clearly observed that impedance of the patient with smaller fat mass decreases in the range of frequency from 1 to 10 kHz (FIG. 11A), however, the curve is almost constant in the same frequency range when patient has a larger fat mass (FIG. 10A). In addition, at the minimum point of phase, the magnitude of the frequency is different between the two patients.
FIG. 12 shows a correlation of muscle (MS) estimation by BIS and by MRI (MSBIS=285*L2/Ri, where L (10 cm) is the distance between measuring electrodes along the patient's calf and Ri is intracellular resistance which was calculated using the model of FIG. 2B). The value of 285 Ω-cm used as muscle resistivity in calculating MSBIS was obtained by a regression analysis based on the MRI data.
In generating FIG. 12, for each patient, Ri was determined by solving Equations 23-29 using the P1, P2, P3, P4, Q1, Q2, and Q3 values obtained for that patient from the measured frequency response data. FIG. 12 thus illustrates the use of a group of parameters to quantify the analysis of the frequency response of the components of an equivalent circuit which can be used to identify the differences in body composition of individual subjects.
Differences in P1, P2, P3, P4, Q1, Q2, and Q3 values were investigated for groups G1, G2, and G3 discussed above. Table 4 shows the significant difference (p<0.05) of P1 between G1 and G2 and the significant difference of P4 between G1 and G3. There were no significant differences except for Pi and P4 in the P and Q parameters.
Example 1 above showed that electrical resistivity of measurement on the skin surface is affected by adipose tissue using the equivalent circuit model of FIG. 2. This is believed to be one of the more important factors which influence the accuracy of bioimpedance technique to estimate body composition such as extracellular and intracellular fluid volumes. Of course, bioimpedance accuracy is not only influenced by the amount of fat mass but also by the interface between electrodes and skin or the degree of skin hydration.
In this example, the results shown in FIG. 9 demonstrate a decrease in impedance related to increase of the current frequency in patients with different fat content. Impedance of the patient with larger fat mass is higher than that of the patient with small fat mass. Simulation results of two patients with similar muscle mass and bone mass show that the frequency response curve is different in the frequency range from 0 to 1 kHz. The difference between FIGS. 10 and 11 demonstrate that the change in curve in response to the frequency from 0 to 1 kHz depends on the fat mass and that will largely affect the estimation of extracellular and intracellular resistance by multi-frequency bioimpedance analysis. Moreover, the difference at the minimum point of phase response between the two patients indicates that using 50 kHz single frequency bioimpedance method could produce error when subject has a larger fat mass.
The correlation between muscle estimation by BIS and by MRI shown in FIG. 12 has a number of important implications for clinical and other applications. Thus, as illustrated in this figure, the electrical model of the present invention correlates well with actual muscle mass as measured by MRI. The model therefore allows one to determine the muscle mass of individual subjects (patients) using only simple, inexpensive, and non-invasive electrical measurements.
The particular muscle mass determined using the model will depend on the particular locations chosen for the electrodes used to make the BIS measurement. For example, the data of FIG. 12 is for the entire calf muscle. By using a different electrode placement, the muscle mass of, for example, a portion of the calf muscle can be determined, e.g., the muscle mass of the gastrocnemius portion of the calf. Alternatively, rather than making measurements on the calf muscle, measurements can be made on other muscles, portions of muscles, and/or or muscle groups, e.g., measurements can be made on all or part of a subject's biceps. If only relative measurements are needed, the bioimpedance measurement can be used as is. Alternatively, correlations of the type shown in FIG. 12 can be obtained between bioimpedance measurements and MRI measurements for particular muscles, portions of muscles, or muscle groups, thus allowing the electrical measurements to provide “absolute” values for muscle mass, where “absolute” values are preferably those which correspond to those obtained using MRI, the acknowledged “gold standard” for muscle mass measurements.
The ability to measure muscle mass, discrete muscles or muscle groups with suitable application of electrodes, with or without validation against MRI measurements, has numerous applications. For example, these techniques are applicable in exercise programs in the home, gymnasium, sports and health clubs and in rehabilitation after surgery and injury, where the effect of muscle mass increases is relevant. Thus, by taking a series of bioimpedance measurements over time, subjects (patients) and/or their health care professionals can monitor changes in muscle mass as a result of exercise programs, diet changes, and/or rehabilitation programs.
Example 3
Single Low Frequency Bioimpedance Measurements
FIG. 13 is a plot based on the data of Examples 1 and 2 which shows the lack of a correlation between measured resistance values at 5 kilohertz and muscle mass determined by MRI. This data supports the conclusion that at low frequencies, measured bioimpedance data is not responsive to the muscle content of a segment. Rather, as set forth in the circuit model of FIG. 3, at these frequencies, the components of a segment that primarily determine the measured values obtained with a bioimpedance procedure are the fat and extracellular fluid components.
Example 4
Comparative Example Multi-Frequency, Three Parallel Path Model Versus 50 Kilohertz
FIGS. 14 and 15 compare correlations between bioimpedance and MRI measurements achieved using the multi-frequency, three parallel path model of FIG. 2 with correlations achieved using impedance at a single frequency, namely, the commonly used frequency of 50 kilohertz. The data used in preparing these figures was the data discussed above in connection with Examples 1 and 2.
FIG. 14 shows muscle mass correlations, while FIG. 15 shows fat mass correlations. As illustrated by these figures, the multi-frequency, three parallel path equivalent circuit frequency response model of the invention achieves higher correlations than the single 50 kilohertz approach both for muscle mass (R2 of 0.6 for the invention versus 0.3 for 50 kilohertz) and for fat mass (R2 of 0.7 for the invention versus 0.2 for 50 kilohertz).
Example 5
Comparative Example Multi-Frequency, Three Parallel Path Model Versus Cole-Cole Model
FIGS. 16 and 17 compare correlations between BIS and MRI measurements achieved using the multi-frequency, three parallel path model of FIG. 2 with correlations achieved using the Cole-Cole model. The data used in preparing these figures was the data discussed above in connection with Examples 1 and 2.
FIG. 16 shows muscle mass correlations, while FIG. 17 shows fat mass correlations. As illustrated by these figures, the multi-frequency, three parallel path equivalent circuit frequency response model of the invention achieves higher correlations than the Cole-Cole model both for muscle mass (R2 of 0.6 for the invention versus 0.4 for the Cole-Cole model) and for fat mass (R2 of 0.7 for the invention versus 0.1 for the Cole-Cole model).
Although specific embodiments of the invention have been described and illustrated, it is to be understood that a variety of modifications which do not depart from the scope and spirit of the invention will be evident to persons of ordinary skill in the art from the foregoing disclosure. As just one example, although the following claims recite various features of the invention, it is to be understood that the invention encompasses any and all combinations of those features, irrespective of whether such combination is currently set forth in the appended set of claims.
REFERENCES
Citations for the various documents referred to above are set forth below. The contents of these documents are incorporated herein by reference.

    [1] K. S. Cole and R. H. Cole, “Dispersion and absorption in dielectrics. I. Alternating current characteristics” J. chem. Phys. Vol. 9, pp. 341-351, 1941
    [2] H. P. Schwan, K. Li, “A dielectric study of the low-conductance surface membrane,” in E. coli Nature Vol. 177, pp. 134-135, 1956
    [3] J. Nyboer, Electrical impedance plethysmography, 2nd ed., Charles C. Thomas, Springfield, Ill. 1970
    [4] R. P. Patterson, “Fundamentals of impedance cardiography,” IEEE Engineering in Medicine and Biology magazine. Vol. 8, pp. 16-18, 1989
    [5] T. Hanai, Electrical properties of emulsions In: Emulsion Science, edited by P. H. Sherman London: Academic, 1968, pp. 354-477,
    [6] A. De Lorenzo, A. Andreoli, J. R. Matthie, and P. O. Withers, “Predicting body cell mass with bioimpedance by using theoretical methods: a technological review,” J Appl Physiol Vol. 82, pp. 1542-1558, 1997
    [7] R. N. Baumgartner, R. Ross and S. B. Heymsfield, “Does adipose tissue influence bioelectric impedance in obese men and women?” J. Appl. Physiol. Vol. 84, pp. 257-262, 1998.
    [8] K. R. Foster and H. C. Lukaski, “Whole-body impedance-what does it measure?” Am J Clin Nutr Vol. 64 (suppl), pp. 388S-396S, 1996
    [9] F. Zhu, S. Sarkar, C. Kaitwatcharachai, R. Greenwood, C. Ronco, N. W. Levin, “Methods and reproducibility of measurement of resistivity in the calf using regional bioimpedance analysis. Blood Purif,” Vol. 21, pp. 131-136, 2003
    [10] Xitron Technologies, Inc., “4000B Bio-Impedance Spectrum Analyzer System Operating Manual,” preliminary edition, San Diego, Calif., 1995, Appendix A, pages 50-61.












TABLE 1







Permittivity
Resistivity



ε
ρ (Ω-cm)




















Bone
640
104



Fat
3 × 104
1.5~5 × 103



Blood
2.8 × 103  
 1.5 × 102



Muscle (parallel)
8 × 104
   2 × 102

















TABLE 2







SUBJECT INFORMATION












Mean
SD
min
max

















SEX
F9/M19






AGE, year
53.4
10.5
33
69



WEIGHT, kg
80.4
18
43.2
119.9



Height, cm
169.7
9.5
149
184.9



BMI, kg/m2
27.7
4.97
19.18
41.11



FAT, g
344.9
118
149.2
533.5



MUSCLE, g
525.3
110.3
326.5
761.5





















TABLE 3






ρG, Mea
ρ*G, Cal
ρG, Cal
Δρ


VF/VG
(Ω-cm)
(Ω-cm)
(Ω-cm)
(Ω-cm)







 >0.2
430.4 ± 62
359.9 ± 43
457.2 ± 64
70.5 ± 24


=<0.2
369.5 ± 84
322.7 ± 57
373.5 ± 70
46.8 ± 29


p value
<0.05
n.s.
<0.005
<0.05





ρG, Cal and ρ*G, Cal are the resistivity values calculated using Eq. 12 and Eq. 13, respectively;


Δρ = ρG, Mea − ρ*G, Cal; a ρE value of 90 Ω□cm was used in calculating ρ*G, Cal and ρG, Cal.













TABLE 4







SUMMARY OF PARAMETERS FROM CURVE FITTING














Fat
P1
P2
P3
P4
Q1
Q2
Q3





G1
32.54 ± 5.9 * 
1412 ± 262
−15330 ± 3352
42640 ± 12130 +
22.5 ± 6.7
−276.4 ± 76 
649.8 ± 430


G2
27.3 ± 3.7 *
1187 ± 548
−13012 ± 3027
34166 ± 13808  
22.6 ± 7.7
 −262 ± 73
 718 ± 372


G3
26 ± 4.4
1168 ± 259
−10973 ± 3448
24191 ± 16426 +
20.7 ± 8.7
−197.2 ± 102
460.3 ± 460





* and + indicate a significant (p < 0.05) difference between groups.

## Claims

1. A method for determining the circumference of a portion of a body segment covered by skin comprising:
(a) applying a series of electrodes around said portion, said series of electrodes having a first electrode and a last electrode, the circumferential distances between all electrodes in the series being known, except for the distance between the first and last electrodes;
(b) determining the resistance between at least two electrodes of the series, other than the first and last electrode, by applying a low frequency current which does not substantially penetrate the skin;
(c) determining a resistivity value per unit length for the skin from the resistance determined in step (b) and the known circumferential distance between the two electrodes;
(d) determining the resistance between the first and last electrodes of the series by applying a low frequency current which does not substantially penetrate the skin; and
(e) calculating the distance between the first and last electrodes of the series from the resistance measured in step (d) and the resistivity value per unit length determined in step (c).

2. The method of claim 1 wherein the resistances between all adjacent pairs of electrodes, other than the first and last electrodes, are determined in step (b) and used in step (c) to determine the resistivity value.

3. The method of claim 1 wherein, except for the first and last electrodes, the distances between adjacent pairs of electrodes are equal.

4. The method of claim 1 wherein the series of electrodes are carried by a belt.

5. Apparatus comprising a programmed computer for determining the circumference of a portion of a body segment covered by skin, said apparatus being adapted for use with a series of electrodes for application around said portion of a body segment, said series of electrodes having a first electrode and a last electrode, the circumferential distances between all electrodes in the series being known, except for the distance between the first and last electrodes, said programmed computer comprising a non-transitory computer readable storage medium having computer executable code embodied therein for performing the steps of:
(i) determining the resistance between at least two electrodes of the series of electrodes, other than the first and last electrode, by applying a low frequency current which does not substantially penetrate the skin;
(ii) determining a resistivity value per unit length for the skin from the resistance determined in step (i) and the known circumferential distance between the two electrodes;
(iii) determining the resistance between the first and last electrodes of the series by applying a low frequency current which does not substantially penetrate the skin; and
(iv) calculating the distance between the first and last electrodes of the series from the resistance determined in step (iii) and the resistivity value per unit length determined in step (ii).