2.Take the fat data, and use the percentage of body fat, siri, as the response and the other variables, except brozek and density, as potential predictors. Remove every tenth observation from the data for use as a test sample. Use the remaining data as a training sample building the following models: (a) Linear regression with all predictors

require(MASS)
## Loading required package: MASS
data(fat, package="faraway")
head(fat)
##   brozek siri density age weight height adipos  free neck chest abdom   hip
## 1   12.6 12.3  1.0708  23 154.25  67.75   23.7 134.9 36.2  93.1  85.2  94.5
## 2    6.9  6.1  1.0853  22 173.25  72.25   23.4 161.3 38.5  93.6  83.0  98.7
## 3   24.6 25.3  1.0414  22 154.00  66.25   24.7 116.0 34.0  95.8  87.9  99.2
## 4   10.9 10.4  1.0751  26 184.75  72.25   24.9 164.7 37.4 101.8  86.4 101.2
## 5   27.8 28.7  1.0340  24 184.25  71.25   25.6 133.1 34.4  97.3 100.0 101.9
## 6   20.6 20.9  1.0502  24 210.25  74.75   26.5 167.0 39.0 104.5  94.4 107.8
##   thigh knee ankle biceps forearm wrist
## 1  59.0 37.3  21.9   32.0    27.4  17.1
## 2  58.7 37.3  23.4   30.5    28.9  18.2
## 3  59.6 38.9  24.0   28.8    25.2  16.6
## 4  60.1 37.3  22.8   32.4    29.4  18.2
## 5  63.2 42.2  24.0   32.2    27.7  17.7
## 6  66.0 42.0  25.6   35.7    30.6  18.8
dim(fat)
## [1] 252  18
testdata = seq(10,252,by=10)
fat
##     brozek siri density age weight height adipos  free neck chest abdom   hip
## 1     12.6 12.3  1.0708  23 154.25  67.75   23.7 134.9 36.2  93.1  85.2  94.5
## 2      6.9  6.1  1.0853  22 173.25  72.25   23.4 161.3 38.5  93.6  83.0  98.7
## 3     24.6 25.3  1.0414  22 154.00  66.25   24.7 116.0 34.0  95.8  87.9  99.2
## 4     10.9 10.4  1.0751  26 184.75  72.25   24.9 164.7 37.4 101.8  86.4 101.2
## 5     27.8 28.7  1.0340  24 184.25  71.25   25.6 133.1 34.4  97.3 100.0 101.9
## 6     20.6 20.9  1.0502  24 210.25  74.75   26.5 167.0 39.0 104.5  94.4 107.8
## 7     19.0 19.2  1.0549  26 181.00  69.75   26.2 146.6 36.4 105.1  90.7 100.3
## 8     12.8 12.4  1.0704  25 176.00  72.50   23.6 153.6 37.8  99.6  88.5  97.1
## 9      5.1  4.1  1.0900  25 191.00  74.00   24.6 181.3 38.1 100.9  82.5  99.9
## 10    12.0 11.7  1.0722  23 198.25  73.50   25.8 174.4 42.1  99.6  88.6 104.1
## 11     7.5  7.1  1.0830  26 186.25  74.50   23.6 172.3 38.5 101.5  83.6  98.2
## 12     8.5  7.8  1.0812  27 216.00  76.00   26.3 197.7 39.4 103.6  90.9 107.7
## 13    20.5 20.8  1.0513  32 180.50  69.50   26.3 143.5 38.4 102.0  91.6 103.9
## 14    20.8 21.2  1.0505  30 205.25  71.25   28.5 162.5 39.4 104.1 101.8 108.6
## 15    21.7 22.1  1.0484  35 187.75  69.50   27.4 147.0 40.5 101.3  96.4 100.1
## 16    20.5 20.9  1.0512  35 162.75  66.00   26.3 129.3 36.4  99.1  92.8  99.2
## 17    28.1 29.0  1.0333  34 195.75  71.00   27.3 140.8 38.9 101.9  96.4 105.2
## 18    22.4 22.9  1.0468  32 209.25  71.00   29.2 162.5 42.1 107.6  97.5 107.0
## 19    16.1 16.0  1.0622  28 183.75  67.75   28.2 154.3 38.0 106.8  89.6 102.4
## 20    16.5 16.5  1.0610  33 211.75  73.50   27.6 176.8 40.0 106.2 100.5 109.0
## 21    19.0 19.1  1.0551  28 179.00  68.00   27.3 145.1 39.1 103.3  95.9 104.9
## 22    15.3 15.2  1.0640  28 200.50  69.75   29.1 169.8 41.3 111.4  98.8 104.8
## 23    15.7 15.6  1.0631  31 140.25  68.25   21.2 118.2 33.9  86.0  76.4  94.6
## 24    17.6 17.7  1.0584  32 148.75  70.00   21.4 122.6 35.5  86.7  80.0  93.4
## 25    14.2 14.0  1.0668  28 151.25  67.75   23.2 129.8 34.5  90.2  76.3  95.8
## 26     4.6  3.7  1.0911  27 159.25  71.50   21.9 151.9 35.7  89.6  79.7  96.5
## 27     8.5  7.9  1.0811  34 131.50  67.50   20.3 120.3 36.2  88.6  74.6  85.3
## 28    22.4 22.9  1.0468  31 148.00  67.50   22.9 114.9 38.8  97.4  88.7  94.7
## 29     4.7  3.7  1.0910  27 133.25  64.75   22.4 127.0 36.4  93.5  73.9  88.5
## 30     9.4  8.8  1.0790  29 160.75  69.00   23.8 145.7 36.7  97.4  83.5  98.7
## 31    12.3 11.9  1.0716  32 182.00  73.75   23.6 159.7 38.7 100.5  88.7  99.8
## 32     6.5  5.7  1.0862  29 160.25  71.25   22.2 149.8 37.3  93.5  84.5 100.6
## 33    13.4 11.8  1.0719  27 168.00  71.25   23.3 142.5 38.1  93.0  79.1  94.5
## 34    20.9 21.3  1.0502  41 218.50  71.00   30.5 172.7 39.8 111.7 100.5 108.3
## 35    31.1 32.3  1.0263  41 247.25  73.50   32.2 170.4 42.1 117.0 115.6 116.1
## 36    38.2 40.1  1.0101  49 191.75  65.00   32.0 118.4 38.4 118.5 113.1 113.8
## 37    23.6 24.2  1.0438  40 202.25  70.00   29.1 154.5 38.5 106.5 100.9 106.2
## 38    27.5 28.4  1.0346  50 196.75  68.25   29.7 142.6 42.1 105.6  98.8 104.8
## 39    33.8 35.2  1.0202  46 363.15  72.25   48.9 240.5 51.2 136.2 148.1 147.7
## 40    31.3 32.6  1.0258  50 203.00  67.00   31.8 139.4 40.2 114.8 108.1 102.5
## 41    33.1 34.5  1.0217  45 262.75  68.75   39.1 175.8 43.2 128.3 126.2 125.6
## 42    31.7 32.9  1.0250  44 205.00  29.50   29.9 140.1 36.6 106.0 104.3 115.5
## 43    30.4 31.6  1.0279  48 217.00  70.00   31.2 151.1 37.3 113.3 111.2 114.1
## 44    30.8 32.0  1.0269  41 212.00  71.50   29.2 146.7 41.5 106.6 104.3 106.0
## 45     8.4  7.7  1.0814  39 125.25  68.00   19.1 114.7 31.5  85.1  76.0  88.2
## 46    14.1 13.9  1.0670  43 164.25  73.25   21.3 141.1 35.7  96.6  81.5  97.2
## 47    11.2 10.8  1.0742  40 133.50  67.50   20.6 118.5 33.6  88.2  73.7  88.5
## 48     6.4  5.6  1.0665  39 148.50  71.25   20.6 139.0 34.6  89.8  79.5  92.7
## 49    13.4 13.6  1.0678  45 135.75  68.50   20.4 117.6 32.8  92.3  83.4  90.4
## 50     5.0  4.0  1.0903  47 127.50  66.75   20.2 121.2 34.0  83.4  70.4  87.2
## 51    10.7 10.2  1.0756  47 158.25  72.25   21.3 141.4 34.9  90.2  86.7  98.3
## 52     7.4  6.6  1.0840  40 139.25  69.00   20.6 129.0 34.3  89.2  77.9  91.0
## 53     8.7  8.0  1.0807  51 137.25  67.75   21.1 125.3 36.5  89.7  82.0  89.1
## 54     7.1  6.3  1.0848  49 152.75  73.50   19.9 142.0 35.1  93.3  79.6  91.6
## 55     4.9  3.9  1.0906  42 136.25  67.50   21.1 129.6 37.8  87.6  77.6  88.6
## 56    22.2 22.6  1.0473  54 198.00  72.00   26.9 154.1 39.9 107.6 100.0  99.6
## 57    20.1 20.4  1.0524  58 181.50  68.00   27.6 145.1 39.1 100.0  99.8 102.5
## 58    27.1 28.0  1.0356  62 201.25  69.50   29.3 146.7 40.5 111.5 104.2 105.8
## 59    30.4 31.5  1.0280  54 202.50  70.75   28.4 141.0 40.5 115.4 105.3  97.0
## 60    24.0 24.6  1.0430  61 179.75  65.75   29.2 136.7 38.4 104.8  98.3  99.6
## 61    25.4 26.1  1.0396  62 216.00  73.25   28.2 161.2 41.4 112.3 104.8 103.1
## 62    28.8 29.8  1.0317  56 178.75  68.50   26.8 127.4 35.6 102.9  94.7 100.8
## 63    29.6 30.7  1.0298  54 193.25  70.25   27.6 136.1 38.0 107.6 102.4  99.4
## 64    25.1 25.8  1.0403  61 178.00  67.00   27.9 133.3 37.4 105.3  99.7  99.7
## 65    31.0 32.3  1.0264  57 205.50  70.00   29.5 141.7 40.1 105.3 105.5 108.3
## 66    28.9 30.0  1.0313  55 183.50  67.50   28.3 130.4 40.9 103.0 100.3 104.2
## 67    21.1 21.5  1.0499  54 151.50  70.75   21.3 119.6 35.6  90.0  83.9  93.9
## 68    14.0 13.8  1.0673  55 154.75  71.50   21.3 133.1 36.9  95.4  86.6  91.8
## 69     7.1  6.3  1.0847  54 155.25  69.25   22.8 144.2 37.5  89.3  78.4  96.1
## 70    13.2 12.9  1.0693  55 156.75  71.50   21.6 136.1 36.3  94.4  84.6  94.3
## 71    23.7 24.3  1.0439  62 167.50  71.50   23.1 127.8 35.5  97.6  91.5  98.5
## 72     9.4  8.8  1.0788  55 146.75  68.75   21.9 132.9 38.7  88.5  82.8  95.5
## 73     9.1  8.5  1.0796  56 160.75  73.75   20.8 146.1 36.4  93.6  82.9  96.3
## 74    13.7 13.5  1.0680  55 125.00  64.00   21.5 107.9 33.2  87.7  76.0  88.6
## 75    12.0 11.8  1.0720  61 143.00  65.75   23.3 125.9 36.5  93.4  83.3  93.0
## 76    18.3 18.5  1.0666  61 148.25  67.50   22.9 121.1 36.0  91.6  81.8  94.8
## 77     9.2  8.8  1.0790  57 162.50  69.50   23.7 147.5 38.7  91.6  78.8  94.3
## 78    21.7 22.2  1.0483  69 177.75  68.50   26.7 139.1 38.7 102.0  95.0  98.3
## 79    21.1 21.5  1.0498  81 161.25  70.25   23.0 127.2 37.8  96.4  95.4  99.3
## 80    18.6 18.8  1.0560  66 171.25  69.25   25.1 139.5 37.4 102.7  98.6 100.2
## 81    30.2 31.4  1.0283  67 163.75  67.75   25.1 114.3 38.4  97.7  95.8  97.1
## 82    26.0 26.8  1.0382  64 150.25  67.25   23.4 111.2 38.1  97.1  89.0  96.9
## 83    18.2 18.4  1.0568  64 190.25  72.75   25.3 155.6 39.3 103.1  97.8  99.6
## 84    26.2 27.0  1.0377  70 170.75  70.00   24.5 126.0 38.7 101.8  94.9  95.0
## 85    26.1 27.0  1.0378  72 168.00  69.25   24.7 124.1 38.5 101.4  99.8  96.2
## 86    25.8 26.6  1.0386  67 167.00  67.50   26.0 123.9 36.5  98.9  89.7  96.2
## 87    15.0 14.9  1.0648  72 157.75  67.25   24.6 134.1 37.7  97.5  88.1  96.9
## 88    22.6 23.1  1.0462  64 160.00  65.75   26.0 123.8 36.5 104.3  90.9  93.8
## 89     8.8  8.3  1.0800  46 176.75  72.50   23.7 161.1 38.0  97.3  86.0  99.3
## 90    14.3 14.1  1.0666  48 176.00  73.00   23.3 150.9 36.7  96.7  86.5  98.3
## 91    20.2 20.5  1.0520  46 177.00  70.00   25.4 141.3 37.2  99.7  95.6 102.2
## 92    18.1 18.2  1.0573  44 179.75  69.50   26.2 147.3 39.2 101.9  93.2 100.6
## 93     9.2  8.5  1.0795  47 165.25  70.50   23.4 150.1 37.5  97.2  83.1  95.4
## 94    24.2 24.9  1.0424  46 192.50  71.75   26.3 145.9 38.0 106.6  97.5 100.6
## 95     9.6  9.0  1.0785  47 184.25  74.50   23.4 166.6 37.3  99.6  88.8 101.4
## 96    17.3 17.4  1.0991  53 224.50  77.75   26.1 185.7 41.1 113.2  99.2 107.5
## 97    10.1  9.6  1.0770  38 188.75  73.25   24.8 169.6 37.5  99.1  91.6 102.4
## 98    11.1 11.3  1.0730  50 162.50  66.50   25.9 143.5 38.7  99.4  86.7  96.2
## 99    17.7 17.8  1.0582  46 156.50  68.25   23.7 128.8 35.9  95.1  88.2  92.8
## 100   21.7 22.2  1.0484  47 197.00  72.00   26.7 154.2 40.0 107.5  94.0 103.7
## 101   20.8 21.2  1.0506  49 198.50  73.50   25.9 157.2 40.1 106.5  95.0 101.7
## 102   20.1 20.4  1.0524  48 173.75  72.00   23.6 138.9 37.0  99.1  92.0  98.3
## 103   19.8 20.1  1.0530  41 172.75  71.25   24.0 138.6 36.3  96.7  89.2  98.3
## 104   21.9 22.3  1.0480  49 196.75  73.75   25.5 153.7 40.7 103.5  95.5 101.6
## 105   24.7 25.4  1.0412  43 177.00  69.25   26.0 133.2 39.6 104.0  98.6  99.5
## 106   17.8 18.0  1.0578  43 165.50  68.50   24.8 136.0 31.1  93.1  87.3  96.6
## 107   19.1 19.3  1.0547  43 200.25  73.50   26.0 162.0 38.6 105.2 102.8 103.6
## 108   18.2 18.3  1.0569  52 203.25  74.25   26.0 166.3 42.0 110.0 101.6 100.7
## 109   17.2 17.3  1.0593  43 194.00  75.50   24.0 160.6 38.5 110.1  88.7 102.1
## 110   21.0 21.4  1.0500  40 168.50  69.25   24.7 133.1 34.2  97.8  92.3 100.6
## 111   19.5 19.7  1.0538  43 170.75  68.50   25.6 137.5 37.2  96.3  90.6  99.3
## 112   27.1 28.0  1.0355  43 183.25  70.00   26.3 133.5 37.1 108.0 105.0 103.0
## 113   21.6 22.1  1.0486  47 178.25  70.00   25.6 139.7 40.2  99.7  95.0  98.6
## 114   20.9 21.3  1.0503  42 163.00  70.25   23.3 128.9 35.3  93.5  89.6  99.8
## 115   25.9 26.7  1.0384  48 175.25  71.75   24.0 129.9 38.0 100.7  92.4  97.5
## 116   16.7 16.7  1.0607  40 158.00  69.25   23.4 131.7 36.3  97.0  86.6  92.6
## 117   19.8 20.1  1.0529  48 177.25  72.75   23.6 142.1 36.8  96.0  90.0  99.7
## 118   14.1 13.9  1.0671  51 179.00  72.00   24.3 153.8 41.0  99.2  90.0  96.4
## 119   25.1 25.8  1.0404  40 191.00  74.00   24.6 143.1 38.3  95.4  92.4 104.3
## 120   17.9 18.1  1.0575  44 187.50  72.25   25.3 153.8 38.0 101.8  87.5 101.0
## 121   27.0 27.9  1.0358  52 206.50  74.50   26.2 150.7 40.8 104.3  99.2 104.1
## 122   24.6 25.3  1.0414  44 185.25  71.50   25.5 139.6 39.5  99.2  98.1 101.4
## 123   14.8 14.7  1.0652  40 160.25  68.75   23.9 136.5 36.9  99.3  83.3  97.5
## 124   16.0 16.0  1.0623  47 151.50  66.75   23.9 127.3 36.9  94.0  86.1  95.2
## 125   14.0 13.8  1.0674  50 161.00  66.50   25.6 138.5 37.7  98.9  84.1  94.0
## 126   17.4 17.5  1.0587  46 167.00  67.00   26.2 137.9 36.6 101.0  89.9 100.0
## 127   26.4 27.2  1.0373  42 177.50  68.75   26.4 130.7 38.9  98.7  92.1  98.5
## 128   17.4 17.4  1.0590  43 152.25  67.75   23.4 125.8 37.5  95.9  78.0  93.2
## 129   20.4 20.8  1.0515  40 192.25  73.25   25.2 153.0 39.8 103.9  93.5  99.5
## 130   15.0 14.9  1.0648  42 165.25  69.75   23.9 140.5 38.3  96.2  87.0  97.8
## 131   18.0 18.1  1.0575  49 171.75  71.50   23.7 140.9 35.5  97.8  90.1  95.8
## 132   22.2 22.7  1.0472  40 171.25  70.50   24.3 133.3 36.3  94.6  90.3  99.1
## 133   23.1 23.6  1.0452  47 197.00  73.25   25.8 151.2 37.8 103.6  99.8 103.2
## 134   25.3 26.1  1.0398  50 157.00  66.75   24.8 117.2 37.8 100.4  89.4  92.3
## 135   23.8 24.4  1.0435  41 168.25  69.50   24.5 128.3 36.5  98.4  87.2  98.4
## 136   26.3 27.1  1.0374  44 186.00  69.75   26.8 137.1 37.8 104.6 101.1 102.1
## 137   21.4 21.8  1.0491  39 166.75  70.75   23.5 131.0 37.0  92.9  86.1  95.6
## 138   28.4 29.4  1.0325  43 187.75  74.00   24.1 134.4 37.7  97.8  98.6 100.6
## 139   21.8 22.4  1.0481  40 168.25  71.25   23.3 131.6 34.3  98.3  88.5  98.3
## 140   20.1 20.4  1.0522  49 212.75  75.00   26.6 169.9 40.8 104.7 106.6 107.7
## 141   24.3 24.9  1.0422  40 176.75  71.00   24.6 133.8 37.4  98.6  93.1 101.6
## 142   18.1 18.3  1.0571  40 173.25  69.50   25.3 141.8 36.5  99.5  93.0  99.3
## 143   22.7 23.3  1.0459  52 167.00  67.75   25.6 129.0 37.5 102.7  91.0  98.9
## 144    9.9  9.4  1.0775  23 159.75  72.25   21.6 143.9 35.5  92.1  77.1  93.9
## 145   10.8 10.3  1.0754  23 188.15  77.50   22.1 168.4 38.0  96.6  85.3 102.5
## 146   14.4 14.2  1.0664  24 156.00  70.75   21.9 133.6 35.7  92.7  81.9  95.3
## 147   19.0 19.2  1.0550  24 208.50  72.75   27.7 168.9 39.2 102.0  99.1 110.1
## 148   28.6 29.6  1.0322  25 206.50  69.75   29.8 147.5 40.9 110.9 100.5 106.2
## 149    6.1  5.3  1.0873  25 143.75  72.50   19.3 135.0 35.2  92.3  76.5  92.1
## 150   24.5 25.2  1.0416  26 223.00  70.25   31.8 168.3 40.6 114.1 106.8 113.9
## 151    9.9  9.4  1.0776  26 152.25  69.00   22.5 137.2 35.4  92.9  77.6  93.5
## 152   19.1 19.6  1.0542  26 241.75  74.50   30.7 195.1 41.8 108.3 102.9 114.4
## 153   10.6 10.1  1.0758  27 146.00  72.25   19.7 130.5 34.1  88.5  72.8  91.1
## 154   16.5 16.5  1.0610  27 156.75  67.25   24.4 130.9 37.9  94.0  88.2  95.2
## 155   20.5 21.0  1.0510  27 200.25  73.50   26.1 159.3 38.2 101.1 100.1 105.0
## 156   17.2 17.3  1.0594  28 171.50  75.25   21.6 142.0 35.6  92.1  83.5  98.3
## 157   30.1 31.2  1.0287  28 205.75  69.00   30.4 143.9 38.5 105.6 105.0 106.4
## 158   10.5 10.0  1.0761  28 182.50  72.25   24.6 163.4 37.0  98.5  90.8 102.5
## 159   12.8 12.5  1.0704  30 136.50  68.75   20.3 119.1 35.9  88.7  76.6  89.8
## 160   22.0 22.5  1.0477  31 177.25  71.50   24.4 138.3 36.2 101.1  92.4  99.3
## 161    9.9  9.4  1.0775  31 151.25  72.25   20.4 136.2 35.0  94.0  81.2  91.5
## 162   14.8 14.6  1.0653  33 196.00  73.00   25.9 167.0 38.5 103.8  95.6 105.1
## 163   13.3 13.0  1.0690  33 184.25  68.75   24.4 159.8 40.7  98.9  92.1 103.5
## 164   15.2 15.1  1.0644  34 140.00  70.50   19.8 118.8 36.0  89.2  83.4  89.6
## 165   26.5 27.3  1.0370  34 218.75  72.00   29.7 160.8 39.5 111.4 106.0 108.8
## 166   19.0 19.2  1.0549  35 217.00  73.75   28.1 175.8 40.5 107.5  95.1 104.5
## 167   21.4 21.8  1.0492  35 166.25  68.00   25.3 130.7 38.5  99.1  90.4  95.6
## 168   20.0 20.3  1.0525  35 224.75  72.25   30.3 179.7 43.9 108.2 100.4 106.8
## 169   34.7 34.3  1.0180  35 228.25  69.50   33.3 149.3 40.4 114.9 115.9 111.9
## 170   16.5 16.5  1.0610  35 172.75  69.50   25.2 144.2 37.6  99.1  90.8  98.1
## 171    4.1  3.0  1.0926  35 152.25  67.75   23.4 146.1 37.0  92.2  81.9  92.8
## 172    1.9  0.7  1.0983  35 125.75  65.50   20.6 123.4 34.0  90.8  75.0  89.2
## 173   20.2 20.5  1.0521  35 177.25  71.00   24.8 141.7 38.4 100.5  90.3  98.7
## 174   16.8 16.9  1.0603  36 176.25  71.50   24.3 146.6 38.7  98.2  90.3  99.9
## 175   24.6 25.3  1.0414  36 226.75  71.75   31.0 170.9 41.5 115.3 108.8 114.4
## 176   10.4  9.9  1.0763  37 145.25  69.25   21.3 130.2 36.0  96.8  79.4  89.2
## 177   13.4 13.1  1.0689  37 151.00  67.00   23.7 130.8 35.3  92.6  83.2  96.4
## 178   28.8 29.9  1.0316  37 241.25  71.50   33.2 171.7 42.1 119.2 110.3 113.9
## 179   22.0 22.5  1.0477  38 187.25  69.25   27.5 146.1 38.0 102.7  92.7 101.9
## 180   16.8 16.9  1.0603  39 234.75  74.50   29.8 195.3 42.8 109.5 104.5 109.9
## 181   25.8 26.6  1.0387  39 219.25  74.25   28.0 162.7 40.0 108.5 104.6 109.8
## 182    0.0  0.0  1.1089  40 118.50  68.00   18.1 118.5 33.8  79.3  69.4  85.0
## 183   11.9 11.5  1.0725  40 145.75  67.25   22.7 128.4 35.5  95.5  83.6  91.6
## 184   12.4 12.1  1.0713  40 159.25  69.75   23.0 139.5 35.3  92.3  86.8  96.1
## 185   17.4 17.5  1.0587  40 170.50  74.25   21.8 140.8 37.7  98.9  90.4  95.5
## 186    9.2  8.6  1.0794  40 167.50  71.50   23.1 152.1 39.4  89.5  83.7  98.1
## 187   23.0 23.6  1.0453  41 232.75  74.25   29.7 179.2 41.9 117.5 109.3 108.8
## 188   20.1 20.4  1.0524  41 210.50  72.00   28.6 168.3 38.5 107.4  98.9 104.1
## 189   20.2 20.5  1.0520  41 202.25  72.50   27.0 161.4 40.8 109.2  98.0 101.8
## 190   23.8 24.4  1.0434  41 185.00  68.25   28.0 141.0 38.0 103.4 101.2 103.1
## 191   11.8 11.4  1.0728  41 153.00  69.25   22.5 135.0 36.4  91.4  80.6  92.3
## 192   36.5 38.1  1.0140  42 244.25  76.00   29.8 155.2 41.8 115.2 113.7 112.4
## 193   16.0 15.9  1.0624  42 193.50  70.50   27.4 162.6 40.7 104.9  94.1 102.7
## 194   24.0 24.7  1.0429  42 224.75  74.75   28.3 170.8 38.5 106.7 105.7 111.8
## 195   22.3 22.8  1.0470  42 162.75  72.75   21.6 126.5 35.4  92.2  85.6  96.5
## 196   24.8 25.5  1.0411  42 180.00  68.25   27.2 135.4 38.5 101.6  96.6 100.6
## 197   21.5 22.0  1.0488  42 156.25  69.00   23.1 122.6 35.5  97.8  86.0  96.2
## 198   17.6 17.7  1.0583  42 168.00  71.50   23.1 138.4 36.5  92.0  89.7 101.0
## 199    7.3  6.6  1.0841  42 167.25  72.75   22.3 155.1 37.6  94.0  78.0  99.0
## 200   22.6 23.6  1.0462  43 170.75  67.50   26.4 132.1 37.4 103.7  89.7  94.2
## 201   12.5 12.2  1.0709  43 178.25  70.25   25.4 155.9 37.8 102.7  89.2  99.2
## 202   21.7 22.1  1.0484  43 150.00  69.25   22.0 117.5 35.2  91.1  85.7  96.9
## 203   27.7 28.7  1.0340  43 200.50  71.50   27.6 144.9 37.9 107.2 103.1 105.5
## 204    6.8  6.0  1.0854  44 184.00  74.00   23.7 171.4 37.9 100.8  89.1 102.6
## 205   33.4 34.8  1.0209  44 223.00  69.75   32.3 148.5 40.9 121.6 113.9 107.1
## 206   16.6 16.6  1.0610  44 208.75  73.00   27.6 174.2 41.9 105.6  96.3 102.0
## 207   31.7 32.9  1.0250  44 166.00  65.50   27.2 113.5 39.1 100.6  93.9 100.1
## 208   31.5 32.8  1.0254  47 195.00  72.50   26.1 133.6 40.2 102.7 101.3 101.7
## 209   10.1  9.6  1.0771  47 160.50  70.25   22.9 144.3 36.0  99.8  83.9  91.8
## 210   11.3 10.8  1.0742  47 159.75  70.75   22.5 141.8 34.5  92.9  84.4  94.0
## 211    7.8  7.1  1.0829  49 140.50  68.00   21.4 129.5 35.8  91.2  79.4  89.0
## 212   26.4 27.2  1.0373  49 216.25  74.50   27.4 159.3 40.2 115.6 104.0 109.0
## 213   19.3 19.5  1.0543  49 168.25  71.75   23.0 135.9 38.3  98.3  89.7  99.1
## 214   18.5 18.7  1.0561  50 194.75  70.75   27.4 158.7 39.0 103.7  97.6 104.2
## 215   19.3 19.5  1.0543  50 172.75  73.00   22.8 139.4 37.4  98.7  87.6  96.1
## 216   45.1 47.5  0.9950  51 219.00  64.00   37.6 120.2 41.2 119.8 122.1 112.8
## 217   13.8 13.6  1.0678  51 149.25  69.75   21.6 128.7 34.8  92.8  81.1  96.3
## 218    8.2  7.5  1.0819  51 154.50  70.00   22.2 141.9 36.9  93.3  81.5  94.4
## 219   23.9 24.5  1.0433  52 199.25  71.75   27.2 151.7 39.4 106.8 100.0 105.0
## 220   15.1 15.0  1.0646  53 154.50  69.25   22.7 131.2 37.6  93.9  88.7  94.5
## 221   12.7 12.4  1.0706  54 153.25  70.50   24.5 151.3 38.5  99.0  91.8  96.2
## 222   25.3 26.0  1.0399  54 230.00  72.25   31.0 171.9 42.5 119.9 110.4 105.5
## 223   11.9 11.5  1.0726  54 161.75  67.50   25.0 142.6 37.4  94.2  87.6  95.6
## 224    6.1  5.2  1.0874  55 142.25  67.25   22.2 133.6 35.2  92.7  82.8  91.9
## 225   11.3 10.9  1.0740  55 179.75  68.75   26.8 159.5 41.1 106.9  95.3  98.2
## 226   12.8 12.5  1.0703  55 126.50  66.75   20.0 110.3 33.4  88.8  78.2  87.5
## 227   14.9 14.8  1.0650  55 169.50  68.25   25.6 144.2 37.2 101.7  91.1  97.1
## 228   24.5 25.2  1.0418  55 198.50  74.25   25.3 149.9 38.3 105.3  96.7 106.6
## 229   15.0 14.9  1.0647  56 174.50  69.50   25.4 148.3 38.1 104.0  89.4  98.4
## 230   16.9 17.0  1.0601  56 167.75  68.50   25.2 139.4 37.4  98.6  93.0  97.0
## 231   11.1 10.6  1.0745  57 147.75  65.75   24.1 131.4 35.2  99.6  86.4  90.1
## 232   16.1 16.1  1.0620  57 182.25  71.75   24.9 152.9 39.4 103.4  96.7 100.7
## 233   15.5 15.4  1.0636  58 175.50  71.50   24.2 148.4 38.0 100.2  88.1  97.8
## 234   25.9 26.7  1.0384  58 161.75  67.25   25.2 119.9 35.1  94.9  94.9 100.2
## 235   25.5 25.8  1.0403  60 157.75  67.50   24.1 117.5 40.4  97.2  93.3  94.0
## 236   18.4 18.6  1.0563  62 168.75  67.50   26.1 137.6 38.3 104.7  95.6  93.7
## 237   24.0 24.8  1.0424  62 191.50  72.25   25.8 145.2 40.6 104.0  98.2 101.1
## 238   26.4 27.3  1.0372  63 219.15  69.50   31.9 161.2 40.2 117.6 113.8 111.8
## 239   12.7 12.4  1.0705  64 155.25  69.50   22.6 135.5 37.9  95.8  82.8  94.5
## 240   28.8 29.9  1.0316  65 189.75  65.75   30.9 135.1 40.8 106.4 100.5 100.5
## 241   17.0 17.0  1.0599  65 127.50  65.75   20.8 105.9 34.7  93.0  79.7  87.6
## 242   33.6 35.0  1.0207  65 224.50  68.25   33.9 149.2 38.8 119.6 118.0 114.3
## 243   29.3 30.4  1.0304  66 234.25  72.00   31.8 165.6 41.4 119.7 109.0 109.1
## 244   31.4 32.6  1.0256  67 227.75  72.75   30.3 156.3 41.3 115.8 113.4 109.8
## 245   28.1 29.0  1.0334  67 199.50  68.50   29.9 143.6 40.7 118.3 106.1 101.6
## 246   15.3 15.2  1.0641  68 155.50  69.25   22.8 131.8 36.3  97.4  84.3  94.4
## 247   29.1 30.2  1.0308  69 215.50  70.50   30.5 152.7 40.8 113.7 107.6 110.0
## 248   11.5 11.0  1.0736  70 134.25  67.00   21.1 118.9 34.9  89.2  83.6  88.8
## 249   32.3 33.6  1.0236  72 201.00  69.75   29.1 136.1 40.9 108.5 105.0 104.5
## 250   28.3 29.3  1.0328  72 186.75  66.00   30.2 133.9 38.9 111.1 111.5 101.7
## 251   25.3 26.0  1.0399  72 190.75  70.50   27.0 142.6 38.9 108.3 101.3  97.8
## 252   30.7 31.9  1.0271  74 207.50  70.00   29.8 143.7 40.8 112.4 108.5 107.1
##     thigh knee ankle biceps forearm wrist
## 1    59.0 37.3  21.9   32.0    27.4  17.1
## 2    58.7 37.3  23.4   30.5    28.9  18.2
## 3    59.6 38.9  24.0   28.8    25.2  16.6
## 4    60.1 37.3  22.8   32.4    29.4  18.2
## 5    63.2 42.2  24.0   32.2    27.7  17.7
## 6    66.0 42.0  25.6   35.7    30.6  18.8
## 7    58.4 38.3  22.9   31.9    27.8  17.7
## 8    60.0 39.4  23.2   30.5    29.0  18.8
## 9    62.9 38.3  23.8   35.9    31.1  18.2
## 10   63.1 41.7  25.0   35.6    30.0  19.2
## 11   59.7 39.7  25.2   32.8    29.4  18.5
## 12   66.2 39.2  25.9   37.2    30.2  19.0
## 13   63.4 38.3  21.5   32.5    28.6  17.7
## 14   66.0 41.5  23.7   36.9    31.6  18.8
## 15   69.0 39.0  23.1   36.1    30.5  18.2
## 16   63.1 38.7  21.7   31.1    26.4  16.9
## 17   64.8 40.8  23.1   36.2    30.8  17.3
## 18   66.9 40.0  24.4   38.2    31.6  19.3
## 19   64.2 38.7  22.9   37.2    30.5  18.5
## 20   65.8 40.6  24.0   37.1    30.1  18.2
## 21   63.5 38.0  22.1   32.5    30.3  18.4
## 22   63.4 40.6  24.6   33.0    32.8  19.9
## 23   57.4 35.3  22.2   27.9    25.9  16.7
## 24   54.9 36.2  22.1   29.8    26.7  17.1
## 25   58.4 35.5  22.9   31.1    28.0  17.6
## 26   55.0 36.7  22.5   29.9    28.2  17.7
## 27   51.7 34.7  21.4   28.7    27.0  16.5
## 28   57.5 36.0  21.0   29.2    26.6  17.0
## 29   50.1 34.5  21.3   30.5    27.9  17.2
## 30   58.9 35.3  22.6   30.1    26.7  17.6
## 31   57.5 38.7  33.9   32.5    27.7  18.4
## 32   58.5 38.8  21.5   30.1    26.4  17.9
## 33   57.3 36.2  24.5   29.0    30.0  18.8
## 34   67.1 44.2  25.2   37.5    31.5  18.7
## 35   71.2 43.3  26.3   37.3    31.7  19.7
## 36   61.9 38.3  21.9   32.0    29.8  17.0
## 37   63.5 39.9  22.6   35.1    30.6  19.0
## 38   66.0 41.5  24.7   33.2    30.5  19.4
## 39   87.3 49.1  29.6   45.0    29.0  21.4
## 40   61.3 41.1  24.7   34.1    31.0  18.3
## 41   72.5 39.6  26.6   36.4    32.7  21.4
## 42   70.6 42.5  23.7   33.6    28.7  17.4
## 43   67.7 40.9  25.0   36.7    29.8  18.4
## 44   65.0 40.2  23.0   35.8    31.5  18.8
## 45   50.0 34.7  21.0   26.1    23.1  16.1
## 46   58.4 38.2  23.4   29.7    27.4  18.3
## 47   53.3 34.5  22.5   27.9    26.2  17.3
## 48   52.7 37.5  21.9   28.8    26.8  17.9
## 49   52.0 35.8  20.6   28.8    25.5  16.3
## 50   50.6 34.4  21.9   26.8    25.8  16.8
## 51   52.6 37.2  22.4   26.0    25.8  17.3
## 52   51.4 34.9  21.0   26.7    26.1  17.2
## 53   49.3 33.7  21.4   29.6    26.0  16.9
## 54   52.6 37.6  22.6   38.5    27.4  18.5
## 55   51.9 34.9  22.5   27.7    27.5  18.5
## 56   57.2 38.0  22.0   35.9    30.2  18.9
## 57   62.1 39.6  22.5   33.1    28.3  18.5
## 58   61.8 39.8  22.7   37.7    30.9  19.2
## 59   59.1 38.0  22.5   31.6    28.8  18.2
## 60   60.6 37.7  22.9   34.5    29.6  18.5
## 61   61.6 40.9  23.1   36.2    31.8  20.2
## 62   60.9 38.0  22.1   32.5    29.8  18.3
## 63   61.0 39.4  23.6   32.7    29.9  19.1
## 64   60.8 40.1  22.7   33.6    29.0  18.8
## 65   65.0 41.2  24.7   35.3    31.1  18.4
## 66   64.8 40.2  22.7   34.8    30.1  18.7
## 67   55.0 36.1  21.7   29.6    27.4  17.4
## 68   54.3 35.4  21.5   32.8    27.4  18.7
## 69   56.0 37.4  22.4   32.6    28.1  18.1
## 70   51.2 37.4  21.6   27.3    27.1  17.3
## 71   56.6 38.6  22.4   31.5    27.3  18.6
## 72   58.9 37.6  21.6   30.3    27.3  18.3
## 73   52.9 37.5  23.1   29.7    27.3  18.2
## 74   50.9 35.4  19.1   29.3    25.7  16.9
## 75   55.5 35.2  20.9   29.4    27.0  16.8
## 76   54.5 37.0  21.4   29.3    27.0  18.3
## 77   56.7 39.7  24.2   30.2    29.2  18.1
## 78   55.0 38.3  21.8   30.8    25.7  18.8
## 79   53.5 37.5  21.5   31.4    26.8  18.3
## 80   56.5 39.3  22.7   30.3    28.7  19.0
## 81   54.8 38.2  23.7   29.4    27.2  19.0
## 82   54.8 38.0  22.0   29.9    25.2  17.7
## 83   58.9 39.0  23.0   34.3    29.6  19.0
## 84   56.0 36.5  24.1   31.2    27.3  19.2
## 85   56.3 36.6  22.0   29.7    26.3  18.0
## 86   54.7 37.8  33.7   32.4    27.7  18.2
## 87   57.2 37.7  21.8   32.6    28.0  18.8
## 88   57.8 39.5  23.3   29.2    28.4  18.1
## 89   61.0 38.4  23.8   30.2    29.3  18.8
## 90   60.4 39.9  24.4   28.8    29.6  18.7
## 91   58.3 38.2  22.5   29.1    27.7  17.7
## 92   58.9 39.7  23.1   31.4    28.4  18.8
## 93   56.9 38.3  22.1   30.1    28.2  18.4
## 94   58.9 40.5  24.5   33.3    29.6  19.1
## 95   57.4 39.6  24.6   30.3    27.9  17.8
## 96   61.7 42.3  23.2   32.9    30.8  20.4
## 97   60.6 39.4  22.9   31.6    30.1  18.5
## 98   62.1 39.3  23.3   30.6    27.8  18.2
## 99   54.7 37.3  21.9   31.6    27.5  18.2
## 100  62.7 39.0  22.3   35.3    30.9  18.3
## 101  59.0 39.4  22.3   32.2    31.0  18.6
## 102  59.3 38.4  22.4   27.9    26.2  17.0
## 103  60.0 38.4  23.2   31.0    29.2  18.4
## 104  59.1 39.8  25.4   31.0    30.3  19.7
## 105  59.5 36.1  22.0   30.1    27.2  17.7
## 106  54.7 39.0  24.8   31.0    29.4  18.8
## 107  61.2 39.3  23.5   30.5    28.5  18.1
## 108  55.8 38.7  23.4   35.1    29.6  19.1
## 109  57.5 40.0  24.8   35.1    30.7  19.2
## 110  57.5 36.8  22.8   32.1    26.0  17.3
## 111  61.9 38.0  22.3   33.3    28.2  18.1
## 112  63.7 40.0  23.6   33.5    27.8  17.4
## 113  62.3 38.1  23.9   35.3    31.1  19.8
## 114  61.5 37.8  21.9   30.7    27.6  17.4
## 115  59.3 38.1  21.8   31.8    27.3  17.5
## 116  55.9 36.3  22.1   29.8    26.3  17.3
## 117  58.8 38.4  22.8   29.9    28.0  18.1
## 118  56.8 38.8  23.3   33.4    29.8  19.5
## 119  64.6 41.1  24.8   33.6    29.5  18.5
## 120  58.5 39.2  24.5   32.1    28.6  18.0
## 121  58.5 39.3  24.6   33.9    31.2  19.5
## 122  57.1 40.5  23.2   33.0    29.6  18.4
## 123  60.5 38.7  22.6   34.4    28.0  17.6
## 124  58.1 36.5  22.1   30.6    27.5  17.6
## 125  58.5 36.6  23.5   34.4    29.2  18.0
## 126  60.7 36.0  21.9   35.6    30.2  17.6
## 127  60.7 36.8  22.2   33.8    30.3  17.2
## 128  53.5 35.8  20.8   33.9    28.2  17.4
## 129  61.7 39.0  21.8   33.3    29.6  18.1
## 130  57.4 36.9  22.2   31.6    27.8  17.7
## 131  57.0 38.7  23.2   27.5    26.5  17.6
## 132  60.3 38.5  23.0   31.2    28.4  17.1
## 133  61.2 38.1  22.6   33.5    28.6  17.9
## 134  56.1 35.6  20.5   33.6    29.3  17.3
## 135  56.0 36.9  23.0   34.0    29.8  18.1
## 136  58.9 37.9  22.7   30.9    28.8  17.6
## 137  58.8 36.1  22.4   32.7    28.3  17.1
## 138  63.6 39.2  23.8   34.3    28.4  17.7
## 139  58.1 38.4  22.5   31.7    27.4  17.6
## 140  66.5 42.5  24.5   35.5    29.8  18.7
## 141  59.1 39.6  21.6   30.8    27.9  16.6
## 142  60.4 38.2  22.0   32.0    28.5  17.8
## 143  57.1 36.7  22.3   31.6    27.5  17.9
## 144  56.1 36.1  22.7   30.5    27.2  18.2
## 145  59.1 37.6  23.2   31.8    29.7  18.3
## 146  56.4 36.5  22.0   33.5    28.3  17.3
## 147  71.2 43.5  25.2   36.1    30.3  18.7
## 148  68.4 40.8  24.6   33.3    29.7  18.4
## 149  51.9 35.7  22.0   25.8    25.2  16.9
## 150  67.6 42.7  24.7   36.0    30.4  18.4
## 151  56.9 35.9  20.4   31.6    29.0  17.8
## 152  72.9 43.5  25.1   38.5    33.8  19.6
## 153  53.6 36.8  23.8   27.8    26.3  17.4
## 154  56.8 37.4  22.8   30.6    28.3  17.9
## 155  62.1 40.0  24.9   33.7    29.2  19.4
## 156  57.3 37.8  21.7   32.2    27.7  17.7
## 157  68.6 40.0  25.2   35.2    30.7  19.1
## 158  60.8 38.5  25.0   31.6    28.0  18.6
## 159  50.1 34.8  21.8   27.0    34.9  16.9
## 160  59.4 39.0  24.6   30.1    28.2  18.2
## 161  52.5 36.6  21.0   27.0    26.3  16.5
## 162  61.4 40.6  25.0   31.3    29.2  19.1
## 163  64.0 37.3  23.5   33.5    30.6  19.7
## 164  52.4 35.6  20.4   28.3    26.2  16.5
## 165  63.8 42.0  23.4   34.0    31.2  18.5
## 166  64.8 41.3  25.6   36.4    33.7  19.4
## 167  55.5 34.2  21.9   30.2    28.7  17.7
## 168  63.3 41.7  24.6   37.2    33.1  19.8
## 169  74.4 40.6  24.0   36.1    31.8  18.8
## 170  60.1 39.1  23.4   32.5    29.8  17.4
## 171  54.7 36.2  22.1   30.4    27.4  17.7
## 172  50.0 34.8  22.0   24.8    25.9  16.9
## 173  57.8 37.3  22.4   31.0    28.7  17.7
## 174  59.2 37.7  21.5   32.4    28.4  17.8
## 175  69.2 42.4  24.0   35.4    21.0  20.1
## 176  50.3 34.8  22.2   31.0    26.9  16.9
## 177  60.0 38.1  22.0   31.5    26.6  16.7
## 178  69.8 42.6  24.8   34.4    29.5  18.4
## 179  64.7 39.5  24.7   34.8    30.3  18.1
## 180  69.5 43.1  25.8   39.1    32.5  19.9
## 181  68.1 42.8  24.1   35.6    29.0  19.0
## 182  47.2 33.5  20.2   27.7    24.6  16.5
## 183  54.1 36.2  21.8   31.4    28.3  17.2
## 184  58.0 39.4  22.7   30.0    26.4  17.4
## 185  55.4 38.9  22.4   30.5    28.9  17.7
## 186  57.3 39.7  22.6   32.9    29.3  18.2
## 187  67.7 41.3  24.7   37.2    31.8  20.0
## 188  63.5 39.8  23.5   36.4    30.4  19.1
## 189  62.8 41.3  24.8   36.6    32.4  18.8
## 190  61.5 40.4  22.9   33.4    29.2  18.5
## 191  54.3 36.3  21.8   29.6    27.3  17.9
## 192  68.5 45.0  25.5   37.1    31.2  19.9
## 193  60.6 38.6  24.7   34.0    30.1  18.7
## 194  65.3 43.3  26.0   33.7    29.9  18.5
## 195  60.2 38.9  22.4   31.7    27.1  17.1
## 196  61.1 38.4  24.1   32.9    29.8  18.8
## 197  57.7 38.6  24.0   31.2    27.3  17.4
## 198  62.3 38.0  22.3   30.8    27.8  16.9
## 199  57.5 40.0  22.5   30.6    30.0  18.5
## 200  58.5 39.0  24.1   33.8    28.8  18.8
## 201  60.2 39.2  23.8   31.7    28.4  18.6
## 202  55.5 35.7  22.0   29.4    26.6  17.4
## 203  68.8 38.3  23.7   32.1    28.9  18.7
## 204  60.6 39.0  24.0   32.9    29.2  18.4
## 205  63.5 40.3  21.8   34.8    30.7  17.4
## 206  63.3 39.8  24.1   37.3    23.1  19.4
## 207  58.9 37.6  21.4   33.1    29.5  17.3
## 208  60.7 39.4  23.3   36.7    31.6  18.4
## 209  53.0 36.2  22.5   31.4    27.5  17.7
## 210  56.0 38.2  22.6   29.0    26.2  17.6
## 211  51.1 35.0  21.7   30.9    28.8  17.4
## 212  63.7 40.3  23.2   36.8    31.0  18.9
## 213  56.3 38.8  23.0   29.5    27.9  18.6
## 214  60.0 40.9  25.5   32.7    30.0  19.0
## 215  57.1 38.1  21.8   28.6    26.7  18.0
## 216  62.5 36.9  23.6   34.7    29.1  18.4
## 217  53.8 36.5  21.5   31.3    26.3  17.8
## 218  54.7 39.0  22.6   27.5    25.9  18.6
## 219  63.9 39.2  22.9   35.7    30.4  19.2
## 220  53.7 36.2  22.0   28.5    25.7  17.1
## 221  57.7 38.1  23.9   31.4    29.9  18.9
## 222  64.2 42.7  27.0   38.4    32.0  19.6
## 223  59.7 40.2  23.4   27.9    27.0  17.8
## 224  54.4 35.2  22.5   29.4    26.8  17.0
## 225  57.4 37.1  21.8   34.1    31.1  19.2
## 226  50.8 33.0  19.7   25.3    22.0  15.8
## 227  56.6 38.5  22.6   33.4    29.3  18.8
## 228  64.0 42.6  23.4   33.2    30.0  18.4
## 229  58.4 37.4  22.5   34.6    30.1  18.8
## 230  55.4 38.8  23.2   32.4    29.7  19.0
## 231  53.0 35.0  21.3   31.7    27.3  16.9
## 232  59.3 38.6  22.8   31.8    29.1  19.0
## 233  57.1 38.9  23.6   30.9    29.6  18.0
## 234  56.8 35.9  21.0   27.8    26.1  17.6
## 235  54.3 35.7  21.0   31.3    28.7  18.3
## 236  54.4 37.1  22.7   30.3    26.3  18.3
## 237  59.3 40.3  23.0   32.6    28.5  19.0
## 238  63.4 41.1  22.3   35.1    29.6  18.5
## 239  61.2 39.1  22.3   29.8    28.9  18.3
## 240  59.2 38.1  24.0   35.9    30.5  19.1
## 241  50.7 33.4  20.1   28.5    24.8  16.5
## 242  61.3 42.1  23.4   34.9    30.1  19.4
## 243  63.7 42.4  24.6   35.6    30.7  19.5
## 244  65.6 46.0  25.4   35.3    29.8  19.5
## 245  58.2 38.8  24.1   32.1    29.3  18.5
## 246  54.3 37.5  22.6   29.2    27.3  18.5
## 247  63.3 44.0  22.6   37.5    32.6  18.8
## 248  49.6 34.8  21.5   25.6    25.7  18.5
## 249  59.6 40.8  23.2   35.2    28.6  20.1
## 250  60.3 37.3  21.5   31.3    27.2  18.0
## 251  56.0 41.6  22.7   30.5    29.4  19.8
## 252  59.3 42.2  24.6   33.7    30.0  20.9
testdata
##  [1]  10  20  30  40  50  60  70  80  90 100 110 120 130 140 150 160 170 180 190
## [20] 200 210 220 230 240 250
train = fat[-testdata,]
subset=fat[testdata,] 
train 
##     brozek siri density age weight height adipos  free neck chest abdom   hip
## 1     12.6 12.3  1.0708  23 154.25  67.75   23.7 134.9 36.2  93.1  85.2  94.5
## 2      6.9  6.1  1.0853  22 173.25  72.25   23.4 161.3 38.5  93.6  83.0  98.7
## 3     24.6 25.3  1.0414  22 154.00  66.25   24.7 116.0 34.0  95.8  87.9  99.2
## 4     10.9 10.4  1.0751  26 184.75  72.25   24.9 164.7 37.4 101.8  86.4 101.2
## 5     27.8 28.7  1.0340  24 184.25  71.25   25.6 133.1 34.4  97.3 100.0 101.9
## 6     20.6 20.9  1.0502  24 210.25  74.75   26.5 167.0 39.0 104.5  94.4 107.8
## 7     19.0 19.2  1.0549  26 181.00  69.75   26.2 146.6 36.4 105.1  90.7 100.3
## 8     12.8 12.4  1.0704  25 176.00  72.50   23.6 153.6 37.8  99.6  88.5  97.1
## 9      5.1  4.1  1.0900  25 191.00  74.00   24.6 181.3 38.1 100.9  82.5  99.9
## 11     7.5  7.1  1.0830  26 186.25  74.50   23.6 172.3 38.5 101.5  83.6  98.2
## 12     8.5  7.8  1.0812  27 216.00  76.00   26.3 197.7 39.4 103.6  90.9 107.7
## 13    20.5 20.8  1.0513  32 180.50  69.50   26.3 143.5 38.4 102.0  91.6 103.9
## 14    20.8 21.2  1.0505  30 205.25  71.25   28.5 162.5 39.4 104.1 101.8 108.6
## 15    21.7 22.1  1.0484  35 187.75  69.50   27.4 147.0 40.5 101.3  96.4 100.1
## 16    20.5 20.9  1.0512  35 162.75  66.00   26.3 129.3 36.4  99.1  92.8  99.2
## 17    28.1 29.0  1.0333  34 195.75  71.00   27.3 140.8 38.9 101.9  96.4 105.2
## 18    22.4 22.9  1.0468  32 209.25  71.00   29.2 162.5 42.1 107.6  97.5 107.0
## 19    16.1 16.0  1.0622  28 183.75  67.75   28.2 154.3 38.0 106.8  89.6 102.4
## 21    19.0 19.1  1.0551  28 179.00  68.00   27.3 145.1 39.1 103.3  95.9 104.9
## 22    15.3 15.2  1.0640  28 200.50  69.75   29.1 169.8 41.3 111.4  98.8 104.8
## 23    15.7 15.6  1.0631  31 140.25  68.25   21.2 118.2 33.9  86.0  76.4  94.6
## 24    17.6 17.7  1.0584  32 148.75  70.00   21.4 122.6 35.5  86.7  80.0  93.4
## 25    14.2 14.0  1.0668  28 151.25  67.75   23.2 129.8 34.5  90.2  76.3  95.8
## 26     4.6  3.7  1.0911  27 159.25  71.50   21.9 151.9 35.7  89.6  79.7  96.5
## 27     8.5  7.9  1.0811  34 131.50  67.50   20.3 120.3 36.2  88.6  74.6  85.3
## 28    22.4 22.9  1.0468  31 148.00  67.50   22.9 114.9 38.8  97.4  88.7  94.7
## 29     4.7  3.7  1.0910  27 133.25  64.75   22.4 127.0 36.4  93.5  73.9  88.5
## 31    12.3 11.9  1.0716  32 182.00  73.75   23.6 159.7 38.7 100.5  88.7  99.8
## 32     6.5  5.7  1.0862  29 160.25  71.25   22.2 149.8 37.3  93.5  84.5 100.6
## 33    13.4 11.8  1.0719  27 168.00  71.25   23.3 142.5 38.1  93.0  79.1  94.5
## 34    20.9 21.3  1.0502  41 218.50  71.00   30.5 172.7 39.8 111.7 100.5 108.3
## 35    31.1 32.3  1.0263  41 247.25  73.50   32.2 170.4 42.1 117.0 115.6 116.1
## 36    38.2 40.1  1.0101  49 191.75  65.00   32.0 118.4 38.4 118.5 113.1 113.8
## 37    23.6 24.2  1.0438  40 202.25  70.00   29.1 154.5 38.5 106.5 100.9 106.2
## 38    27.5 28.4  1.0346  50 196.75  68.25   29.7 142.6 42.1 105.6  98.8 104.8
## 39    33.8 35.2  1.0202  46 363.15  72.25   48.9 240.5 51.2 136.2 148.1 147.7
## 41    33.1 34.5  1.0217  45 262.75  68.75   39.1 175.8 43.2 128.3 126.2 125.6
## 42    31.7 32.9  1.0250  44 205.00  29.50   29.9 140.1 36.6 106.0 104.3 115.5
## 43    30.4 31.6  1.0279  48 217.00  70.00   31.2 151.1 37.3 113.3 111.2 114.1
## 44    30.8 32.0  1.0269  41 212.00  71.50   29.2 146.7 41.5 106.6 104.3 106.0
## 45     8.4  7.7  1.0814  39 125.25  68.00   19.1 114.7 31.5  85.1  76.0  88.2
## 46    14.1 13.9  1.0670  43 164.25  73.25   21.3 141.1 35.7  96.6  81.5  97.2
## 47    11.2 10.8  1.0742  40 133.50  67.50   20.6 118.5 33.6  88.2  73.7  88.5
## 48     6.4  5.6  1.0665  39 148.50  71.25   20.6 139.0 34.6  89.8  79.5  92.7
## 49    13.4 13.6  1.0678  45 135.75  68.50   20.4 117.6 32.8  92.3  83.4  90.4
## 51    10.7 10.2  1.0756  47 158.25  72.25   21.3 141.4 34.9  90.2  86.7  98.3
## 52     7.4  6.6  1.0840  40 139.25  69.00   20.6 129.0 34.3  89.2  77.9  91.0
## 53     8.7  8.0  1.0807  51 137.25  67.75   21.1 125.3 36.5  89.7  82.0  89.1
## 54     7.1  6.3  1.0848  49 152.75  73.50   19.9 142.0 35.1  93.3  79.6  91.6
## 55     4.9  3.9  1.0906  42 136.25  67.50   21.1 129.6 37.8  87.6  77.6  88.6
## 56    22.2 22.6  1.0473  54 198.00  72.00   26.9 154.1 39.9 107.6 100.0  99.6
## 57    20.1 20.4  1.0524  58 181.50  68.00   27.6 145.1 39.1 100.0  99.8 102.5
## 58    27.1 28.0  1.0356  62 201.25  69.50   29.3 146.7 40.5 111.5 104.2 105.8
## 59    30.4 31.5  1.0280  54 202.50  70.75   28.4 141.0 40.5 115.4 105.3  97.0
## 61    25.4 26.1  1.0396  62 216.00  73.25   28.2 161.2 41.4 112.3 104.8 103.1
## 62    28.8 29.8  1.0317  56 178.75  68.50   26.8 127.4 35.6 102.9  94.7 100.8
## 63    29.6 30.7  1.0298  54 193.25  70.25   27.6 136.1 38.0 107.6 102.4  99.4
## 64    25.1 25.8  1.0403  61 178.00  67.00   27.9 133.3 37.4 105.3  99.7  99.7
## 65    31.0 32.3  1.0264  57 205.50  70.00   29.5 141.7 40.1 105.3 105.5 108.3
## 66    28.9 30.0  1.0313  55 183.50  67.50   28.3 130.4 40.9 103.0 100.3 104.2
## 67    21.1 21.5  1.0499  54 151.50  70.75   21.3 119.6 35.6  90.0  83.9  93.9
## 68    14.0 13.8  1.0673  55 154.75  71.50   21.3 133.1 36.9  95.4  86.6  91.8
## 69     7.1  6.3  1.0847  54 155.25  69.25   22.8 144.2 37.5  89.3  78.4  96.1
## 71    23.7 24.3  1.0439  62 167.50  71.50   23.1 127.8 35.5  97.6  91.5  98.5
## 72     9.4  8.8  1.0788  55 146.75  68.75   21.9 132.9 38.7  88.5  82.8  95.5
## 73     9.1  8.5  1.0796  56 160.75  73.75   20.8 146.1 36.4  93.6  82.9  96.3
## 74    13.7 13.5  1.0680  55 125.00  64.00   21.5 107.9 33.2  87.7  76.0  88.6
## 75    12.0 11.8  1.0720  61 143.00  65.75   23.3 125.9 36.5  93.4  83.3  93.0
## 76    18.3 18.5  1.0666  61 148.25  67.50   22.9 121.1 36.0  91.6  81.8  94.8
## 77     9.2  8.8  1.0790  57 162.50  69.50   23.7 147.5 38.7  91.6  78.8  94.3
## 78    21.7 22.2  1.0483  69 177.75  68.50   26.7 139.1 38.7 102.0  95.0  98.3
## 79    21.1 21.5  1.0498  81 161.25  70.25   23.0 127.2 37.8  96.4  95.4  99.3
## 81    30.2 31.4  1.0283  67 163.75  67.75   25.1 114.3 38.4  97.7  95.8  97.1
## 82    26.0 26.8  1.0382  64 150.25  67.25   23.4 111.2 38.1  97.1  89.0  96.9
## 83    18.2 18.4  1.0568  64 190.25  72.75   25.3 155.6 39.3 103.1  97.8  99.6
## 84    26.2 27.0  1.0377  70 170.75  70.00   24.5 126.0 38.7 101.8  94.9  95.0
## 85    26.1 27.0  1.0378  72 168.00  69.25   24.7 124.1 38.5 101.4  99.8  96.2
## 86    25.8 26.6  1.0386  67 167.00  67.50   26.0 123.9 36.5  98.9  89.7  96.2
## 87    15.0 14.9  1.0648  72 157.75  67.25   24.6 134.1 37.7  97.5  88.1  96.9
## 88    22.6 23.1  1.0462  64 160.00  65.75   26.0 123.8 36.5 104.3  90.9  93.8
## 89     8.8  8.3  1.0800  46 176.75  72.50   23.7 161.1 38.0  97.3  86.0  99.3
## 91    20.2 20.5  1.0520  46 177.00  70.00   25.4 141.3 37.2  99.7  95.6 102.2
## 92    18.1 18.2  1.0573  44 179.75  69.50   26.2 147.3 39.2 101.9  93.2 100.6
## 93     9.2  8.5  1.0795  47 165.25  70.50   23.4 150.1 37.5  97.2  83.1  95.4
## 94    24.2 24.9  1.0424  46 192.50  71.75   26.3 145.9 38.0 106.6  97.5 100.6
## 95     9.6  9.0  1.0785  47 184.25  74.50   23.4 166.6 37.3  99.6  88.8 101.4
## 96    17.3 17.4  1.0991  53 224.50  77.75   26.1 185.7 41.1 113.2  99.2 107.5
## 97    10.1  9.6  1.0770  38 188.75  73.25   24.8 169.6 37.5  99.1  91.6 102.4
## 98    11.1 11.3  1.0730  50 162.50  66.50   25.9 143.5 38.7  99.4  86.7  96.2
## 99    17.7 17.8  1.0582  46 156.50  68.25   23.7 128.8 35.9  95.1  88.2  92.8
## 101   20.8 21.2  1.0506  49 198.50  73.50   25.9 157.2 40.1 106.5  95.0 101.7
## 102   20.1 20.4  1.0524  48 173.75  72.00   23.6 138.9 37.0  99.1  92.0  98.3
## 103   19.8 20.1  1.0530  41 172.75  71.25   24.0 138.6 36.3  96.7  89.2  98.3
## 104   21.9 22.3  1.0480  49 196.75  73.75   25.5 153.7 40.7 103.5  95.5 101.6
## 105   24.7 25.4  1.0412  43 177.00  69.25   26.0 133.2 39.6 104.0  98.6  99.5
## 106   17.8 18.0  1.0578  43 165.50  68.50   24.8 136.0 31.1  93.1  87.3  96.6
## 107   19.1 19.3  1.0547  43 200.25  73.50   26.0 162.0 38.6 105.2 102.8 103.6
## 108   18.2 18.3  1.0569  52 203.25  74.25   26.0 166.3 42.0 110.0 101.6 100.7
## 109   17.2 17.3  1.0593  43 194.00  75.50   24.0 160.6 38.5 110.1  88.7 102.1
## 111   19.5 19.7  1.0538  43 170.75  68.50   25.6 137.5 37.2  96.3  90.6  99.3
## 112   27.1 28.0  1.0355  43 183.25  70.00   26.3 133.5 37.1 108.0 105.0 103.0
## 113   21.6 22.1  1.0486  47 178.25  70.00   25.6 139.7 40.2  99.7  95.0  98.6
## 114   20.9 21.3  1.0503  42 163.00  70.25   23.3 128.9 35.3  93.5  89.6  99.8
## 115   25.9 26.7  1.0384  48 175.25  71.75   24.0 129.9 38.0 100.7  92.4  97.5
## 116   16.7 16.7  1.0607  40 158.00  69.25   23.4 131.7 36.3  97.0  86.6  92.6
## 117   19.8 20.1  1.0529  48 177.25  72.75   23.6 142.1 36.8  96.0  90.0  99.7
## 118   14.1 13.9  1.0671  51 179.00  72.00   24.3 153.8 41.0  99.2  90.0  96.4
## 119   25.1 25.8  1.0404  40 191.00  74.00   24.6 143.1 38.3  95.4  92.4 104.3
## 121   27.0 27.9  1.0358  52 206.50  74.50   26.2 150.7 40.8 104.3  99.2 104.1
## 122   24.6 25.3  1.0414  44 185.25  71.50   25.5 139.6 39.5  99.2  98.1 101.4
## 123   14.8 14.7  1.0652  40 160.25  68.75   23.9 136.5 36.9  99.3  83.3  97.5
## 124   16.0 16.0  1.0623  47 151.50  66.75   23.9 127.3 36.9  94.0  86.1  95.2
## 125   14.0 13.8  1.0674  50 161.00  66.50   25.6 138.5 37.7  98.9  84.1  94.0
## 126   17.4 17.5  1.0587  46 167.00  67.00   26.2 137.9 36.6 101.0  89.9 100.0
## 127   26.4 27.2  1.0373  42 177.50  68.75   26.4 130.7 38.9  98.7  92.1  98.5
## 128   17.4 17.4  1.0590  43 152.25  67.75   23.4 125.8 37.5  95.9  78.0  93.2
## 129   20.4 20.8  1.0515  40 192.25  73.25   25.2 153.0 39.8 103.9  93.5  99.5
## 131   18.0 18.1  1.0575  49 171.75  71.50   23.7 140.9 35.5  97.8  90.1  95.8
## 132   22.2 22.7  1.0472  40 171.25  70.50   24.3 133.3 36.3  94.6  90.3  99.1
## 133   23.1 23.6  1.0452  47 197.00  73.25   25.8 151.2 37.8 103.6  99.8 103.2
## 134   25.3 26.1  1.0398  50 157.00  66.75   24.8 117.2 37.8 100.4  89.4  92.3
## 135   23.8 24.4  1.0435  41 168.25  69.50   24.5 128.3 36.5  98.4  87.2  98.4
## 136   26.3 27.1  1.0374  44 186.00  69.75   26.8 137.1 37.8 104.6 101.1 102.1
## 137   21.4 21.8  1.0491  39 166.75  70.75   23.5 131.0 37.0  92.9  86.1  95.6
## 138   28.4 29.4  1.0325  43 187.75  74.00   24.1 134.4 37.7  97.8  98.6 100.6
## 139   21.8 22.4  1.0481  40 168.25  71.25   23.3 131.6 34.3  98.3  88.5  98.3
## 141   24.3 24.9  1.0422  40 176.75  71.00   24.6 133.8 37.4  98.6  93.1 101.6
## 142   18.1 18.3  1.0571  40 173.25  69.50   25.3 141.8 36.5  99.5  93.0  99.3
## 143   22.7 23.3  1.0459  52 167.00  67.75   25.6 129.0 37.5 102.7  91.0  98.9
## 144    9.9  9.4  1.0775  23 159.75  72.25   21.6 143.9 35.5  92.1  77.1  93.9
## 145   10.8 10.3  1.0754  23 188.15  77.50   22.1 168.4 38.0  96.6  85.3 102.5
## 146   14.4 14.2  1.0664  24 156.00  70.75   21.9 133.6 35.7  92.7  81.9  95.3
## 147   19.0 19.2  1.0550  24 208.50  72.75   27.7 168.9 39.2 102.0  99.1 110.1
## 148   28.6 29.6  1.0322  25 206.50  69.75   29.8 147.5 40.9 110.9 100.5 106.2
## 149    6.1  5.3  1.0873  25 143.75  72.50   19.3 135.0 35.2  92.3  76.5  92.1
## 151    9.9  9.4  1.0776  26 152.25  69.00   22.5 137.2 35.4  92.9  77.6  93.5
## 152   19.1 19.6  1.0542  26 241.75  74.50   30.7 195.1 41.8 108.3 102.9 114.4
## 153   10.6 10.1  1.0758  27 146.00  72.25   19.7 130.5 34.1  88.5  72.8  91.1
## 154   16.5 16.5  1.0610  27 156.75  67.25   24.4 130.9 37.9  94.0  88.2  95.2
## 155   20.5 21.0  1.0510  27 200.25  73.50   26.1 159.3 38.2 101.1 100.1 105.0
## 156   17.2 17.3  1.0594  28 171.50  75.25   21.6 142.0 35.6  92.1  83.5  98.3
## 157   30.1 31.2  1.0287  28 205.75  69.00   30.4 143.9 38.5 105.6 105.0 106.4
## 158   10.5 10.0  1.0761  28 182.50  72.25   24.6 163.4 37.0  98.5  90.8 102.5
## 159   12.8 12.5  1.0704  30 136.50  68.75   20.3 119.1 35.9  88.7  76.6  89.8
## 161    9.9  9.4  1.0775  31 151.25  72.25   20.4 136.2 35.0  94.0  81.2  91.5
## 162   14.8 14.6  1.0653  33 196.00  73.00   25.9 167.0 38.5 103.8  95.6 105.1
## 163   13.3 13.0  1.0690  33 184.25  68.75   24.4 159.8 40.7  98.9  92.1 103.5
## 164   15.2 15.1  1.0644  34 140.00  70.50   19.8 118.8 36.0  89.2  83.4  89.6
## 165   26.5 27.3  1.0370  34 218.75  72.00   29.7 160.8 39.5 111.4 106.0 108.8
## 166   19.0 19.2  1.0549  35 217.00  73.75   28.1 175.8 40.5 107.5  95.1 104.5
## 167   21.4 21.8  1.0492  35 166.25  68.00   25.3 130.7 38.5  99.1  90.4  95.6
## 168   20.0 20.3  1.0525  35 224.75  72.25   30.3 179.7 43.9 108.2 100.4 106.8
## 169   34.7 34.3  1.0180  35 228.25  69.50   33.3 149.3 40.4 114.9 115.9 111.9
## 171    4.1  3.0  1.0926  35 152.25  67.75   23.4 146.1 37.0  92.2  81.9  92.8
## 172    1.9  0.7  1.0983  35 125.75  65.50   20.6 123.4 34.0  90.8  75.0  89.2
## 173   20.2 20.5  1.0521  35 177.25  71.00   24.8 141.7 38.4 100.5  90.3  98.7
## 174   16.8 16.9  1.0603  36 176.25  71.50   24.3 146.6 38.7  98.2  90.3  99.9
## 175   24.6 25.3  1.0414  36 226.75  71.75   31.0 170.9 41.5 115.3 108.8 114.4
## 176   10.4  9.9  1.0763  37 145.25  69.25   21.3 130.2 36.0  96.8  79.4  89.2
## 177   13.4 13.1  1.0689  37 151.00  67.00   23.7 130.8 35.3  92.6  83.2  96.4
## 178   28.8 29.9  1.0316  37 241.25  71.50   33.2 171.7 42.1 119.2 110.3 113.9
## 179   22.0 22.5  1.0477  38 187.25  69.25   27.5 146.1 38.0 102.7  92.7 101.9
## 181   25.8 26.6  1.0387  39 219.25  74.25   28.0 162.7 40.0 108.5 104.6 109.8
## 182    0.0  0.0  1.1089  40 118.50  68.00   18.1 118.5 33.8  79.3  69.4  85.0
## 183   11.9 11.5  1.0725  40 145.75  67.25   22.7 128.4 35.5  95.5  83.6  91.6
## 184   12.4 12.1  1.0713  40 159.25  69.75   23.0 139.5 35.3  92.3  86.8  96.1
## 185   17.4 17.5  1.0587  40 170.50  74.25   21.8 140.8 37.7  98.9  90.4  95.5
## 186    9.2  8.6  1.0794  40 167.50  71.50   23.1 152.1 39.4  89.5  83.7  98.1
## 187   23.0 23.6  1.0453  41 232.75  74.25   29.7 179.2 41.9 117.5 109.3 108.8
## 188   20.1 20.4  1.0524  41 210.50  72.00   28.6 168.3 38.5 107.4  98.9 104.1
## 189   20.2 20.5  1.0520  41 202.25  72.50   27.0 161.4 40.8 109.2  98.0 101.8
## 191   11.8 11.4  1.0728  41 153.00  69.25   22.5 135.0 36.4  91.4  80.6  92.3
## 192   36.5 38.1  1.0140  42 244.25  76.00   29.8 155.2 41.8 115.2 113.7 112.4
## 193   16.0 15.9  1.0624  42 193.50  70.50   27.4 162.6 40.7 104.9  94.1 102.7
## 194   24.0 24.7  1.0429  42 224.75  74.75   28.3 170.8 38.5 106.7 105.7 111.8
## 195   22.3 22.8  1.0470  42 162.75  72.75   21.6 126.5 35.4  92.2  85.6  96.5
## 196   24.8 25.5  1.0411  42 180.00  68.25   27.2 135.4 38.5 101.6  96.6 100.6
## 197   21.5 22.0  1.0488  42 156.25  69.00   23.1 122.6 35.5  97.8  86.0  96.2
## 198   17.6 17.7  1.0583  42 168.00  71.50   23.1 138.4 36.5  92.0  89.7 101.0
## 199    7.3  6.6  1.0841  42 167.25  72.75   22.3 155.1 37.6  94.0  78.0  99.0
## 201   12.5 12.2  1.0709  43 178.25  70.25   25.4 155.9 37.8 102.7  89.2  99.2
## 202   21.7 22.1  1.0484  43 150.00  69.25   22.0 117.5 35.2  91.1  85.7  96.9
## 203   27.7 28.7  1.0340  43 200.50  71.50   27.6 144.9 37.9 107.2 103.1 105.5
## 204    6.8  6.0  1.0854  44 184.00  74.00   23.7 171.4 37.9 100.8  89.1 102.6
## 205   33.4 34.8  1.0209  44 223.00  69.75   32.3 148.5 40.9 121.6 113.9 107.1
## 206   16.6 16.6  1.0610  44 208.75  73.00   27.6 174.2 41.9 105.6  96.3 102.0
## 207   31.7 32.9  1.0250  44 166.00  65.50   27.2 113.5 39.1 100.6  93.9 100.1
## 208   31.5 32.8  1.0254  47 195.00  72.50   26.1 133.6 40.2 102.7 101.3 101.7
## 209   10.1  9.6  1.0771  47 160.50  70.25   22.9 144.3 36.0  99.8  83.9  91.8
## 211    7.8  7.1  1.0829  49 140.50  68.00   21.4 129.5 35.8  91.2  79.4  89.0
## 212   26.4 27.2  1.0373  49 216.25  74.50   27.4 159.3 40.2 115.6 104.0 109.0
## 213   19.3 19.5  1.0543  49 168.25  71.75   23.0 135.9 38.3  98.3  89.7  99.1
## 214   18.5 18.7  1.0561  50 194.75  70.75   27.4 158.7 39.0 103.7  97.6 104.2
## 215   19.3 19.5  1.0543  50 172.75  73.00   22.8 139.4 37.4  98.7  87.6  96.1
## 216   45.1 47.5  0.9950  51 219.00  64.00   37.6 120.2 41.2 119.8 122.1 112.8
## 217   13.8 13.6  1.0678  51 149.25  69.75   21.6 128.7 34.8  92.8  81.1  96.3
## 218    8.2  7.5  1.0819  51 154.50  70.00   22.2 141.9 36.9  93.3  81.5  94.4
## 219   23.9 24.5  1.0433  52 199.25  71.75   27.2 151.7 39.4 106.8 100.0 105.0
## 221   12.7 12.4  1.0706  54 153.25  70.50   24.5 151.3 38.5  99.0  91.8  96.2
## 222   25.3 26.0  1.0399  54 230.00  72.25   31.0 171.9 42.5 119.9 110.4 105.5
## 223   11.9 11.5  1.0726  54 161.75  67.50   25.0 142.6 37.4  94.2  87.6  95.6
## 224    6.1  5.2  1.0874  55 142.25  67.25   22.2 133.6 35.2  92.7  82.8  91.9
## 225   11.3 10.9  1.0740  55 179.75  68.75   26.8 159.5 41.1 106.9  95.3  98.2
## 226   12.8 12.5  1.0703  55 126.50  66.75   20.0 110.3 33.4  88.8  78.2  87.5
## 227   14.9 14.8  1.0650  55 169.50  68.25   25.6 144.2 37.2 101.7  91.1  97.1
## 228   24.5 25.2  1.0418  55 198.50  74.25   25.3 149.9 38.3 105.3  96.7 106.6
## 229   15.0 14.9  1.0647  56 174.50  69.50   25.4 148.3 38.1 104.0  89.4  98.4
## 231   11.1 10.6  1.0745  57 147.75  65.75   24.1 131.4 35.2  99.6  86.4  90.1
## 232   16.1 16.1  1.0620  57 182.25  71.75   24.9 152.9 39.4 103.4  96.7 100.7
## 233   15.5 15.4  1.0636  58 175.50  71.50   24.2 148.4 38.0 100.2  88.1  97.8
## 234   25.9 26.7  1.0384  58 161.75  67.25   25.2 119.9 35.1  94.9  94.9 100.2
## 235   25.5 25.8  1.0403  60 157.75  67.50   24.1 117.5 40.4  97.2  93.3  94.0
## 236   18.4 18.6  1.0563  62 168.75  67.50   26.1 137.6 38.3 104.7  95.6  93.7
## 237   24.0 24.8  1.0424  62 191.50  72.25   25.8 145.2 40.6 104.0  98.2 101.1
## 238   26.4 27.3  1.0372  63 219.15  69.50   31.9 161.2 40.2 117.6 113.8 111.8
## 239   12.7 12.4  1.0705  64 155.25  69.50   22.6 135.5 37.9  95.8  82.8  94.5
## 241   17.0 17.0  1.0599  65 127.50  65.75   20.8 105.9 34.7  93.0  79.7  87.6
## 242   33.6 35.0  1.0207  65 224.50  68.25   33.9 149.2 38.8 119.6 118.0 114.3
## 243   29.3 30.4  1.0304  66 234.25  72.00   31.8 165.6 41.4 119.7 109.0 109.1
## 244   31.4 32.6  1.0256  67 227.75  72.75   30.3 156.3 41.3 115.8 113.4 109.8
## 245   28.1 29.0  1.0334  67 199.50  68.50   29.9 143.6 40.7 118.3 106.1 101.6
## 246   15.3 15.2  1.0641  68 155.50  69.25   22.8 131.8 36.3  97.4  84.3  94.4
## 247   29.1 30.2  1.0308  69 215.50  70.50   30.5 152.7 40.8 113.7 107.6 110.0
## 248   11.5 11.0  1.0736  70 134.25  67.00   21.1 118.9 34.9  89.2  83.6  88.8
## 249   32.3 33.6  1.0236  72 201.00  69.75   29.1 136.1 40.9 108.5 105.0 104.5
## 251   25.3 26.0  1.0399  72 190.75  70.50   27.0 142.6 38.9 108.3 101.3  97.8
## 252   30.7 31.9  1.0271  74 207.50  70.00   29.8 143.7 40.8 112.4 108.5 107.1
##     thigh knee ankle biceps forearm wrist
## 1    59.0 37.3  21.9   32.0    27.4  17.1
## 2    58.7 37.3  23.4   30.5    28.9  18.2
## 3    59.6 38.9  24.0   28.8    25.2  16.6
## 4    60.1 37.3  22.8   32.4    29.4  18.2
## 5    63.2 42.2  24.0   32.2    27.7  17.7
## 6    66.0 42.0  25.6   35.7    30.6  18.8
## 7    58.4 38.3  22.9   31.9    27.8  17.7
## 8    60.0 39.4  23.2   30.5    29.0  18.8
## 9    62.9 38.3  23.8   35.9    31.1  18.2
## 11   59.7 39.7  25.2   32.8    29.4  18.5
## 12   66.2 39.2  25.9   37.2    30.2  19.0
## 13   63.4 38.3  21.5   32.5    28.6  17.7
## 14   66.0 41.5  23.7   36.9    31.6  18.8
## 15   69.0 39.0  23.1   36.1    30.5  18.2
## 16   63.1 38.7  21.7   31.1    26.4  16.9
## 17   64.8 40.8  23.1   36.2    30.8  17.3
## 18   66.9 40.0  24.4   38.2    31.6  19.3
## 19   64.2 38.7  22.9   37.2    30.5  18.5
## 21   63.5 38.0  22.1   32.5    30.3  18.4
## 22   63.4 40.6  24.6   33.0    32.8  19.9
## 23   57.4 35.3  22.2   27.9    25.9  16.7
## 24   54.9 36.2  22.1   29.8    26.7  17.1
## 25   58.4 35.5  22.9   31.1    28.0  17.6
## 26   55.0 36.7  22.5   29.9    28.2  17.7
## 27   51.7 34.7  21.4   28.7    27.0  16.5
## 28   57.5 36.0  21.0   29.2    26.6  17.0
## 29   50.1 34.5  21.3   30.5    27.9  17.2
## 31   57.5 38.7  33.9   32.5    27.7  18.4
## 32   58.5 38.8  21.5   30.1    26.4  17.9
## 33   57.3 36.2  24.5   29.0    30.0  18.8
## 34   67.1 44.2  25.2   37.5    31.5  18.7
## 35   71.2 43.3  26.3   37.3    31.7  19.7
## 36   61.9 38.3  21.9   32.0    29.8  17.0
## 37   63.5 39.9  22.6   35.1    30.6  19.0
## 38   66.0 41.5  24.7   33.2    30.5  19.4
## 39   87.3 49.1  29.6   45.0    29.0  21.4
## 41   72.5 39.6  26.6   36.4    32.7  21.4
## 42   70.6 42.5  23.7   33.6    28.7  17.4
## 43   67.7 40.9  25.0   36.7    29.8  18.4
## 44   65.0 40.2  23.0   35.8    31.5  18.8
## 45   50.0 34.7  21.0   26.1    23.1  16.1
## 46   58.4 38.2  23.4   29.7    27.4  18.3
## 47   53.3 34.5  22.5   27.9    26.2  17.3
## 48   52.7 37.5  21.9   28.8    26.8  17.9
## 49   52.0 35.8  20.6   28.8    25.5  16.3
## 51   52.6 37.2  22.4   26.0    25.8  17.3
## 52   51.4 34.9  21.0   26.7    26.1  17.2
## 53   49.3 33.7  21.4   29.6    26.0  16.9
## 54   52.6 37.6  22.6   38.5    27.4  18.5
## 55   51.9 34.9  22.5   27.7    27.5  18.5
## 56   57.2 38.0  22.0   35.9    30.2  18.9
## 57   62.1 39.6  22.5   33.1    28.3  18.5
## 58   61.8 39.8  22.7   37.7    30.9  19.2
## 59   59.1 38.0  22.5   31.6    28.8  18.2
## 61   61.6 40.9  23.1   36.2    31.8  20.2
## 62   60.9 38.0  22.1   32.5    29.8  18.3
## 63   61.0 39.4  23.6   32.7    29.9  19.1
## 64   60.8 40.1  22.7   33.6    29.0  18.8
## 65   65.0 41.2  24.7   35.3    31.1  18.4
## 66   64.8 40.2  22.7   34.8    30.1  18.7
## 67   55.0 36.1  21.7   29.6    27.4  17.4
## 68   54.3 35.4  21.5   32.8    27.4  18.7
## 69   56.0 37.4  22.4   32.6    28.1  18.1
## 71   56.6 38.6  22.4   31.5    27.3  18.6
## 72   58.9 37.6  21.6   30.3    27.3  18.3
## 73   52.9 37.5  23.1   29.7    27.3  18.2
## 74   50.9 35.4  19.1   29.3    25.7  16.9
## 75   55.5 35.2  20.9   29.4    27.0  16.8
## 76   54.5 37.0  21.4   29.3    27.0  18.3
## 77   56.7 39.7  24.2   30.2    29.2  18.1
## 78   55.0 38.3  21.8   30.8    25.7  18.8
## 79   53.5 37.5  21.5   31.4    26.8  18.3
## 81   54.8 38.2  23.7   29.4    27.2  19.0
## 82   54.8 38.0  22.0   29.9    25.2  17.7
## 83   58.9 39.0  23.0   34.3    29.6  19.0
## 84   56.0 36.5  24.1   31.2    27.3  19.2
## 85   56.3 36.6  22.0   29.7    26.3  18.0
## 86   54.7 37.8  33.7   32.4    27.7  18.2
## 87   57.2 37.7  21.8   32.6    28.0  18.8
## 88   57.8 39.5  23.3   29.2    28.4  18.1
## 89   61.0 38.4  23.8   30.2    29.3  18.8
## 91   58.3 38.2  22.5   29.1    27.7  17.7
## 92   58.9 39.7  23.1   31.4    28.4  18.8
## 93   56.9 38.3  22.1   30.1    28.2  18.4
## 94   58.9 40.5  24.5   33.3    29.6  19.1
## 95   57.4 39.6  24.6   30.3    27.9  17.8
## 96   61.7 42.3  23.2   32.9    30.8  20.4
## 97   60.6 39.4  22.9   31.6    30.1  18.5
## 98   62.1 39.3  23.3   30.6    27.8  18.2
## 99   54.7 37.3  21.9   31.6    27.5  18.2
## 101  59.0 39.4  22.3   32.2    31.0  18.6
## 102  59.3 38.4  22.4   27.9    26.2  17.0
## 103  60.0 38.4  23.2   31.0    29.2  18.4
## 104  59.1 39.8  25.4   31.0    30.3  19.7
## 105  59.5 36.1  22.0   30.1    27.2  17.7
## 106  54.7 39.0  24.8   31.0    29.4  18.8
## 107  61.2 39.3  23.5   30.5    28.5  18.1
## 108  55.8 38.7  23.4   35.1    29.6  19.1
## 109  57.5 40.0  24.8   35.1    30.7  19.2
## 111  61.9 38.0  22.3   33.3    28.2  18.1
## 112  63.7 40.0  23.6   33.5    27.8  17.4
## 113  62.3 38.1  23.9   35.3    31.1  19.8
## 114  61.5 37.8  21.9   30.7    27.6  17.4
## 115  59.3 38.1  21.8   31.8    27.3  17.5
## 116  55.9 36.3  22.1   29.8    26.3  17.3
## 117  58.8 38.4  22.8   29.9    28.0  18.1
## 118  56.8 38.8  23.3   33.4    29.8  19.5
## 119  64.6 41.1  24.8   33.6    29.5  18.5
## 121  58.5 39.3  24.6   33.9    31.2  19.5
## 122  57.1 40.5  23.2   33.0    29.6  18.4
## 123  60.5 38.7  22.6   34.4    28.0  17.6
## 124  58.1 36.5  22.1   30.6    27.5  17.6
## 125  58.5 36.6  23.5   34.4    29.2  18.0
## 126  60.7 36.0  21.9   35.6    30.2  17.6
## 127  60.7 36.8  22.2   33.8    30.3  17.2
## 128  53.5 35.8  20.8   33.9    28.2  17.4
## 129  61.7 39.0  21.8   33.3    29.6  18.1
## 131  57.0 38.7  23.2   27.5    26.5  17.6
## 132  60.3 38.5  23.0   31.2    28.4  17.1
## 133  61.2 38.1  22.6   33.5    28.6  17.9
## 134  56.1 35.6  20.5   33.6    29.3  17.3
## 135  56.0 36.9  23.0   34.0    29.8  18.1
## 136  58.9 37.9  22.7   30.9    28.8  17.6
## 137  58.8 36.1  22.4   32.7    28.3  17.1
## 138  63.6 39.2  23.8   34.3    28.4  17.7
## 139  58.1 38.4  22.5   31.7    27.4  17.6
## 141  59.1 39.6  21.6   30.8    27.9  16.6
## 142  60.4 38.2  22.0   32.0    28.5  17.8
## 143  57.1 36.7  22.3   31.6    27.5  17.9
## 144  56.1 36.1  22.7   30.5    27.2  18.2
## 145  59.1 37.6  23.2   31.8    29.7  18.3
## 146  56.4 36.5  22.0   33.5    28.3  17.3
## 147  71.2 43.5  25.2   36.1    30.3  18.7
## 148  68.4 40.8  24.6   33.3    29.7  18.4
## 149  51.9 35.7  22.0   25.8    25.2  16.9
## 151  56.9 35.9  20.4   31.6    29.0  17.8
## 152  72.9 43.5  25.1   38.5    33.8  19.6
## 153  53.6 36.8  23.8   27.8    26.3  17.4
## 154  56.8 37.4  22.8   30.6    28.3  17.9
## 155  62.1 40.0  24.9   33.7    29.2  19.4
## 156  57.3 37.8  21.7   32.2    27.7  17.7
## 157  68.6 40.0  25.2   35.2    30.7  19.1
## 158  60.8 38.5  25.0   31.6    28.0  18.6
## 159  50.1 34.8  21.8   27.0    34.9  16.9
## 161  52.5 36.6  21.0   27.0    26.3  16.5
## 162  61.4 40.6  25.0   31.3    29.2  19.1
## 163  64.0 37.3  23.5   33.5    30.6  19.7
## 164  52.4 35.6  20.4   28.3    26.2  16.5
## 165  63.8 42.0  23.4   34.0    31.2  18.5
## 166  64.8 41.3  25.6   36.4    33.7  19.4
## 167  55.5 34.2  21.9   30.2    28.7  17.7
## 168  63.3 41.7  24.6   37.2    33.1  19.8
## 169  74.4 40.6  24.0   36.1    31.8  18.8
## 171  54.7 36.2  22.1   30.4    27.4  17.7
## 172  50.0 34.8  22.0   24.8    25.9  16.9
## 173  57.8 37.3  22.4   31.0    28.7  17.7
## 174  59.2 37.7  21.5   32.4    28.4  17.8
## 175  69.2 42.4  24.0   35.4    21.0  20.1
## 176  50.3 34.8  22.2   31.0    26.9  16.9
## 177  60.0 38.1  22.0   31.5    26.6  16.7
## 178  69.8 42.6  24.8   34.4    29.5  18.4
## 179  64.7 39.5  24.7   34.8    30.3  18.1
## 181  68.1 42.8  24.1   35.6    29.0  19.0
## 182  47.2 33.5  20.2   27.7    24.6  16.5
## 183  54.1 36.2  21.8   31.4    28.3  17.2
## 184  58.0 39.4  22.7   30.0    26.4  17.4
## 185  55.4 38.9  22.4   30.5    28.9  17.7
## 186  57.3 39.7  22.6   32.9    29.3  18.2
## 187  67.7 41.3  24.7   37.2    31.8  20.0
## 188  63.5 39.8  23.5   36.4    30.4  19.1
## 189  62.8 41.3  24.8   36.6    32.4  18.8
## 191  54.3 36.3  21.8   29.6    27.3  17.9
## 192  68.5 45.0  25.5   37.1    31.2  19.9
## 193  60.6 38.6  24.7   34.0    30.1  18.7
## 194  65.3 43.3  26.0   33.7    29.9  18.5
## 195  60.2 38.9  22.4   31.7    27.1  17.1
## 196  61.1 38.4  24.1   32.9    29.8  18.8
## 197  57.7 38.6  24.0   31.2    27.3  17.4
## 198  62.3 38.0  22.3   30.8    27.8  16.9
## 199  57.5 40.0  22.5   30.6    30.0  18.5
## 201  60.2 39.2  23.8   31.7    28.4  18.6
## 202  55.5 35.7  22.0   29.4    26.6  17.4
## 203  68.8 38.3  23.7   32.1    28.9  18.7
## 204  60.6 39.0  24.0   32.9    29.2  18.4
## 205  63.5 40.3  21.8   34.8    30.7  17.4
## 206  63.3 39.8  24.1   37.3    23.1  19.4
## 207  58.9 37.6  21.4   33.1    29.5  17.3
## 208  60.7 39.4  23.3   36.7    31.6  18.4
## 209  53.0 36.2  22.5   31.4    27.5  17.7
## 211  51.1 35.0  21.7   30.9    28.8  17.4
## 212  63.7 40.3  23.2   36.8    31.0  18.9
## 213  56.3 38.8  23.0   29.5    27.9  18.6
## 214  60.0 40.9  25.5   32.7    30.0  19.0
## 215  57.1 38.1  21.8   28.6    26.7  18.0
## 216  62.5 36.9  23.6   34.7    29.1  18.4
## 217  53.8 36.5  21.5   31.3    26.3  17.8
## 218  54.7 39.0  22.6   27.5    25.9  18.6
## 219  63.9 39.2  22.9   35.7    30.4  19.2
## 221  57.7 38.1  23.9   31.4    29.9  18.9
## 222  64.2 42.7  27.0   38.4    32.0  19.6
## 223  59.7 40.2  23.4   27.9    27.0  17.8
## 224  54.4 35.2  22.5   29.4    26.8  17.0
## 225  57.4 37.1  21.8   34.1    31.1  19.2
## 226  50.8 33.0  19.7   25.3    22.0  15.8
## 227  56.6 38.5  22.6   33.4    29.3  18.8
## 228  64.0 42.6  23.4   33.2    30.0  18.4
## 229  58.4 37.4  22.5   34.6    30.1  18.8
## 231  53.0 35.0  21.3   31.7    27.3  16.9
## 232  59.3 38.6  22.8   31.8    29.1  19.0
## 233  57.1 38.9  23.6   30.9    29.6  18.0
## 234  56.8 35.9  21.0   27.8    26.1  17.6
## 235  54.3 35.7  21.0   31.3    28.7  18.3
## 236  54.4 37.1  22.7   30.3    26.3  18.3
## 237  59.3 40.3  23.0   32.6    28.5  19.0
## 238  63.4 41.1  22.3   35.1    29.6  18.5
## 239  61.2 39.1  22.3   29.8    28.9  18.3
## 241  50.7 33.4  20.1   28.5    24.8  16.5
## 242  61.3 42.1  23.4   34.9    30.1  19.4
## 243  63.7 42.4  24.6   35.6    30.7  19.5
## 244  65.6 46.0  25.4   35.3    29.8  19.5
## 245  58.2 38.8  24.1   32.1    29.3  18.5
## 246  54.3 37.5  22.6   29.2    27.3  18.5
## 247  63.3 44.0  22.6   37.5    32.6  18.8
## 248  49.6 34.8  21.5   25.6    25.7  18.5
## 249  59.6 40.8  23.2   35.2    28.6  20.1
## 251  56.0 41.6  22.7   30.5    29.4  19.8
## 252  59.3 42.2  24.6   33.7    30.0  20.9
subset
##     brozek siri density age weight height adipos  free neck chest abdom   hip
## 10    12.0 11.7  1.0722  23 198.25  73.50   25.8 174.4 42.1  99.6  88.6 104.1
## 20    16.5 16.5  1.0610  33 211.75  73.50   27.6 176.8 40.0 106.2 100.5 109.0
## 30     9.4  8.8  1.0790  29 160.75  69.00   23.8 145.7 36.7  97.4  83.5  98.7
## 40    31.3 32.6  1.0258  50 203.00  67.00   31.8 139.4 40.2 114.8 108.1 102.5
## 50     5.0  4.0  1.0903  47 127.50  66.75   20.2 121.2 34.0  83.4  70.4  87.2
## 60    24.0 24.6  1.0430  61 179.75  65.75   29.2 136.7 38.4 104.8  98.3  99.6
## 70    13.2 12.9  1.0693  55 156.75  71.50   21.6 136.1 36.3  94.4  84.6  94.3
## 80    18.6 18.8  1.0560  66 171.25  69.25   25.1 139.5 37.4 102.7  98.6 100.2
## 90    14.3 14.1  1.0666  48 176.00  73.00   23.3 150.9 36.7  96.7  86.5  98.3
## 100   21.7 22.2  1.0484  47 197.00  72.00   26.7 154.2 40.0 107.5  94.0 103.7
## 110   21.0 21.4  1.0500  40 168.50  69.25   24.7 133.1 34.2  97.8  92.3 100.6
## 120   17.9 18.1  1.0575  44 187.50  72.25   25.3 153.8 38.0 101.8  87.5 101.0
## 130   15.0 14.9  1.0648  42 165.25  69.75   23.9 140.5 38.3  96.2  87.0  97.8
## 140   20.1 20.4  1.0522  49 212.75  75.00   26.6 169.9 40.8 104.7 106.6 107.7
## 150   24.5 25.2  1.0416  26 223.00  70.25   31.8 168.3 40.6 114.1 106.8 113.9
## 160   22.0 22.5  1.0477  31 177.25  71.50   24.4 138.3 36.2 101.1  92.4  99.3
## 170   16.5 16.5  1.0610  35 172.75  69.50   25.2 144.2 37.6  99.1  90.8  98.1
## 180   16.8 16.9  1.0603  39 234.75  74.50   29.8 195.3 42.8 109.5 104.5 109.9
## 190   23.8 24.4  1.0434  41 185.00  68.25   28.0 141.0 38.0 103.4 101.2 103.1
## 200   22.6 23.6  1.0462  43 170.75  67.50   26.4 132.1 37.4 103.7  89.7  94.2
## 210   11.3 10.8  1.0742  47 159.75  70.75   22.5 141.8 34.5  92.9  84.4  94.0
## 220   15.1 15.0  1.0646  53 154.50  69.25   22.7 131.2 37.6  93.9  88.7  94.5
## 230   16.9 17.0  1.0601  56 167.75  68.50   25.2 139.4 37.4  98.6  93.0  97.0
## 240   28.8 29.9  1.0316  65 189.75  65.75   30.9 135.1 40.8 106.4 100.5 100.5
## 250   28.3 29.3  1.0328  72 186.75  66.00   30.2 133.9 38.9 111.1 111.5 101.7
##     thigh knee ankle biceps forearm wrist
## 10   63.1 41.7  25.0   35.6    30.0  19.2
## 20   65.8 40.6  24.0   37.1    30.1  18.2
## 30   58.9 35.3  22.6   30.1    26.7  17.6
## 40   61.3 41.1  24.7   34.1    31.0  18.3
## 50   50.6 34.4  21.9   26.8    25.8  16.8
## 60   60.6 37.7  22.9   34.5    29.6  18.5
## 70   51.2 37.4  21.6   27.3    27.1  17.3
## 80   56.5 39.3  22.7   30.3    28.7  19.0
## 90   60.4 39.9  24.4   28.8    29.6  18.7
## 100  62.7 39.0  22.3   35.3    30.9  18.3
## 110  57.5 36.8  22.8   32.1    26.0  17.3
## 120  58.5 39.2  24.5   32.1    28.6  18.0
## 130  57.4 36.9  22.2   31.6    27.8  17.7
## 140  66.5 42.5  24.5   35.5    29.8  18.7
## 150  67.6 42.7  24.7   36.0    30.4  18.4
## 160  59.4 39.0  24.6   30.1    28.2  18.2
## 170  60.1 39.1  23.4   32.5    29.8  17.4
## 180  69.5 43.1  25.8   39.1    32.5  19.9
## 190  61.5 40.4  22.9   33.4    29.2  18.5
## 200  58.5 39.0  24.1   33.8    28.8  18.8
## 210  56.0 38.2  22.6   29.0    26.2  17.6
## 220  53.7 36.2  22.0   28.5    25.7  17.1
## 230  55.4 38.8  23.2   32.4    29.7  19.0
## 240  59.2 38.1  24.0   35.9    30.5  19.1
## 250  60.3 37.3  21.5   31.3    27.2  18.0
rmse <- function(x,y) sqrt(mean((x-y)^2))
rmse
## function(x,y) sqrt(mean((x-y)^2))
lmodtr1 <-lm(siri~.-brozek -density, train)
summary(lmodtr1)
## 
## Call:
## lm(formula = siri ~ . - brozek - density, data = train)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -5.8314 -0.6722  0.1828  0.9150  6.6619 
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)    
## (Intercept) -12.591885   6.448868  -1.953 0.052193 .  
## age           0.007978   0.012320   0.648 0.517983    
## weight        0.362999   0.023314  15.570  < 2e-16 ***
## height        0.049026   0.040315   1.216 0.225315    
## adipos       -0.514032   0.114074  -4.506 1.09e-05 ***
## free         -0.564773   0.014889 -37.933  < 2e-16 ***
## neck          0.016525   0.089863   0.184 0.854272    
## chest         0.120219   0.039590   3.037 0.002694 ** 
## abdom         0.140108   0.042186   3.321 0.001056 ** 
## hip           0.006197   0.056101   0.110 0.912148    
## thigh         0.195057   0.054460   3.582 0.000424 ***
## knee          0.106637   0.093534   1.140 0.255542    
## ankle         0.125118   0.081303   1.539 0.125325    
## biceps        0.096199   0.064656   1.488 0.138278    
## forearm       0.230775   0.073332   3.147 0.001888 ** 
## wrist         0.139279   0.206804   0.673 0.501378    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.55 on 211 degrees of freedom
## Multiple R-squared:  0.9692, Adjusted R-squared:  0.967 
## F-statistic: 442.5 on 15 and 211 DF,  p-value: < 2.2e-16
rmse(lmodtr1$fit, train$siri)
## [1] 1.494315
pred1 <- predict(lmodtr1, subset)
y10 <- subset$siri
rmse(pred1,y10)
## [1] 1.131529

