Method for determining and representing food products based on nutrient density rating and predicted satiating effect

• US 7,620,531 B1
• Filed: 12/12/2005
• Issued: 11/17/2009
• Est. Priority Date: 12/12/2005
• Status: Active Grant

• Current Assignee: Advance Magazine Publishers Incorporated (Advance Publications Incorporated)
• Original Assignee: Condenet Incorporated (Advance Publications Incorporated)
• Inventors: Johnson, Ronald B.
• Primary Examiner(s): Moran; Marjorie
• Assistant Examiner(s): Skowronek; Karlheinz R
• Application Number:
  US11/302,080
• Time in Patent Office: 1,436 Days
• Field of Search: None
• US Class Current: 703/2
• CPC Class Codes: G01N 33/02 Food
• CPC Class Codes: G09B 7/02 of the type wherein the student is expected to construct an answer to the question which is presented or wherein the machine gives an answer to the question presented by a student

Abstract

A method of assisting individuals in making food choices by providing two nutritional indices. One is a numerical expression of a food'"'"'s overall nutrient density and the other is a separate numerical expression that relates to a food'"'"'s caloric density that is a prediction of the satiating effect of a food. A visual aid in the form of a graph or chart is provided so individuals using the indices can refer to the chart to determine food selections in accordance with the individual'"'"'s dietary goals.

Claims

1. A method of rating food, comprising the steps of:

(a) calculating a first factor FF that represents an estimate of the satiating effect of a food according to the relationship;

FF=max(FFMIN,min(FFMAX,C00

+C01*max(CALMIN,CAL)^E01

+C02*min(PRMAX,PR)^E02

+C03*min(DFMAX,DF)^E03

+C04*min(TFMAX,TF)^E04))

where;

max(CALMIN,CAL) is a function that returns the maximum of either CALMIN or CAL,

min(PRMAX,PR) is a function that returns the minimum of either PRMAX or PR,

min(DFMAX,DF) is a function that returns the minimum of either DFMAX or DF,

min(TFMAX,TF) is a function that returns the minimum of either TFMAX or TF,

max(FFMIN,min(FFMAX,C00+C01*max(CALMIN,CAL)^E1+C02*min(PRMAX,PR)^E02+C03*min(DFMAX,DF)^E03+C04*min(TFMAX,TF)^E04)) is a function that returns the maximum of either (i) FFMIN or (ii) the function that returns the minimum of either FFMAX or the quantity C00+C01*max(CALMIN,CAL)^E01+C02*min(PRMAX,PR)^E02+C03*min(DFMAX,DF)^E03+C04*min(TFMAX,TF)^04

CAL is total Calories per 100 gram portion of the food,

PR is grams of Protein per 100 g,

DF is grams Dietary Fiber per 100 g, and

TF is grams total Fat per 100 g;

FFMIN, FFMAX, CALMIN, PRMAX, DFMAX, TFMAX, and C00 are constants, C01 through C04 are coefficients, and E01 through E04 are exponential powers having the approximate following values;

FFMIN=0.5

FFMAX=5.0

CALMIN=30

DFMAX=12

C03=1/1620

C04=−

1/138,000

E02=1.0

E04=3.0;

(b) calculating a second factor NDR that represents the nutritional benefit of said food according to the relationship;

NDR=max(RMIN,min(RMAX,C00

+C01*(ln(max(PMIN,min(PMAX,100*DF0/DF1)))−

ln(PAVG))

+C02*(ln(max(PMIN,min(PMAX,100*PR0/PR1)))−

ln(PAVG))

+C03*(ln(max(PMIN,min(PMAX,100*VA0/VA1)))−

ln(PAVG))

+C04*(ln(max(PMIN,min(PMAX,100*VC0/VC1)))−

ln(PAVG))

+C05*(ln(max(PMIN,min(PMAX,100*CA0/CA1)))−

ln(PAVG))

+C06*(ln(max(PMIN,min(PMAX,100*FE0/FE1)))−

ln(PAVG))

+C07*(ln(max(PMIN,min(PMAX,100*VD0/VD1)))−

ln(PAVG))

+C08*(ln(max(PMIN,min(PMAX,100*VE0/VE1)))−

ln(PAVG))

+C09*(ln(max(PMIN,min(PMAX,100*VK0/VK1)))−

ln(PAVG))

+C10*(ln(max(PMIN,min(PMAX,100*TH0/TH1)))−

ln(PAVG))

+C11*(ln(max(PMIN,min(PMAX,100*RI0/RI1)))−

ln(PAVG))

+C12*(ln(max(PMIN,min(PMAX,100*NI0/NI1)))−

ln(PAVG))

+C13*(ln(max(PMIN,min(PMAX,100*VB60/VB61)))−

ln(PAVG))

+C14*(ln(max(PMIN,min(PMAX,100*FO0/FO1)))−

ln(PAVG))

+C15*(ln(max(PMIN,min(PMAX,100*VB120/VB121)))−

ln(PAVG))

+C16*(ln(max(PMIN,min(PMAX,100*PA0/PA1)))−

ln(PAVG))

+C17*(ln(max(PMIN,min(PMAX,100*MG0/MG1)))−

ln(PAVG))

+C18*(ln(max(PMIN,min(PMAX,100*P0/P1)))−

ln(PAVG))

+C19*(ln(max(PMIN,min(PMAX,100*K0/K1)))−

ln(PAVG))

+C20*(ln(max(PMIN,min(PMAX,100*ZN0/ZN1)))−

ln(PAVG))

+C21*(ln(max(PMIN,min(PMAX,100*CU0/CU1)))−

ln(PAVG))

+C22*(ln(max(PMIN,min(PMAX,100*MN0/MN1)))−

ln(PAVG))

+C23*(ln(max(PMIN,min(PMAX,100*SE0/SE1)))−

ln(PAVG))

+C24*(ln(max(PMIN,min(PMAX,100*SF0/SF1)))−

ln(PAVG))

+C25*(ln(max(PMIN,min(PMAX,100*CH0/CH1)))−

ln(PAVG))

+C26*(ln(max(PMIN,min(PMAX,100*NA0NA1)))−

ln(PAVG))

+C27*(ln(max(PMIN,min(PMAX,100*SA0/SA1)))−

ln(PAVG))))

where;

ln is the natural (base e) logarithm,

max(PMIN,min(PMAX,100*DF0/DF1)) is a function that returns the maximum of either (i) PMIN or (ii) the function that returns the minimum of either PMAX or 100*DF0/DF1,

max(PMIN min(PMAX,100*PR0/PR1)) is a function that returns the maximum of either (i) PMIN or (ii) the function that returns the minimum of either PMAX or 100*PR0/PR1,

max(PMIN min(PMAX,100*VA0/VA1)) is a function that returns the maximum of either (i) PMIN or (ii) the function that returns the minimum of either PMAX or 100*VA0/VA1,

max(PMIN,min(PMAX,100*VC0/VC1)) is a function that returns the maximum of either (i) PMIN or (ii) the function that returns the minimum of either PMAX or 100*VC0/VC1,

max(PMIN,min(PMAX,100*CA0/CA1)) is a function that returns the maximum of either (i) PMIN or (ii) the function that returns the minimum of either PMAX or 100*CA0/CA1,

max(PMIN,min(PMAX,100*FE0/FE1)) is a function that returns the maximum of either (i) PMIN or (ii) the function that returns the minimum of either PMAX or 100*FE0/FE1,

max(PMIN min(PMAX,100*VD0/VD1)) is a function that returns the maximum of either (i) PMIN or (ii) the function that returns the minimum of either PMAX or 100*VD0/VD1,

max(PMIN,min(PMAX,100*VE0/VE1)) is a function that returns the maximum of either (i) PMIN or (ii) the function that returns the minimum of either PMAX or 100*VE0/VE1,

max(PMIN,min(PMAX,100*VK0/VK1)) is a function that returns the maximum of either (i) PMIN or (ii) the function that returns the minimum of either PMAX or 100*VK0/VK1,

max(PMIN,min(PMAX,100*TH0/TH1)) is a function that returns the maximum of either (i) PMIN or (ii) the function that returns the minimum of either PMAX or 100*TH0/TH1,

max(PMIN,min(PMAX,100*RI0/RI1)) is a function that returns the maximum of either (i) PMIN or (ii) the function that returns the minimum of either PMAX or 100*RI0/RI1,

max(PMIN,min(PMAX,100*NI0/NI1)) is a function that returns the maximum of either (i) PMIN or (ii) the function that returns the minimum of either PMAX or 100*NI0/NI1,

