Paul Bernal
2019-Mar-20 14:13 UTC
[R] How To Extracting Summary Stats for Sepwise Regression
Hello dear friends, What I want is to, somehow data.frame the summary of a stepwise regression (the estimated coefficients and p-values) to show them as a report in table format. However I have not been able to do that. Any help and/or guidance will be greatly appreciated, below are all the details of what I am doing: I am working with R version 3.5.2 in Power BI. The Power BI version is: 2.67.5404.581 64-bit (March 2019) The error I am getting is the following: Error in step(basis_model, direction = "both", scope = formula(FitAll)) : number of rows in use has changed: remove missing values? In addition: Warning message: In add1.lm(fit, scope$add, scale = scale, trace = trace, k = k, : using the 106/107 rows from a combined fit Execution halted BasisModel <- lm(datos$Transits ~ 1) FitAllModel_2 <- lm(datos$Transits ~ datos$CPI + datos$ConCharter + datos$NewBuildingPrice + datos$HFO) StepWiseModel_3 <- step(BasisModel, direction="both", scope=formula(FitAllModel_2)) sum <- summary(StepWiseModel_3) grid.table(sum) I have also tried doing grid.table(sum$coef) and grid.table(sum$coefficients) but none of that worked.> dput(datos)structure(list(Transits = c(23L, 86L, 29L, 85L, 36L, 80L, 36L, 68L, 26L, 97L, 24L, 85L, 27L, 87L, 26L, 66L, 28L, 76L, 23L, 68L, 67L, 38L, 30L, 68L, 65L, 37L, 31L, 85L, 27L, 85L, 26L, 87L, 96L, 30L, 24L, 85L, 22L, 82L, 23L, 80L, 81L, 26L, 25L, 81L, 29L, 94L, 23L, 83L, 90L, 33L, 27L, 94L, 25L, 90L, 17L, 86L, 40L, 89L, 31L, 81L, 87L, 26L, 18L, 76L, 85L, 17L, 23L, 87L, 88L, 23L, 37L, 82L, 38L, 102L, 39L, 90L, 40L, 93L, 33L, 84L, 41L, 84L, 27L, 71L, 37L, 63L, 26L, 76L, 28L, 63L, 34L, 74L, 37L, 68L, 37L, 76L, 42L, 79L, 41L, 81L, 40L, 77L, 34L, 74L, 39L, 58L, 35L, 55L, 36L, 51L, 34L, 57L, 39L, 61L, 35L, 68L, 43L, 70L, 40L, 59L, 46L, 60L, 34L, 63L, 42L, 67L, 43L, 58L, 44L, 59L, 53L, 26L, 32L, 58L, 34L, 41L, 42L, 39L, 30L, 42L, 30L, 48L, 37L, 42L, 36L, 60L, 44L, 65L, 43L, 55L, 22L, 53L, 31L, 8L, 9L, 25L, 12L, 33L, 25L, 45L, 22L, 47L, 25L, 43L, 28L, 37L, 28L, 48L, 19L, 35L, 11L, 28L, 9L, 22L, 12L, 25L, 9L, 14L, 9L, 13L, 13L, 12L, 13L, 25L, 10L, 24L, 13L, 24L, 10L, 22L, 10L, 19L, 11L, 15L, 8L, 16L, 12L, 9L, 8L, 17L, 9L, 4L, 8L, 11L, 8L, 12L, 6L, 11L, 8L, 11L, 21L, 8L, 14L, 11L, 12L, 11L, 10L, 7L, 13L, 6L, 8L), CargoTons = c(327994L, 1606616L, 392457L, 1712810L, 572909L, 1693921L, 696160L, 1633810L, 427511L, 1990180L, 434473L, 1838645L, 501521L, 1970717L, 548017L, 1588291L, 488718L, 1739769L, 457908L, 1623617L, 1677868L, 709159L, 467997L, 1491590L, 1329897L, 540330L, 454764L, 1696188L, 412585L, 1771492L, 514509L, 1843676L, 2033267L, 499482L, 393591L, 1610270L, 333572L, 1546048L, 417399L, 1901424L, 1686199L, 432638L, 428805L, 1861493L, 511944L, 2011693L, 421151L, 1837227L, 1912533L, 593826L, 453480L, 1921809L, 460072L, 1675982L, 341377L, 1774207L, 473299L, 1526080L, 499347L, 1464667L, 1684398L, 437215L, 309439L, 1739618L, 1764891L, 344394L, 393377L, 1671126L, 1663948L, 419590L, 480516L, 1593703L, 527984L, 1730980L, 505520L, 1537961L, 527088L, 1600723L, 484147L, 1541796L, 662549L, 1506180L, 433022L, 1424827L, 643535L, 1337935L, 424964L, 1622296L, 462553L, 1194666L, 514447L, 1840988L, 641739L, 1382457L, 563509L, 1378973L, 667025L, 1345523L, 649433L, 1338687L, 673892L, 1348285L, 584263L, 1528680L, 600019L, 944145L, 649869L, 912873L, 624110L, 850102L, 607789L, 1208901L, 615976L, 1182815L, 615512L, 1399101L, 697432L, 1325899L, 608629L, 872243L, 750331L, 848187L, 524618L, 872075L, 706920L, 986303L, 687968L, 896847L, 636234L, 724704L, 723088L, 350679L, 460177L, 722250L, 622978L, 723952L, 660514L, 689362L, 483923L, 698217L, 569948L, 813910L, 602531L, 671289L, 494437L, 778160L, 540130L, 717010L, 590081L, 664101L, 368966L, 633134L, 447660L, 159223L, 160920L, 413470L, 271818L, 488672L, 390754L, 640038L, 321733L, 577064L, 443330L, 639336L, 430580L, 451870L, 399601L, 621382L, 256492L, 526333L, 211949L, 467734L, 173952L, 343172L, 263817L, 447959L, 166116L, 257188L, 175630L, 241603L, 239032L, 214663L, 293355L, 421452L, 165754L, 394366L, 254257L, 438731L, 201630L, 340764L, 160402L, 343059L, 203165L, 184462L, 141414L, 247706L, 223107L, 153652L, 153972L, 261942L, 200656L, 59416L, 166912L, 234790L, 177974L, 228005L, 138093L, 219065L, 189138L, 176860L, 391862L, 161972L, 260940L, 152540L, 182772L, 218337L, 225203L, 109038L, 260027L, 118230L, 180304L), RcnstPCUMS = c(814089L, 4018838L, 1039070L, 3962734L, 1307110L, 3715584L, 1272525L, 3159759L, 918793L, 4480777L, 862597L, 3959057L, 982768L, 4027394L, 932933L, 3085639L, 1008364L, 3545084L, 829730L, 3151571L, 3106543L, 1378972L, 1049638L, 3143534L, 3004616L, 1314721L, 1099565L, 3926565L, 972551L, 3950212L, 937604L, 4037464L, 4454625L, 1099265L, 870325L, 3935381L, 796631L, 3803778L, 833356L, 3704866L, 3776594L, 928562L, 903368L, 3740199L, 1024314L, 4351353L, 837034L, 3848355L, 4151682L, 1177413L, 988868L, 4359297L, 914531L, 4162950L, 621160L, 3973881L, 1435287L, 4126620L, 1111529L, 3773231L, 4050280L, 941832L, 650553L, 3528985L, 3955725L, 616642L, 825445L, 4023044L, 4097211L, 834924L, 1311944L, 3796557L, 1364762L, 4737589L, 1375555L, 4176912L, 1448418L, 4329880L, 1186202L, 3896554L, 1473688L, 3918132L, 975482L, 3306258L, 1329400L, 2945480L, 929590L, 3522287L, 1020073L, 2917674L, 1225467L, 3434734L, 1337442L, 3163533L, 1312327L, 3546667L, 1519244L, 3672447L, 1464130L, 3779991L, 1443061L, 3585158L, 1222389L, 3455774L, 1412875L, 2727757L, 1270067L, 2566944L, 1310215L, 2416900L, 1227347L, 2663666L, 1401909L, 2846029L, 1271238L, 3177336L, 1538445L, 3274950L, 1436806L, 2761960L, 1667355L, 2829878L, 1242930L, 2959036L, 1522456L, 3163181L, 1573896L, 2718126L, 1603101L, 2758533L, 2462967L, 946629L, 1151827L, 2691777L, 1246668L, 1925651L, 1527348L, 1835005L, 1087464L, 1969332L, 1093925L, 2234340L, 1315999L, 1949677L, 1300176L, 2785853L, 1588889L, 3028544L, 1563391L, 2559034L, 796993L, 2483341L, 1426101L, 283420L, 320329L, 1154490L, 430564L, 1524589L, 908063L, 2103074L, 806928L, 2190153L, 903527L, 2011725L, 1018823L, 1711333L, 1004757L, 2225944L, 689023L, 1610759L, 388765L, 1295317L, 320579L, 1011523L, 446023L, 1164442L, 323091L, 653104L, 315790L, 594555L, 464906L, 548716L, 470282L, 1165447L, 373934L, 1115930L, 464005L, 1112178L, 355084L, 1016473L, 367861L, 885959L, 402602L, 691014L, 290094L, 746782L, 435042L, 417022L, 283322L, 785941L, 320143L, 178424L, 283567L, 397199L, 283370L, 549200L, 216474L, 499255L, 290743L, 515437L, 963622L, 290782L, 508405L, 514425L, 556650L, 403736L, 362223L, 313357L, 467175L, 278034L, 288877L), TotalToll = c(5658660L, 25045209L, 8409456L, 28334952L, 10684800L, 26256312L, 10232640L, 22154904L, 7385184L, 31794048L, 7105176L, 28139544L, 8081568L, 28620648L, 7618608L, 21787128L, 7995024L, 24984648L, 6678216L, 22065480L, 21767976L, 11140344L, 8342568L, 22126968L, 21198528L, 10703664L, 8950680L, 27573120L, 7887024L, 28178856L, 8583088L, 31479312L, 34815908L, 9713854L, 7785698L, 30481326L, 7141340L, 29531732L, 7506990L, 29051552L, 29541774L, 8399644L, 8077564L, 28938572L, 9193712L, 33842728L, 7619546L, 30069734L, 31965450L, 10622842L, 8939174L, 34113946L, 8367102L, 32325008L, 5509564L, 30988230L, 12738622L, 32011730L, 10079752L, 29351100L, 31599596L, 8405222L, 5867970L, 27722768L, 31299036L, 5551330L, 7411918L, 31425586L, 32109240L, 7644508L, 11586334L, 29576948L, 12308666L, 37031904L, 12060814L, 32598000L, 13147498L, 33810750L, 10585596L, 30252354L, 13395666L, 30832518L, 8939520L, 25772004L, 12196422L, 23161516L, 8370134L, 27664010L, 9265958L, 22574300L, 10874402L, 26518982L, 12227228L, 24745802L, 11917848L, 27524048L, 13781322L, 28628988L, 13258712L, 29704254L, 13062000L, 27882468L, 11171612L, 27395884L, 12708357L, 21812064L, 11455579L, 20552168L, 12002006L, 19359130L, 11233494L, 21034682L, 12731992L, 22272600L, 11642944L, 25170646L, 13997224L, 25843164L, 12775186L, 21979684L, 15149512L, 22494414L, 11309236L, 23744618L, 13475898L, 25170222L, 14266136L, 21695278L, 14599552L, 21839098L, 19512312L, 8684618L, 10421828L, 21498934L, 11384046L, 15646056L, 13728128L, 14582010L, 9858034L, 15779454L, 9967770L, 17877786L, 11890900L, 15472342L, 11782498L, 21993570L, 14439268L, 24119724L, 14170544L, 20208244L, 7208820L, 19774008L, 11551690L, 2555626L, 2774579L, 9384376L, 3853045L, 12234653L, 8380684L, 16982380L, 7210410L, 17477956L, 8221924L, 16287822L, 9037974L, 13420990L, 9030589L, 17503524L, 6018223L, 12638167L, 3503589L, 10144213L, 2776575L, 8044772L, 4032809L, 9474450L, 2732395L, 5228934L, 2765976L, 4672118L, 4066749L, 4249184L, 4253170L, 9525707L, 3201832L, 8718645L, 4056891L, 8829727L, 3089209L, 7862348L, 3264051L, 6812088L, 3466098L, 5082556L, 2458062L, 5979517L, 3894452L, 3116222L, 2486857L, 5945299L, 2809445L, 1390359L, 2503865L, 3582466L, 2505757L, 4432436L, 1930800L, 3794537L, 2623668L, 3959920L, 7499419L, 2472087L, 4411731L, 3827331L, 4234585L, 3521434L, 3207605L, 2375325L, 4079966L, 2137927L, 2540208L ), CPI = c(212.193, 212.193, 213.856, 213.856, 215.834, 215.834, 216.687, 216.687, 218.711, 218.711, 218.803, 218.803, 219.179, 219.179, 216.741, 216.741, 217.631, 217.631, 218.009, 218.009, 218.178, 218.178, 217.965, 217.965, 218.011, 218.011, 218.312, 218.312, 218.439, 218.439, 220.223, 220.223, 226.421, 226.421, 226.23, 226.23, 225.672, 225.672, 221.309, 221.309, 223.467, 223.467, 224.906, 224.906, 225.964, 225.964, 225.722, 225.722, 225.922, 225.922, 226.545, 226.545, 226.889, 226.889, 226.665, 226.665, 231.317, 231.317, 230.221, 230.221, 229.601, 229.601, 227.663, 227.663, 229.392, 229.392, 230.085, 230.085, 229.815, 229.815, 229.478, 229.478, 229.104, 229.104, 230.379, 230.379, 231.407, 231.407, 230.28, 230.28, 233.546, 233.546, 233.069, 233.069, 233.049, 233.049, 232.166, 232.166, 232.773, 232.773, 232.531, 232.531, 232.945, 232.945, 233.504, 233.504, 233.596, 233.596, 233.877, 233.877, 234.149, 234.149, 233.916, 233.916, 237.433, 237.433, 236.151, 236.151, 234.812, 234.812, 234.781, 234.781, 236.293, 236.293, 237.072, 237.072, 237.9, 237.9, 238.343, 238.343, 238.25, 238.25, 237.852, 237.852, 238.031, 238.031, 233.707, 233.707, 237.838, 237.838, 237.336, 237.336, 236.525, 236.525, 234.722, 234.722, 236.119, 236.119, 236.599, 236.599, 237.805, 237.805, 238.638, 238.638, 238.654, 238.654, 238.316, 238.316, 237.945, 237.945, 236.916, 236.916, 241.729, 241.729, 241.353, 241.353, 241.432, 241.432, 237.111, 237.111, 238.132, 238.132, 239.261, 239.261, 240.229, 240.229, 241.018, 241.018, 240.628, 240.628, 240.849, 240.849, 241.428, 241.428, 242.839, 242.839, 246.663, 246.663, 246.669, 246.669, 246.524, 246.524, 243.603, 243.603, 243.801, 243.801, 244.524, 244.524, 244.733, 244.733, 244.955, 244.955, 244.786, 244.786, 245.519, 245.519, 246.819, 246.819, 247.867, 247.867, 252.885, 252.885, 252.038, 251.233, 248.991, 248.991, 249.554, 249.554, 250.546, 250.546, 251.588, 251.588, 251.989, 251.989, 252.006, 252.006, 252.146, 252.146, 252.439, 252.439, 251.712), ConCharter = c(9185, 11832.34663, 5650, 7278.47125, 5650, 7278.47125, 4713, 6071, 18567, 23918, 14800, 19065, 13800, 17777, 5117, 6591, 6450, 8309, 7500, 9661, 14170, 11000, 14000, 18035, 23240, 18041, 16492, 21245, 16492, 21245, 16400, 21126, 10305, 8000, 7700, 9919, 7000, 9017, 16700, 21513, 24347, 18900, 20500, 26408, 19350, 24927, 19200, 24733, 21899, 17000, 15500, 19967, 12500, 16102, 7000, 9017, 7000, 9017, 6500, 8373, 8373, 6500, 6800, 8759, 9339, 7250, 7750, 9983, 10305, 8000, 7800, 10048, 7750, 9983, 7500, 9661, 7250, 9339, 6700, 8631, 7750, 9983, 7600, 9790, 7600, 9790, 6900, 8888, 7000, 9017, 7300, 9404, 7200, 9275, 7450, 9597, 7500, 9661, 7700, 9919, 7900, 10176, 7600, 9790, 8400, 10821, 8400, 10821, 8500, 10949, 7600, 9790, 7500, 9661, 7800, 10048, 8000, 10305, 8100, 10434, 8200, 10563, 8200, 10563, 8300, 10692, 8500, 10949, 10100, 13011, 12238, 9500, 8300, 10692, 8500, 10949, 8700, 11207, 9100, 11722, 9500, 12238, 10000, 12882, 10900, 14041, 10600, 13655, 10200, 13139, 7000, 9017, 6698, 5200, 5150, 6634, 5200, 6698, 6700, 8631, 6300, 8115, 6100, 7858, 5700, 7342, 5700, 7342, 5700, 7342, 5700, 7342, 5300, 6827, 5200, 6698, 8200, 10563, 8200, 10563, 8700, 11207, 5500, 7085, 5700, 7342, 7500, 9661, 9400, 12109, 9500, 12238, 9100, 11722, 8000, 10305, 8200, 10563, 8900, 11465, 10800, 13912, 9500, 9500, 9000, 11594, 9900, 12753, 10500, 13526, 14170, 11000, 11100, 14299, 15973, 12400, 12300, 15845, 12000, 15458, 8500), NewBuildingPrice = c(53700000, 68360100, 48000000, 61104000, 41700000, 53084100, 35705882.35, 45453588.24, 40411764.71, 51444176.47, 41764705.88, 53166470.59, 43117647.06, 54888764.71, 35705882.35, 45453588.24, 35876470.59, 45670747.06, 36132352.94, 45996485.29, 46539382.35, 36558823.53, 36558823.53, 46539382.35, 47625176.47, 37411764.71, 38264705.88, 48710970.59, 39411764.71, 50171176.47, 44000000, 56012000, 57285000, 45000000, 45000000, 57285000, 45000000, 57285000, 44500000, 56648500, 59831000, 47000000, 49000000, 62377000, 50000000, 63650000, 46000000, 58558000, 58558000, 46000000, 46000000, 58558000, 45300000, 57666900, 44500000, 56648500, 41500000, 52829500, 41000000, 52193000, 52193000, 41000000, 44500000, 56648500, 56648500, 44500000, 44000000, 56012000, 56012000, 44000000, 44000000, 56012000, 43500000, 55375500, 42000000, 53466000, 42000000, 53466000, 41000000, 52193000, 40032258.06, 50961064.52, 40000000, 50920000, 40000000, 50920000, 38176470.59, 48598647.06, 38134328.36, 48545000, 41000000, 52193000, 41000000, 52193000, 38093750, 48493343.75, 40031250, 50959781.25, 38125000, 48533125, 38142857.14, 48555857.14, 40000000, 50920000, 41000000, 52193000, 41000000, 52193000, 41000000, 52193000, 40000000, 50920000, 40000000, 50920000, 41000000, 52193000, 41000000, 52193000, 41000000, 52193000, 41000000, 52193000, 41000000, 52193000, 41000000, 52193000, 41000000, 52193000, 41000000, 52193000, 52193000, 41000000, 41000000, 52193000, 41000000, 52193000, 41000000, 52193000, 41000000, 52193000, 41000000, 52193000, 41000000, 52193000, 41000000, 52193000, 41000000, 52193000, 41000000, 52193000, 41000000, 52193000, 50920000, 40000000, 40000000, 50920000, 40000000, 50920000, 41000000, 52193000, 41000000, 52193000, 40000000, 50920000, 40000000, 50920000, 40000000, 50920000, 40000000, 50920000, 40000000, 50920000, 40000000, 50920000, 40000000, 50920000, 41000000, 52193000, 41000000, 52193000, 42000000, 53466000, 40000000, 50920000, 40000000, 50920000, 40000000, 50920000, 40000000, 50920000, 41000000, 52193000, 41000000, 52193000, 41000000, 52193000, 41000000, 52193000, 42000000, 53466000, 42000000, 53466000, 42000000, 42000000, 42000000, 53466000, 42000000, 53466000, 43000000, 54739000, 54739000, 43000000, 43000000, 54739000, 53466000, 42000000, 42000000, 53466000, 42000000, 53466000, 42000000 ), HFO = c(265.5382699, 265.5382699, 345.6949275, 345.6949275, 445.9030797, 445.9030797, 481.5938406, 481.5938406, 474.6855046, 474.6855046, 494.0818807, 494.0818807, 508.073211, 508.073211, 466.3808877, 466.3808877, 469.3727355, 469.3727355, 485.209058, 485.209058, 462.8575634, 462.8575634, 447.1970109, 447.1970109, 446.4384058, 446.4384058, 454.4567482, 454.4567482, 447.6084692, 447.6084692, 535.3327982, 535.3327982, 661.2970183, 661.2970183, 681.5729358, 681.5729358, 663.2363303, 663.2363303, 600.8077982, 600.8077982, 642.783945, 642.783945, 673.2849541, 673.2849541, 653.8043578, 653.8043578, 662.915367, 662.915367, 671.2282569, 671.2282569, 660.6488532, 660.6488532, 660.7256881, 660.7256881, 708.3672018, 708.3672018, 648.4905963, 648.4905963, 624.0271186, 624.0271186, 616.5951036, 616.5951036, 725.7183486, 725.7183486, 739.9016514, 739.9016514, 726.0190367, 726.0190367, 682.0949541, 682.0949541, 600.1110092, 600.1110092, 611.887844, 611.887844, 661.6696798, 661.6696798, 667.6007533, 667.6007533, 634.8754708, 634.8754708, 615.2807895, 615.2807895, 608.9688034, 608.9688034, 615.4396368, 615.4396368, 658.9545669, 658.9545669, 634.4610169, 634.4610169, 614.9274809, 614.9274809, 606.5840954, 606.5840954, 606.4353877, 606.4353877, 609.7552923, 609.7552923, 616.1560484, 616.1560484, 618.7678942, 618.7678942, 614.3444685, 614.3444685, 511.3463592, 511.3463592, 454.7987654, 454.7987654, 365.7330402, 365.7330402, 616.6238987, 616.6238987, 610.7062225, 610.7062225, 603.4524609, 603.4524609, 604.0654362, 604.0654362, 609.8914773, 609.8914773, 600.3314088, 600.3314088, 592.171831, 592.171831, 576.8988663, 576.8988663, 292.1725441, 292.1725441, 246.4201183, 246.4201183, 229.0865385, 229.0865385, 184.9615385, 184.9615385, 342.7136076, 342.7136076, 333.5895619, 333.5895619, 336.7940292, 336.7940292, 370.07888, 370.07888, 360.3186667, 360.3186667, 316.5268617, 316.5268617, 255.4885216, 255.4885216, 249.9428571, 249.9428571, 155.6958042, 155.6958042, 285.6822917, 285.6822917, 276.8942308, 276.8942308, 323.1230769, 323.1230769, 161.2596154, 161.2596154, 176.7019231, 176.7019231, 191.8435897, 191.8435897, 228.8461538, 228.8461538, 245.875, 245.875, 249.4307692, 249.4307692, 241.2403846, 241.2403846, 258.3384615, 258.3384615, 337.5432692, 337.5432692, 340.9807692, 340.9807692, 373.3173077, 373.3173077, 371.4692308, 371.4692308, 326.7980769, 326.7980769, 309.1347633, 309.1347633, 314.2120726, 314.2120726, 309.9615385, 309.9615385, 299.1931818, 299.1931818, 308.4807692, 308.4807692, 316.7596154, 316.7596154, 333.7692308, 333.7692308, 395.0576923, 395.0576923, 495.5480769, 495.5480769, 469.5985183, 388.7115385, 379.6442308, 379.6442308, 374.0428571, 374.0428571, 393.9711538, 393.9711538, 440.5865385, 440.5865385, 449.6538462, 449.6538462, 455.9711538, 455.9711538, 453.4538462, 453.4538462, 465.5096154, 465.5096154, 386.3028846), IFO = c(286.9981846, 286.9981846, 361.859751, 361.859751, 460.747666, 460.747666, 496.6584232, 496.6584232, 489.6275789, 489.6275789, 507.4778947, 507.4778947, 523.4875789, 523.4875789, 483.5749481, 483.5749481, 490.6869813, 490.6869813, 502.8487552, 502.8487552, 480.3807054, 480.3807054, 469.1151452, 469.1151452, 465.9790456, 465.9790456, 474.4701763, 474.4701763, 465.2258817, 465.2258817, 556.2789474, 556.2789474, 684.9715789, 684.9715789, 705.7531579, 705.7531579, 687.4496421, 687.4496421, 620.2205263, 620.2205263, 667.03, 667.03, 697.0783158, 697.0783158, 677.82, 677.82, 688.0405263, 688.0405263, 696.2004211, 696.2004211, 685.1592105, 685.1592105, 683.9185263, 683.9185263, 729.9213158, 729.9213158, 678.1105263, 678.1105263, 660.056833, 660.056833, 655.233731, 655.233731, 748.6018421, 748.6018421, 763.4444211, 763.4444211, 747.6252632, 747.6252632, 706.8092105, 706.8092105, 626.1515789, 626.1515789, 635.915, 635.915, 683.818872, 683.818872, 704.5911063, 704.5911063, 670.9761388, 670.9761388, 650.054321, 650.054321, 643.8965432, 643.8965432, 652.5866834, 652.5866834, 693.9902386, 693.9902386, 673.625705, 673.625705, 649.965859, 649.965859, 644.3718822, 644.3718822, 644.1045035, 644.1045035, 645.8142606, 645.8142606, 654.8107981, 654.8107981, 657.1025485, 657.1025485, 651.8107417, 651.8107417, 563.7994152, 563.7994152, 506.9328358, 506.9328358, 430.4456301, 430.4456301, 649.5345052, 649.5345052, 645.6389974, 645.6389974, 640.6578249, 640.6578249, 646.264191, 646.264191, 646.4189189, 646.4189189, 638.6136364, 638.6136364, 635.7407303, 635.7407303, 624.032235, 624.032235, 348.2040248, 348.2040248, 284.7455808, 284.7455808, 262.2386364, 262.2386364, 219.9318182, 219.9318182, 377.4450464, 377.4450464, 373.4734968, 373.4734968, 374.8380802, 374.8380802, 408.3593443, 408.3593443, 403.8184426, 403.8184426, 364.935339, 364.935339, 304.8829225, 304.8829225, 289.5515511, 289.5515511, 181.3295455, 181.3295455, 317.329072, 317.329072, 309.0113636, 309.0113636, 351.1181818, 351.1181818, 186.4318182, 186.4318182, 205.5681818, 205.5681818, 221.92625, 221.92625, 256.3977273, 256.3977273, 279, 279, 286.2818182, 286.2818182, 278.2045455, 278.2045455, 293.3545455, 293.3545455, 362.599026, 362.599026, 368.8636364, 368.8636364, 406.1136364, 406.1136364, 407.1727273, 407.1727273, 355.2727273, 355.2727273, 339.004021, 339.004021, 343.7635732, 343.7635732, 340.4659091, 340.4659091, 335.6105372, 335.6105372, 345.3181818, 345.3181818, 356.8863636, 356.8863636, 367.1727273, 367.1727273, 430.4431818, 430.4431818, 521.25, 521.25, 492.6477273, 419.2954545, 411.1704545, 411.1704545, 410.2459416, 410.2459416, 422.9318182, 422.9318182, 465.5340909, 465.5340909, 478.2545455, 478.2545455, 486.5909091, 486.5909091, 483.9090909, 483.9090909, 496.2045455, 496.2045455, 407.7159091 ), MGO = c(453.2946429, 453.2946429, 477.9595238, 477.9595238, 591.0833333, 591.0833333, 623.2161905, 623.2161905, 692.458624, 692.458624, 711.9907946, 711.9907946, 742.2348837, 742.2348837, 614.5327381, 614.5327381, 643.2946429, 643.2946429, 686.0642857, 686.0642857, 657.7142857, 657.7142857, 651.0505952, 651.0505952, 643.5, 643.5, 656.4315476, 656.4315476, 657.985119, 657.985119, 786.6656977, 786.6656977, 922.6119186, 922.6119186, 956.1061047, 956.1061047, 930.755814, 930.755814, 866.1065891, 866.1065891, 938.9481589, 938.9481589, 971.4492248, 971.4492248, 937.7955426, 937.7955426, 941.8498062, 941.8498062, 948.1639535, 948.1639535, 937.3231589, 937.3231589, 933.9395349, 933.9395349, 959.095057, 959.095057, 1019.252193, 1019.252193, 991.6284211, 991.6284211, 985.2127193, 985.2127193, 979.1858365, 979.1858365, 1000.876426, 1000.876426, 985.2276616, 985.2276616, 948.5893536, 948.5893536, 872.8977186, 872.8977186, 896.7889734, 896.7889734, 998.4410526, 998.4410526, 1019.827632, 1019.827632, 994.0855263, 994.0855263, 979.0090201, 979.0090201, 967.2960674, 967.2960674, 983.7078652, 983.7078652, 1027.866228, 1027.866228, 977.4957895, 977.4957895, 942.8925439, 942.8925439, 930.0765464, 930.0765464, 937.4226804, 937.4226804, 963.5524055, 963.5524055, 980.5042955, 980.5042955, 985.0053763, 985.0053763, 968.5701149, 968.5701149, 837.2136986, 837.2136986, 788.818662, 788.818662, 685.173913, 685.173913, 972.9114379, 972.9114379, 962.829902, 962.829902, 964.2464859, 964.2464859, 961.3871486, 961.3871486, 957.4074074, 957.4074074, 942.0369198, 942.0369198, 931.2835498, 931.2835498, 908.5083333, 908.5083333, 585.1183168, 585.1183168, 491.1208333, 491.1208333, 476.0892857, 476.0892857, 408.6428571, 408.6428571, 621.1639851, 621.1639851, 605.6785714, 605.6785714, 598.0431548, 598.0431548, 642.6857143, 642.6857143, 630.0462963, 630.0462963, 571.7781746, 571.7781746, 507.6554487, 507.6554487, 505.9642857, 505.9642857, 351.8231293, 351.8231293, 497.1383929, 497.1383929, 469.3392857, 469.3392857, 508.4142857, 508.4142857, 340.9642857, 340.9642857, 371.6964286, 371.6964286, 386.1060714, 386.1060714, 449.6785714, 449.6785714, 475.9464286, 475.9464286, 453.5428571, 453.5428571, 441.3214286, 441.3214286, 452.0428571, 452.0428571, 517.7104592, 517.7104592, 547.625, 547.625, 577.0892857, 577.0892857, 586.4142857, 586.4142857, 520.0535714, 520.0535714, 500.5695055, 500.5695055, 510.5292659, 510.5292659, 491.9642857, 491.9642857, 470.5665584, 470.5665584, 480.5892857, 480.5892857, 501, 501, 541, 541, 627.3214286, 627.3214286, 746.9642857, 746.9642857, 699.3928571, 600.1607143, 613.4285714, 613.4285714, 611.1628401, 611.1628401, 648.0535714, 648.0535714, 694.8214286, 694.8214286, 686.0571429, 686.0571429, 678.6428571, 678.6428571, 683.2142857, 683.2142857, 720.8214286, 720.8214286, 597.1875)), class "data.frame", row.names = c(NA, -221L)) [[alternative HTML version deleted]]