How To Extracting Summary Stats for Sepwise Regression

classic Classic list List threaded Threaded
1 message Options
Reply | Threaded
Open this post in threaded view
|

How To Extracting Summary Stats for Sepwise Regression

PaulJr
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]]

______________________________________________
[hidden email] mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.