2-food-optimization-nl.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "%run imports.py"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Preprocessing"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 105,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "# Read in the food nutrients and recommended daily allowance DataFrames.\n",
    "nutrients_df = pd.read_msgpack('nutrients-nl.msg')\n",
    "rda_df = pd.read_msgpack('rda.msg')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 106,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "# Convert to numeric.\n",
    "nutrients_df = nutrients_df \\\n",
    "    .replace(r'([0-9.]+)[^0-9]+', r'\\1', regex=True) \\\n",
    "    .replace('--', np.nan) \\\n",
    "    .replace('spoor', 0)\n",
    "nutrients_df.iloc[:, 3:] = nutrients_df.iloc[:, 3:].astype(np.float64)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 107,
   "metadata": {
    "collapsed": true,
    "scrolled": false
   },
   "outputs": [],
   "source": [
    "# Convert all measurements to grams.\n",
    "def to_grams(df):\n",
    "    def col_to_grams(col):\n",
    "        si = re.findall(r'\\(([^A]?)(?:g|J|L|cal|IU|%)', col.name)\n",
    "        if len(si) != 1:\n",
    "            raise ValueError(f'Bad column format: {si} ({col.name})')\n",
    "        si = si[0]\n",
    "        if not si:\n",
    "            return col\n",
    "        elif si == 'k':\n",
    "            return 1e3 * col\n",
    "        elif si == 'm':\n",
    "            return 1e-3 * col\n",
    "        elif si in ['µ', 'μ', 'u']:  # These two mus are different chars.\n",
    "            return 1e-6 * col\n",
    "        else:\n",
    "            raise ValueError(f'Unknown unit: {si} ({col.name})')\n",
    "    for col in df.select_dtypes(exclude=['object']).columns:\n",
    "        df[col] = col_to_grams(df[col])\n",
    "    return df\n",
    "\n",
    "nutrients_df = to_grams(nutrients_df)\n",
    "rda_df = to_grams(rda_df)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 108,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "# Remove RDA columns that are not measured in the food nutrients DataFrame.\n",
    "verified_missing = [\n",
    "    # Elements\n",
    "    'Chromium(μg/d)',\n",
    "    'Fluoride(mg/d)',\n",
    "    'Manganese(mg/d)',\n",
    "    'Molybdenum(μg/d)',\n",
    "    'Chloride(g/d)',\n",
    "    # Vitamins\n",
    "    'Pantothenic Acid(mg/d)',\n",
    "    'Biotin(μg/d)',\n",
    "    'Choline(mg/d)g',\n",
    "    # Macronutrients\n",
    "    'Total Watera(L/d)'\n",
    "]\n",
    "rda_df = rda_df.drop(verified_missing, axis=1)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 109,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "# Rename RDA columns to match the food nutrients column names.\n",
    "exceptions = {\n",
    "    # Elements\n",
    "    'Copper(μg/d)': 'Koper (mg)',\n",
    "    'Iodine(μg/d)': 'Jodium (ug)',\n",
    "    'Iron(mg/d)': 'IJzer totaal (mg)',\n",
    "    'Phosphorus(mg/d)': 'Fosfor (mg)',\n",
    "    'Zinc(mg/d)': 'Zink (mg)',\n",
    "    'Potassium(g/d)': 'Kalium (mg)',\n",
    "    'Sodium(g/d)': 'Natrium (mg)',\n",
    "    # Vitamins\n",
    "    'Vitamin A(μg/d)a': 'Retinol activiteit equiv.(RAE) (ug)',\n",
    "    'Vitamin C(mg/d)': 'Vitamine C (mg)',\n",
    "    'Vitamin D(μg/d)b,c': 'Cholecalciferol - vit D3 (ug)',  # 'Vitamine D totaal (ug)'\n",
    "    'Vitamin E(mg/d)d': 'Alfa-tocoferol (mg)',  # 'Vitamine E totaal (mg)',\n",
    "    'Vitamin K(μg/d)': 'Vitamine K totaal (ug)',\n",
    "    'Thiamin(mg/d)': 'Vitamine B1 (mg)',\n",
    "    'Riboflavin(mg/d)': 'Vitamine B2 (mg)',\n",
    "    'Niacin(mg/d)e': 'Nicotinezuur (mg)',\n",
    "    'Vitamin B6(mg/d)': 'Vitamine B6 (mg)',\n",
    "    'Folate(μg/d)f': 'Foliumzuur (ug)',\n",
    "    'Vitamin B12(μg/d)': 'Vitamine B12 (ug)',\n",
    "    # Macronutrients\n",
    "    'Carbohydrate(g/d)': 'Koolhydraten totaal (g)',\n",
    "    'Total Fiber(g/d)': 'Voedingsvezel totaal (g)',\n",
    "    'Fat(g/d)': 'Vet totaal (g)',\n",
    "    'Linoleic Acid(g/d)': 'C18:2 n-6 cis (linolzuur) (g)',\n",
    "    'α-Linolenic Acid(g/d)': 'C18:3 n-3 cis (ALA) (g)',\n",
    "    'Proteinb(g/d)': 'Eiwit totaal (g)'\n",
    "}\n",
    "columns = {**exceptions}\n",
    "for col in rda_df.columns:\n",
    "    if col in exceptions:\n",
    "        continue\n",
    "    query = re.findall(r'^([a-z\\s]+)[^a-z]', col.lower())[0]\n",
    "    results = [c for c in nutrients_df.columns if query in c.lower()]\n",
    "    if len(results) != 1:\n",
    "        raise ValueError(f'Found {results} for \"{query}\" in column {col}')\n",
    "    columns[col] = results[0]\n",
    "rda_df = rda_df.rename(columns=columns)\n",
    "assert all([col in nutrients_df.columns for col in rda_df.columns])"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# MILP model that optimizes health, time, cost\n",
    "\n",
    "Possible losses:\n",
    "- Time (Preparation): Waste\n",
    "- Time (Shopping): Number of items\n",
    "- Time (Eating): Meal weight\n",
    "- Health: Recommended Dietary Allowances (Adequate Intake)\n",
    "- Money: Cost"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 234,
   "metadata": {
    "scrolled": false
   },
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style>\n",
       "    .dataframe thead tr:only-child th {\n",
       "        text-align: right;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: left;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>Voedselgroep</th>\n",
       "      <th>Voedingsmiddel</th>\n",
       "      <th>Calcium (mg)</th>\n",
       "      <th>Koper (mg)</th>\n",
       "      <th>Jodium (ug)</th>\n",
       "      <th>IJzer totaal (mg)</th>\n",
       "      <th>Fosfor (mg)</th>\n",
       "      <th>Zink (mg)</th>\n",
       "      <th>Kalium (mg)</th>\n",
       "      <th>Natrium (mg)</th>\n",
       "      <th>Retinol activiteit equiv.(RAE) (ug)</th>\n",
       "      <th>Vitamine C (mg)</th>\n",
       "      <th>Cholecalciferol - vit D3 (ug)</th>\n",
       "      <th>Alfa-tocoferol (mg)</th>\n",
       "      <th>Vitamine K totaal (ug)</th>\n",
       "      <th>Vitamine B1 (mg)</th>\n",
       "      <th>Vitamine B2 (mg)</th>\n",
       "      <th>Nicotinezuur (mg)</th>\n",
       "      <th>Vitamine B6 (mg)</th>\n",
       "      <th>Vitamine B12 (ug)</th>\n",
       "      <th>Koolhydraten totaal (g)</th>\n",
       "      <th>Voedingsvezel totaal (g)</th>\n",
       "      <th>C18:2 n-6 cis (linolzuur) (g)</th>\n",
       "      <th>C18:3 n-3 cis (ALA) (g)</th>\n",
       "      <th>Eiwit totaal (g)</th>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>NEVO-code</th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>2783</th>\n",
       "      <td>Brood</td>\n",
       "      <td>Brood meergranen- m div zaden obv bloem</td>\n",
       "      <td>0.046</td>\n",
       "      <td>26.666667</td>\n",
       "      <td>4.706667</td>\n",
       "      <td>2.1250</td>\n",
       "      <td>0.225714</td>\n",
       "      <td>10.727273</td>\n",
       "      <td>0.040851</td>\n",
       "      <td>0.296000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.066667</td>\n",
       "      <td>0.000</td>\n",
       "      <td>9.166667</td>\n",
       "      <td>5.384615</td>\n",
       "      <td>1.0000</td>\n",
       "      <td>72.307692</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.310769</td>\n",
       "      <td>1.421053</td>\n",
       "      <td>1.058824</td>\n",
       "      <td>2.500</td>\n",
       "      <td>0.223214</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2423</th>\n",
       "      <td>Vetten, oliën en hartige sauzen</td>\n",
       "      <td>Halvarineproduct Blue Band Idee</td>\n",
       "      <td>0.600</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.100000</td>\n",
       "      <td>0.0125</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.088000</td>\n",
       "      <td>0.888889</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000</td>\n",
       "      <td>1.416667</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.0000</td>\n",
       "      <td>1.615385</td>\n",
       "      <td>0.158333</td>\n",
       "      <td>0.011538</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.705882</td>\n",
       "      <td>19.375</td>\n",
       "      <td>0.053571</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3270</th>\n",
       "      <td>Preparaten</td>\n",
       "      <td>Duocal ss poeder Milupa</td>\n",
       "      <td>0.005</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.0000</td>\n",
       "      <td>0.007143</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.001064</td>\n",
       "      <td>0.013333</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.0000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>5.592308</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>305</th>\n",
       "      <td>Melk en melkproducten</td>\n",
       "      <td>Kwark magere</td>\n",
       "      <td>0.128</td>\n",
       "      <td>2.222222</td>\n",
       "      <td>0.993333</td>\n",
       "      <td>0.0000</td>\n",
       "      <td>0.142857</td>\n",
       "      <td>5.181818</td>\n",
       "      <td>0.034681</td>\n",
       "      <td>0.030000</td>\n",
       "      <td>0.006667</td>\n",
       "      <td>0.011111</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>2.125</td>\n",
       "      <td>3.333333</td>\n",
       "      <td>2.307692</td>\n",
       "      <td>0.0625</td>\n",
       "      <td>0.769231</td>\n",
       "      <td>0.291667</td>\n",
       "      <td>0.037692</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000</td>\n",
       "      <td>0.151786</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1255</th>\n",
       "      <td>Vis</td>\n",
       "      <td>Makreel in water blik</td>\n",
       "      <td>0.020</td>\n",
       "      <td>31.111111</td>\n",
       "      <td>0.326667</td>\n",
       "      <td>0.2500</td>\n",
       "      <td>0.128571</td>\n",
       "      <td>2.454545</td>\n",
       "      <td>0.048511</td>\n",
       "      <td>0.277333</td>\n",
       "      <td>0.041111</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>2.866667</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000</td>\n",
       "      <td>2.500000</td>\n",
       "      <td>20.769231</td>\n",
       "      <td>3.3125</td>\n",
       "      <td>10.769231</td>\n",
       "      <td>28.916667</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.176471</td>\n",
       "      <td>0.175</td>\n",
       "      <td>0.332143</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "                              Voedselgroep  \\\n",
       "NEVO-code                                    \n",
       "2783                                 Brood   \n",
       "2423       Vetten, oliën en hartige sauzen   \n",
       "3270                            Preparaten   \n",
       "305                  Melk en melkproducten   \n",
       "1255                                   Vis   \n",
       "\n",
       "                                    Voedingsmiddel  Calcium (mg)  Koper (mg)  \\\n",
       "NEVO-code                                                                      \n",
       "2783       Brood meergranen- m div zaden obv bloem         0.046   26.666667   \n",
       "2423               Halvarineproduct Blue Band Idee         0.600    0.000000   \n",
       "3270                       Duocal ss poeder Milupa         0.005    0.000000   \n",
       "305                                   Kwark magere         0.128    2.222222   \n",
       "1255                         Makreel in water blik         0.020   31.111111   \n",
       "\n",
       "           Jodium (ug)  IJzer totaal (mg)  Fosfor (mg)  Zink (mg)  \\\n",
       "NEVO-code                                                           \n",
       "2783          4.706667             2.1250     0.225714  10.727273   \n",
       "2423          0.100000             0.0125     0.000000   0.000000   \n",
       "3270          0.000000             0.0000     0.007143   0.000000   \n",
       "305           0.993333             0.0000     0.142857   5.181818   \n",
       "1255          0.326667             0.2500     0.128571   2.454545   \n",
       "\n",
       "           Kalium (mg)  Natrium (mg)  Retinol activiteit equiv.(RAE) (ug)  \\\n",
       "NEVO-code                                                                   \n",
       "2783          0.040851      0.296000                             0.000000   \n",
       "2423          0.000000      0.088000                             0.888889   \n",
       "3270          0.001064      0.013333                             0.000000   \n",
       "305           0.034681      0.030000                             0.006667   \n",
       "1255          0.048511      0.277333                             0.041111   \n",
       "\n",
       "           Vitamine C (mg)  Cholecalciferol - vit D3 (ug)  \\\n",
       "NEVO-code                                                   \n",
       "2783              0.000000                       0.000000   \n",
       "2423              0.000000                       0.000000   \n",
       "3270              0.000000                       0.000000   \n",
       "305               0.011111                       0.000000   \n",
       "1255              0.000000                       2.866667   \n",
       "\n",
       "           Alfa-tocoferol (mg)  Vitamine K totaal (ug)  Vitamine B1 (mg)  \\\n",
       "NEVO-code                                                                  \n",
       "2783                  0.066667                   0.000          9.166667   \n",
       "2423                  0.000000                   0.000          1.416667   \n",
       "3270                  0.000000                   0.000          0.000000   \n",
       "305                   0.000000                   2.125          3.333333   \n",
       "1255                  0.000000                   0.000          2.500000   \n",
       "\n",
       "           Vitamine B2 (mg)  Nicotinezuur (mg)  Vitamine B6 (mg)  \\\n",
       "NEVO-code                                                          \n",
       "2783               5.384615             1.0000         72.307692   \n",
       "2423               0.000000             0.0000          1.615385   \n",
       "3270               0.000000             0.0000          0.000000   \n",
       "305                2.307692             0.0625          0.769231   \n",
       "1255              20.769231             3.3125         10.769231   \n",
       "\n",
       "           Vitamine B12 (ug)  Koolhydraten totaal (g)  \\\n",
       "NEVO-code                                               \n",
       "2783                0.000000                 0.310769   \n",
       "2423                0.158333                 0.011538   \n",
       "3270                0.000000                 5.592308   \n",
       "305                 0.291667                 0.037692   \n",
       "1255               28.916667                 0.000000   \n",
       "\n",
       "           Voedingsvezel totaal (g)  C18:2 n-6 cis (linolzuur) (g)  \\\n",
       "NEVO-code                                                            \n",
       "2783                       1.421053                       1.058824   \n",
       "2423                       0.000000                       0.705882   \n",
       "3270                       0.000000                       0.000000   \n",
       "305                        0.000000                       0.000000   \n",
       "1255                       0.000000                       0.176471   \n",
       "\n",
       "           C18:3 n-3 cis (ALA) (g)  Eiwit totaal (g)  \n",
       "NEVO-code                                             \n",
       "2783                         2.500          0.223214  \n",
       "2423                        19.375          0.053571  \n",
       "3270                         0.000          0.000000  \n",
       "305                          0.000          0.151786  \n",
       "1255                         0.175          0.332143  "
      ]
     },
     "execution_count": 234,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# Extract the RDA for the user.\n",
    "me = 'Males: 19–30 y'\n",
    "my_rda = rda_df.loc[me].dropna()\n",
    "# Drop Magnesium because it is very difficult to find in foods.\n",
    "my_rda.drop(['Magnesium (mg)', 'Foliumzuur (ug)', 'Selenium (ug)'], axis=0, inplace=True) \n",
    "# Keep only the nutrients we will be evaluating on and\n",
    "# normalize the nutrients to (% of RDA) / (household measure).\n",
    "pct_rda_df = nutrients_df.fillna(0.)[['Voedselgroep', 'Voedingsmiddel'] + my_rda.index.tolist()].copy()\n",
    "pct_rda_df.loc[:, my_rda.index] = pct_rda_df.loc[:, my_rda.index].divide(my_rda, axis=1)\n",
    "pct_rda_df.sample(5)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 235,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "# Drop exotic or difficult to find foods.\n",
    "drop = [\n",
    "    #3131,  # Graanproducten en bindmiddelen\tOntbijtproduct Pyjamapapje fijne tarwe granen,\n",
    "    439,   # Graanproducten en bindmiddelen\tGelatine\n",
    "    83,    # Eieren\tEi kippen- rauw gem\n",
    "    #1233,  # Kruiden en specerijen\tSambal oelek Na-\n",
    "    #1232,  # Sambal oelek normaal\n",
    "]\n",
    "pct_rda_df.drop(drop, axis=0, inplace=True)\n",
    "is_preparaat = \\\n",
    "    pct_rda_df.Voedingsmiddel.str.contains('iroop') | \\\n",
    "    pct_rda_df.Voedingsmiddel.str.contains('nergy') | \\\n",
    "    pct_rda_df.Voedingsmiddel.str.contains('ortdrank') | \\\n",
    "    pct_rda_df.Voedingsmiddel.str.contains('oeder') | \\\n",
    "    pct_rda_df.Voedingsmiddel.str.contains('iksap') | \\\n",
    "    pct_rda_df.Voedingsmiddel.str.contains('portvoed') | \\\n",
    "    pct_rda_df.Voedingsmiddel.str.contains('hicken') | \\\n",
    "    pct_rda_df.Voedingsmiddel.str.contains('owder')\n",
    "pct_rda_df = pct_rda_df.loc[~is_preparaat, :]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 238,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "# Variables.\n",
    "purchases = pulp.LpVariable.dicts(\n",
    "    'purchases',\n",
    "    indexs=range(len(pct_rda_df)),\n",
    "    lowBound=0,\n",
    "    upBound=1,\n",
    "    cat='Binary')\n",
    "portions = pulp.LpVariable.dicts(\n",
    "    'portions',\n",
    "    indexs=range(len(pct_rda_df)),\n",
    "    lowBound=0,\n",
    "    upBound=50,\n",
    "    cat='Continuous')\n",
    "rda_errors = pulp.LpVariable.dicts(\n",
    "    'rda_errors',\n",
    "    indexs=range(len(my_rda)),\n",
    "    lowBound=0,\n",
    "    upBound=2,\n",
    "    cat='Continuous')\n",
    "max_rda_error = pulp.LpVariable(\n",
    "    'max_rda_error',\n",
    "    lowBound=0,\n",
    "    upBound=2,\n",
    "    cat='Continuous')\n",
    "total_weight = pulp.LpVariable(\n",
    "    'total_weight',\n",
    "    lowBound=0,\n",
    "    upBound=2000,\n",
    "    cat='Continuous')\n",
    "\n",
    "# Cost function.\n",
    "rda_vs_weight = 0.95\n",
    "max_foods = 14\n",
    "max_weight_per_food = 200\n",
    "max_spices_weight = 3\n",
    "prob = pulp.LpProblem('food', pulp.LpMinimize)\n",
    "# prob += rda_vs_weight * max_rda_error + (1 - rda_vs_weight) * total_weight / 1e3\n",
    "prob += rda_vs_weight * pulp.lpSum(rda_errors) / len(rda_errors) + 0.25 * max_rda_error + (1 - rda_vs_weight) * total_weight / 1e3\n",
    "\n",
    "# Constraints:\n",
    "# - Total weight.\n",
    "prob += total_weight == pulp.lpSum(100 * portions[portion] for portion in portions)\n",
    "# - Max RDA error.\n",
    "for rda_error in rda_errors:\n",
    "    prob += max_rda_error >= rda_errors[rda_error]\n",
    "# - Not too much spices.\n",
    "for portion in portions:\n",
    "    if pct_rda_df.iat[portion, 0] == 'Kruiden en specerijen':\n",
    "        prob += 100 * portions[portion] <= max_spices_weight\n",
    "# - No more than max_weight per food.\n",
    "for portion in portions:\n",
    "    prob += 100 * portions[portion] <= max_weight_per_food\n",
    "# - No more than max_foods.\n",
    "prob += pulp.lpSum(purchases[purchase] for purchase in purchases) <= max_foods\n",
    "# - Purchase is binary.\n",
    "for portion, purchase in zip(portions, purchases):\n",
    "    prob += portions[portion] <= 10000 * purchases[purchase]\n",
    "# - RDA error is absolute error per nutrient.\n",
    "for rda_error in rda_errors:\n",
    "    raw_error = pulp.lpSum(pct_rda_df.iat[portion, rda_error + 2] * portions[portion] for portion in portions) - 1.\n",
    "    prob += -rda_errors[rda_error] <= raw_error\n",
    "    prob += raw_error <= rda_errors[rda_error]\n",
    "    # prob += raw_error == rda_errors[rda_error]\n",
    "prob.writeLP('food.lp')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 239,
   "metadata": {
    "scrolled": false
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "0"
      ]
     },
     "execution_count": 239,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Med RDA error is 0.049559443\n",
      "Avg RDA error is 0.12606019669565213\n",
      "Max RDA error is 0.33224947\n",
      "Total weight is 537.14015g\n"
     ]
    },
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style>\n",
       "    .dataframe thead tr:only-child th {\n",
       "        text-align: right;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: left;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>Voedselgroep</th>\n",
       "      <th>Voedingsmiddel</th>\n",
       "      <th>Gewicht</th>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>NEVO-code</th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>85</th>\n",
       "      <td>Eieren</td>\n",
       "      <td>Eidooier kippen- rauw</td>\n",
       "      <td>29.692627</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>126</th>\n",
       "      <td>Groenten</td>\n",
       "      <td>Sterkers rauw</td>\n",
       "      <td>2.175914</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>317</th>\n",
       "      <td>Vetten, oliën en hartige sauzen</td>\n",
       "      <td>Olie zonnebloem-</td>\n",
       "      <td>2.056726</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1140</th>\n",
       "      <td>Groenten</td>\n",
       "      <td>Groente erwtensoep- diepvries onbereid</td>\n",
       "      <td>200.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1212</th>\n",
       "      <td>Diversen</td>\n",
       "      <td>Zout mineraal- Na- LoSalt</td>\n",
       "      <td>7.824516</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1214</th>\n",
       "      <td>Sojaproducten en vegetarische producten</td>\n",
       "      <td>Ketjap zoet Na-</td>\n",
       "      <td>1.519350</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1612</th>\n",
       "      <td>Vis</td>\n",
       "      <td>Forel regenboog- bereid in magnetron z toev</td>\n",
       "      <td>2.214068</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1897</th>\n",
       "      <td>Diversen</td>\n",
       "      <td>Zeewier kelp rauw</td>\n",
       "      <td>1.313881</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2005</th>\n",
       "      <td>Graanproducten en bindmiddelen</td>\n",
       "      <td>Ontbijtproduct Special K Original Kellogg's</td>\n",
       "      <td>21.304192</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2065</th>\n",
       "      <td>Vetten, oliën en hartige sauzen</td>\n",
       "      <td>Margarineproduct 7017</td>\n",
       "      <td>25.204891</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2132</th>\n",
       "      <td>Preparaten</td>\n",
       "      <td>PreOp p 100</td>\n",
       "      <td>68.175961</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2151</th>\n",
       "      <td>Alcoholische en niet-alcoholische dranken</td>\n",
       "      <td>Sap vruchten- Coolbest Vitaday original</td>\n",
       "      <td>48.882432</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3078</th>\n",
       "      <td>Brood</td>\n",
       "      <td>Brood meergranen- glutenvrij Rustico Schar</td>\n",
       "      <td>46.543716</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3166</th>\n",
       "      <td>Kaas</td>\n",
       "      <td>Kaasproduct m plantaardige vetten Kees gerijpt</td>\n",
       "      <td>80.231874</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "                                        Voedselgroep  \\\n",
       "NEVO-code                                              \n",
       "85                                            Eieren   \n",
       "126                                         Groenten   \n",
       "317                  Vetten, oliën en hartige sauzen   \n",
       "1140                                        Groenten   \n",
       "1212                                        Diversen   \n",
       "1214         Sojaproducten en vegetarische producten   \n",
       "1612                                             Vis   \n",
       "1897                                        Diversen   \n",
       "2005                  Graanproducten en bindmiddelen   \n",
       "2065                 Vetten, oliën en hartige sauzen   \n",
       "2132                                      Preparaten   \n",
       "2151       Alcoholische en niet-alcoholische dranken   \n",
       "3078                                           Brood   \n",
       "3166                                            Kaas   \n",
       "\n",
       "                                           Voedingsmiddel     Gewicht  \n",
       "NEVO-code                                                              \n",
       "85                                  Eidooier kippen- rauw   29.692627  \n",
       "126                                         Sterkers rauw    2.175914  \n",
       "317                                      Olie zonnebloem-    2.056726  \n",
       "1140               Groente erwtensoep- diepvries onbereid  200.000000  \n",
       "1212                            Zout mineraal- Na- LoSalt    7.824516  \n",
       "1214                                      Ketjap zoet Na-    1.519350  \n",
       "1612          Forel regenboog- bereid in magnetron z toev    2.214068  \n",
       "1897                                    Zeewier kelp rauw    1.313881  \n",
       "2005          Ontbijtproduct Special K Original Kellogg's   21.304192  \n",
       "2065                                Margarineproduct 7017   25.204891  \n",
       "2132                                          PreOp p 100   68.175961  \n",
       "2151              Sap vruchten- Coolbest Vitaday original   48.882432  \n",
       "3078           Brood meergranen- glutenvrij Rustico Schar   46.543716  \n",
       "3166       Kaasproduct m plantaardige vetten Kees gerijpt   80.231874  "
      ]
     },
     "execution_count": 239,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style>\n",
       "    .dataframe thead tr:only-child th {\n",
       "        text-align: right;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: left;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>Calcium (mg)</th>\n",
       "      <th>Koper (mg)</th>\n",
       "      <th>Jodium (ug)</th>\n",
       "      <th>IJzer totaal (mg)</th>\n",
       "      <th>Fosfor (mg)</th>\n",
       "      <th>Zink (mg)</th>\n",
       "      <th>Kalium (mg)</th>\n",
       "      <th>Natrium (mg)</th>\n",
       "      <th>Retinol activiteit equiv.(RAE) (ug)</th>\n",
       "      <th>Vitamine C (mg)</th>\n",
       "      <th>Cholecalciferol - vit D3 (ug)</th>\n",
       "      <th>Alfa-tocoferol (mg)</th>\n",
       "      <th>Vitamine K totaal (ug)</th>\n",
       "      <th>Vitamine B1 (mg)</th>\n",
       "      <th>Vitamine B2 (mg)</th>\n",
       "      <th>Nicotinezuur (mg)</th>\n",
       "      <th>Vitamine B6 (mg)</th>\n",
       "      <th>Vitamine B12 (ug)</th>\n",
       "      <th>Koolhydraten totaal (g)</th>\n",
       "      <th>Voedingsvezel totaal (g)</th>\n",
       "      <th>C18:2 n-6 cis (linolzuur) (g)</th>\n",
       "      <th>C18:3 n-3 cis (ALA) (g)</th>\n",
       "      <th>Eiwit totaal (g)</th>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>NEVO-code</th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>85</th>\n",
       "      <td>0.044539</td>\n",
       "      <td>0.049488</td>\n",
       "      <td>0.136982</td>\n",
       "      <td>0.218983</td>\n",
       "      <td>0.262568</td>\n",
       "      <td>0.116071</td>\n",
       "      <td>0.006444</td>\n",
       "      <td>0.007918</td>\n",
       "      <td>0.177166</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.890779</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.049488</td>\n",
       "      <td>0.114202</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.036545</td>\n",
       "      <td>0.047013</td>\n",
       "      <td>0.000457</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.036679</td>\n",
       "      <td>0.005567</td>\n",
       "      <td>0.088548</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>126</th>\n",
       "      <td>0.005440</td>\n",
       "      <td>0.217591</td>\n",
       "      <td>0.001306</td>\n",
       "      <td>0.006800</td>\n",
       "      <td>0.002487</td>\n",
       "      <td>0.090993</td>\n",
       "      <td>0.002315</td>\n",
       "      <td>0.001160</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.012088</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.982606</td>\n",
       "      <td>0.018133</td>\n",
       "      <td>0.016738</td>\n",
       "      <td>0.001360</td>\n",
       "      <td>0.284543</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.001890</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000777</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>317</th>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.667751</td>\n",
       "      <td>0.009769</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000032</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.751310</td>\n",
       "      <td>0.012855</td>\n",
       "      <td>0.001102</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1140</th>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.250000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.033333</td>\n",
       "      <td>0.335556</td>\n",
       "      <td>0.222222</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.015385</td>\n",
       "      <td>0.178947</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.071429</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1212</th>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.576018</td>\n",
       "      <td>0.683341</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1214</th>\n",
       "      <td>0.000912</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.100657</td>\n",
       "      <td>0.001085</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.020366</td>\n",
       "      <td>0.000810</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.006428</td>\n",
       "      <td>0.000400</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.001221</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1612</th>\n",
       "      <td>0.000443</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.003690</td>\n",
       "      <td>0.002768</td>\n",
       "      <td>0.006958</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.001578</td>\n",
       "      <td>0.001299</td>\n",
       "      <td>0.000344</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.313660</td>\n",
       "      <td>0.051094</td>\n",
       "      <td>0.156369</td>\n",
       "      <td>0.068125</td>\n",
       "      <td>0.004613</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.001302</td>\n",
       "      <td>0.000138</td>\n",
       "      <td>0.009370</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1897</th>\n",
       "      <td>0.002207</td>\n",
       "      <td>0.189783</td>\n",
       "      <td>0.661583</td>\n",
       "      <td>0.047628</td>\n",
       "      <td>0.000788</td>\n",
       "      <td>0.146916</td>\n",
       "      <td>0.000249</td>\n",
       "      <td>0.002041</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000438</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.054745</td>\n",
       "      <td>0.151602</td>\n",
       "      <td>0.004106</td>\n",
       "      <td>0.020214</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000839</td>\n",
       "      <td>0.004495</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000399</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2005</th>\n",
       "      <td>0.012783</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.034087</td>\n",
       "      <td>0.308911</td>\n",
       "      <td>0.051739</td>\n",
       "      <td>0.290512</td>\n",
       "      <td>0.009972</td>\n",
       "      <td>0.063913</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.317196</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.003551</td>\n",
       "      <td>0.319563</td>\n",
       "      <td>0.376920</td>\n",
       "      <td>0.400785</td>\n",
       "      <td>0.376920</td>\n",
       "      <td>0.037282</td>\n",
       "      <td>0.129464</td>\n",
       "      <td>0.025229</td>\n",
       "      <td>0.075191</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.034239</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2065</th>\n",
       "      <td>0.002520</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.025205</td>\n",
       "      <td>0.003151</td>\n",
       "      <td>0.008282</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.002252</td>\n",
       "      <td>0.038984</td>\n",
       "      <td>0.224043</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.109221</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.326181</td>\n",
       "      <td>0.472592</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2132</th>\n",
       "      <td>0.004091</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000974</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.017697</td>\n",
       "      <td>0.022725</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.660782</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2151</th>\n",
       "      <td>0.004399</td>\n",
       "      <td>0.543138</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.061103</td>\n",
       "      <td>0.007682</td>\n",
       "      <td>0.355509</td>\n",
       "      <td>0.015497</td>\n",
       "      <td>0.000652</td>\n",
       "      <td>0.108084</td>\n",
       "      <td>0.217255</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.004074</td>\n",
       "      <td>0.244412</td>\n",
       "      <td>0.263213</td>\n",
       "      <td>0.244412</td>\n",
       "      <td>0.263213</td>\n",
       "      <td>0.026478</td>\n",
       "      <td>0.039106</td>\n",
       "      <td>0.005146</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.026187</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3078</th>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.024823</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.155146</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.147508</td>\n",
       "      <td>0.783894</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.036570</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3166</th>\n",
       "      <td>0.922667</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.112325</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.767934</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.015364</td>\n",
       "      <td>0.320927</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.601739</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.141586</td>\n",
       "      <td>0.300870</td>\n",
       "      <td>0.415486</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "           Calcium (mg)  Koper (mg)  Jodium (ug)  IJzer totaal (mg)  \\\n",
       "NEVO-code                                                             \n",
       "85             0.044539    0.049488     0.136982           0.218983   \n",
       "126            0.005440    0.217591     0.001306           0.006800   \n",
       "317            0.000000    0.000000     0.000000           0.000000   \n",
       "1140           0.000000    0.000000     0.000000           0.250000   \n",
       "1212           0.000000    0.000000     0.000000           0.000000   \n",
       "1214           0.000912    0.000000     0.000000           0.100657   \n",
       "1612           0.000443    0.000000     0.003690           0.002768   \n",
       "1897           0.002207    0.189783     0.661583           0.047628   \n",
       "2005           0.012783    0.000000     0.034087           0.308911   \n",
       "2065           0.002520    0.000000     0.025205           0.003151   \n",
       "2132           0.004091    0.000000     0.000000           0.000000   \n",
       "2151           0.004399    0.543138     0.000000           0.061103   \n",
       "3078           0.000000    0.000000     0.024823           0.000000   \n",
       "3166           0.922667    0.000000     0.112325           0.000000   \n",
       "\n",
       "           Fosfor (mg)  Zink (mg)  Kalium (mg)  Natrium (mg)  \\\n",
       "NEVO-code                                                      \n",
       "85            0.262568   0.116071     0.006444      0.007918   \n",
       "126           0.002487   0.090993     0.002315      0.001160   \n",
       "317           0.000000   0.000000     0.000000      0.000000   \n",
       "1140          0.000000   0.000000     0.000000      0.033333   \n",
       "1212          0.000000   0.000000     0.576018      0.683341   \n",
       "1214          0.001085   0.000000     0.020366      0.000810   \n",
       "1612          0.006958   0.000000     0.001578      0.001299   \n",
       "1897          0.000788   0.146916     0.000249      0.002041   \n",
       "2005          0.051739   0.290512     0.009972      0.063913   \n",
       "2065          0.008282   0.000000     0.002252      0.038984   \n",
       "2132          0.000974   0.000000     0.017697      0.022725   \n",
       "2151          0.007682   0.355509     0.015497      0.000652   \n",
       "3078          0.000000   0.000000     0.000000      0.155146   \n",
       "3166          0.767934   0.000000     0.015364      0.320927   \n",
       "\n",
       "           Retinol activiteit equiv.(RAE) (ug)  Vitamine C (mg)  \\\n",
       "NEVO-code                                                         \n",
       "85                                    0.177166         0.000000   \n",
       "126                                   0.000000         0.012088   \n",
       "317                                   0.000000         0.000000   \n",
       "1140                                  0.335556         0.222222   \n",
       "1212                                  0.000000         0.000000   \n",
       "1214                                  0.000000         0.000000   \n",
       "1612                                  0.000344         0.000000   \n",
       "1897                                  0.000000         0.000438   \n",
       "2005                                  0.000000         0.317196   \n",
       "2065                                  0.224043         0.000000   \n",
       "2132                                  0.000000         0.000000   \n",
       "2151                                  0.108084         0.217255   \n",
       "3078                                  0.000000         0.000000   \n",
       "3166                                  0.000000         0.000000   \n",
       "\n",
       "           Cholecalciferol - vit D3 (ug)  Alfa-tocoferol (mg)  \\\n",
       "NEVO-code                                                       \n",
       "85                              0.890779             0.000000   \n",
       "126                             0.000000             0.000000   \n",
       "317                             0.000000             0.667751   \n",
       "1140                            0.000000             0.000000   \n",
       "1212                            0.000000             0.000000   \n",
       "1214                            0.000000             0.000000   \n",
       "1612                            0.000000             0.000000   \n",
       "1897                            0.000000             0.000000   \n",
       "2005                            0.000000             0.000000   \n",
       "2065                            0.109221             0.000000   \n",
       "2132                            0.000000             0.000000   \n",
       "2151                            0.000000             0.000000   \n",
       "3078                            0.000000             0.000000   \n",
       "3166                            0.000000             0.000000   \n",
       "\n",
       "           Vitamine K totaal (ug)  Vitamine B1 (mg)  Vitamine B2 (mg)  \\\n",
       "NEVO-code                                                               \n",
       "85                       0.000000          0.049488          0.114202   \n",
       "126                      0.982606          0.018133          0.016738   \n",
       "317                      0.009769          0.000000          0.000000   \n",
       "1140                     0.000000          0.000000          0.000000   \n",
       "1212                     0.000000          0.000000          0.000000   \n",
       "1214                     0.000000          0.000000          0.000000   \n",
       "1612                     0.000000          0.313660          0.051094   \n",
       "1897                     0.000000          0.054745          0.151602   \n",
       "2005                     0.003551          0.319563          0.376920   \n",
       "2065                     0.000000          0.000000          0.000000   \n",
       "2132                     0.000000          0.000000          0.000000   \n",
       "2151                     0.004074          0.244412          0.263213   \n",
       "3078                     0.000000          0.000000          0.000000   \n",
       "3166                     0.000000          0.000000          0.000000   \n",
       "\n",
       "           Nicotinezuur (mg)  Vitamine B6 (mg)  Vitamine B12 (ug)  \\\n",
       "NEVO-code                                                           \n",
       "85                  0.000000          0.036545           0.047013   \n",
       "126                 0.001360          0.284543           0.000000   \n",
       "317                 0.000000          0.000000           0.000000   \n",
       "1140                0.000000          0.000000           0.000000   \n",
       "1212                0.000000          0.000000           0.000000   \n",
       "1214                0.000000          0.000000           0.000000   \n",
       "1612                0.156369          0.068125           0.004613   \n",
       "1897                0.004106          0.020214           0.000000   \n",
       "2005                0.400785          0.376920           0.037282   \n",
       "2065                0.000000          0.000000           0.000000   \n",
       "2132                0.000000          0.000000           0.000000   \n",
       "2151                0.244412          0.263213           0.026478   \n",
       "3078                0.000000          0.000000           0.000000   \n",
       "3166                0.000000          0.000000           0.601739   \n",
       "\n",
       "           Koolhydraten totaal (g)  Voedingsvezel totaal (g)  \\\n",
       "NEVO-code                                                      \n",
       "85                        0.000457                  0.000000   \n",
       "126                       0.000000                  0.001890   \n",
       "317                       0.000032                  0.000000   \n",
       "1140                      0.015385                  0.178947   \n",
       "1212                      0.000000                  0.000000   \n",
       "1214                      0.006428                  0.000400   \n",
       "1612                      0.000000                  0.000000   \n",
       "1897                      0.000839                  0.004495   \n",
       "2005                      0.129464                  0.025229   \n",
       "2065                      0.000000                  0.000000   \n",
       "2132                      0.660782                  0.000000   \n",
       "2151                      0.039106                  0.005146   \n",
       "3078                      0.147508                  0.783894   \n",
       "3166                      0.000000                  0.000000   \n",
       "\n",
       "           C18:2 n-6 cis (linolzuur) (g)  C18:3 n-3 cis (ALA) (g)  \\\n",
       "NEVO-code                                                           \n",
       "85                              0.036679                 0.005567   \n",
       "126                             0.000000                 0.000000   \n",
       "317                             0.751310                 0.012855   \n",
       "1140                            0.000000                 0.000000   \n",
       "1212                            0.000000                 0.000000   \n",
       "1214                            0.000000                 0.000000   \n",
       "1612                            0.001302                 0.000138   \n",
       "1897                            0.000000                 0.000000   \n",
       "2005                            0.075191                 0.000000   \n",
       "2065                            0.326181                 0.472592   \n",
       "2132                            0.000000                 0.000000   \n",
       "2151                            0.000000                 0.000000   \n",
       "3078                            0.000000                 0.000000   \n",
       "3166                            0.141586                 0.300870   \n",
       "\n",
       "           Eiwit totaal (g)  \n",
       "NEVO-code                    \n",
       "85                 0.088548  \n",
       "126                0.000777  \n",
       "317                0.001102  \n",
       "1140               0.071429  \n",
       "1212               0.000000  \n",
       "1214               0.001221  \n",
       "1612               0.009370  \n",
       "1897               0.000399  \n",
       "2005               0.034239  \n",
       "2065               0.000000  \n",
       "2132               0.000000  \n",
       "2151               0.026187  \n",
       "3078               0.036570  \n",
       "3166               0.415486  "
      ]
     },
     "execution_count": 239,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "text/plain": [
       "Calcium (mg)                           1.000000\n",
       "Koper (mg)                             1.000000\n",
       "Jodium (ug)                            1.000000\n",
       "IJzer totaal (mg)                      1.000000\n",
       "Fosfor (mg)                            1.110496\n",
       "Zink (mg)                              1.000000\n",
       "Kalium (mg)                            0.667751\n",
       "Natrium (mg)                           1.332249\n",
       "Retinol activiteit equiv.(RAE) (ug)    0.845194\n",
       "Vitamine C (mg)                        0.769200\n",
       "Cholecalciferol - vit D3 (ug)          1.000000\n",
       "Alfa-tocoferol (mg)                    0.667751\n",
       "Vitamine K totaal (ug)                 1.000000\n",
       "Vitamine B1 (mg)                       1.000000\n",
       "Vitamine B2 (mg)                       0.973769\n",
       "Nicotinezuur (mg)                      0.807032\n",
       "Vitamine B6 (mg)                       1.049559\n",
       "Vitamine B12 (ug)                      0.717125\n",
       "Koolhydraten totaal (g)                1.000000\n",
       "Voedingsvezel totaal (g)               1.000000\n",
       "C18:2 n-6 cis (linolzuur) (g)          1.332249\n",
       "C18:3 n-3 cis (ALA) (g)                0.792022\n",
       "Eiwit totaal (g)                       0.685328\n",
       "dtype: float64"
      ]
     },
     "execution_count": 239,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "solver = pulp.solvers.COIN_CMD(\n",
    "    maxSeconds=30,\n",
    "    options=['feasibilityPump off', 'probing off', 'knapsack off', 'clique off', 'flowcover off'])\n",
    "prob.solve(solver=solver)\n",
    "\n",
    "print(f'Med RDA error is {np.median([rda_errors[i].value() for i in rda_errors])}')\n",
    "print(f'Avg RDA error is {np.mean([rda_errors[i].value() for i in rda_errors])}')\n",
    "print(f'Max RDA error is {np.max([rda_errors[i].value() for i in rda_errors])}')\n",
    "print(f'Total weight is {total_weight.value()}g')\n",
    "\n",
    "x = np.array([portions[p].value() for p in portions])\n",
    "tmp = pct_rda_df.iloc[:, :2].loc[x > 0].copy()\n",
    "tmp['Gewicht'] = 100 * x[x > 0]\n",
    "tmp\n",
    "pct_rda_df.iloc[:, 2:].loc[x > 0].multiply(x[x > 0], axis=0)\n",
    "pct_rda_df.iloc[:, 2:].loc[x > 0].multiply(x[x > 0], axis=0).sum(axis=0)"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.6.2"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}