max(PMIN,min(PMAX,100*VB60/VB61)) is a function that returns the maximum of either (i) PMIN or (ii) the function that returns the minimum of either PMAX or 100*VB60/VB61,

max(PMIN,min(PMAX,100*FO0/FO1)) is a function that returns the maximum of either (i) PMIN or (ii) the function that returns the minimum of either PMAX or 100*FO0/FO1,

max(PMIN,min(PMAX,100*VB120/VB121)) is a function that returns the maximum of either (i) PMIN or (ii) the function that returns the minimum of either PMAX or 100*VB120/VB121,

max(PMIN,min(PMAX,100*PA0/PA1)) is a function that returns the maximum of either (i) PMIN or (ii) the function that returns the minimum of either PMAX or 100*PA0/PA1,

max(PMIN,min(PMAX,100*MG0/MG1)) is a function that returns the maximum of either (i) PMIN or (ii) the function that returns the minimum of either PMAX or 100*MG0/MG1,

max(PMIN,min(PMAX,100*P0/P1)) is a function that returns the maximum of either (i) PMIN or (ii) the function that returns the minimum of either PMAX or 100*P0/P1,

max(PMIN,min(PMAX,100*K0/K1)) is a function that returns the maximum of either (i) PMIN or (ii) the function that returns the minimum of either PMAX or 100*K0/K1,

max(PMIN,min(PMAX,100*ZN0/ZN1)) is a function that returns the maximum of either (i) PMIN or (ii) the function that returns the minimum of either PMAX or 100*ZN0/ZN1,

max(PMIN,min(PMAX,100*CU0/CU1)) is a function that returns the maximum of either (i) PMIN or (ii) the function that returns the minimum of either PMAX or 100*CU0/CU1,

max(PMIN,min(PMAX,100*MN0/MN1)) is a function that returns the maximum of either (i) PMIN or (ii) the function that returns the minimum of either PMAX or 100*MN0/MN1,

max(PMIN,min(PMAX,100*SE0/SE1)) is a function that returns the maximum of either (i) PMIN or (ii) the function that returns the minimum of either PMAX or 100*SE0/SE1,

max(PMIN,min(PMAX,100*SF0/SF1)) is a function that returns the maximum of either (i) PMIN or (ii) the function that returns the minimum of either PMAX or 100*SF0/SF1,

max(PMIN,min(PMAX,100*CH0/CH1)) is a function that returns the maximum of either (i) PMIN or (ii) the function that returns the minimum of either PMAX or 100*CH0/CH1,

max(PMIN,min(PMAX,100*NA0/NA1)) is a function that returns the maximum of either (i) PMIN or (ii) the function that returns the minimum of either PMAX or 100*NA0/NA1,

max(PMIN,min(PMAX,100*SA0/SA1)) is a function that returns the maximum of either (i) PMIN or (ii) the function that returns the minimum of either PMAX or 100*SA0/SA1,

max(RMIN,min(RMAX,C00+C01*(ln(max(PMIN,min(PMAX,100*DF0/DF1)))−

ln(PAVG))+C02*(ln(max(PMIN,min(PMAX,100*PR0/PR1)))−

ln(PAVG))+C03*(ln(max(PMIN,min(PMAX,100*VA0/VA1)))−

ln(PAVG))+C04*(ln(max(PMIN,min(PMAX,100*VC0/VC1)))−

ln(PAVG))+C05*(ln(max(PMIN,min(PMAX,100*CA0/CA1)))−

ln(PAVG))+C06*(ln(max(PMIN,min(PMAX,100*FE0/FE1)))−

ln(PAVG))+C07*(ln(max(PMIN,min(PMAX,100*VD0/VD1)))−

ln(PAVG))+C08*(ln(max(PMIN,min(PMAX,100*VE0/VE1)))−

ln(PAVG))+C09*(ln(max(PMIN,min(PMAX,100*VK0/VK1)))−

ln(PAVG))+C10*(ln(max(PMIN,min(PMAX,100*TH0/TH1)))−

ln(PAVG))+C11*(ln(max(PMIN,min(PMAX,100*RI0/RI1)))−

ln(PAVG))+C12*(ln(max(PMIN,min(PMAX,100*NI0/NI1)))−

ln(PAVG))+C13*(ln(max(PMIN,min(PMAX,100*VB60/VB61)))−

ln(PAVG))+C14*(ln(max(PMIN,min(PMAX,100*FO0/FO1)))−

ln(PAVG))+C15*(ln(max(PMIN,min(PMAX,100*VB120/VB121)))−

ln(PAVG))+C16*(ln(max(PMIN,min(PMAX,100*PA0/PA1)))−

ln(PAVG))+C17*(ln(max(PMIN,min(PMAX,100*MG0/MG1)))−

ln(PAVG))+C18*(ln(max(PMIN,min(PMAX,100*P0/P1)))−

ln(PAVG))+C19*(ln(max(PMIN,min(PMAX,100*K0/K1)))−

ln(PAVG))+C20*(ln(max(PMIN,min(PMAX,100*ZN0/ZN1)))−

ln(PAVG))+C21*(ln(max(PMIN,min(PMAX,100*CU0/CU1)))−

ln(PAVG))+C22*(ln(max(PMIN,min(PMAX,100*MN0/MN1)))−

ln(PAVG))+C23*(ln(max(PMIN,min(PMAX,100*SE0/SE1)))−

ln(PAVG))+C24*(ln(max(PMIN,min(PMAX,100*SF0/SF11)))−

ln(PAVG))+C25*(ln(max(PMIN,min(PMAX,100*CH0/CH1)))−

ln(PAVG))+C26*(ln(max(PMIN,min(PMAX,100*NA0/NA1)))−

ln(PAVG))+C27*(ln(max(PMIN,min(PMAX,100*SA0/SA1)))−

ln(PAVG)))) is a function that returns the maximum of either (i) RMIN or (ii) the function that returns the minimum of either RMAX or the quantity C00+C01*(ln(max(PMIN,min(PMAX,100*DF0/DF1)))−

ln(PAVG))+C02*(ln(max(PMIN,min(PMAX,100*PR0/PR1)))−

ln(PAVG))+C03*(ln(max(PMIN,min(PMAX,100*VA0/VA1)))−

ln(PAVG))+C04*(ln(max(PMIN,min(PMAX,100*VC0/VC1)))−

ln(PAVG))+C05*(ln(max(PMIN,min(PMAX,100*CA0/CA1)))−

ln(PAVG))+C06*(ln(max(PMIN,min(PMAX,100*FE0/FE1)))−

ln(PAVG))+C07*(ln(max(PMIN,min(PMAX,100*VD0/VD1)))−

ln(PAVG))+C08*(ln(max(PMIN,min(PMAX,100*VE0/VE1)))−

ln(PAVG))+C09*(ln(max(PMIN,min(PMAX,100*VK0/VK1)))−

ln(PAVG))+C10*(ln(max(PMIN,min(PMAX,100*TH0/TH1)))−

ln(PAVG))+C11*(ln(max(PMIN,min(PMAX,100*RI0/RI1)))−

ln(PAVG))+C12*(ln(max(PMIN,min(PMAX,100*NI0/NI1)))−

ln(PAVG))+C13*(ln(max(PMIN,min(PMAX,100*VB60/VB61)))−

ln(PAVG))+C14*(ln(max(PMIN,min(PMAX,100*FO0/FO1)))−

ln(PAVG))+C15*(ln(max(PMIN,min(PMAX,100*VB120/VB121)))−

ln(PAVG))+C16*(ln(max(PMIN,min(PMAX,100*PA0/PA1)))−

ln(PAVG))+C17*(ln(max(PMIN,min(PMAX,100*MG0/MG1)))−

ln(PAVG))+C18*(In(max(PMIN,min(PMAX,100*P0/P1)))−

ln(PAVG))+C19*(ln(max(PMIN,min(PMAX,100*K0/K1)))−

ln(PAVG))+C20*(ln(max(PMIN,min(PMAX,100*ZN0/ZN1)))−

ln(PAVG))+C21*(ln(max(PMIN,min(PMAX,100*CU0/CU1)))−

ln(PAVG))+C22*(ln(max(PMIN,min(PMAX,100*MN0/MN1)))−

ln(PAVG))+C23*(ln(max(PMIN,min(PMAX,100*SE0/SE1)))−

ln(PAVG))+C24*(ln(max(PMIN,min(PMAX,100*SF0/SF1)))−

ln(PAVG))+C25*(ln(max(PMIN,min(PMAX,100*CH0/CH1)))−

ln(PAVG))+C26*(ln(max(PMIN,min(PMAX,100*NA0/NA1)))−

ln(PAVG))+C27*(ln(max(PMIN,min(PMAX,100*SA0/SA1)))−

ln(PAVG))