The predicted result seems to be pretty good (with RMSE 1.1315288), even better than the overall fitted model (RMSE 1.4943146) .

  1. Linear regression with stepwise variable selection
lmodtr2 <- step(lmodtr1)
## Start:  AIC=214.36
## siri ~ (brozek + density + age + weight + height + adipos + free + 
##     neck + chest + abdom + hip + thigh + knee + ankle + biceps + 
##     forearm + wrist) - brozek - density
## 
##           Df Sum of Sq    RSS    AIC
## - hip      1       0.0  506.9 212.37
## - neck     1       0.1  507.0 212.39
## - age      1       1.0  507.9 212.81
## - wrist    1       1.1  508.0 212.84
## - knee     1       3.1  510.0 213.75
## - height   1       3.6  510.4 213.94
## <none>                  506.9 214.36
## - biceps   1       5.3  512.2 214.73
## - ankle    1       5.7  512.6 214.89
## - chest    1      22.2  529.0 222.07
## - forearm  1      23.8  530.7 222.77
## - abdom    1      26.5  533.4 223.92
## - thigh    1      30.8  537.7 225.75
## - adipos   1      48.8  555.7 233.21
## - weight   1     582.4 1089.3 386.01
## - free     1    3456.8 3963.7 679.21
## 
## Step:  AIC=212.37
## siri ~ age + weight + height + adipos + free + neck + chest + 
##     abdom + thigh + knee + ankle + biceps + forearm + wrist
## 
##           Df Sum of Sq    RSS    AIC
## - neck     1       0.1  507.0 210.40
## - age      1       1.0  507.9 210.81
## - wrist    1       1.1  508.0 210.86
## - knee     1       3.2  510.1 211.80
## - height   1       3.5  510.4 211.95
## <none>                  506.9 212.37
## - biceps   1       5.3  512.2 212.73
## - ankle    1       5.7  512.6 212.89
## - chest    1      23.1  530.0 220.50
## - forearm  1      23.8  530.7 220.78
## - abdom    1      27.9  534.9 222.55
## - thigh    1      34.2  541.2 225.21
## - adipos   1      50.3  557.2 231.85
## - weight   1     683.9 1190.8 404.23
## - free     1    3488.9 3995.8 679.05
## 
## Step:  AIC=210.4
## siri ~ age + weight + height + adipos + free + chest + abdom + 
##     thigh + knee + ankle + biceps + forearm + wrist
## 
##           Df Sum of Sq    RSS    AIC
## - age      1       1.1  508.1 208.88
## - wrist    1       1.3  508.3 208.99
## - knee     1       3.1  510.1 209.80
## - height   1       3.6  510.6 210.02
## <none>                  507.0 210.40
## - biceps   1       5.4  512.4 210.80
## - ankle    1       5.6  512.6 210.89
## - chest    1      23.2  530.2 218.55
## - forearm  1      24.6  531.6 219.15
## - abdom    1      28.0  535.0 220.60
## - thigh    1      34.4  541.4 223.29
## - adipos   1      50.8  557.8 230.07
## - weight   1     689.6 1196.6 403.34
## - free     1    3532.0 4039.0 679.49
## 
## Step:  AIC=208.88
## siri ~ weight + height + adipos + free + chest + abdom + thigh + 
##     knee + ankle + biceps + forearm + wrist
## 
##           Df Sum of Sq    RSS    AIC
## - wrist    1       2.9  511.0 208.19
## - height   1       3.3  511.4 208.35
## - knee     1       4.5  512.5 208.87
## <none>                  508.1 208.88
## - ankle    1       5.2  513.2 209.18
## - biceps   1       6.0  514.0 209.53
## - forearm  1      23.6  531.6 217.18
## - chest    1      24.2  532.3 217.46
## - abdom    1      33.7  541.8 221.48
## - thigh    1      35.3  543.3 222.12
## - adipos   1      51.1  559.1 228.63
## - weight   1     699.1 1207.2 403.34
## - free     1    3598.0 4106.0 681.23
## 
## Step:  AIC=208.19
## siri ~ weight + height + adipos + free + chest + abdom + thigh + 
##     knee + ankle + biceps + forearm
## 
##           Df Sum of Sq    RSS    AIC
## - height   1       3.8  514.8 207.89
## <none>                  511.0 208.19
## - knee     1       5.7  516.7 208.72
## - ankle    1       6.9  517.9 209.24
## - biceps   1       7.0  518.0 209.30
## - chest    1      23.8  534.8 216.53
## - forearm  1      27.7  538.7 218.16
## - thigh    1      32.4  543.4 220.13
## - abdom    1      37.3  548.3 222.19
## - adipos   1      49.3  560.3 227.11
## - weight   1     696.5 1207.5 401.40
## - free     1    3798.4 4309.4 690.20
## 
## Step:  AIC=207.89
## siri ~ weight + adipos + free + chest + abdom + thigh + knee + 
##     ankle + biceps + forearm
## 
##           Df Sum of Sq    RSS    AIC
## <none>                  514.8 207.89
## - knee     1       5.1  519.9 208.12
## - ankle    1       7.4  522.2 209.11
## - biceps   1       7.5  522.4 209.18
## - chest    1      24.0  538.9 216.25
## - forearm  1      28.8  543.6 218.23
## - thigh    1      30.0  544.8 218.73
## - abdom    1      39.1  553.9 222.49
## - adipos   1      86.6  601.4 241.18
## - weight   1     819.8 1334.7 422.13
## - free     1    3809.4 4324.2 688.98
summary(lmodtr2)
## 
## Call:
## lm(formula = siri ~ weight + adipos + free + chest + abdom + 
##     thigh + knee + ankle + biceps + forearm, data = train)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -5.7926 -0.6839  0.2371  0.8824  6.8655 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept) -7.24766    3.94171  -1.839 0.067331 .  
## weight       0.36973    0.01994  18.546  < 2e-16 ***
## adipos      -0.57022    0.09459  -6.028 7.09e-09 ***
## free        -0.55965    0.01400 -39.978  < 2e-16 ***
## chest        0.12099    0.03809   3.176 0.001709 ** 
## abdom        0.15824    0.03908   4.049 7.18e-05 ***
## thigh        0.16140    0.04551   3.546 0.000479 ***
## knee         0.12767    0.08749   1.459 0.145947    
## ankle        0.13817    0.07858   1.758 0.080106 .  
## biceps       0.11222    0.06317   1.777 0.077055 .  
## forearm      0.24281    0.06990   3.474 0.000620 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.544 on 216 degrees of freedom
## Multiple R-squared:  0.9687, Adjusted R-squared:  0.9673 
## F-statistic: 668.6 on 10 and 216 DF,  p-value: < 2.2e-16
rmse(lmodtr2$fit, train$siri)
## [1] 1.505982
pred2 <- predict(lmodtr2, subset)
rmse(pred2,y10)
## [1] 1.12202

