MedDietCalc (version 0.1.1)

computeMDS12: computeMDS12

Description

Computes a 2012 update of the widely used Mediterranean Diet Score.

Usage

computeMDS12(data, Vegetables, Legumes, FruitAndNuts, Cereals, Potatoes = NULL,
             Fish, Dairy, Meat, Alcohol,
             OOprincipal, Sex, men = "male", women = "female",
             frequency = NULL, output = "percent", rm.na = FALSE)

Arguments

data

Your data set with nutritional information about food or nutrient consumption. Each row is expected to be a person, and food or nutrient intake are in columns

Vegetables

Numeric variable with Vegetables consumption as grams

Legumes

Numeric variable with Legumes consumption as grams

FruitAndNuts

Numeric variable with consumption of Fruits and Nuts as grams

Cereals

Numeric variable with Cereals consumption as grams

Potatoes

Numeric variable with Potatoes consumption as grams

Fish

Numeric variable with Fish consumption as grams

Dairy

Numeric variable with Dairy consumption as grams

Meat

Numeric variable with Meat consumption as grams

Alcohol

Numeric variable with Alcohol consumption as etanol grams from any beberage origin

OOprincipal

Integer. This item scores wether olive oil is the main dietary fat as a dichotomous variable (1-yes, 0-no).

Sex

Vector with gender, it can be numeric, factor or character, as long as its values are provided by 'men' and 'women' arguments.

men

A character string with the value of male gender, default is "male"

women

A character string with the value of female gender, default is "female"

frequency

A character string. Allowed values are 'daily', 'weekly' and 'monthly'. It informs about the frequency which food or nutrient consumption refers to (i.e. wether the rest of arguments are 'grams per day' or 'grams per week' or 'grams per month')

output

A character string to set which output should the formula give, allowed values are 'data.frame', 'score' and 'percent' (default).

rm.na

Logical. If set to FALSE (default), a diet score will be computed only if a person has all score components informed. If set to TRUE, NA values in score components will be drop off and a value of available components will be returned, but percent of score adherence will be computed with basis of the whole score range (see Details)

Value

Computed MDS score according to 2012 version. Depending on 'output' argument, value can be a data.frame, or a vector:

if output = 'data.frame'

A data frame with a row corresponding to each person in data. Columns are the score of each component, as well as the global score as natural sum ('absolute' column) and as percentage ('percent' column)

if output = 'score'

Instead of the full data.frame, just the integer vector corresponding to the absolute points of adherence to Mediterranean Diet for each person, from 0 (min.) to 18 (max.)

if output = 'percent'

Instead of the full data.frame, just the numeric vector corresponding to the percent of adherence to Mediterranean Diet for each person, from 0 (min. adherence) to 100 percent (max. adhrence)

Details

This score is an update of the widely used MDS (Mediterranean Diet Score), with some modifications, the most relevant are the following: First, it uses fixed ('a priori') cut-offs, instead of using sample derived medians. Second, instead of scoring all variables dichotomously (0-1), it scores from 0 (minimum) to 2 (maximum), with items which can receive 1 point. As another difference, it stops evaluating Mono and Poli Unsaturated fats, but instead scores Olive Oil consumption. Olive Oil is considered dichotomously.

References

Sofi, Francesco, Rosanna Abbate, Gian Franco Gensini, Alessandro Casini, Antonia Trichopoulou, and Christina Bamia. 2012. <U+2018>Identification of Change-Points in the Relationship between Food Groups in the Mediterranean Diet and Overall Mortality: An <U+201C>a Posteriori<U+201D> Approach<U+2019>. European Journal of Nutrition 51 (2): 167<U+2013>72. doi:10.1007/s00394-011-0202-7.

See Also

computeMDS95 computeMDS03 computeMDS05

Examples

# NOT RUN {
data(nutriSample)
MedDiet <- computeMDS12(data = nutriSample,
        Vegetables = P41grCom + P42grCom,
        Legumes =  P46grCom,
        FruitAndNuts =  P50grCom + P52grCom + P53grCom,
        Cereals = P55grCom + P56grCom + P57grCom + P59grCom +
                  P60grCom + P61grCom + P62grCom,
        Fish = P35grCom + P36grCom + P37grCom + P38grCom,
        Meat = P29grCom + P30grCom + P31grCom + P32grCom,
        Dairy = P19grCom + P20grCom + P22grCom + P23grCom +
                P24grCom + P25grCom + P26grCom + P27grCom,
        Alcohol =  12 * (P94rac + P96rac + P97rac + P98rac + P99rac),
        Potatoes =  NULL,
        OOprincipal = ifelse(nutriSample$AceiteTipo == 1, 1, 0),
        Sex =  SEXO, men = "Hombre", women = "Mujer", frequency = "daily",
        output = "percent", rm.na = FALSE)
hist(MedDiet)
# }