DF0 is the amount of Dietary Fiber present in a serving,

DF1 is the amount of Dietary Fiber specified as the Daily Value by the US FDA;

PR0 is the amount of Protein present in a serving;

PR1 is the amount of Protein specified as the Daily Value by the US FDA;

VA0 is the amount of Vitamin A present in a serving;

VA1 is the amount of Vitamin A specified as the Daily Value by the US FDA;

VC0 is the amount of Vitamin C present in a serving;

VC1 is the amount of Vitamin C specified as the Daily Value by the US FDA;

CA0 is the amount of Calcium present in a serving;

CA1 is the amount of Calcium specified as the Daily Value by the US FDA;

FE0 is the amount of Iron present in a serving;

FE1 is the amount of Iron specified as the Daily Value by the US FDA;

VD0 is the amount of Vitamin D present in a serving;

VD1 is the amount of Vitamin D specified as the Daily Value by the US FDA;

VE0 is the amount of Vitamin E present in a serving;

VE1 is the amount of Vitamin E specified as the Daily Value by the US FDA;

VK0 is the amount of Vitamin K present in a serving;

VK1 is the amount of Vitamin K specified as the Daily Value by the US FDA;

TH0 is the amount of Thiamin present in a serving;

TH1 is the amount of Thiamin specified as the Daily Value by the US FDA;

RI0 is the amount of Riboflavin present in a serving;

RI1 is the amount of Riboflavin specified as the Daily Value by the US FDA;

NI0 is the amount of Niacin present in a serving;

NI1 is the amount of Niacin specified as the Daily Value by the US FDA;

VB60 is the amount of Vitamin B6 present in a serving;

VB61 is the amount of Vitamin B6 specified as the Daily Value by the US FDA;

FO0 is the amount of Folate present in a serving;

FO1 is the amount of Folate specified as the Daily Value by the US FDA;

VB120 is the amount of Vitamin B12 present in a serving;

VB121 is the amount of Vitamin B12 specified as the Daily Value by the US FDA;

PA0 is the amount of Pantothenic Acid present in a serving;

PA1 is the amount of Pantothenic Acid specified as the Daily Value by the US FDA;

MG0 is the amount of Magnesium present in a serving;

MG1 is the amount of Magnesium specified as the Daily Value by the US FDA;

P0 is the amount of Phosphorus present in a serving;

P1 is the amount of Phosphorus specified as the Daily Value by the US FDA;

K0 is the amount of Potassium present in a serving;

K1 is the amount of Potassium specified as the Daily Value by the US FDA;

ZN0 is the amount of Zinc present in a serving;

ZN1 is the amount of Zinc specified as the Daily Value by the US FDA;

CU0 is the amount of Copper present in a serving;

CU1 is the amount of Copper specified as the Daily Value by the US FDA;

MN0 is the amount of Manganese present in a serving;

MN1 is the amount of Manganese specified as the Daily Value by the US FDA;

SE0 is the amount of Selenium present in a serving;

SE1 is the amount of Selenium specified as the Daily Value by the US FDA;

SF0 is the amount of Saturated Fat present in a serving;

SF1 is the amount of Saturated Fat specified as the Daily Value by the US FDA;

CH0 is the amount of Cholesterol present in a serving;

CH1 is the amount of Cholesterol specified as the Daily Value by the US FDA;

NA0 is the amount of Sodium present in a serving;

NA1 is the amount of Sodium specified as the Daily Value by the US FDA;

SA0 is the number of Calories in a serving from sugars (as reported on Nutrition Facts labels) and alcohol; and

SA1 is the total number of Calories in a serving;

wherein;

RMIN, RMAX, PMIN, PMAX, PAVG, and C00 are constants,

C01 thru C27 are coefficients that weight the effects of the individual nutrients on the NDR; and

wherein these constants and coefficients have the following approximate values;

RMIN=0.0

RMAX=5.0

PMAX=40

C00=3.0

C01=0.238

C02=0.238

C03=0.238

C04=0.238

C05=0.238

C06=0.238

C07=0.0238

C08=0.0238

C09=0.0238

C10=0.0238

C11=0.0238

C12=0.0238

C13=0.0238

C14=0.0238

C15=0.0238

C16=0.0238

C17=0.0238

C18=0.0238

C19=0.0238

C20=0.0238

C21=0.0238

C22=0.0238

C23=0.0238

C24=−

0.238

C25=−

0.238

C26=−

0.238

C27=−

0.238; and

(c) providing to a user a guideline using said first and second factors to enable food selection in accordance with at least one dietary goal;

wherein at least one of steps (a) and (b) is effected using a suitably programmed computer.

2. The method of claim 1, wherein said step of providing to a user a guideline includes plotting said factors NDR and FF on a rectilinear graph.

3. The method of claim 1, wherein said step of providing to a user a guideline includes providing said first and second factors in a visual format.

4. The method of claim 3, wherein said visual format includes regions based on nutritional characteristics of said food.

5. The method of claim 1, wherein said step of providing to a user a guideline includes providing a graphical representation including said first and second factors as coordinates.

6. The method of claim 5, wherein said graphical representation includes regions associated with different ones of said at least one dietary goal.

7. The method of claim 6, further comprising the step of combining said first and second factors into a single indicator corresponding to one of said regions.

8. The method of claim 1, further comprising the step of combining said first and second factors into a single indicator of suitability of said food for at least one dietary goal.

9. A food rating system, comprising:

(a) a data processor adapted to calculate;

(i) a first factor that represents an estimate of the satiating effect of a food according to the relationship;

FF=max(FFMIN,min(FFMAX,C00

+C01*max(CALMIN,CAL)^E01

+C02*min(PRMAX,PR)^E02

+C03*min(DFMAX,DF)^E03

+C04*min(TFMAX,TF)^E04))

where;

FF refers to “

fullness factor”

max(CALMIN,CAL) is a function that returns the maximum of either CALMIN or CAL,

min(PRMAX,PR) is a function that returns the minimum of either PRMAX or PR,

min(DFMAX,DF) is a function that returns the minimum of either DFMAX or DF,

min(TFMAX,TF) is a function that returns the minimum of either TFMAX or TF,

max(FFMIN,min(FFMAX,C00+C01*max(CALMIN,CAL)^E01+C02*min(PRMAX,PR)^E02+C03*min(DFMAX,DF)^E03+C04*min(TFMAX,TF)^E04)) is a function that returns the maximum of either (i) FFMIN or (ii) the function that returns the minimum of either FFMAX or the quantity C00+C01*max(CALMIN,CAL)^E01+C02*min(PRMAX PR)^E02+C03*min(DFMAX,DF)^E03+C04*min(TFMAX,TF)^E04,

CAL is total Calories per 100 gram portion of the food,

PR is grams of Protein per 100 g,

DF is grams Dietary Fiber per 100 g, and

TF is grams total Fat per 100 g;

FFMIN, FFMAX, CALMIN, PRMAX, DFMAX, TFMAX, and C00 are constants, C01 through C04 are coefficients, and E01 through E04 are exponential powers having the approximate following values;

FFMIN=0.5

FFMAX=5.0

CALMIN=30

PRMAX=30

DFMAX=12

TFMAX=50

C00=37/60

C01=2500/60

C02=3/60

C03=1/1620

C04=−

1/138,000

E01=−

0.7

E02=1.0

E03=3.0

E04=3.0; and

(ii) a second factor that represents the nutritional benefit of said food according to the relationship;