Using AIC, we can get rid of 5 variables, and the predicted value (RMSE 1.1220197) is better than using the full model (RMSE 1.1315288)

  1. Principal component regression
pca <- prcomp(train[,4:18])
#the square root of eigenvalues:
round(pca$sdev,3)
##  [1] 37.047 14.540  9.511  3.634  3.462  2.544  2.200  1.854  1.577  1.412
## [11]  1.324  1.249  1.037  0.828  0.491
plot(pca$sdev, type="l",ylab="the square root of eigenvalues", xlab="PC number")

#We can see that the contribution drops very quickly, most variance can be explained by the first 5 PCs.
matplot(1:15,pca$rot[,1:5],type="l", xlab="Frequency")

#Seems that the 1st PC is a linear combination of most the 2nd and 5th variables (weight and free), which agrees with results from previous steps, where the regression coefficients of weight and free are most significant.
# Do regression on PCs:
lmodtr3 <- lm(train$siri ~ pca$x[,1:6])
summary(lmodtr3)
## 
## Call:
## lm(formula = train$siri ~ pca$x[, 1:6])
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -10.3999  -0.5813   0.3739   0.9771   7.9487 
## 
## Coefficients:
##                  Estimate Std. Error t value Pr(>|t|)    
## (Intercept)     19.180176   0.110619 173.389  < 2e-16 ***
## pca$x[, 1:6]PC1 -0.126730   0.002992 -42.350  < 2e-16 ***
## pca$x[, 1:6]PC2 -0.356175   0.007625 -46.714  < 2e-16 ***
## pca$x[, 1:6]PC3  0.464783   0.011656  39.875  < 2e-16 ***
## pca$x[, 1:6]PC4  0.309275   0.030508  10.138  < 2e-16 ***
## pca$x[, 1:6]PC5 -0.091998   0.032023  -2.873  0.00447 ** 
## pca$x[, 1:6]PC6 -0.222500   0.043574  -5.106 7.11e-07 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.667 on 220 degrees of freedom
## Multiple R-squared:  0.9629, Adjusted R-squared:  0.9618 
## F-statistic: 950.5 on 6 and 220 DF,  p-value: < 2.2e-16
#Now do prediction. We need to first center the predictors.

