Creates polynomial variables, to be used in regression. Will create polynomials of degree less than
or equal to the degree
specified, and will mean center variables by default.
Usage
polynomial(x, degree = 2, center = mean(x, na.rm = TRUE))
Arguments
- x
variable used to create the polynomials.
- degree
the maximum degree polynomial to be returned. Polynomials of degree <=
degree
will be returned.- center
the value to center the polynomials at.
Examples
# Reading in a dataset
data(mri)
# Create a polynomial on ldl
polynomial(mri$ldl, degree=3)
#> Linear(ctr) Square(ctr) Cube(ctr)
#> [1,] 9.1972414 8.458925e+01 7.779877e+02
#> [2,] -41.8027586 1.747471e+03 -7.304909e+04
#> [3,] -10.8027586 1.166996e+02 -1.260678e+03
#> [4,] -64.8027586 4.199398e+03 -2.721325e+05
#> [5,] 22.1972414 4.927175e+02 1.093697e+04
#> [6,] 37.1972414 1.383635e+03 5.146740e+04
#> [7,] -24.8027586 6.151768e+02 -1.525808e+04
#> [8,] -9.8027586 9.609408e+01 -9.419870e+02
#> [9,] -1.8027586 3.249939e+00 -5.858855e+00
#> [10,] -15.8027586 2.497272e+02 -3.946378e+03
#> [11,] 10.1972414 1.039837e+02 1.060347e+03
#> [12,] -17.8027586 3.169382e+02 -5.642375e+03
#> [13,] 80.1972414 6.431598e+03 5.157964e+05
#> [14,] 4.1972414 1.761684e+01 7.394211e+01
#> [15,] -4.8027586 2.306649e+01 -1.107828e+02
#> [16,] -30.8027586 9.488099e+02 -2.922596e+04
#> [17,] -47.8027586 2.285104e+03 -1.092343e+05
#> [18,] -15.8027586 2.497272e+02 -3.946378e+03
#> [19,] 10.1972414 1.039837e+02 1.060347e+03
#> [20,] 34.1972414 1.169451e+03 3.999201e+04
#> [21,] 4.1972414 1.761684e+01 7.394211e+01
#> [22,] -42.8027586 1.832076e+03 -7.841791e+04
#> [23,] 28.1972414 7.950844e+02 2.241919e+04
#> [24,] -55.8027586 3.113948e+03 -1.737669e+05
#> [25,] -50.8027586 2.580920e+03 -1.311179e+05
#> [26,] -11.8027586 1.393051e+02 -1.644185e+03
#> [27,] 13.1972414 1.741672e+02 2.298526e+03
#> [28,] -44.8027586 2.007287e+03 -8.993200e+04
#> [29,] -4.8027586 2.306649e+01 -1.107828e+02
#> [30,] -20.8027586 4.327548e+02 -9.002493e+03
#> [31,] -3.8027586 1.446097e+01 -5.499159e+01
#> [32,] 9.1972414 8.458925e+01 7.779877e+02
#> [33,] 10.1972414 1.039837e+02 1.060347e+03
#> [34,] 62.1972414 3.868497e+03 2.406098e+05
#> [35,] 19.1972414 3.685341e+02 7.074838e+03
#> [36,] -26.8027586 7.183879e+02 -1.925478e+04
#> [37,] 45.1972414 2.042791e+03 9.232850e+04
#> [38,] -42.8027586 1.832076e+03 -7.841791e+04
#> [39,] 29.1972414 8.524789e+02 2.489003e+04
#> [40,] 17.1972414 2.957451e+02 5.086000e+03
#> [41,] 22.1972414 4.927175e+02 1.093697e+04
#> [42,] 36.1972414 1.310240e+03 4.742708e+04
#> [43,] 56.1972414 3.158130e+03 1.774782e+05
#> [44,] -43.8027586 1.918682e+03 -8.404355e+04
#> [45,] 25.1972414 6.349010e+02 1.599775e+04
#> [46,] -13.8027586 1.905161e+02 -2.629648e+03
#> [47,] 33.1972414 1.102057e+03 3.658525e+04
#> [48,] -50.8027586 2.580920e+03 -1.311179e+05
#> [49,] -27.8027586 7.729934e+02 -2.149135e+04
#> [50,] 63.1972414 3.993891e+03 2.524029e+05
#> [51,] -2.8027586 7.855456e+00 -2.201695e+01
#> [52,] 101.1972414 1.024088e+04 1.036349e+06
#> [53,] -52.8027586 2.788131e+03 -1.472210e+05
#> [54,] -26.8027586 7.183879e+02 -1.925478e+04
#> [55,] 90.1972414 8.135542e+03 7.338035e+05
#> [56,] -21.8027586 4.753603e+02 -1.036417e+04
#> [57,] -45.8027586 2.097893e+03 -9.608927e+04
#> [58,] -36.8027586 1.354443e+03 -4.984724e+04
#> [59,] -17.8027586 3.169382e+02 -5.642375e+03
#> [60,] 21.1972414 4.493230e+02 9.524409e+03
#> [61,] 17.1972414 2.957451e+02 5.086000e+03
#> [62,] 25.1972414 6.349010e+02 1.599775e+04
#> [63,] 25.1972414 6.349010e+02 1.599775e+04
#> [64,] 12.1972414 1.487727e+02 1.814616e+03
#> [65,] -12.8027586 1.639106e+02 -2.098508e+03
#> [66,] 20.1972414 4.079286e+02 8.239032e+03
#> [67,] 23.1972414 5.381120e+02 1.248271e+04
#> [68,] -53.8027586 2.894737e+03 -1.557448e+05
#> [69,] 6.1972414 3.840580e+01 2.380100e+02
#> [70,] -31.8027586 1.011415e+03 -3.216580e+04
#> [71,] -40.8027586 1.664865e+03 -6.793109e+04
#> [72,] -46.8027586 2.190498e+03 -1.025214e+05
#> [73,] 9.1972414 8.458925e+01 7.779877e+02
#> [74,] 13.1972414 1.741672e+02 2.298526e+03
#> [75,] -40.8027586 1.664865e+03 -6.793109e+04
#> [76,] -22.8027586 5.199658e+02 -1.185665e+04
#> [77,] 1.1972414 1.433387e+00 1.716110e+00
#> [78,] 14.1972414 2.015617e+02 2.861620e+03
#> [79,] -24.8027586 6.151768e+02 -1.525808e+04
#> [80,] -4.8027586 2.306649e+01 -1.107828e+02
#> [81,] 20.1972414 4.079286e+02 8.239032e+03
#> [82,] -38.8027586 1.505654e+03 -5.842353e+04
#> [83,] 24.1972414 5.855065e+02 1.416764e+04
#> [84,] 8.1972414 6.719477e+01 5.508117e+02
#> [85,] -45.8027586 2.097893e+03 -9.608927e+04
#> [86,] 3.1972414 1.022235e+01 3.268333e+01
#> [87,] 28.1972414 7.950844e+02 2.241919e+04
#> [88,] -26.8027586 7.183879e+02 -1.925478e+04
#> [89,] -8.8027586 7.748856e+01 -6.821131e+02
#> [90,] 47.1972414 2.227580e+03 1.051356e+05
#> [91,] -5.8027586 3.367201e+01 -1.953905e+02
#> [92,] -0.8027586 6.444214e-01 -5.173148e-01
#> [93,] -15.8027586 2.497272e+02 -3.946378e+03
#> [94,] -40.8027586 1.664865e+03 -6.793109e+04
#> [95,] -55.8027586 3.113948e+03 -1.737669e+05
#> [96,] 41.1972414 1.697213e+03 6.992048e+04
#> [97,] -3.8027586 1.446097e+01 -5.499159e+01
#> [98,] 37.1972414 1.383635e+03 5.146740e+04
#> [99,] -25.8027586 6.657824e+02 -1.717902e+04
#> [100,] 32.1972414 1.036662e+03 3.337767e+04
#> [101,] -10.8027586 1.166996e+02 -1.260678e+03
#> [102,] -28.8027586 8.295989e+02 -2.389474e+04
#> [103,] -41.8027586 1.747471e+03 -7.304909e+04
#> [104,] -44.8027586 2.007287e+03 -8.993200e+04
#> [105,] 25.1972414 6.349010e+02 1.599775e+04
#> [106,] 65.1972414 4.250680e+03 2.771326e+05
#> [107,] 41.1972414 1.697213e+03 6.992048e+04
#> [108,] -8.8027586 7.748856e+01 -6.821131e+02
#> [109,] -21.8027586 4.753603e+02 -1.036417e+04
#> [110,] -12.8027586 1.639106e+02 -2.098508e+03
#> [111,] -7.8027586 6.088304e+01 -4.750557e+02
#> [112,] -5.8027586 3.367201e+01 -1.953905e+02
#> [113,] -27.8027586 7.729934e+02 -2.149135e+04
#> [114,] 0.1972414 3.890416e-02 7.673511e-03
#> [115,] 16.1972414 2.623506e+02 4.249356e+03
#> [116,] 20.1972414 4.079286e+02 8.239032e+03
#> [117,] -42.8027586 1.832076e+03 -7.841791e+04
#> [118,] 7.1972414 5.180028e+01 3.728191e+02
#> [119,] 2.1972414 4.827870e+00 1.060800e+01
#> [120,] 4.1972414 1.761684e+01 7.394211e+01
#> [121,] -35.8027586 1.281838e+03 -4.589332e+04
#> [122,] 10.1972414 1.039837e+02 1.060347e+03
#> [123,] -18.8027586 3.535437e+02 -6.647597e+03
#> [124,] -2.8027586 7.855456e+00 -2.201695e+01
#> [125,] -21.8027586 4.753603e+02 -1.036417e+04
#> [126,] -23.8027586 5.665713e+02 -1.348596e+04
#> [127,] -25.8027586 6.657824e+02 -1.717902e+04
#> [128,] 51.1972414 2.621158e+03 1.341960e+05
#> [129,] -39.8027586 1.584260e+03 -6.305790e+04
#> [130,] NA NA NA
#> [131,] 38.1972414 1.459029e+03 5.573089e+04
#> [132,] 30.1972414 9.118734e+02 2.753606e+04
#> [133,] -17.8027586 3.169382e+02 -5.642375e+03
#> [134,] -72.8027586 5.300242e+03 -3.858722e+05
#> [135,] 24.1972414 5.855065e+02 1.416764e+04
#> [136,] 69.1972414 4.788258e+03 3.313343e+05
#> [137,] 8.1972414 6.719477e+01 5.508117e+02
#> [138,] 16.1972414 2.623506e+02 4.249356e+03
#> [139,] -37.8027586 1.429049e+03 -5.402198e+04
#> [140,] 4.1972414 1.761684e+01 7.394211e+01
#> [141,] 48.1972414 2.322974e+03 1.119609e+05
#> [142,] 25.1972414 6.349010e+02 1.599775e+04
#> [143,] 22.1972414 4.927175e+02 1.093697e+04
#> [144,] -6.8027586 4.627752e+01 -3.148148e+02
#> [145,] 35.1972414 1.238846e+03 4.360395e+04
#> [146,] 10.1972414 1.039837e+02 1.060347e+03
#> [147,] 30.1972414 9.118734e+02 2.753606e+04
#> [148,] -26.8027586 7.183879e+02 -1.925478e+04
#> [149,] -11.8027586 1.393051e+02 -1.644185e+03
#> [150,] 12.1972414 1.487727e+02 1.814616e+03
#> [151,] -21.8027586 4.753603e+02 -1.036417e+04
#> [152,] -12.8027586 1.639106e+02 -2.098508e+03
#> [153,] 7.1972414 5.180028e+01 3.728191e+02
#> [154,] -44.8027586 2.007287e+03 -8.993200e+04
#> [155,] -47.8027586 2.285104e+03 -1.092343e+05
#> [156,] -10.8027586 1.166996e+02 -1.260678e+03
#> [157,] 15.1972414 2.309561e+02 3.509896e+03
#> [158,] NA NA NA
#> [159,] -47.8027586 2.285104e+03 -1.092343e+05
#> [160,] -24.8027586 6.151768e+02 -1.525808e+04
#> [161,] 29.1972414 8.524789e+02 2.489003e+04
#> [162,] 10.1972414 1.039837e+02 1.060347e+03
#> [163,] 5.1972414 2.701132e+01 1.403843e+02
#> [164,] 12.1972414 1.487727e+02 1.814616e+03
#> [165,] 14.1972414 2.015617e+02 2.861620e+03
#> [166,] -10.8027586 1.166996e+02 -1.260678e+03
#> [167,] 18.1972414 3.311396e+02 6.025827e+03
#> [168,] 0.1972414 3.890416e-02 7.673511e-03
#> [169,] 91.1972414 8.316937e+03 7.584817e+05
#> [170,] 19.1972414 3.685341e+02 7.074838e+03
#> [171,] 1.1972414 1.433387e+00 1.716110e+00
#> [172,] -49.8027586 2.480315e+03 -1.235265e+05
#> [173,] NA NA NA
#> [174,] -74.8027586 5.595453e+03 -4.185553e+05
#> [175,] 6.1972414 3.840580e+01 2.380100e+02
#> [176,] -11.8027586 1.393051e+02 -1.644185e+03
#> [177,] -64.8027586 4.199398e+03 -2.721325e+05
#> [178,] 23.1972414 5.381120e+02 1.248271e+04
#> [179,] 11.1972414 1.253782e+02 1.403890e+03
#> [180,] 1.1972414 1.433387e+00 1.716110e+00
#> [181,] 7.1972414 5.180028e+01 3.728191e+02
#> [182,] 7.1972414 5.180028e+01 3.728191e+02
#> [183,] 18.1972414 3.311396e+02 6.025827e+03
#> [184,] -16.8027586 2.823327e+02 -4.743968e+03
#> [185,] 26.1972414 6.862955e+02 1.797905e+04
#> [186,] -30.8027586 9.488099e+02 -2.922596e+04
#> [187,] 3.1972414 1.022235e+01 3.268333e+01
#> [188,] NA NA NA
#> [189,] 2.1972414 4.827870e+00 1.060800e+01
#> [190,] -26.8027586 7.183879e+02 -1.925478e+04
#> [191,] 13.1972414 1.741672e+02 2.298526e+03
#> [192,] 16.1972414 2.623506e+02 4.249356e+03
#> [193,] -22.8027586 5.199658e+02 -1.185665e+04
#> [194,] -11.8027586 1.393051e+02 -1.644185e+03
#> [195,] 59.1972414 3.504313e+03 2.074457e+05
#> [196,] -1.8027586 3.249939e+00 -5.858855e+00
#> [197,] -27.8027586 7.729934e+02 -2.149135e+04
#> [198,] -24.8027586 6.151768e+02 -1.525808e+04
#> [199,] 65.1972414 4.250680e+03 2.771326e+05
#> [200,] -7.8027586 6.088304e+01 -4.750557e+02
#> [201,] -28.8027586 8.295989e+02 -2.389474e+04
#> [202,] -8.8027586 7.748856e+01 -6.821131e+02
#> [203,] -30.8027586 9.488099e+02 -2.922596e+04
#> [204,] 9.1972414 8.458925e+01 7.779877e+02
#> [205,] 80.1972414 6.431598e+03 5.157964e+05
#> [206,] 24.1972414 5.855065e+02 1.416764e+04
#> [207,] -3.8027586 1.446097e+01 -5.499159e+01
#> [208,] 14.1972414 2.015617e+02 2.861620e+03
#> [209,] -34.8027586 1.211232e+03 -4.215422e+04
#> [210,] -4.8027586 2.306649e+01 -1.107828e+02
#> [211,] 34.1972414 1.169451e+03 3.999201e+04
#> [212,] 36.1972414 1.310240e+03 4.742708e+04
#> [213,] 62.1972414 3.868497e+03 2.406098e+05
#> [214,] -30.8027586 9.488099e+02 -2.922596e+04
#> [215,] 38.1972414 1.459029e+03 5.573089e+04
#> [216,] 4.1972414 1.761684e+01 7.394211e+01
#> [217,] 28.1972414 7.950844e+02 2.241919e+04
#> [218,] NA NA NA
#> [219,] 10.1972414 1.039837e+02 1.060347e+03
#> [220,] -30.8027586 9.488099e+02 -2.922596e+04
#> [221,] 121.1972414 1.468877e+04 1.780239e+06
#> [222,] 18.1972414 3.311396e+02 6.025827e+03
#> [223,] -42.8027586 1.832076e+03 -7.841791e+04
#> [224,] -3.8027586 1.446097e+01 -5.499159e+01
#> [225,] 31.1972414 9.732679e+02 3.036327e+04
#> [226,] -18.8027586 3.535437e+02 -6.647597e+03
#> [227,] 5.1972414 2.701132e+01 1.403843e+02
#> [228,] 14.1972414 2.015617e+02 2.861620e+03
#> [229,] 74.1972414 5.505231e+03 4.084729e+05
#> [230,] -30.8027586 9.488099e+02 -2.922596e+04
#> [231,] 20.1972414 4.079286e+02 8.239032e+03
#> [232,] -23.8027586 5.665713e+02 -1.348596e+04
#> [233,] 11.1972414 1.253782e+02 1.403890e+03
#> [234,] 24.1972414 5.855065e+02 1.416764e+04
#> [235,] 22.1972414 4.927175e+02 1.093697e+04
#> [236,] -53.8027586 2.894737e+03 -1.557448e+05
#> [237,] 16.1972414 2.623506e+02 4.249356e+03
#> [238,] -34.8027586 1.211232e+03 -4.215422e+04
#> [239,] 27.1972414 7.396899e+02 2.011753e+04
#> [240,] 82.1972414 6.756386e+03 5.553563e+05
#> [241,] -7.8027586 6.088304e+01 -4.750557e+02
#> [242,] 24.1972414 5.855065e+02 1.416764e+04
#> [243,] 25.1972414 6.349010e+02 1.599775e+04
#> [244,] 1.1972414 1.433387e+00 1.716110e+00
#> [245,] -25.8027586 6.657824e+02 -1.717902e+04
#> [246,] -11.8027586 1.393051e+02 -1.644185e+03
#> [247,] -13.8027586 1.905161e+02 -2.629648e+03
#> [248,] 3.1972414 1.022235e+01 3.268333e+01
#> [249,] 53.1972414 2.829946e+03 1.505453e+05
#> [250,] 8.1972414 6.719477e+01 5.508117e+02
#> [251,] -10.8027586 1.166996e+02 -1.260678e+03
#> [252,] -17.8027586 3.169382e+02 -5.642375e+03
#> [253,] -26.8027586 7.183879e+02 -1.925478e+04
#> [254,] 81.1972414 6.592992e+03 5.353328e+05
#> [255,] -48.8027586 2.381709e+03 -1.162340e+05
#> [256,] 15.1972414 2.309561e+02 3.509896e+03
#> [257,] 4.1972414 1.761684e+01 7.394211e+01
#> [258,] 36.1972414 1.310240e+03 4.742708e+04
#> [259,] 1.1972414 1.433387e+00 1.716110e+00
#> [260,] -33.8027586 1.142626e+03 -3.862393e+04
#> [261,] NA NA NA
#> [262,] -25.8027586 6.657824e+02 -1.717902e+04
#> [263,] 61.1972414 3.745102e+03 2.291899e+05
#> [264,] -12.8027586 1.639106e+02 -2.098508e+03
#> [265,] 2.1972414 4.827870e+00 1.060800e+01
#> [266,] 57.1972414 3.271524e+03 1.871222e+05
#> [267,] NA NA NA
#> [268,] -15.8027586 2.497272e+02 -3.946378e+03
#> [269,] -22.8027586 5.199658e+02 -1.185665e+04
#> [270,] 30.1972414 9.118734e+02 2.753606e+04
#> [271,] -26.8027586 7.183879e+02 -1.925478e+04
#> [272,] 29.1972414 8.524789e+02 2.489003e+04
#> [273,] 55.1972414 3.046735e+03 1.681714e+05
#> [274,] -7.8027586 6.088304e+01 -4.750557e+02
#> [275,] 6.1972414 3.840580e+01 2.380100e+02
#> [276,] 22.1972414 4.927175e+02 1.093697e+04
#> [277,] -35.8027586 1.281838e+03 -4.589332e+04
#> [278,] 5.1972414 2.701132e+01 1.403843e+02
#> [279,] 12.1972414 1.487727e+02 1.814616e+03
#> [280,] -16.8027586 2.823327e+02 -4.743968e+03
#> [281,] 10.1972414 1.039837e+02 1.060347e+03
#> [282,] 13.1972414 1.741672e+02 2.298526e+03
#> [283,] -34.8027586 1.211232e+03 -4.215422e+04
#> [284,] 50.1972414 2.519763e+03 1.264852e+05
#> [285,] 22.1972414 4.927175e+02 1.093697e+04
#> [286,] -38.8027586 1.505654e+03 -5.842353e+04
#> [287,] -12.8027586 1.639106e+02 -2.098508e+03
#> [288,] 29.1972414 8.524789e+02 2.489003e+04
#> [289,] 59.1972414 3.504313e+03 2.074457e+05
#> [290,] 11.1972414 1.253782e+02 1.403890e+03
#> [291,] 11.1972414 1.253782e+02 1.403890e+03
#> [292,] 43.1972414 1.866002e+03 8.060612e+04
#> [293,] 80.1972414 6.431598e+03 5.157964e+05
#> [294,] -37.8027586 1.429049e+03 -5.402198e+04
#> [295,] 13.1972414 1.741672e+02 2.298526e+03
#> [296,] 56.1972414 3.158130e+03 1.774782e+05
#> [297,] 11.1972414 1.253782e+02 1.403890e+03
#> [298,] 59.1972414 3.504313e+03 2.074457e+05
#> [299,] -9.8027586 9.609408e+01 -9.419870e+02
#> [300,] -11.8027586 1.393051e+02 -1.644185e+03
#> [301,] -56.8027586 3.226553e+03 -1.832771e+05
#> [302,] -6.8027586 4.627752e+01 -3.148148e+02
#> [303,] -43.8027586 1.918682e+03 -8.404355e+04
#> [304,] 4.1972414 1.761684e+01 7.394211e+01
#> [305,] -15.8027586 2.497272e+02 -3.946378e+03
#> [306,] 36.1972414 1.310240e+03 4.742708e+04
#> [307,] -12.8027586 1.639106e+02 -2.098508e+03
#> [308,] -32.8027586 1.076021e+03 -3.529646e+04
#> [309,] -7.8027586 6.088304e+01 -4.750557e+02
#> [310,] 70.1972414 4.927653e+03 3.459076e+05
#> [311,] -47.8027586 2.285104e+03 -1.092343e+05
#> [312,] -14.8027586 2.191217e+02 -3.243605e+03
#> [313,] -54.8027586 3.003342e+03 -1.645914e+05
#> [314,] 15.1972414 2.309561e+02 3.509896e+03
#> [315,] -13.8027586 1.905161e+02 -2.629648e+03
#> [316,] 69.1972414 4.788258e+03 3.313343e+05
#> [317,] -9.8027586 9.609408e+01 -9.419870e+02
#> [318,] 18.1972414 3.311396e+02 6.025827e+03
#> [319,] 47.1972414 2.227580e+03 1.051356e+05
#> [320,] 27.1972414 7.396899e+02 2.011753e+04
#> [321,] 10.1972414 1.039837e+02 1.060347e+03
#> [322,] 56.1972414 3.158130e+03 1.774782e+05
#> [323,] -35.8027586 1.281838e+03 -4.589332e+04
#> [324,] 2.1972414 4.827870e+00 1.060800e+01
#> [325,] 62.1972414 3.868497e+03 2.406098e+05
#> [326,] 58.1972414 3.386919e+03 1.971093e+05
#> [327,] -4.8027586 2.306649e+01 -1.107828e+02
#> [328,] -0.8027586 6.444214e-01 -5.173148e-01
#> [329,] -12.8027586 1.639106e+02 -2.098508e+03
#> [330,] 8.1972414 6.719477e+01 5.508117e+02
#> [331,] 14.1972414 2.015617e+02 2.861620e+03
#> [332,] -45.8027586 2.097893e+03 -9.608927e+04
#> [333,] 45.1972414 2.042791e+03 9.232850e+04
#> [334,] -44.8027586 2.007287e+03 -8.993200e+04
#> [335,] -13.8027586 1.905161e+02 -2.629648e+03
#> [336,] -9.8027586 9.609408e+01 -9.419870e+02
#> [337,] 32.1972414 1.036662e+03 3.337767e+04
#> [338,] -19.8027586 3.921492e+02 -7.765637e+03
#> [339,] 32.1972414 1.036662e+03 3.337767e+04
#> [340,] -34.8027586 1.211232e+03 -4.215422e+04
#> [341,] 6.1972414 3.840580e+01 2.380100e+02
#> [342,] -37.8027586 1.429049e+03 -5.402198e+04
#> [343,] -12.8027586 1.639106e+02 -2.098508e+03
#> [344,] 21.1972414 4.493230e+02 9.524409e+03
#> [345,] 16.1972414 2.623506e+02 4.249356e+03
#> [346,] -58.8027586 3.457764e+03 -2.033261e+05
#> [347,] 35.1972414 1.238846e+03 4.360395e+04
#> [348,] 52.1972414 2.724552e+03 1.422141e+05
#> [349,] -17.8027586 3.169382e+02 -5.642375e+03
#> [350,] -38.8027586 1.505654e+03 -5.842353e+04
#> [351,] 11.1972414 1.253782e+02 1.403890e+03
#> [352,] 34.1972414 1.169451e+03 3.999201e+04
#> [353,] 1.1972414 1.433387e+00 1.716110e+00
#> [354,] -29.8027586 8.882044e+02 -2.647094e+04
#> [355,] -15.8027586 2.497272e+02 -3.946378e+03
#> [356,] 49.1972414 2.420369e+03 1.190755e+05
#> [357,] 45.1972414 2.042791e+03 9.232850e+04
#> [358,] 5.1972414 2.701132e+01 1.403843e+02
#> [359,] -32.8027586 1.076021e+03 -3.529646e+04
#> [360,] -17.8027586 3.169382e+02 -5.642375e+03
#> [361,] -7.8027586 6.088304e+01 -4.750557e+02
#> [362,] -9.8027586 9.609408e+01 -9.419870e+02
#> [363,] -68.8027586 4.733820e+03 -3.256998e+05
#> [364,] 1.1972414 1.433387e+00 1.716110e+00
#> [365,] 55.1972414 3.046735e+03 1.681714e+05
#> [366,] -17.8027586 3.169382e+02 -5.642375e+03
#> [367,] 7.1972414 5.180028e+01 3.728191e+02
#> [368,] 16.1972414 2.623506e+02 4.249356e+03
#> [369,] 8.1972414 6.719477e+01 5.508117e+02
#> [370,] 0.1972414 3.890416e-02 7.673511e-03
#> [371,] 38.1972414 1.459029e+03 5.573089e+04
#> [372,] -48.8027586 2.381709e+03 -1.162340e+05
#> [373,] -28.8027586 8.295989e+02 -2.389474e+04
#> [374,] -70.8027586 5.013031e+03 -3.549364e+05
#> [375,] -32.8027586 1.076021e+03 -3.529646e+04
#> [376,] -57.8027586 3.341159e+03 -1.931282e+05
#> [377,] 10.1972414 1.039837e+02 1.060347e+03
#> [378,] -1.8027586 3.249939e+00 -5.858855e+00
#> [379,] 7.1972414 5.180028e+01 3.728191e+02
#> [380,] -4.8027586 2.306649e+01 -1.107828e+02
#> [381,] -45.8027586 2.097893e+03 -9.608927e+04
#> [382,] 17.1972414 2.957451e+02 5.086000e+03
#> [383,] 16.1972414 2.623506e+02 4.249356e+03
#> [384,] -48.8027586 2.381709e+03 -1.162340e+05
#> [385,] -53.8027586 2.894737e+03 -1.557448e+05
#> [386,] -36.8027586 1.354443e+03 -4.984724e+04
#> [387,] -48.8027586 2.381709e+03 -1.162340e+05
#> [388,] -32.8027586 1.076021e+03 -3.529646e+04
#> [389,] 75.1972414 5.654625e+03 4.252122e+05
#> [390,] 55.1972414 3.046735e+03 1.681714e+05
#> [391,] -55.8027586 3.113948e+03 -1.737669e+05
#> [392,] -52.8027586 2.788131e+03 -1.472210e+05
#> [393,] -20.8027586 4.327548e+02 -9.002493e+03
#> [394,] 8.1972414 6.719477e+01 5.508117e+02
#> [395,] -18.8027586 3.535437e+02 -6.647597e+03
#> [396,] 94.1972414 8.873120e+03 8.358235e+05
#> [397,] 40.1972414 1.615818e+03 6.495143e+04
#> [398,] 7.1972414 5.180028e+01 3.728191e+02
#> [399,] 99.1972414 9.840093e+03 9.761101e+05
#> [400,] 6.1972414 3.840580e+01 2.380100e+02
#> [401,] 17.1972414 2.957451e+02 5.086000e+03
#> [402,] -41.8027586 1.747471e+03 -7.304909e+04
#> [403,] 33.1972414 1.102057e+03 3.658525e+04
#> [404,] -11.8027586 1.393051e+02 -1.644185e+03
#> [405,] -50.8027586 2.580920e+03 -1.311179e+05
#> [406,] -14.8027586 2.191217e+02 -3.243605e+03
#> [407,] -11.8027586 1.393051e+02 -1.644185e+03
#> [408,] 44.1972414 1.953396e+03 8.633472e+04
#> [409,] -16.8027586 2.823327e+02 -4.743968e+03
#> [410,] 32.1972414 1.036662e+03 3.337767e+04
#> [411,] 47.1972414 2.227580e+03 1.051356e+05
#> [412,] -31.8027586 1.011415e+03 -3.216580e+04
#> [413,] 27.1972414 7.396899e+02 2.011753e+04
#> [414,] 0.1972414 3.890416e-02 7.673511e-03
#> [415,] 6.1972414 3.840580e+01 2.380100e+02
#> [416,] -10.8027586 1.166996e+02 -1.260678e+03
#> [417,] -35.8027586 1.281838e+03 -4.589332e+04
#> [418,] 2.1972414 4.827870e+00 1.060800e+01
#> [419,] 19.1972414 3.685341e+02 7.074838e+03
#> [420,] 18.1972414 3.311396e+02 6.025827e+03
#> [421,] 13.1972414 1.741672e+02 2.298526e+03
#> [422,] 21.1972414 4.493230e+02 9.524409e+03
#> [423,] -44.8027586 2.007287e+03 -8.993200e+04
#> [424,] 69.1972414 4.788258e+03 3.313343e+05
#> [425,] 35.1972414 1.238846e+03 4.360395e+04
#> [426,] -28.8027586 8.295989e+02 -2.389474e+04
#> [427,] -19.8027586 3.921492e+02 -7.765637e+03
#> [428,] -73.8027586 5.446847e+03 -4.019923e+05
#> [429,] 20.1972414 4.079286e+02 8.239032e+03
#> [430,] 2.1972414 4.827870e+00 1.060800e+01
#> [431,] -3.8027586 1.446097e+01 -5.499159e+01
#> [432,] 7.1972414 5.180028e+01 3.728191e+02
#> [433,] -14.8027586 2.191217e+02 -3.243605e+03
#> [434,] 30.1972414 9.118734e+02 2.753606e+04
#> [435,] -16.8027586 2.823327e+02 -4.743968e+03
#> [436,] 13.1972414 1.741672e+02 2.298526e+03
#> [437,] -11.8027586 1.393051e+02 -1.644185e+03
#> [438,] -114.8027586 1.317967e+04 -1.513063e+06
#> [439,] 0.1972414 3.890416e-02 7.673511e-03
#> [440,] -4.8027586 2.306649e+01 -1.107828e+02
#> [441,] 7.1972414 5.180028e+01 3.728191e+02
#> [442,] -86.8027586 7.534719e+03 -6.540344e+05
#> [443,] -19.8027586 3.921492e+02 -7.765637e+03
#> [444,] 11.1972414 1.253782e+02 1.403890e+03
#> [445,] 30.1972414 9.118734e+02 2.753606e+04
#> [446,] 36.1972414 1.310240e+03 4.742708e+04
#> [447,] -15.8027586 2.497272e+02 -3.946378e+03
#> [448,] -3.8027586 1.446097e+01 -5.499159e+01
#> [449,] 21.1972414 4.493230e+02 9.524409e+03
#> [450,] 43.1972414 1.866002e+03 8.060612e+04
#> [451,] 40.1972414 1.615818e+03 6.495143e+04
#> [452,] -31.8027586 1.011415e+03 -3.216580e+04
#> [453,] -33.8027586 1.142626e+03 -3.862393e+04
#> [454,] -9.8027586 9.609408e+01 -9.419870e+02
#> [455,] -29.8027586 8.882044e+02 -2.647094e+04
#> [456,] -36.8027586 1.354443e+03 -4.984724e+04
#> [457,] 13.1972414 1.741672e+02 2.298526e+03
#> [458,] 61.1972414 3.745102e+03 2.291899e+05
#> [459,] 47.1972414 2.227580e+03 1.051356e+05
#> [460,] -35.8027586 1.281838e+03 -4.589332e+04
#> [461,] -23.8027586 5.665713e+02 -1.348596e+04
#> [462,] 18.1972414 3.311396e+02 6.025827e+03
#> [463,] -35.8027586 1.281838e+03 -4.589332e+04
#> [464,] 42.1972414 1.780607e+03 7.513671e+04
#> [465,] -14.8027586 2.191217e+02 -3.243605e+03
#> [466,] -24.8027586 6.151768e+02 -1.525808e+04
#> [467,] -23.8027586 5.665713e+02 -1.348596e+04
#> [468,] -55.8027586 3.113948e+03 -1.737669e+05
#> [469,] 63.1972414 3.993891e+03 2.524029e+05
#> [470,] -21.8027586 4.753603e+02 -1.036417e+04
#> [471,] -65.8027586 4.330003e+03 -2.849261e+05
#> [472,] -42.8027586 1.832076e+03 -7.841791e+04
#> [473,] 0.1972414 3.890416e-02 7.673511e-03
#> [474,] -63.8027586 4.070792e+03 -2.597278e+05
#> [475,] -4.8027586 2.306649e+01 -1.107828e+02
#> [476,] -12.8027586 1.639106e+02 -2.098508e+03
#> [477,] -23.8027586 5.665713e+02 -1.348596e+04
#> [478,] 12.1972414 1.487727e+02 1.814616e+03
#> [479,] -1.8027586 3.249939e+00 -5.858855e+00
#> [480,] 78.1972414 6.114809e+03 4.781612e+05
#> [481,] -5.8027586 3.367201e+01 -1.953905e+02
#> [482,] -11.8027586 1.393051e+02 -1.644185e+03
#> [483,] 25.1972414 6.349010e+02 1.599775e+04
#> [484,] -43.8027586 1.918682e+03 -8.404355e+04
#> [485,] -31.8027586 1.011415e+03 -3.216580e+04
#> [486,] 7.1972414 5.180028e+01 3.728191e+02
#> [487,] -88.8027586 7.885930e+03 -7.002923e+05
#> [488,] 30.1972414 9.118734e+02 2.753606e+04
#> [489,] 45.1972414 2.042791e+03 9.232850e+04
#> [490,] -23.8027586 5.665713e+02 -1.348596e+04
#> [491,] -57.8027586 3.341159e+03 -1.931282e+05
#> [492,] -9.8027586 9.609408e+01 -9.419870e+02
#> [493,] 39.1972414 1.536424e+03 6.022357e+04
#> [494,] 22.1972414 4.927175e+02 1.093697e+04
#> [495,] 15.1972414 2.309561e+02 3.509896e+03
#> [496,] -17.8027586 3.169382e+02 -5.642375e+03
#> [497,] 48.1972414 2.322974e+03 1.119609e+05
#> [498,] -5.8027586 3.367201e+01 -1.953905e+02
#> [499,] -38.8027586 1.505654e+03 -5.842353e+04
#> [500,] 3.1972414 1.022235e+01 3.268333e+01
#> [501,] 23.1972414 5.381120e+02 1.248271e+04
#> [502,] 24.1972414 5.855065e+02 1.416764e+04
#> [503,] -13.8027586 1.905161e+02 -2.629648e+03
#> [504,] 32.1972414 1.036662e+03 3.337767e+04
#> [505,] -2.8027586 7.855456e+00 -2.201695e+01
#> [506,] NA NA NA
#> [507,] 90.1972414 8.135542e+03 7.338035e+05
#> [508,] 14.1972414 2.015617e+02 2.861620e+03
#> [509,] 35.1972414 1.238846e+03 4.360395e+04
#> [510,] 25.1972414 6.349010e+02 1.599775e+04
#> [511,] 6.1972414 3.840580e+01 2.380100e+02
#> [512,] 6.1972414 3.840580e+01 2.380100e+02
#> [513,] 0.1972414 3.890416e-02 7.673511e-03
#> [514,] 26.1972414 6.862955e+02 1.797905e+04
#> [515,] -35.8027586 1.281838e+03 -4.589332e+04
#> [516,] -30.8027586 9.488099e+02 -2.922596e+04
#> [517,] 10.1972414 1.039837e+02 1.060347e+03
#> [518,] 17.1972414 2.957451e+02 5.086000e+03
#> [519,] -12.8027586 1.639106e+02 -2.098508e+03
#> [520,] -21.8027586 4.753603e+02 -1.036417e+04
#> [521,] -48.8027586 2.381709e+03 -1.162340e+05
#> [522,] -26.8027586 7.183879e+02 -1.925478e+04
#> [523,] -20.8027586 4.327548e+02 -9.002493e+03
#> [524,] 1.1972414 1.433387e+00 1.716110e+00
#> [525,] 15.1972414 2.309561e+02 3.509896e+03
#> [526,] -57.8027586 3.341159e+03 -1.931282e+05
#> [527,] -4.8027586 2.306649e+01 -1.107828e+02
#> [528,] 36.1972414 1.310240e+03 4.742708e+04
#> [529,] -17.8027586 3.169382e+02 -5.642375e+03
#> [530,] -12.8027586 1.639106e+02 -2.098508e+03
#> [531,] -16.8027586 2.823327e+02 -4.743968e+03
#> [532,] -27.8027586 7.729934e+02 -2.149135e+04
#> [533,] 60.1972414 3.623708e+03 2.181372e+05
#> [534,] -4.8027586 2.306649e+01 -1.107828e+02
#> [535,] 41.1972414 1.697213e+03 6.992048e+04
#> [536,] -64.8027586 4.199398e+03 -2.721325e+05
#> [537,] 7.1972414 5.180028e+01 3.728191e+02
#> [538,] 22.1972414 4.927175e+02 1.093697e+04
#> [539,] 18.1972414 3.311396e+02 6.025827e+03
#> [540,] 21.1972414 4.493230e+02 9.524409e+03
#> [541,] -15.8027586 2.497272e+02 -3.946378e+03
#> [542,] -21.8027586 4.753603e+02 -1.036417e+04
#> [543,] 31.1972414 9.732679e+02 3.036327e+04
#> [544,] 4.1972414 1.761684e+01 7.394211e+01
#> [545,] NA NA NA
#> [546,] 7.1972414 5.180028e+01 3.728191e+02
#> [547,] -22.8027586 5.199658e+02 -1.185665e+04
#> [548,] 16.1972414 2.623506e+02 4.249356e+03
#> [549,] 4.1972414 1.761684e+01 7.394211e+01
#> [550,] 26.1972414 6.862955e+02 1.797905e+04
#> [551,] -30.8027586 9.488099e+02 -2.922596e+04
#> [552,] 5.1972414 2.701132e+01 1.403843e+02
#> [553,] 23.1972414 5.381120e+02 1.248271e+04
#> [554,] -38.8027586 1.505654e+03 -5.842353e+04
#> [555,] 29.1972414 8.524789e+02 2.489003e+04
#> [556,] -33.8027586 1.142626e+03 -3.862393e+04
#> [557,] 28.1972414 7.950844e+02 2.241919e+04
#> [558,] NA NA NA
#> [559,] -0.8027586 6.444214e-01 -5.173148e-01
#> [560,] 12.1972414 1.487727e+02 1.814616e+03
#> [561,] -4.8027586 2.306649e+01 -1.107828e+02
#> [562,] 17.1972414 2.957451e+02 5.086000e+03
#> [563,] -13.8027586 1.905161e+02 -2.629648e+03
#> [564,] -7.8027586 6.088304e+01 -4.750557e+02
#> [565,] -37.8027586 1.429049e+03 -5.402198e+04
#> [566,] -2.8027586 7.855456e+00 -2.201695e+01
#> [567,] 29.1972414 8.524789e+02 2.489003e+04
#> [568,] -79.8027586 6.368480e+03 -5.082223e+05
#> [569,] 14.1972414 2.015617e+02 2.861620e+03
#> [570,] -13.8027586 1.905161e+02 -2.629648e+03
#> [571,] 6.1972414 3.840580e+01 2.380100e+02
#> [572,] -1.8027586 3.249939e+00 -5.858855e+00
#> [573,] 25.1972414 6.349010e+02 1.599775e+04
#> [574,] -48.8027586 2.381709e+03 -1.162340e+05
#> [575,] 33.1972414 1.102057e+03 3.658525e+04
#> [576,] -14.8027586 2.191217e+02 -3.243605e+03
#> [577,] 5.1972414 2.701132e+01 1.403843e+02
#> [578,] -14.8027586 2.191217e+02 -3.243605e+03
#> [579,] -2.8027586 7.855456e+00 -2.201695e+01
#> [580,] 18.1972414 3.311396e+02 6.025827e+03
#> [581,] -50.8027586 2.580920e+03 -1.311179e+05
#> [582,] -39.8027586 1.584260e+03 -6.305790e+04
#> [583,] -16.8027586 2.823327e+02 -4.743968e+03
#> [584,] 90.1972414 8.135542e+03 7.338035e+05
#> [585,] -4.8027586 2.306649e+01 -1.107828e+02
#> [586,] -42.8027586 1.832076e+03 -7.841791e+04
#> [587,] 16.1972414 2.623506e+02 4.249356e+03
#> [588,] 2.1972414 4.827870e+00 1.060800e+01
#> [589,] -3.8027586 1.446097e+01 -5.499159e+01
#> [590,] 30.1972414 9.118734e+02 2.753606e+04
#> [591,] -4.8027586 2.306649e+01 -1.107828e+02
#> [592,] 42.1972414 1.780607e+03 7.513671e+04
#> [593,] -61.8027586 3.819581e+03 -2.360606e+05
#> [594,] -16.8027586 2.823327e+02 -4.743968e+03
#> [595,] -15.8027586 2.497272e+02 -3.946378e+03
#> [596,] -47.8027586 2.285104e+03 -1.092343e+05
#> [597,] 39.1972414 1.536424e+03 6.022357e+04
#> [598,] 33.1972414 1.102057e+03 3.658525e+04
#> [599,] -3.8027586 1.446097e+01 -5.499159e+01
#> [600,] 27.1972414 7.396899e+02 2.011753e+04
#> [601,] -36.8027586 1.354443e+03 -4.984724e+04
#> [602,] 38.1972414 1.459029e+03 5.573089e+04
#> [603,] -25.8027586 6.657824e+02 -1.717902e+04
#> [604,] -31.8027586 1.011415e+03 -3.216580e+04
#> [605,] -47.8027586 2.285104e+03 -1.092343e+05
#> [606,] -4.8027586 2.306649e+01 -1.107828e+02
#> [607,] 41.1972414 1.697213e+03 6.992048e+04
#> [608,] -20.8027586 4.327548e+02 -9.002493e+03
#> [609,] 9.1972414 8.458925e+01 7.779877e+02
#> [610,] -3.8027586 1.446097e+01 -5.499159e+01
#> [611,] -42.8027586 1.832076e+03 -7.841791e+04
#> [612,] -10.8027586 1.166996e+02 -1.260678e+03
#> [613,] -23.8027586 5.665713e+02 -1.348596e+04
#> [614,] 43.1972414 1.866002e+03 8.060612e+04
#> [615,] -19.8027586 3.921492e+02 -7.765637e+03
#> [616,] -33.8027586 1.142626e+03 -3.862393e+04
#> [617,] 35.1972414 1.238846e+03 4.360395e+04
#> [618,] 21.1972414 4.493230e+02 9.524409e+03
#> [619,] 30.1972414 9.118734e+02 2.753606e+04
#> [620,] 22.1972414 4.927175e+02 1.093697e+04
#> [621,] -10.8027586 1.166996e+02 -1.260678e+03
#> [622,] -25.8027586 6.657824e+02 -1.717902e+04
#> [623,] 56.1972414 3.158130e+03 1.774782e+05
#> [624,] -31.8027586 1.011415e+03 -3.216580e+04
#> [625,] 20.1972414 4.079286e+02 8.239032e+03
#> [626,] 71.1972414 5.069047e+03 3.609022e+05
#> [627,] -4.8027586 2.306649e+01 -1.107828e+02
#> [628,] 5.1972414 2.701132e+01 1.403843e+02
#> [629,] 7.1972414 5.180028e+01 3.728191e+02
#> [630,] -1.8027586 3.249939e+00 -5.858855e+00
#> [631,] -33.8027586 1.142626e+03 -3.862393e+04
#> [632,] 30.1972414 9.118734e+02 2.753606e+04
#> [633,] -26.8027586 7.183879e+02 -1.925478e+04
#> [634,] -25.8027586 6.657824e+02 -1.717902e+04
#> [635,] -39.8027586 1.584260e+03 -6.305790e+04
#> [636,] 4.1972414 1.761684e+01 7.394211e+01
#> [637,] 46.1972414 2.134185e+03 9.859346e+04
#> [638,] -9.8027586 9.609408e+01 -9.419870e+02
#> [639,] 24.1972414 5.855065e+02 1.416764e+04
#> [640,] 8.1972414 6.719477e+01 5.508117e+02
#> [641,] -25.8027586 6.657824e+02 -1.717902e+04
#> [642,] 92.1972414 8.500331e+03 7.837071e+05
#> [643,] 30.1972414 9.118734e+02 2.753606e+04
#> [644,] -1.8027586 3.249939e+00 -5.858855e+00
#> [645,] -13.8027586 1.905161e+02 -2.629648e+03
#> [646,] 35.1972414 1.238846e+03 4.360395e+04
#> [647,] -17.8027586 3.169382e+02 -5.642375e+03
#> [648,] -22.8027586 5.199658e+02 -1.185665e+04
#> [649,] -50.8027586 2.580920e+03 -1.311179e+05
#> [650,] -9.8027586 9.609408e+01 -9.419870e+02
#> [651,] 12.1972414 1.487727e+02 1.814616e+03
#> [652,] 43.1972414 1.866002e+03 8.060612e+04
#> [653,] 26.1972414 6.862955e+02 1.797905e+04
#> [654,] -1.8027586 3.249939e+00 -5.858855e+00
#> [655,] -26.8027586 7.183879e+02 -1.925478e+04
#> [656,] -42.8027586 1.832076e+03 -7.841791e+04
#> [657,] -32.8027586 1.076021e+03 -3.529646e+04
#> [658,] 42.1972414 1.780607e+03 7.513671e+04
#> [659,] 15.1972414 2.309561e+02 3.509896e+03
#> [660,] -11.8027586 1.393051e+02 -1.644185e+03
#> [661,] -34.8027586 1.211232e+03 -4.215422e+04
#> [662,] -2.8027586 7.855456e+00 -2.201695e+01
#> [663,] -48.8027586 2.381709e+03 -1.162340e+05
#> [664,] 26.1972414 6.862955e+02 1.797905e+04
#> [665,] 6.1972414 3.840580e+01 2.380100e+02
#> [666,] -20.8027586 4.327548e+02 -9.002493e+03
#> [667,] -29.8027586 8.882044e+02 -2.647094e+04
#> [668,] 7.1972414 5.180028e+01 3.728191e+02
#> [669,] 45.1972414 2.042791e+03 9.232850e+04
#> [670,] 13.1972414 1.741672e+02 2.298526e+03
#> [671,] -39.8027586 1.584260e+03 -6.305790e+04
#> [672,] 2.1972414 4.827870e+00 1.060800e+01
#> [673,] -23.8027586 5.665713e+02 -1.348596e+04
#> [674,] 10.1972414 1.039837e+02 1.060347e+03
#> [675,] 62.1972414 3.868497e+03 2.406098e+05
#> [676,] 2.1972414 4.827870e+00 1.060800e+01
#> [677,] 57.1972414 3.271524e+03 1.871222e+05
#> [678,] -36.8027586 1.354443e+03 -4.984724e+04
#> [679,] -46.8027586 2.190498e+03 -1.025214e+05
#> [680,] -19.8027586 3.921492e+02 -7.765637e+03
#> [681,] 45.1972414 2.042791e+03 9.232850e+04
#> [682,] -13.8027586 1.905161e+02 -2.629648e+03
#> [683,] 10.1972414 1.039837e+02 1.060347e+03
#> [684,] -29.8027586 8.882044e+02 -2.647094e+04
#> [685,] -13.8027586 1.905161e+02 -2.629648e+03
#> [686,] -48.8027586 2.381709e+03 -1.162340e+05
#> [687,] -16.8027586 2.823327e+02 -4.743968e+03
#> [688,] 5.1972414 2.701132e+01 1.403843e+02
#> [689,] -34.8027586 1.211232e+03 -4.215422e+04
#> [690,] -9.8027586 9.609408e+01 -9.419870e+02
#> [691,] -41.8027586 1.747471e+03 -7.304909e+04
#> [692,] 6.1972414 3.840580e+01 2.380100e+02
#> [693,] 52.1972414 2.724552e+03 1.422141e+05
#> [694,] 42.1972414 1.780607e+03 7.513671e+04
#> [695,] -56.8027586 3.226553e+03 -1.832771e+05
#> [696,] 20.1972414 4.079286e+02 8.239032e+03
#> [697,] -33.8027586 1.142626e+03 -3.862393e+04
#> [698,] 12.1972414 1.487727e+02 1.814616e+03
#> [699,] -5.8027586 3.367201e+01 -1.953905e+02
#> [700,] -15.8027586 2.497272e+02 -3.946378e+03
#> [701,] -8.8027586 7.748856e+01 -6.821131e+02
#> [702,] -36.8027586 1.354443e+03 -4.984724e+04
#> [703,] 13.1972414 1.741672e+02 2.298526e+03
#> [704,] -22.8027586 5.199658e+02 -1.185665e+04
#> [705,] 25.1972414 6.349010e+02 1.599775e+04
#> [706,] -40.8027586 1.664865e+03 -6.793109e+04
#> [707,] -2.8027586 7.855456e+00 -2.201695e+01
#> [708,] 34.1972414 1.169451e+03 3.999201e+04
#> [709,] -13.8027586 1.905161e+02 -2.629648e+03
#> [710,] -4.8027586 2.306649e+01 -1.107828e+02
#> [711,] 35.1972414 1.238846e+03 4.360395e+04
#> [712,] -32.8027586 1.076021e+03 -3.529646e+04
#> [713,] 27.1972414 7.396899e+02 2.011753e+04
#> [714,] 30.1972414 9.118734e+02 2.753606e+04
#> [715,] -5.8027586 3.367201e+01 -1.953905e+02
#> [716,] -54.8027586 3.003342e+03 -1.645914e+05
#> [717,] 2.1972414 4.827870e+00 1.060800e+01
#> [718,] 19.1972414 3.685341e+02 7.074838e+03
#> [719,] 15.1972414 2.309561e+02 3.509896e+03
#> [720,] -25.8027586 6.657824e+02 -1.717902e+04
#> [721,] 55.1972414 3.046735e+03 1.681714e+05
#> [722,] -68.8027586 4.733820e+03 -3.256998e+05
#> [723,] 56.1972414 3.158130e+03 1.774782e+05
#> [724,] 7.1972414 5.180028e+01 3.728191e+02
#> [725,] -32.8027586 1.076021e+03 -3.529646e+04
#> [726,] 12.1972414 1.487727e+02 1.814616e+03
#> [727,] -15.8027586 2.497272e+02 -3.946378e+03
#> [728,] -39.8027586 1.584260e+03 -6.305790e+04
#> [729,] 22.1972414 4.927175e+02 1.093697e+04
#> [730,] 1.1972414 1.433387e+00 1.716110e+00
#> [731,] -1.8027586 3.249939e+00 -5.858855e+00
#> [732,] 18.1972414 3.311396e+02 6.025827e+03
#> [733,] -3.8027586 1.446097e+01 -5.499159e+01
#> [734,] 29.1972414 8.524789e+02 2.489003e+04
#> [735,] 21.1972414 4.493230e+02 9.524409e+03
#> attr(,"transformation")
#> [1] "polynomial"
#> attr(,"name")
#> [1] "polynomial(mri$ldl)"
#> attr(,"prnm")
#> [1] "mri$ldl"
#> attr(,"degree")
#> [1] 1 2 3
#> attr(,"center")
#> [1] 125.8028
#> attr(,"original")
#> [1] 9.1972414 -41.8027586 -10.8027586 -64.8027586 22.1972414
#> [6] 37.1972414 -24.8027586 -9.8027586 -1.8027586 -15.8027586
#> [11] 10.1972414 -17.8027586 80.1972414 4.1972414 -4.8027586
#> [16] -30.8027586 -47.8027586 -15.8027586 10.1972414 34.1972414
#> [21] 4.1972414 -42.8027586 28.1972414 -55.8027586 -50.8027586
#> [26] -11.8027586 13.1972414 -44.8027586 -4.8027586 -20.8027586
#> [31] -3.8027586 9.1972414 10.1972414 62.1972414 19.1972414
#> [36] -26.8027586 45.1972414 -42.8027586 29.1972414 17.1972414
#> [41] 22.1972414 36.1972414 56.1972414 -43.8027586 25.1972414
#> [46] -13.8027586 33.1972414 -50.8027586 -27.8027586 63.1972414
#> [51] -2.8027586 101.1972414 -52.8027586 -26.8027586 90.1972414
#> [56] -21.8027586 -45.8027586 -36.8027586 -17.8027586 21.1972414
#> [61] 17.1972414 25.1972414 25.1972414 12.1972414 -12.8027586
#> [66] 20.1972414 23.1972414 -53.8027586 6.1972414 -31.8027586
#> [71] -40.8027586 -46.8027586 9.1972414 13.1972414 -40.8027586
#> [76] -22.8027586 1.1972414 14.1972414 -24.8027586 -4.8027586
#> [81] 20.1972414 -38.8027586 24.1972414 8.1972414 -45.8027586
#> [86] 3.1972414 28.1972414 -26.8027586 -8.8027586 47.1972414
#> [91] -5.8027586 -0.8027586 -15.8027586 -40.8027586 -55.8027586
#> [96] 41.1972414 -3.8027586 37.1972414 -25.8027586 32.1972414
#> [101] -10.8027586 -28.8027586 -41.8027586 -44.8027586 25.1972414
#> [106] 65.1972414 41.1972414 -8.8027586 -21.8027586 -12.8027586
#> [111] -7.8027586 -5.8027586 -27.8027586 0.1972414 16.1972414
#> [116] 20.1972414 -42.8027586 7.1972414 2.1972414 4.1972414
#> [121] -35.8027586 10.1972414 -18.8027586 -2.8027586 -21.8027586
#> [126] -23.8027586 -25.8027586 51.1972414 -39.8027586 NA
#> [131] 38.1972414 30.1972414 -17.8027586 -72.8027586 24.1972414
#> [136] 69.1972414 8.1972414 16.1972414 -37.8027586 4.1972414
#> [141] 48.1972414 25.1972414 22.1972414 -6.8027586 35.1972414
#> [146] 10.1972414 30.1972414 -26.8027586 -11.8027586 12.1972414
#> [151] -21.8027586 -12.8027586 7.1972414 -44.8027586 -47.8027586
#> [156] -10.8027586 15.1972414 NA -47.8027586 -24.8027586
#> [161] 29.1972414 10.1972414 5.1972414 12.1972414 14.1972414
#> [166] -10.8027586 18.1972414 0.1972414 91.1972414 19.1972414
#> [171] 1.1972414 -49.8027586 NA -74.8027586 6.1972414
#> [176] -11.8027586 -64.8027586 23.1972414 11.1972414 1.1972414
#> [181] 7.1972414 7.1972414 18.1972414 -16.8027586 26.1972414
#> [186] -30.8027586 3.1972414 NA 2.1972414 -26.8027586
#> [191] 13.1972414 16.1972414 -22.8027586 -11.8027586 59.1972414
#> [196] -1.8027586 -27.8027586 -24.8027586 65.1972414 -7.8027586
#> [201] -28.8027586 -8.8027586 -30.8027586 9.1972414 80.1972414
#> [206] 24.1972414 -3.8027586 14.1972414 -34.8027586 -4.8027586
#> [211] 34.1972414 36.1972414 62.1972414 -30.8027586 38.1972414
#> [216] 4.1972414 28.1972414 NA 10.1972414 -30.8027586
#> [221] 121.1972414 18.1972414 -42.8027586 -3.8027586 31.1972414
#> [226] -18.8027586 5.1972414 14.1972414 74.1972414 -30.8027586
#> [231] 20.1972414 -23.8027586 11.1972414 24.1972414 22.1972414
#> [236] -53.8027586 16.1972414 -34.8027586 27.1972414 82.1972414
#> [241] -7.8027586 24.1972414 25.1972414 1.1972414 -25.8027586
#> [246] -11.8027586 -13.8027586 3.1972414 53.1972414 8.1972414
#> [251] -10.8027586 -17.8027586 -26.8027586 81.1972414 -48.8027586
#> [256] 15.1972414 4.1972414 36.1972414 1.1972414 -33.8027586
#> [261] NA -25.8027586 61.1972414 -12.8027586 2.1972414
#> [266] 57.1972414 NA -15.8027586 -22.8027586 30.1972414
#> [271] -26.8027586 29.1972414 55.1972414 -7.8027586 6.1972414
#> [276] 22.1972414 -35.8027586 5.1972414 12.1972414 -16.8027586
#> [281] 10.1972414 13.1972414 -34.8027586 50.1972414 22.1972414
#> [286] -38.8027586 -12.8027586 29.1972414 59.1972414 11.1972414
#> [291] 11.1972414 43.1972414 80.1972414 -37.8027586 13.1972414
#> [296] 56.1972414 11.1972414 59.1972414 -9.8027586 -11.8027586
#> [301] -56.8027586 -6.8027586 -43.8027586 4.1972414 -15.8027586
#> [306] 36.1972414 -12.8027586 -32.8027586 -7.8027586 70.1972414
#> [311] -47.8027586 -14.8027586 -54.8027586 15.1972414 -13.8027586
#> [316] 69.1972414 -9.8027586 18.1972414 47.1972414 27.1972414
#> [321] 10.1972414 56.1972414 -35.8027586 2.1972414 62.1972414
#> [326] 58.1972414 -4.8027586 -0.8027586 -12.8027586 8.1972414
#> [331] 14.1972414 -45.8027586 45.1972414 -44.8027586 -13.8027586
#> [336] -9.8027586 32.1972414 -19.8027586 32.1972414 -34.8027586
#> [341] 6.1972414 -37.8027586 -12.8027586 21.1972414 16.1972414
#> [346] -58.8027586 35.1972414 52.1972414 -17.8027586 -38.8027586
#> [351] 11.1972414 34.1972414 1.1972414 -29.8027586 -15.8027586
#> [356] 49.1972414 45.1972414 5.1972414 -32.8027586 -17.8027586
#> [361] -7.8027586 -9.8027586 -68.8027586 1.1972414 55.1972414
#> [366] -17.8027586 7.1972414 16.1972414 8.1972414 0.1972414
#> [371] 38.1972414 -48.8027586 -28.8027586 -70.8027586 -32.8027586
#> [376] -57.8027586 10.1972414 -1.8027586 7.1972414 -4.8027586
#> [381] -45.8027586 17.1972414 16.1972414 -48.8027586 -53.8027586
#> [386] -36.8027586 -48.8027586 -32.8027586 75.1972414 55.1972414
#> [391] -55.8027586 -52.8027586 -20.8027586 8.1972414 -18.8027586
#> [396] 94.1972414 40.1972414 7.1972414 99.1972414 6.1972414
#> [401] 17.1972414 -41.8027586 33.1972414 -11.8027586 -50.8027586
#> [406] -14.8027586 -11.8027586 44.1972414 -16.8027586 32.1972414
#> [411] 47.1972414 -31.8027586 27.1972414 0.1972414 6.1972414
#> [416] -10.8027586 -35.8027586 2.1972414 19.1972414 18.1972414
#> [421] 13.1972414 21.1972414 -44.8027586 69.1972414 35.1972414
#> [426] -28.8027586 -19.8027586 -73.8027586 20.1972414 2.1972414
#> [431] -3.8027586 7.1972414 -14.8027586 30.1972414 -16.8027586
#> [436] 13.1972414 -11.8027586 -114.8027586 0.1972414 -4.8027586
#> [441] 7.1972414 -86.8027586 -19.8027586 11.1972414 30.1972414
#> [446] 36.1972414 -15.8027586 -3.8027586 21.1972414 43.1972414
#> [451] 40.1972414 -31.8027586 -33.8027586 -9.8027586 -29.8027586
#> [456] -36.8027586 13.1972414 61.1972414 47.1972414 -35.8027586
#> [461] -23.8027586 18.1972414 -35.8027586 42.1972414 -14.8027586
#> [466] -24.8027586 -23.8027586 -55.8027586 63.1972414 -21.8027586
#> [471] -65.8027586 -42.8027586 0.1972414 -63.8027586 -4.8027586
#> [476] -12.8027586 -23.8027586 12.1972414 -1.8027586 78.1972414
#> [481] -5.8027586 -11.8027586 25.1972414 -43.8027586 -31.8027586
#> [486] 7.1972414 -88.8027586 30.1972414 45.1972414 -23.8027586
#> [491] -57.8027586 -9.8027586 39.1972414 22.1972414 15.1972414
#> [496] -17.8027586 48.1972414 -5.8027586 -38.8027586 3.1972414
#> [501] 23.1972414 24.1972414 -13.8027586 32.1972414 -2.8027586
#> [506] NA 90.1972414 14.1972414 35.1972414 25.1972414
#> [511] 6.1972414 6.1972414 0.1972414 26.1972414 -35.8027586
#> [516] -30.8027586 10.1972414 17.1972414 -12.8027586 -21.8027586
#> [521] -48.8027586 -26.8027586 -20.8027586 1.1972414 15.1972414
#> [526] -57.8027586 -4.8027586 36.1972414 -17.8027586 -12.8027586
#> [531] -16.8027586 -27.8027586 60.1972414 -4.8027586 41.1972414
#> [536] -64.8027586 7.1972414 22.1972414 18.1972414 21.1972414
#> [541] -15.8027586 -21.8027586 31.1972414 4.1972414 NA
#> [546] 7.1972414 -22.8027586 16.1972414 4.1972414 26.1972414
#> [551] -30.8027586 5.1972414 23.1972414 -38.8027586 29.1972414
#> [556] -33.8027586 28.1972414 NA -0.8027586 12.1972414
#> [561] -4.8027586 17.1972414 -13.8027586 -7.8027586 -37.8027586
#> [566] -2.8027586 29.1972414 -79.8027586 14.1972414 -13.8027586
#> [571] 6.1972414 -1.8027586 25.1972414 -48.8027586 33.1972414
#> [576] -14.8027586 5.1972414 -14.8027586 -2.8027586 18.1972414
#> [581] -50.8027586 -39.8027586 -16.8027586 90.1972414 -4.8027586
#> [586] -42.8027586 16.1972414 2.1972414 -3.8027586 30.1972414
#> [591] -4.8027586 42.1972414 -61.8027586 -16.8027586 -15.8027586
#> [596] -47.8027586 39.1972414 33.1972414 -3.8027586 27.1972414
#> [601] -36.8027586 38.1972414 -25.8027586 -31.8027586 -47.8027586
#> [606] -4.8027586 41.1972414 -20.8027586 9.1972414 -3.8027586
#> [611] -42.8027586 -10.8027586 -23.8027586 43.1972414 -19.8027586
#> [616] -33.8027586 35.1972414 21.1972414 30.1972414 22.1972414
#> [621] -10.8027586 -25.8027586 56.1972414 -31.8027586 20.1972414
#> [626] 71.1972414 -4.8027586 5.1972414 7.1972414 -1.8027586
#> [631] -33.8027586 30.1972414 -26.8027586 -25.8027586 -39.8027586
#> [636] 4.1972414 46.1972414 -9.8027586 24.1972414 8.1972414
#> [641] -25.8027586 92.1972414 30.1972414 -1.8027586 -13.8027586
#> [646] 35.1972414 -17.8027586 -22.8027586 -50.8027586 -9.8027586
#> [651] 12.1972414 43.1972414 26.1972414 -1.8027586 -26.8027586
#> [656] -42.8027586 -32.8027586 42.1972414 15.1972414 -11.8027586
#> [661] -34.8027586 -2.8027586 -48.8027586 26.1972414 6.1972414
#> [666] -20.8027586 -29.8027586 7.1972414 45.1972414 13.1972414
#> [671] -39.8027586 2.1972414 -23.8027586 10.1972414 62.1972414
#> [676] 2.1972414 57.1972414 -36.8027586 -46.8027586 -19.8027586
#> [681] 45.1972414 -13.8027586 10.1972414 -29.8027586 -13.8027586
#> [686] -48.8027586 -16.8027586 5.1972414 -34.8027586 -9.8027586
#> [691] -41.8027586 6.1972414 52.1972414 42.1972414 -56.8027586
#> [696] 20.1972414 -33.8027586 12.1972414 -5.8027586 -15.8027586
#> [701] -8.8027586 -36.8027586 13.1972414 -22.8027586 25.1972414
#> [706] -40.8027586 -2.8027586 34.1972414 -13.8027586 -4.8027586
#> [711] 35.1972414 -32.8027586 27.1972414 30.1972414 -5.8027586
#> [716] -54.8027586 2.1972414 19.1972414 15.1972414 -25.8027586
#> [721] 55.1972414 -68.8027586 56.1972414 7.1972414 -32.8027586
#> [726] 12.1972414 -15.8027586 -39.8027586 22.1972414 1.1972414
#> [731] -1.8027586 18.1972414 -3.8027586 29.1972414 21.1972414
# Use a polynomial in regress
regress("mean", atrophy ~ polynomial(age, degree = 2), data = mri)
#>
#> Call:
#> regress(fnctl = "mean", formula = atrophy ~ polynomial(age, degree = 2),
#> data = mri)
#>
#> Residuals:
#> Min 1Q Median 3Q Max
#> -36.828 -9.063 -1.018 7.816 50.887
#>
#> Coefficients:
#> Estimate Naive SE Robust SE 95%L
#> [1] Intercept 35.34 0.5598 0.5622 34.24
#> polynomial(age, degree = 2)
#> [2] age^1 0.5869 0.1009 0.09808 0.3943
#> [3] age^2 0.02159 0.01099 0.01150 -9.897e-04
#> 95%H F stat df Pr(>F)
#> [1] Intercept 36.45 3952.13 1 < 0.00005
#> polynomial(age, degree = 2) 32.91 2 < 0.00005
#> [2] age^1 0.7794 35.80 1 < 0.00005
#> [3] age^2 0.04417 3.52 1 0.0609
#>
#> Polynomial terms calculated from age, centered at 74.566
#>
#> Residual standard error: 12.33 on 732 degrees of freedom
#> Multiple R-squared: 0.09148, Adjusted R-squared: 0.089
#> F-statistic: 32.91 on 2 and 732 DF, p-value: 2.06e-14
#>