NDR=max(RMIN,min(RMAX,C00

+C01*(ln(max(PMIN,min(PMAX,100*DF0/DF1)))−

ln(PAVG))

+C02*(ln(max(PMIN,min(PMAX,100*PR0/PR1)))−

ln(PAVG))

+C03*(ln(max(PMIN,min(PMAX,100*VA0/VA1)))−

ln(PAVG))

+C04*(ln(max(PMIN,min(PMAX,100*VC0/VC1)))−

ln(PAVG))

+C05*(ln(max(PMIN,min(PMAX,100*CA0/CA1)))−

ln(PAVG))

+C06*(ln(max(PMIN,min(PMAX,100*FE0/FE1)))−

ln(PAVG))

+C07*(ln(max(PMIN,min(PMAX,100*VD0/VD1)))−

ln(PAVG))

+C08*(ln(max(PMIN,min(PMAX,100*VE0/VE1)))−

ln(PAVG))

+C09*(ln(max(PMIN,min(PMAX,100*VK0/VK1)))−

ln(PAVG))

+C10*(ln(max(PMIN,min(PMAX,100*TH0/TH1)))−

ln(PAVG))

+C11*(ln(max(PMIN,min(PMAX,100*RI0/RI1)))−

ln(PAVG))

+C12*(ln(max(PMIN,min(PMAX,100*NI0/NI1)))−

ln(PAVG))

+C13*(ln(max(PMIN,min(PMAX,100*VB60/VB61)))−

ln(PAVG))

+C14*(ln(max(PMIN,min(PMAX,100*FO00/FO1)))−

ln(PAVG))

+C15*(ln(max(PMIN,min(PMAX,100*VB120/VB121)))−

ln(PAVG))

+C16*(ln(max(PMIN,min(PMAX,100*PA0/PA1)))−

ln(PAVG))

+C17*(ln(max(PMIN,min(PMAX,100*MG0/MG1)))−

ln(PAVG))

+C18*(ln(max(PMIN,min(PMAX,100*P0/P1)))−

ln(PAVG))

+C19*(ln(max(PMIN,min(PMAX,00*K0/K1)))−

ln(PAVG))

+C20*(ln(max(PMIN,min(PMAX,100*ZN0/ZN1)))−

ln(PAVG))

+C21*(ln(max(PMIN,min(PMAX,100*CU0/CU1)))−

ln(PAVG))

+C22*(ln(max(PMIN,min(PMAX,100*MN0/MN1)))−

ln(PAVG))

+C23*(ln(max(PMIN,min(PMAX,100*SE0/SE1)))−

ln(PAVG))

+C24*(ln(max(PMIN,min(PMAX,100*SF0/SF1)))−

ln(PAVG))

+C25*(ln(max(PMIN,min(PMAX,100*CH0/CH1)))−

ln(PAVG))

+C26*(ln(max(PMIN,min(PMAX,100*NA0/NA1)))−

ln(PAVG))

+C27*(ln(max(PMIN,min(PMAX,100*SA0/SA1)))−

ln(PAVG))))

where;

NDR refers to “

nutrient density”

ln is the natural (base e) logarithm,

