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) .
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)
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.
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
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