mm = apply(train[,4:18],2,mean) #the mean of each X variables
tx = as.matrix(sweep(subset[,4:18],2,mm)) #centering the predicting data.

nx = tx %*% pca$rot[,1:6]
pred3 = cbind(1,nx) %*% lmodtr3$coef

rmse(pred3,y10)
## [1] 1.141839
#So the RMSE of PCA using the first 6 PCs is about the same with AIC.

a <- numeric(15)
for (i in 1:15) {
  nx = tx %*% pca$rot[,1:i]
  model4 = lm(train$siri ~pca$x[,1:i])
  pred4 = cbind(1,nx) %*% model4$coef
  a[i] = rmse(pred4,y10)
}
plot(a, ylab="Test RMSE",xlab="PC number")

which.min(a)
## [1] 7
#So when using the first 7 PCs gives the best RMSE.
  1. Partial least squares
library(pls)
## Warning: package 'pls' was built under R version 3.6.2
## 
## Attaching package: 'pls'
## The following object is masked from 'package:stats':
## 
##     loadings
par(mfrow=c(1,2))
trainx = as.matrix(sweep(train[,4:18],2,mm))
gpls <- plsr(train$siri~trainx, ncomp=15, validation="CV",grpsize=10)
plot(gpls$validat$adj[1,], xlab="PC number",ylab="CV adj")
plot(MSEP(gpls)$val[1,,], xlab="PC number", ylab="CV RMS")