max(PMIN,min(PMAX,100*DF0/DF1) is a function that returns the maximum of either (i) PMIN or (ii) the function that returns the minimum of either PMAZ or 100*DF0/DF1,

max(PMIN,min(PMAZ,100*PR0/PR1)) is a function that returns the maximum of either (i) PMIN or (ii) the function that returns the minimum of either PMAX or 100*PR0/PR1,

max(PMIN,min(PMAX,100*VA0/VA1)) is a function that returns the maximum of either (i) PMIN or (ii) the function that returns the minimum of either PMAX or 100*VA0/VA1,

max(PMIN min(PMAX,100*VC0/VC1)) is a function that returns the maximum of either (i) PMIN or (ii) the function that returns the minimum of either PMAX or 100*VC0/VC1,

max(PMIN,min(PMAX,100*CA0/CA1)) is a function that returns the maximum of either (i) PMIN or (ii) the function that returns the minimum of either PMAX or 100*CA0/CA1,

max(PMIN,min(PMAX,100*FE0/FE1)) is a function that returns the maximum of either (i) PMIN or (ii) the function that returns the minimum of either PMAX or 100*FE0/FE1,

max(PMIN,min(PMAX,100*VD0/VD1)) is a function that returns the maximum of either (i) PMIN or (ii) the function that returns the minimum of either PMAX or 100*VD0/VD1,

max(PMIN,min(PMAX,100*VE0/VE1)) is a function that returns the maximum of either (i) PMIN or (ii) the function that returns the minimum of either PMAX or 100*VE0/VE1,

max(PMIN,min(PMAX,100*VK0/VK1)) is a function that returns the maximum of either (i) PMIN or (ii) the function that returns the minimum of either PMAX or 100*VK0/VK1,

max(PMIN,min(PMAX,100*TH0/TH1)) is a function that returns the maximum of either (i) PMIN or (ii) the function that returns the minimum of either PMAX or 100*TH0/TH1,

max(PMIN,min(PMAX,100*RI0/RI1)) is a function that returns the maximum of either (i) PMIN or (ii) the function that returns the minimum of either PMAX or 100*RI0/RI1,

max(PMIN,min(PMAX,100*NI0/NI1)) is a function that returns the maximum of either (i) PMIN or (ii) the function that returns the minimum of either PMAX or 100*NI0/NI1,

max(PMIN,min(PMAX,100*VB60/VB61)) is a function that returns the maximum of either (i) PMIN or (ii) the function that returns the minimum of either PMAX or 100*VB60/VB61,

max(PMIN,min(PMAX,100*FO0/FO1)) is a function that returns the maximum of either (i) PMIN or (ii) the function that returns the minimum of either PMAX or 100*FO0/FO1,

max(PMIN,min(PMAX,100*VB120/VB121)) is a function that returns the maximum of either (i) PMIN or (ii) the function that returns the minimum of either PMAX or 100*VB120/VB121,

max(PMIN,min(PMAX,100*PA0/PA1)) is a function that returns the maximum of either (i) PMIN or (ii) the function that returns the minimum of either PMAX or 100*PA0/PA1,

max(PMIN,min(PMAX,100*MG0/MG1)) is a function that returns the maximum of either (i) PMIN or (ii) the function that returns the minimum of either PMAX or 100*MG0/MG1,

max(PMIN,min(PMAX,100*P0/P1)) is a function that returns the maximum of either (i) PMIN or (ii) the function that returns the minimum of either PMAX or 100*P0/P1,

max(PMIN,min(PMAX,100*K0/K1)) is a function that returns the maximum of either (i) PMIN or (ii) the function that returns the minimum of either PMAX or 100*K0/K1,

max(PMIN,min(PMAX,100*ZN0/ZN1)) is a function that returns the maximum of either (i) PMIN or (ii) the function that returns the minimum of either PMAX or 100*ZN0/ZN1,

max(PMIN,min(PMAX,100*CU0/CU1)) is a function that returns the maximum of either (i) PMIN or (ii) the function that returns the minimum of either PMAX or 100*CU0/CU1,

max(PMIN,min(PMAX,100*MN0/MN1)) is a function that returns the maximum of either (i) PMIN or (ii) the function that returns the minimum of either PMAX or 100*MN0/MN1,

max(PMIN,min(PMAX,100*SE0/SE1)) is a function that returns the maximum of either (i) PMIN or (ii) the function that returns the minimum of either PMAX or 100*SE0/SE1,

max(PMIN,min(PMAX,100*SF0/SF1)) is a function that returns the maximum of either (i) PMIN or (ii) the function that returns the minimum of either PMAX or 100*SF0/SF1,

max(PMIN,min(PMAX,100*CH0/CH1)) is a function that returns the maximum of either (i) PMIN or (ii) the function that returns the minimum of either PMAX or 100*CH0/CH1,

max(PMIN,min(PMAX,100*NA0/NA1)) is a function that returns the maximum of either (i) PMIN or (ii) the function that returns the minimum of either PMAX or 100*NA0/NA1,

max(PMIN,min(PMAX,100*SA0/SA1)) is a function that returns the maximum of either (i) PMIN or (ii) the function that returns the minimum of either PMAX or 100*SA0/SA1,

max(RMIN,min(RMAX,C00+C01*(ln(max(PMIN,min(PMAX,100*DF0/DF1)))−

ln(PAVG))+C02*(ln(max(PMIN,min(PMAX,100*PR0/PR1)))−

ln(PAVG))+C03*(ln(max(PMIN,min(PMAX,100*VA0/VA1)))−

ln(PAVG))+C04*(ln(max(PMIN,min(PMAX,100*VC0/VC1)))−

ln(PAVG))+C05*(ln(max(PMIN,min(PMAX,100*CA0/CA1)))−

ln(PAVG))+C06*(ln(max(PMIN,min(PMAX,100*FE0/FE1)))−

ln(PAVG))+C07*(ln(max(PMIN,min(PMAX,100*VD0/VD1)))−

ln(PAVG))+C08*(ln(max(PMIN,min(PMAX,100*VE0/VE1)))−

ln(PAVG))+C09*(ln(max(PMIN,min(PMAX,100*VK0/VK1)))−

ln(PAVG))+C10*(ln(max(PMIN,min(PMAX,100*TH0/TH1)))−

ln(PAVG))+C11*(ln(max(PMIN,min(PMAX,100*RI0/RI1)))−

ln(PAVG)+C12*(ln(max(PMIN,min(PMAX,100*NI0/NI1)))−

ln(PAVG))+C13*(ln(max(PMIN,min(PMAX,100*VB60/VB61)))−

ln(PAVG))+C14*(ln(max(PMIN,min(PMAX,100*FO0/FO1)))−

ln(PAVG))+C15*(ln(max(PMIN,min(PMAX,100*VB120/VB121)))−

ln(PAVG))+C16*(ln(max(PMIN,min(PMAX,100*PA0/PA1)))−

ln(PAVG))+C17*(ln(max(PMIN,min(PMAX,100*MG0/MG1)))−

ln(PAVG))+C18*(ln(max(PMIN,min(PMAX,100*P0/P1)))−

ln(PAVG))+C19*(ln(max(PMIN,min(PMAX,100*K0/K1)))−

ln(PAVG))+C20*(ln(max(PMIN,min(PMAX,100*ZN0/ZN1)))−

ln(PAVG))+C21*(ln(max(PMIN,min(PMAX,100*CU0/CU1)))−

ln(PAVG))+C22*(ln(max(PMIN,min(PMAX,100*MN0/MN1)))−

ln(PAVG))+C23*(ln(max(PMIN,min(PMAX,100*SE0/SE1)))−

ln(PAVG))+C24*(ln(max(PMIN,min(PMAX,100*SF0/SF1)))−

ln(PAVG))+C25*(ln(max(PMIN,min(PMAX,100*CH0/CH1)))−

ln(PAVG))+C26*(ln(max(PMIN,min(PMAX,100*NA0/NA1)))−

ln(PAVG))+C27*(ln(max(PMIN,min(PMAX,100*SA0/SA1)))−

ln(PAVG)))) is a function that returns the maximum of either (i) RMIN or (ii) the function that returns the minimum of either RMAX or the quantity C00+C01*(ln(max(PMIN,min(PMAX,100*DF0/DF1)))−

ln(PAVG))+C02*(ln(max(PMIN,min(PMAX,100*PR0/PR1)))−

ln(PAVG))+C03*(ln(max(PMIN,min(PMAX,100*VA0/VA1)))−

ln(PAVG))+C04*(ln(max(PMIN,min(PMAX,100*VC0/VC1)))−

ln(PAVG))+C05*(ln(max(PMIN,min(PMAX,100*CA0/CA1)))−

ln(PAVG))+C06*(ln(max(PMIN,min(PMAX,100*FE0/FE1)))−

ln(PAVG))+C07*(ln(max(PMIN,min(PMAX,100*VD0/VD1)))−

ln(PAVG))+C08*(ln(max(PMIN,min(PMAX,100*VE0/VE1)))−

ln(PAVG))+C09*(ln(max(PMIN,min(PMAX,100*VK0/VK1)))−

ln(PAVG))+C10*(ln(max(PMIN,min(PMAX,100*TH0/TH1)))−

ln(PAVG))+C11*(ln(max(PMIN,min(PMAX,100*RI0/RI1)))−

ln(PAVG))+C12*(ln(max(PMIN,min(PMAX,100*NI0/NI1)))−

ln(PAVG))+C13*(ln(max(PMIN,min(PMAX,100*VB60/VB61)))−

ln(PAVG))+C14*(ln(max(PMIN,min(PMAX,100*FO0/FO1)))−

ln(PAVG))+C15*(ln(max(PMIN,min(PMAX,100*VB120/VB121)))−

ln(PAVG))+C16*(ln(max(PMIN,min(PMAX,100*PA0/PA1)))−

ln(PAVG))+C17*(ln(max(PMIN,min(PMAX,100*MG0/MG1)))−

ln(PAVG))+C18*(ln(max(PMIN,min(PMAX,100*P0/P1)))−

ln(PAVG))+C19*(ln(max(PMIN,min(PMAX,100*K0/K1)))−

ln(PAVG))+C20*(ln(max(PMIN,min(PMAX,100*ZN0/ZN1)))−

ln(PAVG))+C21*(ln(max(PMIN,min(PMAX,100*CU0/CU1)))−

ln(PAVG))+C22*(ln(max(PMIN,min(PMAX,100*MN0/MN1)))−

ln(PAVG))+C23*(ln(max(PMIN,min(PMAX,100*SE0/SE1)))−

ln(PAVG))+C24*(ln(max(PMIN,min(PMAX,100*SF0/SF1)))−

ln(PAVG))+C25*(ln(max(PMIN,min(PMAX,100*CH0/CH1)))−

ln(PAVG))+C26*(ln(max(PMIN,min(PMAX,100*NA0/NA1)))−

ln(PAVG))+C27*(ln(max(PMIN,min(PMAX,100*SA0/SA1)))−

ln(PAVG))

DF0 is the amount of Dietary Fiber present in a serving,

DF1 is the amount of Dietary Fiber specified as the Daily Value by the US FDA;

PR0 is the amount of Protein present in a serving;

PR1 is the amount of Protein specified as the Daily Value by the US FDA;

VA0 is the amount of Vitamin A present in a serving;

VA1 is the amount of Vitamin A specified as the Daily Value by the US FDA;

VC0 is the amount of Vitamin C present in a serving;

VC1 is the amount of Vitamin C specified as the Daily Value by the US FDA;

CA0 is the amount of Calcium present in a serving;

CA1 is the amount of Calcium specified as the Daily Value by the US FDA;

FE0 is the amount of Iron present in a serving;

FE1 is the amount of Iron specified as the Daily Value by the US FDA;

VD0 is the amount of Vitamin D present in a serving;

VD1 is the amount of Vitamin D specified as the Daily Value by the US FDA;

VE0 is the amount of Vitamin E present in a serving;

VE1 is the amount of Vitamin E specified as the Daily Value by the US FDA;

VK0 is the amount of Vitamin K present in a serving;

VK1 is the amount of Vitamin K specified as the Daily Value by the US FDA;

TH0 is the amount of Thiamin present in a serving;

TH1 is the amount of Thiamin specified as the Daily Value by the US FDA;

RI0 is the amount of Riboflavin present in a serving;

RI1 is the amount of Riboflavin specified as the Daily Value by the US FDA;

NI0 is the amount of Niacin present in a serving;

NI1 is the amount of Niacin specified as the Daily Value by the US FDA;

VB60 is the amount of Vitamin B6 present in a serving;

VB61 is the amount of Vitamin B6 specified as the Daily Value by the US FDA;

FO0 is the amount of Folate present in a serving;

FO1 is the amount of Folate specified as the Daily Value by the US FDA;

VB120 is the amount of Vitamin B12 present in a serving;

VB121 is the amount of Vitamin B12 specified as the Daily Value by the US FDA;

PA0 is the amount of Pantothenic Acid present in a serving;

PA1 is the amount of Pantothenic Acid specified as the Daily Value by the US FDA;

MG0 is the amount of Magnesium present in a serving;

MG1 is the amount of Magnesium specified as the Daily Value by the US FDA;

P0 is the amount of Phosphorus present in a serving;

P1 is the amount of Phosphorus specified as the Daily Value by the US FDA;

K0 is the amount of Potassium present in a serving;

K1 is the amount of Potassium specified as the Daily Value by the US FDA;

ZN0 is the amount of Zinc present in a serving;

ZN1 is the amount of Zinc specified as the Daily Value by the US FDA;

CU0 is the amount of Copper present in a serving;

CU1 is the amount of Copper specified as the Daily Value by the US FDA;

MN0 is the amount of Manganese present in a serving;

MN1 is the amount of Manganese specified as the Daily Value by the US FDA;

SE0 is the amount of Selenium present in a serving;

SE1 is the amount of Selenium specified as the Daily Value by the US FDA;

SF0 is the amount of Saturated Fat present in a serving;

SF1 is the amount of Saturated Fat specified as the Daily Value by the US FDA;

CH0 is the amount of Cholesterol present in a serving;

CH1 is the amount of Cholesterol specified as the Daily Value by the US FDA;

NA0 is the amount of Sodium present in a serving;

NA1 is the amount of Sodium specified as the Daily Value by the US FDA;

SA0 is the number of Calories in a serving from sugars (as reported on Nutrition Facts labels) and alcohol; and

SA1 is the total number of Calories in a serving;

wherein;

RMIN, RMAX, PMIN, PMAX, PAVG, and C00 are constants,

C01 thru C27 are coefficients that weight the effects of the individual nutrients on the NDR; and

wherein these constants and coefficients have the following approximate values;

C00=3.0

C01=0.238

C02=0.238

C03=0.238

C04=0.238

C05=0.238

C06=0.238

C07=0.0238

C08=0.0238

C09=0.0238

C10=0.0238

C11=0.0238

C12=0.0238

C13=0.0238

C14=0.0238

C15=0.0238

C16=0.0238

C17=0.0238

C18=0.0238

C19=0.0238

C20=0.0238

C21=0.0238

C22=0.0238

C23=0.0238

C24=−

0.238

C25=−

0.238

C26=−

0.238

C27=−

0.238; and

(b) output means for providing a guideline to a user, said guideline including said first and second factors to enable food selection in accordance with at least one dietary goal.

10. A computer program product comprising a computer usable physical medium storing a computer executable program to:

(a) compute a first factor that represents an estimate of the satiating effect of a food according to the relationship;

FF=max(FFMIN,min(FFMAX,C00

+C01*max(CALMIN,CAL)^E01

+C02*min(PRMAX,PR)^E02

+C03*min(DFMAX,DF)^E03

+C04*min(TFMAX,TF)^E04))

where;

FF refers to “

fullness factor”

max(CALMIN,CAL) is a function that returns the maximum of either (CALMIN or CAL,

min(PRMAX,PR) is a function that returns the minimum of either PRMAX or PR,

min(DFMAX,DF) is a function that returns the minimum of either DFMAX or DF,

min(TFMAX,TF) is a function that returns the minimum of either TFMAX or TF,

max(FFMIN,min(FFMAX,C00+C01*max(CALMIN,CAL)^E01+C02*min(PRMAX,PR)^E02+C03*min(DFMAX,DF)^E03+C04*min(TFMAX,TF)^E04)) is a function that returns the maximum of either (i) FFMIN or (ii) the function that returns the minimum of either FFMAX or the quantity C00+C01*max(CALMIN,CAL)^E01+C02*min(PRMAX,PR)^E02+C03*min(DFMAX,DF)^E03+C04*min(TFMAX,TF)^E04

CAL is total Calories per 100 gram portion of the food,

PR is grams of Protein per 100 g,

DF is grams Dietary Fiber per 100 g, and

TF is grams total Fat per 100 g;

FFMIN, FFMAX, CALMIN, PRMAX, DFMAX, TFMAX, and C00 are constants, C01 through C04 are coefficients, and E01 through E04 are exponential powers having the approximate following values;

FFMIN=0.5

FFMAX=5.0

CALMIN=30

PRMAX=30

DFMAX=12

TFMAX=50

C00=37/60

C01=2500/60

C02=3/60

C03=1/1620

C04=−

1/138,000

E01=−

0.7

E02=1.0

E03=3.0

E04=3.0; and

(b) compute a second factor that represents the nutritional benefit of said food according to the relationship;

NDR=max(RMIN,min(RMAX,C00

+C01*(ln(max(PMIN,min(PMAX,100*DF0/DF1)))−

ln(PAVG))

+C02*(ln(max(PMIN,min(PMAX,100*PR0/PR1)))−

ln(PAVG))

+C03*(ln(max(PMIN,min(PMAX,100*VA0/VA1)))−

ln(PAVG))

+C04*(ln(max(PMIN,min(PMAX,100*VC0/VC1)))−

ln(PAVG))

+C05*(ln(max(PMIN,min(PMAX,100*CA0/CA1)))−

ln(PAVG))

+C06*(ln(max(PMIN,min(PMAX,100*FE0/FE1)))−

ln(PAVG))

+C07*(ln(max(PMIN,min(PMAX,100*VD0/VD1)))−

ln(PAVG))

+C08*(ln(max(PMIN,min(PMAX,100*VE0/VE1)))−

ln(PAVG))

+C09*(ln(max(PMIN,min(PMAX,100*VK0/VK1)))−

ln(PAVG))

+C10*(ln(max(PMIN,min(PMAX,100*TH0/TH1)))−

ln(PAVG))

+C11*(ln(max(PMIN,min(PMAX,100*RI0/RI1)))−

ln(PAVG))

+C12*(ln(max(PMIN,min(PMAX,100*NI0/NI1)))−

ln(PAVG))

+C13*(ln(max(PMIN,min(PMAX,100*VB60/VB61)))−

ln(PAVG))

+C14*(ln(max(PMIN,min(PMAX,100*FO0/FO1)))−

ln(PAVG))

+C15*(ln(max(PMIN,min(PMAX,100*VB120/VB121)))−

ln(PAVG))

+C16*(ln(max(PMIN,min(PMAX,100*PA0/PA1)))−

ln(PAVG))

+C17*(ln(max(PMIN,min(PMAX,100*MG0/MG1)))−

ln(PAVG))

+C18*(ln(max(PMIN,min(PMAX,100*P0/P1)))−

ln(PAVG))

+C19*(ln(max(PMIN,min(PMAX,100*K0/K1)))−

ln(PAVG))

+C20*(ln(max(PMIN,min(PMAX,100*ZN0/ZN1)))−

ln(PAVG))

+C21*(ln(max(PMIN,min(PMAX,100*CU0/CU1)))−

ln(PAVG))

+C22*(ln(max(PMIN,min(PMAX,100*MN0/MN1)))−

ln(PAVG))

+C23*(ln(max(PMIN,min(PMAX,100*SE0/SE1)))−

ln(PAVG))

+C24*(ln(max(PMIN,min(PMAX,100*SF0/SF1)))−

ln(PAVG))

+C25*(ln(max(PMIN,min(PMAX,100*CH0/CH1)))−

ln(PAVG))

+C26*(ln(max(PMIN,min(PMAX,100*NA0/NA1)))−

ln(PAVG))

+C27*(ln(max(PMIN,min(PMAX,100*SA0/SA1)))−

ln(PAVG))))

where;

NDR refers to “

nutrient density”

ln is the natural (base e) logarithm,

max(PMIN,min(PMAX,100*DF0/DF1) is a function that returns the maximum of either (i) PMIN or (ii) the function that returns the minimum of either PMAX or 100*DF0/DF1,

max(PMIN,min(PMAX,100*PR0/PR1)) is a function that returns the maximum of either (i) PMIN or (ii) the function that returns the minimum of either PMAX or 100*PR0/PR1,

max(PMIN,min(PMAX,100*VA0/VA1)) is a function that returns the maximum of either (i) PMIN or (ii) the function that returns the minimum of either PMAX or 100*VA0/VA1,

max(PMIN min(PMAX,100*VC0/VC1)) is a function that returns the maximum of either (i) PMIN or (ii) the function that returns the minimum of either PMAX or 100*VC0/VC1,

max(PMIN,min(PMAX,100*CA0/CA

     1)) is a function that returns the maximum of either (i) PMIN or (ii) the function that returns the minimum of either PMAX or 100*CA0/CA1,

max(PMIN,min(PMAX,100*FE0/FE1)) is a function that returns the maximum of either (i) PMIN or (ii) the function that returns the minimum of either PMAX or 100*FE0/FE1,

max(PMIN,min(PMAX,100*VD0/VD1)) is a function that returns the maximum of either (i) PMIN or (ii) the function that returns the minimum of either PMAX or 100*VD0/VD1,

max(PMIN,min(PMAX,100*VE0V/E1)) is a function that returns the maximum of either (i) PMIN or (ii) the function that returns the minimum of either PMAX or 100*VE0/VE1,

max(PMIN,min(PMAX,100*VK0/VK1)) is a function that returns the maximum of either (i) PMIN or (ii) the function that returns the minimum of either PMAX or 100*VK0/VK1,

max(PMIN,min(PMAX,100*TH0/TH1)) is a function that returns the maximum of either (i) PMIN or (ii) the function that returns the minimum of either PMAX or 100*TH0/TH1,

max(PMIN,min(PMAX,100*RI0/RI1)) is a function that returns the maximum of either (i) PMIN or (ii) the function that returns the minimum of either PMAX or 100*RI0/RI1,

max(PMIN,min(PMAX,100*NI0/NI1)) is a function that returns the maximum of either (i) PMIN or (ii) the function that returns the minimum of either PMAX or 100*NI0/NI1,

max(PMIN,min(PMAX,100*VB60/VB61)) is a function that returns the maximum of either (i) PMIN or (ii) the function that returns the minimum of either PMAX or 100*VB60/VB61,

max(PMIN,min(PMAX,100*FO0/FO1)) is a function that returns the maximum of either (i) PMIN or (ii) the function that returns the minimum of either PMAX or 100*FO0/FO1,

max(PMIN,min(PMAX,100*VB120/VB121)) is a function that returns the maximum of either (i) PMIN or (ii) the function that returns the minimum of either PMAX or 100*VB120/VB121,

max(PMIN,min(PMAX,100*PA0/PA1)) is a function that returns the maximum of either (i) PMIN or (ii) the function that returns the minimum of either PMAX or 100*PA0/PA1,

max(PMIN,min(PMAX,100*MG0/MG1)) is a function that returns the maximum of either (i) PMIN or (ii) the function that returns the minimum of either PMAX or 100*MG0/MG1,

max(PMIN,min(PMAX,100*P0/P1)) is a function that returns the maximum of either (i) PMIN or (ii) the function that returns the minimum of either PMAX or 100*P0/P1,

max(PMIN,min(PMAX,100*K0/K1)) is a function that returns the maximum of either (i) PMIN or (ii) the function that returns the minimum of either PMAX or 100*K0/K1,

max(PMIN,min(PMAX,100*ZN0/ZN1)) is a function that returns the maximum of either (i) PMIN or (ii) the function that returns the minimum of either PMAX or 100*ZN0/ZN1,

max(PMIN,min(PMAX,100*CU0/CU1)) is a function that returns the maximum of either (i) PMIN or (ii) the function that returns the minimum of either PMAX or 100*CU0/CU1,

max(PMIN,min(PMAX,100*MN0/MN1)) is a function that returns the maximum of either (i) PMIN or (ii) the function that returns the minimum of either PMAX or 100*MN0/MN1,

max(PMIN,min(PMAX,100*SE0/SE1)) is a function that returns the maximum of either (i) PMIN or (ii) the function that returns the minimum of either PMAX or 100*SE0/SE1,

max(PMIN,min(PMAX,100*SF0/SF1)) is a function that returns the maximum of either (i) PMIN or (ii) the function that returns the minimum of either PMAX or 100*SF0/SF1,

max(PMIN,min(PMAX,100*CH0/CH1)) is a function that returns the maximum of either (i) PMIN or (ii) the function that returns the minimum of either PMAX or 100*CH0/CH1,

max(PMIN,min(PMAX,100*NA0/NA1)) is a function that returns the maximum of either (i) PMIN or (ii) the function that returns the minimum of either PMAX or 100*NA0/NA1,

max(PMIN,min(PMAX,100*SA0/SA1)) is a function that returns the maximum of either (i) PMIN or (ii) the function that returns the minimum of either PMAX or 100*SA0/SA1,

max(RMIN,min(RMAX,C00+C01*(ln(max(PMIN,min(PMAX,100*DF0/DF1)))−

ln(PAVG))+C02*(ln(max(PMIN,min(PMAX,100*PR0/PR1)))−

ln(PAVG))+C03*(ln(max(PMIN,min(PMAX,100*VA0/VA1)))−

ln(PAVG))+C04*(ln(max(PMIN,min(PMAX,100*VC0/VC1)))−

ln(PAVG))+C05*(ln(max(PMIN,min(PMAX,100*CA0/CA1)))−

ln(PAVG))+C06*(ln(max(PMIN,min(PMAX,100*FE0/FE1)))−

ln(PAVG))+C07*(ln(max(PMIN,min(PMAX,100*VD0/VD1)))−

ln(PAVG))+C08*(ln(max(PMIN,min(PMAX,100*VE0/VE1)))−

ln(PAVG))+C09*(ln(max(PMIN,min(PMAX,100*VK0/VK1)))−

ln(PAVG))+C10*(ln(max(PMIN,min(PMAX,100*TH0/TH1)))−

ln(PAVG))+C11*(ln(max(PMIN,min(PMAX,100*RI0/RI1)))−

ln(PAVG)+C12*(ln(max(PMIN,min(PMAX,100*NI0/NI1)))−

ln(PAVG))+C13*(ln(max(PMIN,min(PMAX,100*VB60/VB61)))−

ln(PAVG))+C14*(ln(max(PMIN,min(PMAX,100*FO0/FO1)))−

ln(PAVG))+C15*(ln(max(PMIN,min(PMAX,100*VB120/VB121)))−

ln(PAVG))+C16*(ln(max(PMIN,min(PMAX,100*PA0/PA1)))−

ln(PAVG))+C17*(ln(max(PMIN,min(PMAX,100*MG0/MG1)))−

ln(PAVG))+C18*(ln(max(PMIN,min(PMAX,100*P0/P1)))−

ln(PAVG))+C19*(ln(max(PMIN,min(PMAX,100*K0/K1)))−

ln(PAVG))+C20*(ln(max(PMIN,min(PMAX,100*ZN0/ZN1)))−

ln(PAVG))+C21*(ln(max(PMIN,min(PMAX,100*CU0/CU1)))−

ln(PAVG))+C22*(ln(max(PMIN,min(PMAX,100*MN0/MN1)))−

ln(PAVG))+C23*(ln(max(PMIN,min(PMAX,100*SE0/SE1)))−

ln(PAVG))+C24*(ln(max(PMIN,min(PMAX,100*SF0/SF1)))−

ln(PAVG))+C25*(ln(max(PMIN,min(PMAX,100*CH0/CH1)))−

ln(PAVG))+C26*(ln(max(PMIN,min(PMAX,100*NA0/NA1)))−

ln(PAVG))+C27*(ln(max(PMIN,min(PMAX,100*SA0/SA1)))−

ln(PAVG)))) is a function that returns the maximum of either (i) RMIN or (ii) the function that returns the minimum of either RMAX or the quantity C00+C01*(ln(max(PMIN,min(PMAX,100*DF0/DF1)))−

ln(PAVG))+C02*(ln(max(PMIN,min(PMAX,100*PR0/PR1)))−

ln(PAVG))+C03*(ln(max(PMIN,min(PMAX,100*VA0/VA1)))−

ln(PAVG))+C04*(ln(max(PMIN,min(PMAX,100*VC0/VC1)))−

ln(PAVG))+C05*(ln(max(PMIN,min(PMAX,100*CA0/CA1)))−

ln(PAVG))+C06*(ln(max(PMIN,min(PMAX,100*FE0/FE1)))−

ln(PAVG))+C07*(ln(max(PMIN,min(PMAX,100*VD0/VD1)))−

ln(PAVG))+C08*(ln(max(PMIN,min(PMAX,100*VE0/VE1)))−

ln(PAVG))+C09*(ln(max(PMIN,min(PMAX,100*VK0/VK1)))−

ln(PAVG)+C10*(ln(max(PMIN,min(PMAX,100*TH0/TH1)))−

ln(PAVG))+C11*(ln(max(PMIN,min(PMAX,100*RI0/RI1)))−

ln(PAVG))+C12*(ln(max(PMIN,min(PMAX,100*NI0/NI1)))−

ln(PAVG))+C13*(ln(max(PMIN,min(PMAX,100*VB60/VB61)))−

ln(PAVG))+C14*(ln(max(PMIN,min(PMAX,100*FO0/FO1)))−

ln(PAVG))+C15*(ln(max(PMIN,min(PMAX,100*VB120/VB121)))−

ln(PAVG))+C16*(ln(max(PMIN,min(PMAX,100*PA0/PA1)))−

ln(PAVG))+C17*(ln(max(PMIN,min(PMAX,100*MG0/MG1)))−

ln(PAVG))+C18*(ln(max(PMIN,min(PMAX,100*P0/P1)))−

ln(PAVG))+C19*(ln(max(PMIN,min(PMAX,100*K0/K1)))−

ln(PAVG))+C20*(ln(max(PMIN,min(PMAX,100*ZN0/ZN1)))−

ln(PAVG))+C21*(ln(max(PMIN,min(PMAX,100*CU0/CU1)))−

ln(PAVG))+C22*(ln(max(PMIN,min(PMAX,100*MN0/MN1)))−

ln(PAVG))+C23*(ln(max(PMIN,min(PMAX,100*SE0/SE1)))−

ln(PAVG))+C24*(ln(max(PMIN,min(PMAX,100*SF0/SF1)))−

ln(PAVG))+C25*(ln(max(PMIN,min(PMAX,100*CH0/CH1)))−

ln(PAVG))+C26*(ln(max(PMIN,min(PMAX,100*NA0/NA1)))−

ln(PAVG))+C27*(ln(max(PMIN,min(PMAX,100*SA0/SA1)))−

ln(PAVG)),

DF0 is the amount of Dietary Fiber present in a serving,

DF1 is the amount of Dietary Fiber specified as the Daily Value by the US FDA;

PR0 is the amount of Protein present in a serving;

PR1 is the amount of Protein specified as the Daily Value by the US FDA;

VA0 is the amount of Vitamin A present in a serving;

VA1 is the amount of Vitamin A specified as the Daily Value by the US FDA;

VC0 is the amount of Vitamin C present in a serving;

VC1 is the amount of Vitamin C specified as the Daily Value by the US FDA;

CA0 is the amount of Calcium present in a serving;

CA1 is the amount of Calcium specified as the Daily Value by the US FDA;

FE0 is the amount of Iron present in a serving;

FE1 is the amount of Iron specified as the Daily Value by the US FDA;

VD0 is the amount of Vitamin D present in a serving;

VD1 is the amount of Vitamin D specified as the Daily Value by the US FDA;

VE0 is the amount of Vitamin E present in a serving;

VE1 is the amount of Vitamin E specified as the Daily Value by the US FDA;

VK0 is the amount of Vitamin K present in a serving;

VK1 is the amount of Vitamin K specified as the Daily Value by the US FDA;

TH0 is the amount of Thiamin present in a serving;

TH1 is the amount of Thiamin specified as the Daily Value by the US FDA;

RI0 is the amount of Riboflavin present in a serving;

RI1 is the amount of Riboflavin specified as the Daily Value by the US FDA;

NI0 is the amount of Niacin present in a serving;

NI1 is the amount of Niacin specified as the Daily Value by the US FDA;

VB60 is the amount of Vitamin B6 present in a serving;

VB61 is the amount of Vitamin B6 specified as the Daily Value by the US FDA;

FO0 is the amount of Folate present in a serving;

FO1 is the amount of Folate specified as the Daily Value by the US FDA;

VB120 is the amount of Vitamin B12 present in a serving;

VB121 is the amount of Vitamin B12 specified as the Daily Value by the US FDA;

PA0 is the amount of Pantothenic Acid present in a serving;

PA1 is the amount of Pantothenic Acid specified as the Daily Value by the US FDA;

MG0 is the amount of Magnesium present in a serving;

MG1 is the amount of Magnesium specified as the Daily Value by the US FDA;

P0 is the amount of Phosphorus present in a serving;

P1 is the amount of Phosphorus specified as the Daily Value by the US FDA;

K0 is the amount of Potassium present in a serving;

K1 is the amount of Potassium specified as the Daily Value by the US FDA;

ZN0 is the amount of Zinc present in a serving;

ZN1 is the amount of Zinc specified as the Daily Value by the US FDA;

CU0 is the amount of Copper present in a serving;

CU1 is the amount of Copper specified as the Daily Value by the US FDA;

MN0 is the amount of Manganese present in a serving;

MN1 is the amount of Manganese specified as the Daily Value by the US FDA;

SE0 is the amount of Selenium present in a serving;

SE1 is the amount of Selenium specified as the Daily Value by the US FDA;

SF0 is the amount of Saturated Fat present in a serving;

SF1 is the amount of Saturated Fat specified as the Daily Value by the US FDA;

CH0 is the amount of Cholesterol present in a serving;

CH1 is the amount of Cholesterol specified as the Daily Value by the US FDA;

NA0 is the amount of Sodium present in a serving;

NA1 is the amount of Sodium specified as the Daily Value by the US FDA;

SA0 is the number of Calories in a serving from sugars (as reported on Nutrition Facts labels) and alcohol; and

SA1 is the total number of Calories in a serving;

wherein;

RMIN, RMAX, PMIN, PMAX, PAVG, and C00 are constants,

C01 thru C27 are coefficients that weight the effects of the individual nutrients on the NDR; and

wherein these constants and coefficients have the following approximate values;

RMIN=0.0

RMAX=5.0

PMIN=2.5

PMAX=40

PAVG=10

C00=3.0

C01=0.238

C02=0.238

C03=0.238

C04=0.238

C05=0.238

C06=0.238

C07=0.0238

C08=0.0238

C09=0.0238

C10=0.0238

C11=0.0238

C12=0.0238

C13=0.0238

C14=0.0238

C15=0.0238

C16=0.0238

C17=0.0238

C18=1.0238

C19=0.0238

C20=0.0238

C21=0.0238

C22=0.238

C23=0.0238

C24=−

0.238

C25=−

0.238

C26=−

0.238

C27=−

0.238; and

(c) provide a guideline to a user including said first and second factors to enable food selection in accordance with at least one dietary goal.

11. A method of rating food, comprising the steps of:

(a) calculating a first factor that represents an estimate of the satiating effect of a food;

(b) calculating a second factor that represents the nutritional benefit of said food;

(c) combining said first and second factors to yield a single indicator of suitability of said food for at least one dietary goal, said single indicator being representative of a relationship between (i) said first and second factors and (ii) satiety and nutritional values which represent fulfillment of said at least one dietary goal; and

(d) providing to a user a guideline including said single indicator to enable food selection in accordance with said at least one dietary goal;

wherein at least one of steps (a), (b) and (c) is effected using a suitably programmed computer.

12. The method of claim 11, wherein said step of providing to a user a guideline includes providing a graphical representation including regions associated with different ones of said at least one dietary goal.

13. The method of claim 11, wherein said guideline includes regions based on nutritional characteristics for preselected nutritional goals.

14. The method of claim 11, further comprising the step of plotting said first and second factors as a first graphical coordinate pair and said satiety and nutritional values which represent fulfillment of said at least one dietary goal as a second graphical coordinate pair, and wherein said relationship describes a difference between said first graphical coordinate pair and said second graphical coordinate pair.

15. The method of claim 11 wherein said step of calculating a first factor is based on the calorie, protein, dietary fiber and total fat content of said food and said step of calculating a second factor is based on content of nutrients, said content of nutrients including said protein content, selected vitamin and selected mineral content, and Daily Values for said nutrients.

16. The method of claim 15, wherein said Daily Values are specified by the United States Food and Drug Administration.

17. A food rating system, comprising:

a data processor adapted to (i) calculate a first factor that represents an estimate of the satiating effect of a food and a second factor that represents the nutritional benefit of said food and (ii) combine said first and second factors into a single indicator representative of a relationship between (a) said first and second factors and (b) satiety and nutritional values which represent fulfillment of at least one dietary goal; and

output means for providing a guideline to a user, said guideline including said single indicator to enable food selection in accordance with said at least one dietary goal.

18. The system of claim 17, wherein said first factor is a numerical expression of caloric density of said food.

19. The system of claim 17, wherein said second factor is a numerical expression of nutrient density of said food.

20. The system of claim 17, wherein said data processor is adapted to plot said first and second factors as a first graphical coordinate pair and said satiety and nutritional values which represent fulfillment of said at least one dietary goal as a second graphical coordinate pair, and wherein said relationship describes a difference between said first graphical coordinate pair and said second graphical coordinate pair.

21. The system of claim 17, wherein said guideline is divided into regions based on nutritional characteristics of said food.

22. The system of claim 17, wherein said guideline is a graphical representation including regions associated with different ones of said at least one dietary goal.

23. The system of claim 22, wherein said single indicator corresponds to one of said regions.

24. The system of claim 17, wherein said first factor is calculated based on calorie, protein, dietary fiber and total fat content of said food;

and said second factor is calculated based on content of nutrients, said content of nutrients including said dietary fiber and said protein content, and selected vitamin and selected mineral content, and Daily Values for said nutrients.

25. The system of claim 24, wherein said Daily Values are specified by the United States Food and Drug Administration.

26. A computer program product comprising a computer usable physical medium storing a computer executable program to:

(a) compute a first factor that represents an estimate of the satiating effect of a food;

(b) compute a second factor that represents the nutritional benefit of said food;

(c) combine said first and second factors into a single indicator of suitability of said food for at least one dietary goal, said single indicator being representative of a relationship between (i) said first and second factors and (ii) satiety and nutritional values which represent fulfillment of said at least one dietary goal; and

(d) provide a guideline including said single indicator to a user to enable food selection in accordance with said at least one dietary goal.

27. The computer program product of claim 26, wherein said first factor is a numerical expression of caloric density of said food.

28. The computer program product of claim 26, wherein said second factor is a numerical expression of nutrient density of said food.

29. The computer program product of claim 26, wherein said guideline is a graphical representation including regions associated with different ones of said at least one dietary goal.

30. The computer program product of claim 29, wherein said single indicator corresponds to one of said regions.

31. The computer program product of claim 26, wherein said first and second factors are plotted as a first graphical coordinate pair and said satiety and nutritional values which represent fulfillment of said at least one dietary goal are plotted as a second graphical coordinate pair, and wherein said relationship describes a difference between said first graphical coordinate pair and said second graphical coordinate pair.

32. The computer program product of claim 26, wherein said guideline is divided into regions based on nutritional characteristics of said food.

33. The computer program product of claim 26, wherein said first factor is computed based on calorie, protein, dietary fiber and total fat content of said food;

and said second factor is computed based on content of nutrients, said content of nutrients including said dietary fiber and said protein content, and selected vitamin and selected mineral content, and Daily Values for said nutrients.

34. The computer program product of claim 33, wherein said Daily Values are specified by the United States Food and Drug Administration.