ypred.train = predict(gpls,  ncomp=4)
rmse(ypred.train, train$siri)
## [1] 1.612759
testx = as.matrix(sweep(subset[,4:18],2,mm))
ypred.subset = predict(gpls,newdata=testx,ncomp=4)
rmse(ypred.subset,subset$siri)
## [1] 1.12449

model with low RMSE

  1. Ridge regression:
par(mfrow=c(1,1))
trainx = as.matrix(sweep(train[,4:18],2,mm))
testx = as.matrix(sweep(subset[,4:18],2,mm))
yc = train$siri - mean(train$siri)
gridge = lm.ridge(yc ~trainx, lambda=seq(0,10,1e-4))
matplot(gridge$lambda,t(gridge$coef), type="l",lty=1,xlab=expression(lambda),ylab=expression(hat(beta)),main="the Ridge Trace plot")

#Various automatic selections for lambda are available:
select(gridge)
## modified HKB estimator is 0.1340921 
## modified L-W estimator is 0.4446052 
## smallest value of GCV  at 0.0467
yfit = scale(trainx, center=F, scale=gridge$scales) %*% gridge$coef[,468] + mean(train$siri)
rmse(yfit,train$siri )
## [1] 1.494422
ypred = scale(testx, center=F, scale=gridge$scales) %*% gridge$coef[,468] + mean(train$siri)
rmse(ypred,subset$siri )
## [1] 1.128089