I have some time series data and I would like to run "mini" regression models on it but on a rolling basis.
I can run a linear model and plot it using:
mod <- lm(Y1 ~ time(data), data = data)
ggplot(data, aes(x = time(data), y = Y1)) +
geom_line() +
geom_line(aes(x = as.Date(rownames(data.frame(mod$fitted.values))),
y = mod$fitted.values), colour = "blue")
Now I want to use the rollapply function to apply the linear model. I can use the roll_lm function.
library(roll)
rolledLinearModel <- roll_lm(
y = data[, "Y1"],
x = time(data),
width = 200
)
However, this does not store the fitted values. I am trying to use the rollapply function to do so.
library(zoo)
lmFunction <- function(x){
lm(Y1 ~ time(x), data = x)
}
data %>%
zoo() %>%
rollapply(
width = 40,
FUN = lmFunction
#fill = NA,
#by = 30,
#by.column = FALSE,
#align = "right"
)
However, I cannot seem to get it to work. Once I get the results I would like to plot all the (rolling) fitted values as (many different) lines on the plot.
Data:
data <- structure(c(46.423225, 46.59877, 46.010689, 46.265236, 45.975586,
44.685474, 45.480743, 45.312847, 44.950558, 45.454235, 45.683979,
45.957905, 44.570602, 46.320194, 45.904892, 46.408562, 45.675144,
46.072777, 46.691334, 46.373215, 46.470417, 45.065445, 44.641296,
45.330524, 44.323177, 44.181801, 44.641296, 44.199478, 44.349686,
44.552933, 44.031586, 43.916706, 47.884232, 48.379078, 49.138992,
49.810566, 50.455616, 50.013802, 49.969612, 49.050632, 49.545479,
50.38493, 49.775219, 49.465939, 49.412922, 49.377586, 49.943104,
50.367249, 51.940125, 52.523323, 52.134529, 52.152199, 51.507149,
50.968124, 51.339252, 51.374603, 51.825249, 52.797256, 52.97398,
53.159546, 52.850277, 52.947468, 53.804592, 54.84729, 52.761902,
53.000492, 53.327438, 54.529179, 54.405468, 53.733906, 53.680893,
53.954819, 53.323814, 53.172737, 52.58617, 52.639492, 53.137177,
53.066078, 51.404156, 51.706329, 51.937397, 51.466366, 51.688549,
51.848522, 52.381763, 52.959442, 53.01276, 52.523964, 51.679668,
51.990719, 51.555241, 51.830746, 52.195129, 52.390652, 52.008492,
51.750771, 52.097378, 52.666157, 53.057194, 53.341595, 53.483791,
52.515076, 51.475258, 50.693172, 50.293247, 49.991081, 50.168823,
50.355453, 50.311031, 50.51543, 50.444328, 50.959793, 51.125267,
51.644028, 51.259422, 50.463383, 50.919544, 51.098431, 51.474087,
51.590363, 50.588608, 51.00005, 50.427608, 52.234341, 52.449009,
52.610001, 51.500927, 51.679802, 52.350624, 52.547401, 53.710148,
54.121578, 53.754868, 56.295029, 55.874653, 55.677883, 55.320114,
54.175247, 54.139469, 54.139469, 53.844307, 53.826427, 54.282574,
54.202076, 54.56879, 54.112637, 54.989174, 54.550903, 53.683315,
54.41674, 54.479347, 55.311165, 55.498993, 54.604576, 53.397095,
53.128769, 50.901653, 50.454445, 50.508106, 49.917789, 50.373947,
51.169979, 51.769249, 51.903408, 50.812214, 50.436554, 50.186123,
49.694187, 49.747845, 50.257668, 49.676296, 50.284504, 50.758549,
49.676296, 49.381134, 48.334663, 48.397266, 47.395515, 47.189796,
47.064575, 47.359734, 47.589401, 46.832859, 46.706779, 47.121067,
47.742504, 47.742504, 48.012695, 47.94965, 48.120773, 46.895912,
46.949951, 47.598412, 47.98568, 47.940647, 48.354935, 48.039719,
49.435692, 49.976078, 49.940052, 49.994087, 50.399372, 50.642544,
50.021107, 49.624825, 49.093456, 49.273579, 48.706184, 48.652142,
49.38166, 48.571087, 47.589401, 47.283188, 47.27417, 46.82386,
45.346828, 46.310501, 46.490631, 46.832859, 44.085945, 43.923824,
43.923824, 43.464508, 43.662643, 43.356422, 43.806747, 45.26577,
44.004883, 43.833763, 44.004883, 43.050209, 43.221336, 43.536556,
43.284374, 43.059219, 43.050209, 42.617908, 42.816051, 43.140274,
42.942135, 42.825058, 42.419773, 42.221638, 42.527843, 43.203323,
42.194607, 42.509838, 42.978161, 42.383747, 42.392757, 42.500824,
42.338715, 42.662945, 41.366032, 42.014492, 42.284679, 42.671947,
43.175678, 43.130299, 43.529648, 44.391899, 43.575035, 44.07423,
43.69302, 43.883621, 43.529648, 43.783783, 44.428204, 45.009079,
44.745869, 44.60973, 43.756557, 44.065147, 44.428204, 44.273903,
42.522182, 41.941311, 41.551029, 42.041149, 41.995766, 42.168213,
42.313438, 42.268051, 43.130299, 43.584103, 43.647644, 44.927391,
44.028839, 44.900166, 44.373741, 44.292057, 43.420738, 42.839863,
42.966927, 43.38443, 43.856403, 42.694641, 40.979225, 40.643402,
41.305969, 41.614567, 42.930618, 42.549423, 42.594803, 43.084919,
43.656719, 42.885235, 42.440498, 42.367893, 43.275513, 43.021385,
43.130299, 44.700481, 44.455433, 45.562733, 45.408432, 47.02401,
45.962086, 46.179916, 46.452209, 46.461273, 46.098225, 46.30698,
46.297901, 45.726559, 46.128777, 47.216648, 46.796124, 46.357327,
47.362907, 47.390331, 47.326347, 46.878399, 47.143509, 46.512733,
47.591454, 45.799683, 46.595013, 46.595013, 47.116085, 47.783428,
47.673729, 49.026688, 49.739738, 47.518318, 48.002823, 47.847416,
48.450764, 49.191235, 49.428928, 51.147552, 50.004845, 50.416222,
50.544209, 46.476173, 44.300457, 44.501572, 43.422855, 42.417274,
44.44672, 44.318741, 43.459419, 43.724533, 44.19989, 43.527405,
44.163044, 43.895897, 43.4445, 42.965462, 42.652256, 42.993103,
42.937824, 42.956257, 41.952137, 41.804737, 41.556007, 40.717705,
41.66655, 40.837463, 40.441341, 40.846676, 40.855885, 41.28886,
41.316494, 40.036007, 39.999153, 40.662434, 40.109707, 39.612244,
39.870186, 38.810791, 38.17514, 39.059517, 38.709457, 39.317463,
39.501705, 38.737087, 39.26218, 39.326668, 39.252972, 38.866058,
38.64497, 39.041096, 38.737087, 38.820004, 38.875278, 38.727882,
38.681816, 38.598911, 37.73296, 37.806664, 37.77903, 38.700245,
38.995029, 38.525215, 38.525215, 37.889572, 37.52108, 37.051262,
36.867023, 37.769817, 37.539513, 37.134182, 36.94072, 37.428967,
37.557945, 37.659267, 37.889572, 38.985821, 39.216122, 39.151634,
39.317463, 39.271393, 38.746304, 38.543636, 38.976608, 39.022667,
39.022663, 38.176979, 37.805244, 38.269909, 38.539417, 39.905521,
39.543095, 39.636028, 39.747547, 40.286556, 41.169422, 41.076492,
42.693531, 42.349682, 41.550449, 42.451904, 42.042999, 41.085785,
38.632347, 38.400021, 38.706692, 38.93903, 39.106308, 38.976196,
39.292168, 39.217831, 39.394398, 40.32373, 40.602531, 39.914825,
40.18433, 39.89624, 39.905521, 39.636028, 39.143475, 39.143475,
38.67881, 38.520832, 41.476101, 41.903599, 41.448223, 41.243774,
41.020733, 39.859066, 39.645321, 39.924118, 40.946388, 40.277267,
40.546772, 40.574654, 40.658291, 40.193626, 38.455772, 38.660225,
38.102631, 38.204857, 39.273586, 39.440865, 40.221504, 39.812599,
39.719662, 39.821892, 40.444546, 40.072819, 40.175037, 39.450161,
39.608143, 39.859066, 40.955681, 40.955681, 41.234474, 41.448223,
41.262356, 41.336704, 41.40176, 41.141548, 41.582008, 41.619503,
42.369228, 42.584782, 42.500431, 41.863163, 41.97562, 42.013115,
42.172424, 42.453575, 42.781593, 42.697239, 42.125572, 40.907253,
41.197773, 40.9916, 39.820141, 38.573715, 38.358166, 38.348793,
38.151989, 38.629944, 39.051666, 37.823982, 37.083618, 36.296402,
35.49044, 35.36861, 34.96563, 33.728573, 34.225269, 35.387352,
35.818447, 36.821217, 37.299168, 36.746246, 37.02739, 38.386276,
38.292564, 38.170727, 38.217587, 38.986061, 38.911098, 38.779888,
38.779888, 37.664661, 37.561581, 37.5522, 37.383514, 38.798634,
39.023552, 40.110668, 39.820141, 39.63271, 40.663593, 39.229729,
36.971169, 35.846561, 34.646996, 33.841026, 33.756687, 33.074928,
32.876083, 32.137501, 32.03334, 31.484144, 31.256889, 30.518316,
30.215313, 32.014408, 32.099625, 32.156437, 32.307938, 32.270069,
31.379988, 32.8382, 33.813499, 34.495262, 34.542606, 34.656235,
34.608887, 34.883484, 35.016052, 35.101273, 35.300117, 35.442154,
34.400574, 34.135445, 34.561543, 34.921364, 34.665699, 32.913948,
34.542606, 34.902424, 34.703579, 34.627827, 34.665699, 34.466854,
33.709335, 33.37793, 33.794556, 34.06916, 34.561543, 34.731983,
35.319054, 35.139149, 35.290649, 35.442154, 35.139149, 35.593655,
34.902424, 35.356934, 35.148617, 35.868252, 35.546307, 36.739395,
36.152321, 35.470562, 35.659939, 35.356934, 34.85508, 34.898132,
34.859859, 35.012924, 34.553741, 34.764198, 33.463173, 33.635368,
32.372612, 32.525665, 33.319679, 34.008457, 34.084988, 34.113686,
35.357315, 36.141758, 36.313953, 37.366253, 38.562046, 38.542919,
37.146229, 37.624546, 37.71064, 38.475952, 38.332451, 38.992538,
39.260391, 39.116898, 37.997635, 37.232327, 36.141758, 34.9077,
35.500813, 35.912163, 35.615608, 34.917259, 35.46254, 35.596478,
35.051189, 34.725933, 34.649403, 34.391109, 34.343281, 32.525665,
33.233582, 33.444038, 33.577969, 33.06139, 33.300545, 34.113686,
33.654503, 32.774391, 33.003986, 32.439568, 32.183231, 32.318657,
31.651194, 32.415386, 33.121548, 33.247292, 33.208603, 33.51815,
33.914757, 33.943775, 33.421413, 33.759983, 33.508476, 33.431091,
33.837368, 33.769657, 34.833725, 35.443142, 34.146915, 34.040508,
34.514507, 34.427444, 34.659607, 34.514507, 34.021156, 35.075558,
34.843399, 34.582218, 34.292019, 34.475811, 34.92078, 35.433472,
35.201309, 35.404449, 34.688625, 34.572548, 33.711617, 33.160236,
35.288372, 34.524178, 34.408096, 34.727318, 34.417774, 34.050182,
34.611237, 33.701942, 33.86639, 32.212246, 32.308979, 32.270287,
33.227947, 32.734612, 32.067146, 32.686241, 30.809612, 30.316271,
31.206217, 31.49642, 30.838631, 31.022425, 30.92569, 29.716522,
29.88097, 29.397303, 29.803583, 30.383986, 30.857977, 30.016397,
30.451698, 31.515768, 31.486744, 33.276314, 34.233978, 34.533848,
34.872417, 35.462494, 35.320778, 34.128433, 33.845005, 34.001381,
34.470497, 34.607327, 33.248837, 34.245716, 34.421635, 34.900524,
35.154633, 34.705059, 34.020927, 34.275032, 34.56823, 33.336796,
32.418102, 32.506058, 33.033817, 34.304352, 34.411861, 35.604202,
35.731258, 35.868084, 35.741028, 36.239471, 39.132374, 38.663254,
37.871616, 38.497108, 39.435345, 39.562401, 39.220333, 38.027985,
38.438469, 37.969349, 38.555748, 37.666374, 38.360283, 38.497108,
37.646828, 37.490456, 36.952927, 36.816093, 37.119068, 36.698818,
36.024456, 34.714832, 34.959164, 35.643299, 36.288338, 36.327427,
36.131962, 35.55534, 35.47715, 35.095993, 35.154633, 35.066673,
36.073326, 35.88763, 36.346977, 36.718361, 37.38295, 36.864769,
36.430485, 36.943729, 37.318794, 36.874641, 36.943729, 37.427364,
37.397755, 37.348404, 37.269444, 37.062172), index = structure(c(1472688000,
1472774400, 1473120000, 1473206400, 1473292800, 1473379200, 1473638400,
1473724800, 1473811200, 1473897600, 1473984000, 1474243200, 1474329600,
1474416000, 1474502400, 1474588800, 1474848000, 1474934400, 1475020800,
1475107200, 1475193600, 1475452800, 1475539200, 1475625600, 1475712000,
1475798400, 1476057600, 1476144000, 1476230400, 1476316800, 1476403200,
1476662400, 1476748800, 1476835200, 1476921600, 1477008000, 1477267200,
1477353600, 1477440000, 1477526400, 1477612800, 1477872000, 1477958400,
1478044800, 1478131200, 1478217600, 1478476800, 1478563200, 1478649600,
1478736000, 1478822400, 1479081600, 1479168000, 1479254400, 1479340800,
1479427200, 1479686400, 1479772800, 1479859200, 1480032000, 1480291200,
1480377600, 1480464000, 1480550400, 1480636800, 1480896000, 1480982400,
1481068800, 1481155200, 1481241600, 1481500800, 1481587200, 1481673600,
1481760000, 1481846400, 1482105600, 1482192000, 1482278400, 1482364800,
1482451200, 1482796800, 1482883200, 1482969600, 1483056000, 1483401600,
1483488000, 1483574400, 1483660800, 1483920000, 1484006400, 1484092800,
1484179200, 1484265600, 1484611200, 1484697600, 1484784000, 1484870400,
1485129600, 1485216000, 1485302400, 1485388800, 1485475200, 1485734400,
1485820800, 1485907200, 1485993600, 1486080000, 1486339200, 1486425600,
1486512000, 1486598400, 1486684800, 1486944000, 1487030400, 1487116800,
1487203200, 1487289600, 1487635200, 1487721600, 1487808000, 1487894400,
1488153600, 1488240000, 1488326400, 1488412800, 1488499200, 1488758400,
1488844800, 1488931200, 1489017600, 1489104000, 1489363200, 1489449600,
1489536000, 1489622400, 1489708800, 1489968000, 1490054400, 1490140800,
1490227200, 1490313600, 1490572800, 1490659200, 1490745600, 1490832000,
1490918400, 1491177600, 1491264000, 1491350400, 1491436800, 1491523200,
1491782400, 1491868800, 1491955200, 1492041600, 1492387200, 1492473600,
1492560000, 1492646400, 1492732800, 1492992000, 1493078400, 1493164800,
1493251200, 1493337600, 1493596800, 1493683200, 1493769600, 1493856000,
1493942400, 1494201600, 1494288000, 1494374400, 1494460800, 1494547200,
1494806400, 1494892800, 1494979200, 1495065600, 1495152000, 1495411200,
1495497600, 1495584000, 1495670400, 1495756800, 1496102400, 1496188800,
1496275200, 1496361600, 1496620800, 1496707200, 1496793600, 1496880000,
1496966400, 1497225600, 1497312000, 1497398400, 1497484800, 1497571200,
1497830400, 1497916800, 1498003200, 1498089600, 1498176000, 1498435200,
1498521600, 1498608000, 1498694400, 1498780800, 1499040000, 1499212800,
1499299200, 1499385600, 1499644800, 1499731200, 1499817600, 1499904000,
1499990400, 1500249600, 1500336000, 1500422400, 1500508800, 1500595200,
1500854400, 1500940800, 1501027200, 1501113600, 1501200000, 1501459200,
1501545600, 1501632000, 1501718400, 1501804800, 1502064000, 1502150400,
1502236800, 1502323200, 1502409600, 1502668800, 1502755200, 1502841600,
1502928000, 1503014400, 1503273600, 1503360000, 1503446400, 1503532800,
1503619200, 1503878400, 1503964800, 1504051200, 1504137600, 1504224000,
1504569600, 1504656000, 1504742400, 1504828800, 1505088000, 1505174400,
1505260800, 1505347200, 1505433600, 1505692800, 1505779200, 1505865600,
1505952000, 1506038400, 1506297600, 1506384000, 1506470400, 1506556800,
1506643200, 1506902400, 1506988800, 1507075200, 1507161600, 1507248000,
1507507200, 1507593600, 1507680000, 1507766400, 1507852800, 1508112000,
1508198400, 1508284800, 1508371200, 1508457600, 1508716800, 1508803200,
1508889600, 1508976000, 1509062400, 1509321600, 1509408000, 1509494400,
1509580800, 1509667200, 1509926400, 1510012800, 1510099200, 1510185600,
1510272000, 1510531200, 1510617600, 1510704000, 1510790400, 1510876800,
1511136000, 1511222400, 1511308800, 1511481600, 1511740800, 1511827200,
1511913600, 1512000000, 1512086400, 1512345600, 1512432000, 1512518400,
1512604800, 1512691200, 1512950400, 1513036800, 1513123200, 1513209600,
1513296000, 1513555200, 1513641600, 1513728000, 1513814400, 1513900800,
1514246400, 1514332800, 1514419200, 1514505600, 1514851200, 1514937600,
1515024000, 1515110400, 1515369600, 1515456000, 1515542400, 1515628800,
1515715200, 1516060800, 1516147200, 1516233600, 1516320000, 1516579200,
1516665600, 1516752000, 1516838400, 1516924800, 1517184000, 1517270400,
1517356800, 1517443200, 1517529600, 1517788800, 1517875200, 1517961600,
1518048000, 1518134400, 1518393600, 1518480000, 1518566400, 1518652800,
1518739200, 1519084800, 1519171200, 1519257600, 1519344000, 1519603200,
1519689600, 1519776000, 1519862400, 1519948800, 1520208000, 1520294400,
1520380800, 1520467200, 1520553600, 1520812800, 1520899200, 1520985600,
1521072000, 1521158400, 1521417600, 1521504000, 1521590400, 1521676800,
1521763200, 1522022400, 1522108800, 1522195200, 1522281600, 1522627200,
1522713600, 1522800000, 1522886400, 1522972800, 1523232000, 1523318400,
1523404800, 1523491200, 1523577600, 1523836800, 1523923200, 1524009600,
1524096000, 1524182400, 1524441600, 1524528000, 1524614400, 1524700800,
1524787200, 1525046400, 1525132800, 1525219200, 1525305600, 1525392000,
1525651200, 1525737600, 1525824000, 1525910400, 1525996800, 1526256000,
1526342400, 1526428800, 1526515200, 1526601600, 1526860800, 1526947200,
1527033600, 1527120000, 1527206400, 1527552000, 1527638400, 1527724800,
1527811200, 1528070400, 1528156800, 1528243200, 1528329600, 1528416000,
1528675200, 1528761600, 1528848000, 1528934400, 1529020800, 1529280000,
1529366400, 1529452800, 1529539200, 1529625600, 1529884800, 1529971200,
1530057600, 1530144000, 1530230400, 1530489600, 1530576000, 1530748800,
1530835200, 1531094400, 1531180800, 1531267200, 1531353600, 1531440000,
1531699200, 1531785600, 1531872000, 1531958400, 1532044800, 1532304000,
1532390400, 1532476800, 1532563200, 1532649600, 1532908800, 1532995200,
1533081600, 1533168000, 1533254400, 1533513600, 1533600000, 1533686400,
1533772800, 1533859200, 1534118400, 1534204800, 1534291200, 1534377600,
1534464000, 1534723200, 1534809600, 1534896000, 1534982400, 1535068800,
1535328000, 1535414400, 1535500800, 1535587200, 1535673600, 1536019200,
1536105600, 1536192000, 1536278400, 1536537600, 1536624000, 1536710400,
1536796800, 1536883200, 1537142400, 1537228800, 1537315200, 1537401600,
1537488000, 1537747200, 1537833600, 1537920000, 1538006400, 1538092800,
1538352000, 1538438400, 1538524800, 1538611200, 1538697600, 1538956800,
1539043200, 1539129600, 1539216000, 1539302400, 1539561600, 1539648000,
1539734400, 1539820800, 1539907200, 1540166400, 1540252800, 1540339200,
1540425600, 1540512000, 1540771200, 1540857600, 1540944000, 1541030400,
1541116800, 1541376000, 1541462400, 1541548800, 1541635200, 1541721600,
1541980800, 1542067200, 1542153600, 1542240000, 1542326400, 1542585600,
1542672000, 1542758400, 1542931200, 1543190400, 1543276800, 1543363200,
1543449600, 1543536000, 1543795200, 1543881600, 1544054400, 1544140800,
1544400000, 1544486400, 1544572800, 1544659200, 1544745600, 1545004800,
1545091200, 1545177600, 1545264000, 1545350400, 1545609600, 1545782400,
1545868800, 1545955200, 1546214400, 1546387200, 1546473600, 1546560000,
1546819200, 1546905600, 1546992000, 1547078400, 1547164800, 1547424000,
1547510400, 1547596800, 1547683200, 1547769600, 1548115200, 1548201600,
1548288000, 1548374400, 1548633600, 1548720000, 1548806400, 1548892800,
1548979200, 1549238400, 1549324800, 1549411200, 1549497600, 1549584000,
1549843200, 1549929600, 1550016000, 1550102400, 1550188800, 1550534400,
1550620800, 1550707200, 1550793600, 1551052800, 1551139200, 1551225600,
1551312000, 1551398400, 1551657600, 1551744000, 1551830400, 1551916800,
1552003200, 1552262400, 1552348800, 1552435200, 1552521600, 1552608000,
1552867200, 1552953600, 1553040000, 1553126400, 1553212800, 1553472000,
1553558400, 1553644800, 1553731200, 1553817600, 1554076800, 1554163200,
1554249600, 1554336000, 1554422400, 1554681600, 1554768000, 1554854400,
1554940800, 1555027200, 1555286400, 1555372800, 1555459200, 1555545600,
1555891200, 1555977600, 1556064000, 1556150400, 1556236800, 1556496000,
1556582400, 1556668800, 1556755200, 1556841600, 1557100800, 1557187200,
1557273600, 1557360000, 1557446400, 1557705600, 1557792000, 1557878400,
1557964800, 1558051200, 1558310400, 1558396800, 1558483200, 1558569600,
1558656000, 1559001600, 1559088000, 1559174400, 1559260800, 1559520000,
1559606400, 1559692800, 1559779200, 1559865600, 1560124800, 1560211200,
1560297600, 1560384000, 1560470400, 1560729600, 1560816000, 1560902400,
1560988800, 1561075200, 1561334400, 1561420800, 1561507200, 1561593600,
1561680000, 1561939200, 1562025600, 1562112000, 1562284800, 1562544000,
1562630400, 1562716800, 1562803200, 1562889600, 1563148800, 1563235200,
1563321600, 1563408000, 1563494400, 1563753600, 1563840000, 1563926400,
1564012800, 1564099200, 1564358400, 1564444800, 1564531200, 1564617600,
1564704000, 1564963200, 1565049600, 1565136000, 1565222400, 1565308800,
1565568000, 1565654400, 1565740800, 1565827200, 1565913600, 1566172800,
1566259200, 1566345600, 1566432000, 1566518400, 1566777600, 1566864000,
1566950400, 1567036800, 1567123200, 1567468800, 1567555200, 1567641600,
1567728000, 1567987200, 1568073600, 1568160000, 1568246400, 1568332800,
1568592000, 1568678400, 1568764800, 1568851200, 1568937600, 1569196800,
1569283200, 1569369600, 1569456000, 1569542400, 1569801600, 1569888000,
1569974400, 1570060800, 1570147200, 1570406400, 1570492800, 1570579200,
1570665600, 1570752000, 1571011200, 1571097600, 1571184000, 1571270400,
1571356800, 1571616000, 1571702400, 1571788800, 1571875200, 1571961600,
1572220800, 1572307200, 1572393600, 1572480000, 1572566400, 1572825600,
1572912000, 1572998400, 1573084800, 1573171200, 1573430400, 1573516800,
1573603200, 1573689600, 1573776000, 1574035200, 1574121600, 1574208000,
1574294400, 1574380800, 1574640000, 1574726400, 1574812800, 1574985600,
1575244800, 1575331200, 1575417600, 1575504000, 1575590400, 1575849600,
1575936000, 1576022400, 1576108800, 1576195200, 1576454400, 1576540800,
1576627200, 1576713600, 1576800000, 1577059200, 1577145600, 1577318400,
1577404800, 1577664000), tzone = "UTC", tclass = "Date"), class = c("xts",
"zoo"), src = "yahoo", updated = structure(1585412559.29406, class = c("POSIXct",
"POSIXt")), .Dim = c(837L, 1L), .Dimnames = list(NULL, "Y1"))
The function used in rollapply must produce an ordinary vector. An lm object is not an acceptable output.
Also by default it just passes the coredata of the object into the function but if you want a zoo object passed so that you can get its time then use coredata = FALSE.
lm may have problems with the zoo object passed so use its coredata.
In lmFunction shown in the question it is trying to use lm(..., x) as if x were a data frame but x is not a data frame.
The two examples below return fitted values and coefficients respectively.
Plotting a line for each time point is going to look like a real mess with all those lines but we show how to produce it anyways at the end.
library(xts)
Y1 <- as.zoo(data)$Y1
fitfun <- function(x) fitted(lm(coredata(x) ~ time(x)))
r1 <- rollapplyr(Y1, 40, fitfun, coredata = FALSE, fill = NA)
dim(r1)
## [1] 837 40
coeffun <- function(x) coef(lm(coredata(x) ~ time(x)))
r2 <- rollapplyr(Y1, 40, coeffun, coredata = FALSE, fill = NA)
dim(r2)
## [1] 837 2
# plot
r2na <- na.omit(r2)
plot(Y1)
junk <- Map(abline, a = r2na[, 1], b = r2na[, 2])
Related
I used following code to generate the sets:
H=list('ASD'=c("SCN3A", "TSC2", "RASSF8", "XDH", "CNTN5", "CHD7", "DPYD", "BUB1", "DPP4", "AMPD1", "ARID4A", "QSER1", "OCA2", "GPR141", "DNAH9", "LOXHD1", "SETX", "GSTP1", "SLC6A15", "ASTN1", "CACNA1D", "WFS1", "KCNC1", "TPRA1", "NBAS", "LAMB1", "TSNARE1", "GABRB3", "PAG1", "SOS1", "LHFPL4", "KMT2C", "GATM", "NPTX2", "ANK3", "SH3TC2", "CLPP", "TP53BP1", "EBF4", "ANO1", "DNAJC10", "TTC25", "NDUFB7", "PCSK1", "ROBO3", "MYO9A", "RPS23", "ZNF292", "CRBN", "HIVEP2", "CYP1B1", "KCNE2", "NRP2", "ADSL", "XIRP2", "SMARCA4", "HDAC2", "JAK2", "ATP9B", "ESR1", "MDN1", "CTNNA2", "CNTNAP2", "SPTB", "ACY1", "CEP290", "WNT2", "IQGAP3", "CACNG3", "CACNA1H", "FOXN1", "SPAST", "ADCY1", "IGF2", "RNF4", "CTTNBP2", "TMF1", "CPT1A", "TENM3", "RABGGTA", "MBD3", "CCDC18", "SMG1", "SYNE1", "PITX3", "HMGA2", "PDE6C", "CTU2", "PNPLA7", "PIEZO1", "MAP2", "MCF2L", "PRF1", "VWA8", "LRP2", "RYR3", "IPO9", "DIXDC1", "CNTNAP5", "ELP4", "LRRN1", "TRIO", "SACS", "ASIC2", "KIF1B", "AXIN2", "SORL1", "DLG2", "VPS13C", "F2RL1", "NRP1", "LAMA1", "LARP1", "SDK1", "LAMC2", "LSG1", "EP300", "DGCR6L", "TSHB", "APC", "TENM4", "IGSF21"),
'Epilepsy'=c("ADCY1", "ADSL", "ANK3", "APC", "ASIC2", "ASTN1", "ATP9B", "AXIN2", "BUB1", "CACNA1D", "CACNA1H", "CACNG3", "CEP290", "CLPP", "CNTN5", "CNTNAP2", "CNTNAP5", "CTNNA2", "CTU2", "DGCR6L", "DLG2", "DNAJC10", "ELP4", "EP300", "GABRB3", "HDAC2", "HMGA2", "IGSF21", "IPO9", "IQGAP3", "JAK2", "KCNC1", "KCNE2", "KIF1B", "LAMB1", "LAMC2", "LARP1", "LHFPL4", "MCF2L", "MYO9A", "NDUFB7", "NRP2", "RABGGTA", "RYR3", "SACS", "SCN3A", "SETX", "SMARCA4", "SMG1", "SORL1", "SPAST", "SPTB", "TMF1", "TP53BP1", "TPRA1", "TRIO", "TSC2", "TSHB", "VPS13C", "WFS1", "XDH", "ZNF292", "ABCA1", "ABCG5", "ACP2", "ADCY6", "ADGRA3", "ADGRV1", "AGA", "AGRN", "ALG3", "ANK2", "ANKRD44", "AP5Z1", "ARAP1", "ARHGEF11", "ASCC3", "ASPM", "ATRN", "AURKB", "BRINP2", "BRSK2", "CAPN15", "CAPN2", "CC2D2A", "CD36", "CDX1", "CHAT", "CHM", "CHTF18", "CLCN2", "CLCN7", "CLEC12A", "COG5", "COG6", "COL12A1", "COL17A1", "COL28A1", "COL6A2", "COQ3", "CYBB", "DBNL", "DBR1", "DEAF1", "DGKQ", "DHX29", "DMXL1", "DMXL2", "DNAH12", "DOCK7", "DOLK", "DSG4", "DZIP1", "EFHC1", "EPHA7", "EPHB1", "EPHB3", "ESYT1", "EXOC7", "F13A1", "F5", "FAH", "FARP1", "FASN", "FAT1", "FERMT1", "FGFR3", "FLNB", "FLRT1", "FMNL1", "FN1", "FOLR1", "FZD4", "GAA", "GAK", "GALC", "GANAB", "GLB1", "GM2A", "GNA14", "GOSR2", "GPC1", "GSS", "GTPBP4", "HEATR1", "HGS", "HIPK2", "HSPA13", "HTR2B", "INO80E", "KALRN", "KCNH2", "KCNJ16", "KIAA0586", "KRIT1", "LAMB2", "LMAN1", "LTBP3", "MASTL", "MED17", "METTL25", "MGAT5B", "MPL", "MRPL24", "MRPS27", "MSX2", "MUS81", "NALCN", "NEB", "NFE2L1", "NMBR", "NOL11", "NUCB1", "OBSL1", "OGDHL", "OGFOD2", "PARP1", "PCCB", "PDGFRA", "PDYN", "PFAS", "PFKFB3", "PLEC", "PLOD2", "PNKD", "PNPO", "POLD1", "POLE", "POLL", "PPAT", "PRDX3", "PSEN2", "PTCH1", "PTK2B", "RAB5A", "RAD54B", "RELB", "RHOT2", "RIMS2", "RTN4IP1", "SAE1", "SCAPER", "SCN1B", "SCN5A", "SEC23IP", "SERPINE2", "SGCB", "SIGLEC1", "SLC18A3", "SLC26A4", "SLC2A1", "SLC4A1", "SLC4A3", "SLITRK6", "SPG11", "STARD13", "STT3A", "SYF2", "SYNJ1", "SZT2", "TECPR2", "TES", "TMEM67", "TNK2", "TNS1", "TNS3", "TRAP1", "TRAPPC11", "TRAPPC8", "TRIP10", "TRPM2", "TRPM3", "TRPV1", "TUBGCP6", "TYW1", "UNC13B", "USP45", "VCP", "WDR62", "WWOX", "YAP1", "ZBTB48", "ZFHX4", "ZW10"),
'Intellectual'=c("YAP1", "TSC2", "FMNL1", "MADD", "KRIT1", "TNK2", "NUP155", "COG6", "TRAPPC8", "IPO9", "CCT8", "ZFYVE16", "FAT1", "SLC37A1", "KAT6B", "WWOX", "PTK2B", "CNTNAP5", "HEATR1", "TRAPPC11", "ARID4A", "ADGRV1", "GPC1", "AURKB", "TTC39B", "GSS", "GANAB", "RFWD3", "CUL9", "SAE1", "CACNG3", "SMC4", "ASIC2", "SYF2", "PFAS", "GAK", "METTL25", "OGDHL", "RAB5A", "VPS13C", "FGFR3", "KCNC1", "DDOST", "CTNND1", "ATRN", "UNC13B", "FASN", "MYO18A", "TMF1", "CRBN", "FLNB", "LARP1", "RCC2", "MASTL", "ADGRA3", "SCAPER", "SZT2", "LRP4", "SMG1", "SYNE1", "ADCY6", "RNGTT", "EP300", "ARAP1", "ANKRD44", "VPS39", "DGCR6L", "STT3A", "COG5", "NFE2L1", "SLC7A2", "ASCC3", "LMAN1", "OBSL1", "KIF20B", "PNPLA7"),
'Schizophrenia'=c("CRB1", "FBXW4", "CNTN5", "CHM", "DPYD", "BUB1", "DPP4", "ADGRV1", "GPC1", "DLD", "RBM19", "ITIH1", "ADCY7", "KCNC1", "ACBD5", "ELF1", "TSNARE1", "GABRB3", "TMEM182", "STT3A", "KMT2C", "ANK3", "SAMD4A", "NALCN", "PDLIM5", "LIMA1", "DYNC1I2", "TMEM26", "ANK2", "GRIN3A", "EPHB1", "BRINP2", "KALRN", "MICAL2", "TTLL6", "ZNF292", "SERPINB1", "SCAPER", "KCNS3", "VAMP1", "HDAC2", "DLG1", "CFAP58", "JAK2", "ESR1", "EFCAB14", "CNTNAP2", "PENK", "WWOX", "GLT8D1", "SFMBT2", "EEFSEC", "PDZD8", "GALC", "CACNG3", "SPG7", "ADCY1", "CLDN10", "TENM3", "OTOS", "CAPN2", "LRP4", "SYNE1", "HMGA2", "APOL4", "ARAP1", "DOPEY1", "EPN2", "MCF2L", "PLAA", "RYR3", "CNTNAP5", "POLE", "TRIO", "RAB11FIP5", "BNIP3L", "RUNX1", "AXIN2", "KIF1B", "DLG2", "MAP3K9", "CERCAM", "VPS13C", "MYO18B", "NRP1", "IRS1", "MANEA", "EP300", "DGCR6L", "TSHB", "AASDHPPT", "TRAPPC3", "APC", "TENM4"),
'ADHD'=c("TSC2", "CSMD2", "ANK3", "ASPM", "CNTN5", "TP53BP1", "CNTNAP2", "CHD7", "GLT8D1", "CNTNAP5", "SFMBT2", "MLH1", "FLII", "PTCH1", "SPRED1", "ARHGEF12", "CHMP7", "SVEP1", "CREB5", "MSH6", "TTLL8", "DLG2", "MICAL2", "DNAH3", "ADGRA3", "ACBD5", "SYNE1", "NUCB1", "NRP2", "TMEM182", "STT3A", "GNE", "PAX5", "TENM4"),
'Biopolar Disoder'=c("CSMD2", "ANK3", "NUMB", "CNTN5", "CTNNA2", "CNTNAP2", "NALCN", "RDH14", "WWOX", "GLT8D1", "SFMBT2", "SPRED1", "ADAMTS14", "TTC39B", "SVEP1", "COL28A1", "PDZD8", "GALC", "SORL1", "ITIH1", "CLSPN", "SLCO4A1", "VPS13C", "MYO18B", "TTLL6", "CRBN", "NFATC2", "SYNE1", "INTS7", "TMEM182", "STT3A", "WSCD1", "KCNS3", "PDGFRA"),
'circadian'=c("YTHDC2", "LRP1B", "SAMD4A", "VTA1", "NALCN", "F13A1", "SLC4A1", "TMEM26", "ADAMTS14", "COL28A1", "RNF19A", "BRINP2", "TRHR", "CCBE1", "CACNA2D3", "MYO18A", "SLC26A7", "FRAS1", "EP300", "GNE", "LMAN1", "NMBR", "TGFBI", "CHAT", "HIPK2", "SERPINE2", "RELB", "ARHGEF12", "SETX", "CACNA1D", "ADCY7", "NBAS", "GPR83", "LTBP3", "GABRB3", "HMBS", "PSMA4", "MADD", "SIX3", "EPHB1", "GAK", "SPRY4", "ABCA13", "SLC8A3", "HIVEP2", "FAH", "ADCY6", "MYLIP", "NFE2L1", "PRKAG1", "HDAC2", "TNS1", "GSS", "CACNA1H", "POLL", "CPT1A", "CAPN2", "PITX3", "NFE2L3", "PNPLA7", "AGA", "STX16", "TNK2", "PRF1", "RYR3", "DIXDC1", "PSEN2", "SREBF1", "POLE", "RAB11FIP5", "KIF1B", "RAB5A", "NGB", "ARHGEF11", "SIX6", "SDK1", "AATK", "CYP7A1", "CUL9"),
'dementia'=c("LMAN1", "PRDX3", "NPLOC4", "PSEN2", "ZFYVE16", "CHAT", "PENK", "LTBP1", "SNCAIP", "NOL11", "GANAB", "DLD", "SPG11", "IGSF21", "RNF19A", "ASIC2", "ASTN1", "SORL1", "SPRY4", "VCP", "TATDN1", "RAB5A", "VPS13C", "GM2A", "IRS1", "PSMB4", "NGB", "F2RL1", "RAB38", "SMOX", "PARP1", "ABCC9", "RPN2", "CAPN2", "GALNT7", "SMG1", "TBK1", "MASP2", "GNE", "GTPBP4", "CCT8", "CD63", "PSMA4", "TUBA4A", "CUL9")
)
venn(H)
Throws me an error:
Error in drawVennDiagram(data = counts, small = small, showSetLogicLabel = showSetLogicLabel, : Venn diagrams for 8 dimensions are not yet supported.
Here is a solution with package venneuler.
The data must first be restructured. A lapply loop will create data.frames with element name and set as 1st and 2nd columns, respectively. Then they are rbinded and the diagram is created.
library(venneuler)
Hmat <- lapply(seq_along(H), \(i) {
Set <- names(H)[i]
data.frame(Element = H[[i]], Set)
})
Hmat <- do.call(rbind, Hmat)
vd <- venneuler(Hmat)
plot(vd)
A double sapply loop is a way to get the overlap counts.
overlap <- function(x, y) length(intersect(x, y))
sapply(H, function(h1) sapply(H, function(h2) overlap(h1, h2)))
Using ggvenn package,
library(ggvenn)
ggvenn(H)
I developed a dashboard to show the dummy example of the problem that I am running. The shiny dashboard has four boxes. Top two boxes display a plotly chart and forecasted data using non-seasonal arima model. The bottom two boxes demonstrate a plotly chart and forecasted data using seasonal arima model. I also have a radio button for displaying items in either seasonal or non-seasonal boxes. I want the appropriate boxes to be displayed only when I click submit button. For example: if I check seasonal option in radio-button and click predict, I would like to see the seasonal chart and seasonally forecasted data. Currently, the dashboard display the boxes only based on the selected option in the radio-button.
skin <- Sys.getenv("DASHBOARD_SKIN")
skin <- tolower(skin)
if (skin == "")
skin <- "black"
sidebar <- dashboardSidebar(
sidebarSearchForm(label = "Search...", "searchText", "searchButton"),
sidebarMenu(
menuItem("Dashboard", tabName = "dashboard", icon = icon("dashboard"))
)
)
)
body <- dashboardBody(
tabItems(
tabItem("dashboard",
# Boxes with solid headers
fluidRow(
box(
title = "Enter Stock Code", width = 4, solidHeader = TRUE, status = "primary",
radioButtons("seasonal", "Select", c(NonSeasonal = "NonSeasonal", Seasonal = "Seasonal")),
actionButton(inputId = "click", label = "Predict")
)
),
fluidRow(
box(
title = "Auto Arima - Non Seasonal",
status = "primary",
uiOutput("arima_nonseasonal", height = 350),
height = 470
),
box(
title = "Auto Arima - Non Seasonal",
width = "6",
tableOutput("arima_nonseasonal_Data"),
height = 470
)
),
fluidRow(
box(
title = "Auto Arima Seasonal",
status = "primary",
uiOutput("arima.seasonal", height = 350),
height = 470
),
box(
title = "Auto Arima Seasonal",
width = "6",
tableOutput("arima.seasonalData"),
height = 470
)
)
)
)
)
header <- dashboardHeader(
title = "Executive Dashboard"
)
ui <- dashboardPage(header, sidebar, body, skin = skin)
server <- function(input, output, session) {
output$nonseasonalplot <- renderPlotly({
qxts = structure(list(y=c(-0.4120000005, -0.452666665083333, -0.458666667316667, -0.460666666416667, -0.505333344133333, -0.506666670233333, -0.50999999045, -0.498666668933333, -0.502666672066667, -0.505333344133333, -0.50200000405, -0.506666680166667, -0.508000006266667, -0.5040000081, -0.501333341, -0.50200000405, -0.502000009016667, -0.504000003133333, -0.50466667115, -0.505333344133333, -0.501333341, -0.503333340083333, -0.496666664883333, -0.5040000081, -0.49799999595, -0.389333337533333, -0.3640000075, -0.365333338566667, -0.3586666733, -0.3599999994, -0.355333338183333, -0.365333338566667, -0.360666667416667, -0.3640000075, -0.3546666652, -0.369999999783333, -0.3653333336, -0.36733333765, -0.396666670833333, -0.469333335766667, -0.4586666673, -0.4639999916, -0.457333331283333, -0.4646666596, -0.4680000047, -0.4687999964, -0.467333336683333, -0.460666661433333, -0.4653333326, -0.4720000029, -0.453333328166667, -0.466000000616667, -0.4646666646, -0.45933333535, -0.464666659633333, -0.4579999993, -0.464666669566667, -0.473999996983333, -0.417333334666667, -0.46400000155, -0.4599999934, -0.46866666775, -0.47200000288, -0.465333337583333, -0.4719999979, -0.462666660516667, -0.46240000126, -0.465333322666667, -0.4673333317, -0.459999993416667, -0.464666659616667, -0.461333329466667, -0.46599999565, -0.461999997483333, -0.461999992516667, -0.466666663683333, -0.457333336283333, -0.46600000065, -0.463999996566667, -0.46400000155, -0.488666658566667, -0.543333326766667, -0.489333341533333, -0.3600000143, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -0.3573333323, -0.464666662133333, -0.3360000037, -0.445333333066667, -0.49333334465, -0.444666661316667, -0.325333339483333, -0.60199998815, -0.33039999454, -0.394666674216667, -0.406666671233333, -0.4066666588, -0.52533333005, -0.552000006066667, -0.6973333458, -0.526000003, -0.565333336583333, -0.451428575185714, -0.3973333314, -0.591999994871429, -0.49120000004, -0.4313333432, -0.197333325966667, -0.545999998833333, -0.49600000976, -0.4595000035875, -0.45920000072, -0.453714281285714, -0.455333337166667, -0.461333334433333, -0.461999992533333, -0.457333331283333, -0.4606666664, -0.456000005183333, -0.4593333254, -0.4579999993, -0.45679999586, -0.455428570485714, -0.454666664183333, -0.453999991183333, -0.452666665083333, -0.464000001533333, -0.456666668216667, -0.45439999696, -0.45666666825, -0.454666664166667, -0.459999998383333, -0.45840000508, -0.459333330366667, -0.394666661816667, -0.367333332983333, -0.2823999971, -0.3446666648, -0.464666659633333, -0.443333327766667, -0.33920000494, -0.332666660333333, -0.363333340716667, -0.43599999546, -0.518000001716667, -0.15933333585, -0.2493333357, -0.43359999356, -0.462666669166667, -0.410000008833333, -0.357333337883333, -0.3615999937, -0.494666665783333, -0.52200000485, -0.29400000225, -0.4527999997, -0.474000004416667, -0.342666663233333, -0.53680000306, -0.423333331933333, -0.498666671416667, -0.439333344516667, -0.46480000914, -0.368666668733333, -0.2673333399, -0.303333333366667, -0.21279999912, -0.398000004383333, -0.444666676216667, -0.63679999114, -0.337333336466667, -0.2599999905 ),
date = c(1493139120, 1493139180, 1493139240, 1493139300, 1493139360, 1493139420, 1493139480, 1493139540, 1493139600, 1493139660, 1493139720, 1493139780, 1493139840, 1493139900, 1493139960, 1493140020, 1493140080, 1493140140, 1493140200, 1493140260, 1493140320, 1493140380, 1493140440, 1493140500, 1493140560, 1493140620, 1493140680, 1493140740, 1493140800, 1493140860, 1493140920, 1493140980, 1493141040, 1493141100, 1493141160, 1493141220, 1493141280, 1493141340, 1493141400, 1493141460, 1493141520, 1493141580, 1493141640, 1493141700, 1493141760, 1493141820, 1493141880, 1493141940, 1493142000, 1493142060, 1493142120, 1493142180, 1493142240, 1493142300, 1493142360, 1493142420, 1493142480, 1493142540, 1493142600, 1493142660, 1493142720, 1493142780, 1493142840, 1493142900, 1493142960, 1493143020, 1493143080, 1493143140, 1493143200, 1493143260, 1493143320, 1493143380, 1493143440, 1493143500, 1493143560, 1493143620, 1493143680, 1493143740, 1493143800, 1493143860, 1493143920, 1493143980, 1493144040, 1493144100, 1493144160, 1493144220, 1493144280, 1493144340, 1493144400, 1493144460, 1493144520, 1493144580, 1493144640, 1493144700, 1493144760, 1493144820, 1493144880, 1493144940, 1493145000, 1493145060, 1493145120, 1493145180, 1493145240, 1493145300, 1493145360, 1493145420, 1493145480, 1493145540, 1493145600, 1493145660, 1493145720, 1493145780, 1493145840, 1493145900, 1493145960, 1493146020, 1493146080, 1493146140, 1493146200, 1493146260, 1493146320, 1493146380, 1493146440, 1493146500, 1493146560, 1493146620, 1493146680, 1493146740, 1493146800, 1493146860, 1493146920, 1493146980, 1493147040, 1493147100, 1493147160, 1493147220, 1493147280, 1493147340, 1493147400, 1493147460, 1493147520, 1493147580, 1493147640, 1493147700, 1493147760, 1493147820, 1493147880, 1493147940, 1493148000, 1493148060, 1493148120, 1493148180, 1493148240, 1493148300, 1493148360, 1493148420, 1493148480, 1493148540, 1493148600, 1493148660, 1493148720, 1493148780, 1493148840, 1493148900, 1493148960, 1493149020, 1493149080, 1493149140, 1493149200, 1493149260, 1493149320, 1493149380, 1493149440, 1493149500, 1493149560, 1493149620, 1493149680, 1493149740, 1493149800, 1493149860, 1493149920, 1493149980, 1493150040, 1493150100, 1493150160, 1493150220, 1493150280, 1493150340, 1493150400, 1493150460, 1493150520))
)
fit.xts <- auto.arima(qxts$y, seasonal = FALSE)
forecast_length <- 10
fore.xts <- forecast(fit.xts, h=forecast_length)
fore.dates <- seq(as.POSIXct(qxts$date[length(qxts$date)], origin="1970-01-01"),
by=qxts$date[length(qxts$date)] - qxts$date[length(qxts$date)-1], len=forecast_length)
nonseasonal_p <- plot_ly() %>%
add_lines(x = as.POSIXct(qxts$date, origin="1970-01-01"), y = qxts$y,
color = I("black"),
name = "observed",
marker=list(mode="lines")) %>%
add_lines(x = fore.dates, y = fore.xts$mean, color = I("blue"), name = "prediction") %>%
add_ribbons(x = fore.dates,
ymin = fore.xts$lower[, 2],
ymax = fore.xts$upper[, 2],
color = I("gray95"),
name = "95% confidence") %>%
add_ribbons(x = fore.dates,
ymin = fore.xts$lower[, 1],
ymax = fore.xts$upper[, 1],
color = I("gray80"), name = "80% confidence")
return(nonseasonal_p)
})
#Nonseasonal plot Auto.Arima - plot here Tile#4
output$arima_nonseasonal <- renderUI({
if (input$seasonal == "NonSeasonal"){
plotlyOutput("nonseasonalplot")
}
})
nonseasonaltable <- reactive({
qxts = structure(list(y=c(-0.4120000005, -0.452666665083333, -0.458666667316667, -0.460666666416667, -0.505333344133333, -0.506666670233333, -0.50999999045, -0.498666668933333, -0.502666672066667, -0.505333344133333, -0.50200000405, -0.506666680166667, -0.508000006266667, -0.5040000081, -0.501333341, -0.50200000405, -0.502000009016667, -0.504000003133333, -0.50466667115, -0.505333344133333, -0.501333341, -0.503333340083333, -0.496666664883333, -0.5040000081, -0.49799999595, -0.389333337533333, -0.3640000075, -0.365333338566667, -0.3586666733, -0.3599999994, -0.355333338183333, -0.365333338566667, -0.360666667416667, -0.3640000075, -0.3546666652, -0.369999999783333, -0.3653333336, -0.36733333765, -0.396666670833333, -0.469333335766667, -0.4586666673, -0.4639999916, -0.457333331283333, -0.4646666596, -0.4680000047, -0.4687999964, -0.467333336683333, -0.460666661433333, -0.4653333326, -0.4720000029, -0.453333328166667, -0.466000000616667, -0.4646666646, -0.45933333535, -0.464666659633333, -0.4579999993, -0.464666669566667, -0.473999996983333, -0.417333334666667, -0.46400000155, -0.4599999934, -0.46866666775, -0.47200000288, -0.465333337583333, -0.4719999979, -0.462666660516667, -0.46240000126, -0.465333322666667, -0.4673333317, -0.459999993416667, -0.464666659616667, -0.461333329466667, -0.46599999565, -0.461999997483333, -0.461999992516667, -0.466666663683333, -0.457333336283333, -0.46600000065, -0.463999996566667, -0.46400000155, -0.488666658566667, -0.543333326766667, -0.489333341533333, -0.3600000143, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -0.3573333323, -0.464666662133333, -0.3360000037, -0.445333333066667, -0.49333334465, -0.444666661316667, -0.325333339483333, -0.60199998815, -0.33039999454, -0.394666674216667, -0.406666671233333, -0.4066666588, -0.52533333005, -0.552000006066667, -0.6973333458, -0.526000003, -0.565333336583333, -0.451428575185714, -0.3973333314, -0.591999994871429, -0.49120000004, -0.4313333432, -0.197333325966667, -0.545999998833333, -0.49600000976, -0.4595000035875, -0.45920000072, -0.453714281285714, -0.455333337166667, -0.461333334433333, -0.461999992533333, -0.457333331283333, -0.4606666664, -0.456000005183333, -0.4593333254, -0.4579999993, -0.45679999586, -0.455428570485714, -0.454666664183333, -0.453999991183333, -0.452666665083333, -0.464000001533333, -0.456666668216667, -0.45439999696, -0.45666666825, -0.454666664166667, -0.459999998383333, -0.45840000508, -0.459333330366667, -0.394666661816667, -0.367333332983333, -0.2823999971, -0.3446666648, -0.464666659633333, -0.443333327766667, -0.33920000494, -0.332666660333333, -0.363333340716667, -0.43599999546, -0.518000001716667, -0.15933333585, -0.2493333357, -0.43359999356, -0.462666669166667, -0.410000008833333, -0.357333337883333, -0.3615999937, -0.494666665783333, -0.52200000485, -0.29400000225, -0.4527999997, -0.474000004416667, -0.342666663233333, -0.53680000306, -0.423333331933333, -0.498666671416667, -0.439333344516667, -0.46480000914, -0.368666668733333, -0.2673333399, -0.303333333366667, -0.21279999912, -0.398000004383333, -0.444666676216667, -0.63679999114, -0.337333336466667, -0.2599999905 ),
date = c(1493139120, 1493139180, 1493139240, 1493139300, 1493139360, 1493139420, 1493139480, 1493139540, 1493139600, 1493139660, 1493139720, 1493139780, 1493139840, 1493139900, 1493139960, 1493140020, 1493140080, 1493140140, 1493140200, 1493140260, 1493140320, 1493140380, 1493140440, 1493140500, 1493140560, 1493140620, 1493140680, 1493140740, 1493140800, 1493140860, 1493140920, 1493140980, 1493141040, 1493141100, 1493141160, 1493141220, 1493141280, 1493141340, 1493141400, 1493141460, 1493141520, 1493141580, 1493141640, 1493141700, 1493141760, 1493141820, 1493141880, 1493141940, 1493142000, 1493142060, 1493142120, 1493142180, 1493142240, 1493142300, 1493142360, 1493142420, 1493142480, 1493142540, 1493142600, 1493142660, 1493142720, 1493142780, 1493142840, 1493142900, 1493142960, 1493143020, 1493143080, 1493143140, 1493143200, 1493143260, 1493143320, 1493143380, 1493143440, 1493143500, 1493143560, 1493143620, 1493143680, 1493143740, 1493143800, 1493143860, 1493143920, 1493143980, 1493144040, 1493144100, 1493144160, 1493144220, 1493144280, 1493144340, 1493144400, 1493144460, 1493144520, 1493144580, 1493144640, 1493144700, 1493144760, 1493144820, 1493144880, 1493144940, 1493145000, 1493145060, 1493145120, 1493145180, 1493145240, 1493145300, 1493145360, 1493145420, 1493145480, 1493145540, 1493145600, 1493145660, 1493145720, 1493145780, 1493145840, 1493145900, 1493145960, 1493146020, 1493146080, 1493146140, 1493146200, 1493146260, 1493146320, 1493146380, 1493146440, 1493146500, 1493146560, 1493146620, 1493146680, 1493146740, 1493146800, 1493146860, 1493146920, 1493146980, 1493147040, 1493147100, 1493147160, 1493147220, 1493147280, 1493147340, 1493147400, 1493147460, 1493147520, 1493147580, 1493147640, 1493147700, 1493147760, 1493147820, 1493147880, 1493147940, 1493148000, 1493148060, 1493148120, 1493148180, 1493148240, 1493148300, 1493148360, 1493148420, 1493148480, 1493148540, 1493148600, 1493148660, 1493148720, 1493148780, 1493148840, 1493148900, 1493148960, 1493149020, 1493149080, 1493149140, 1493149200, 1493149260, 1493149320, 1493149380, 1493149440, 1493149500, 1493149560, 1493149620, 1493149680, 1493149740, 1493149800, 1493149860, 1493149920, 1493149980, 1493150040, 1493150100, 1493150160, 1493150220, 1493150280, 1493150340, 1493150400, 1493150460, 1493150520))
)
fit.xts <- auto.arima(qxts$y, seasonal = FALSE)
forecast_length <- 10
fore.xts <- forecast(fit.xts, h=forecast_length)
return(fore.xts)
return(nonseasonal_nwforcast)
})
#Non seasonal dataAuto.Arima1 - plot here Tile#5
output$arima_nonseasonal_Data <- renderTable({
if (input$seasonal == "NonSeasonal"){
nonseasonaltable()
}
})
output$seasonalplot<- renderPlotly({
qxts = structure(list(y=c(-0.4120000005, -0.452666665083333, -0.458666667316667, -0.460666666416667, -0.505333344133333, -0.506666670233333, -0.50999999045, -0.498666668933333, -0.502666672066667, -0.505333344133333, -0.50200000405, -0.506666680166667, -0.508000006266667, -0.5040000081, -0.501333341, -0.50200000405, -0.502000009016667, -0.504000003133333, -0.50466667115, -0.505333344133333, -0.501333341, -0.503333340083333, -0.496666664883333, -0.5040000081, -0.49799999595, -0.389333337533333, -0.3640000075, -0.365333338566667, -0.3586666733, -0.3599999994, -0.355333338183333, -0.365333338566667, -0.360666667416667, -0.3640000075, -0.3546666652, -0.369999999783333, -0.3653333336, -0.36733333765, -0.396666670833333, -0.469333335766667, -0.4586666673, -0.4639999916, -0.457333331283333, -0.4646666596, -0.4680000047, -0.4687999964, -0.467333336683333, -0.460666661433333, -0.4653333326, -0.4720000029, -0.453333328166667, -0.466000000616667, -0.4646666646, -0.45933333535, -0.464666659633333, -0.4579999993, -0.464666669566667, -0.473999996983333, -0.417333334666667, -0.46400000155, -0.4599999934, -0.46866666775, -0.47200000288, -0.465333337583333, -0.4719999979, -0.462666660516667, -0.46240000126, -0.465333322666667, -0.4673333317, -0.459999993416667, -0.464666659616667, -0.461333329466667, -0.46599999565, -0.461999997483333, -0.461999992516667, -0.466666663683333, -0.457333336283333, -0.46600000065, -0.463999996566667, -0.46400000155, -0.488666658566667, -0.543333326766667, -0.489333341533333, -0.3600000143, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -0.3573333323, -0.464666662133333, -0.3360000037, -0.445333333066667, -0.49333334465, -0.444666661316667, -0.325333339483333, -0.60199998815, -0.33039999454, -0.394666674216667, -0.406666671233333, -0.4066666588, -0.52533333005, -0.552000006066667, -0.6973333458, -0.526000003, -0.565333336583333, -0.451428575185714, -0.3973333314, -0.591999994871429, -0.49120000004, -0.4313333432, -0.197333325966667, -0.545999998833333, -0.49600000976, -0.4595000035875, -0.45920000072, -0.453714281285714, -0.455333337166667, -0.461333334433333, -0.461999992533333, -0.457333331283333, -0.4606666664, -0.456000005183333, -0.4593333254, -0.4579999993, -0.45679999586, -0.455428570485714, -0.454666664183333, -0.453999991183333, -0.452666665083333, -0.464000001533333, -0.456666668216667, -0.45439999696, -0.45666666825, -0.454666664166667, -0.459999998383333, -0.45840000508, -0.459333330366667, -0.394666661816667, -0.367333332983333, -0.2823999971, -0.3446666648, -0.464666659633333, -0.443333327766667, -0.33920000494, -0.332666660333333, -0.363333340716667, -0.43599999546, -0.518000001716667, -0.15933333585, -0.2493333357, -0.43359999356, -0.462666669166667, -0.410000008833333, -0.357333337883333, -0.3615999937, -0.494666665783333, -0.52200000485, -0.29400000225, -0.4527999997, -0.474000004416667, -0.342666663233333, -0.53680000306, -0.423333331933333, -0.498666671416667, -0.439333344516667, -0.46480000914, -0.368666668733333, -0.2673333399, -0.303333333366667, -0.21279999912, -0.398000004383333, -0.444666676216667, -0.63679999114, -0.337333336466667, -0.2599999905 ),
date = c(1493139120, 1493139180, 1493139240, 1493139300, 1493139360, 1493139420, 1493139480, 1493139540, 1493139600, 1493139660, 1493139720, 1493139780, 1493139840, 1493139900, 1493139960, 1493140020, 1493140080, 1493140140, 1493140200, 1493140260, 1493140320, 1493140380, 1493140440, 1493140500, 1493140560, 1493140620, 1493140680, 1493140740, 1493140800, 1493140860, 1493140920, 1493140980, 1493141040, 1493141100, 1493141160, 1493141220, 1493141280, 1493141340, 1493141400, 1493141460, 1493141520, 1493141580, 1493141640, 1493141700, 1493141760, 1493141820, 1493141880, 1493141940, 1493142000, 1493142060, 1493142120, 1493142180, 1493142240, 1493142300, 1493142360, 1493142420, 1493142480, 1493142540, 1493142600, 1493142660, 1493142720, 1493142780, 1493142840, 1493142900, 1493142960, 1493143020, 1493143080, 1493143140, 1493143200, 1493143260, 1493143320, 1493143380, 1493143440, 1493143500, 1493143560, 1493143620, 1493143680, 1493143740, 1493143800, 1493143860, 1493143920, 1493143980, 1493144040, 1493144100, 1493144160, 1493144220, 1493144280, 1493144340, 1493144400, 1493144460, 1493144520, 1493144580, 1493144640, 1493144700, 1493144760, 1493144820, 1493144880, 1493144940, 1493145000, 1493145060, 1493145120, 1493145180, 1493145240, 1493145300, 1493145360, 1493145420, 1493145480, 1493145540, 1493145600, 1493145660, 1493145720, 1493145780, 1493145840, 1493145900, 1493145960, 1493146020, 1493146080, 1493146140, 1493146200, 1493146260, 1493146320, 1493146380, 1493146440, 1493146500, 1493146560, 1493146620, 1493146680, 1493146740, 1493146800, 1493146860, 1493146920, 1493146980, 1493147040, 1493147100, 1493147160, 1493147220, 1493147280, 1493147340, 1493147400, 1493147460, 1493147520, 1493147580, 1493147640, 1493147700, 1493147760, 1493147820, 1493147880, 1493147940, 1493148000, 1493148060, 1493148120, 1493148180, 1493148240, 1493148300, 1493148360, 1493148420, 1493148480, 1493148540, 1493148600, 1493148660, 1493148720, 1493148780, 1493148840, 1493148900, 1493148960, 1493149020, 1493149080, 1493149140, 1493149200, 1493149260, 1493149320, 1493149380, 1493149440, 1493149500, 1493149560, 1493149620, 1493149680, 1493149740, 1493149800, 1493149860, 1493149920, 1493149980, 1493150040, 1493150100, 1493150160, 1493150220, 1493150280, 1493150340, 1493150400, 1493150460, 1493150520))
)
fit.xts <- auto.arima(qxts$y, seasonal = TRUE)
forecast_length <- 10
fore.xts <- forecast(fit.xts, h=forecast_length)
fore.dates <- seq(as.POSIXct(qxts$date[length(qxts$date)], origin="1970-01-01"),
by=qxts$date[length(qxts$date)] - qxts$date[length(qxts$date)-1], len=forecast_length)
seasonal_p <- plot_ly() %>%
add_lines(x = as.POSIXct(qxts$date, origin="1970-01-01"), y = qxts$y,
color = I("black"),
name = "observed",
marker=list(mode="lines")) %>%
add_lines(x = fore.dates, y = fore.xts$mean, color = I("blue"), name = "prediction") %>%
add_ribbons(x = fore.dates,
ymin = fore.xts$lower[, 2],
ymax = fore.xts$upper[, 2],
color = I("gray95"),
name = "95% confidence") %>%
add_ribbons(x = fore.dates,
ymin = fore.xts$lower[, 1],
ymax = fore.xts$upper[, 1],
color = I("gray80"), name = "80% confidence")
return(seasonal_p)
})
#Auto.Arima Seasonal
output$arima.seasonal <- renderUI({
if (input$seasonal == "Seasonal"){
plotlyOutput("seasonalplot")
}
})
seasonaldata <- reactive({
qxts = structure(list(y=c(-0.4120000005, -0.452666665083333, -0.458666667316667, -0.460666666416667, -0.505333344133333, -0.506666670233333, -0.50999999045, -0.498666668933333, -0.502666672066667, -0.505333344133333, -0.50200000405, -0.506666680166667, -0.508000006266667, -0.5040000081, -0.501333341, -0.50200000405, -0.502000009016667, -0.504000003133333, -0.50466667115, -0.505333344133333, -0.501333341, -0.503333340083333, -0.496666664883333, -0.5040000081, -0.49799999595, -0.389333337533333, -0.3640000075, -0.365333338566667, -0.3586666733, -0.3599999994, -0.355333338183333, -0.365333338566667, -0.360666667416667, -0.3640000075, -0.3546666652, -0.369999999783333, -0.3653333336, -0.36733333765, -0.396666670833333, -0.469333335766667, -0.4586666673, -0.4639999916, -0.457333331283333, -0.4646666596, -0.4680000047, -0.4687999964, -0.467333336683333, -0.460666661433333, -0.4653333326, -0.4720000029, -0.453333328166667, -0.466000000616667, -0.4646666646, -0.45933333535, -0.464666659633333, -0.4579999993, -0.464666669566667, -0.473999996983333, -0.417333334666667, -0.46400000155, -0.4599999934, -0.46866666775, -0.47200000288, -0.465333337583333, -0.4719999979, -0.462666660516667, -0.46240000126, -0.465333322666667, -0.4673333317, -0.459999993416667, -0.464666659616667, -0.461333329466667, -0.46599999565, -0.461999997483333, -0.461999992516667, -0.466666663683333, -0.457333336283333, -0.46600000065, -0.463999996566667, -0.46400000155, -0.488666658566667, -0.543333326766667, -0.489333341533333, -0.3600000143, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -0.3573333323, -0.464666662133333, -0.3360000037, -0.445333333066667, -0.49333334465, -0.444666661316667, -0.325333339483333, -0.60199998815, -0.33039999454, -0.394666674216667, -0.406666671233333, -0.4066666588, -0.52533333005, -0.552000006066667, -0.6973333458, -0.526000003, -0.565333336583333, -0.451428575185714, -0.3973333314, -0.591999994871429, -0.49120000004, -0.4313333432, -0.197333325966667, -0.545999998833333, -0.49600000976, -0.4595000035875, -0.45920000072, -0.453714281285714, -0.455333337166667, -0.461333334433333, -0.461999992533333, -0.457333331283333, -0.4606666664, -0.456000005183333, -0.4593333254, -0.4579999993, -0.45679999586, -0.455428570485714, -0.454666664183333, -0.453999991183333, -0.452666665083333, -0.464000001533333, -0.456666668216667, -0.45439999696, -0.45666666825, -0.454666664166667, -0.459999998383333, -0.45840000508, -0.459333330366667, -0.394666661816667, -0.367333332983333, -0.2823999971, -0.3446666648, -0.464666659633333, -0.443333327766667, -0.33920000494, -0.332666660333333, -0.363333340716667, -0.43599999546, -0.518000001716667, -0.15933333585, -0.2493333357, -0.43359999356, -0.462666669166667, -0.410000008833333, -0.357333337883333, -0.3615999937, -0.494666665783333, -0.52200000485, -0.29400000225, -0.4527999997, -0.474000004416667, -0.342666663233333, -0.53680000306, -0.423333331933333, -0.498666671416667, -0.439333344516667, -0.46480000914, -0.368666668733333, -0.2673333399, -0.303333333366667, -0.21279999912, -0.398000004383333, -0.444666676216667, -0.63679999114, -0.337333336466667, -0.2599999905 ),
date = c(1493139120, 1493139180, 1493139240, 1493139300, 1493139360, 1493139420, 1493139480, 1493139540, 1493139600, 1493139660, 1493139720, 1493139780, 1493139840, 1493139900, 1493139960, 1493140020, 1493140080, 1493140140, 1493140200, 1493140260, 1493140320, 1493140380, 1493140440, 1493140500, 1493140560, 1493140620, 1493140680, 1493140740, 1493140800, 1493140860, 1493140920, 1493140980, 1493141040, 1493141100, 1493141160, 1493141220, 1493141280, 1493141340, 1493141400, 1493141460, 1493141520, 1493141580, 1493141640, 1493141700, 1493141760, 1493141820, 1493141880, 1493141940, 1493142000, 1493142060, 1493142120, 1493142180, 1493142240, 1493142300, 1493142360, 1493142420, 1493142480, 1493142540, 1493142600, 1493142660, 1493142720, 1493142780, 1493142840, 1493142900, 1493142960, 1493143020, 1493143080, 1493143140, 1493143200, 1493143260, 1493143320, 1493143380, 1493143440, 1493143500, 1493143560, 1493143620, 1493143680, 1493143740, 1493143800, 1493143860, 1493143920, 1493143980, 1493144040, 1493144100, 1493144160, 1493144220, 1493144280, 1493144340, 1493144400, 1493144460, 1493144520, 1493144580, 1493144640, 1493144700, 1493144760, 1493144820, 1493144880, 1493144940, 1493145000, 1493145060, 1493145120, 1493145180, 1493145240, 1493145300, 1493145360, 1493145420, 1493145480, 1493145540, 1493145600, 1493145660, 1493145720, 1493145780, 1493145840, 1493145900, 1493145960, 1493146020, 1493146080, 1493146140, 1493146200, 1493146260, 1493146320, 1493146380, 1493146440, 1493146500, 1493146560, 1493146620, 1493146680, 1493146740, 1493146800, 1493146860, 1493146920, 1493146980, 1493147040, 1493147100, 1493147160, 1493147220, 1493147280, 1493147340, 1493147400, 1493147460, 1493147520, 1493147580, 1493147640, 1493147700, 1493147760, 1493147820, 1493147880, 1493147940, 1493148000, 1493148060, 1493148120, 1493148180, 1493148240, 1493148300, 1493148360, 1493148420, 1493148480, 1493148540, 1493148600, 1493148660, 1493148720, 1493148780, 1493148840, 1493148900, 1493148960, 1493149020, 1493149080, 1493149140, 1493149200, 1493149260, 1493149320, 1493149380, 1493149440, 1493149500, 1493149560, 1493149620, 1493149680, 1493149740, 1493149800, 1493149860, 1493149920, 1493149980, 1493150040, 1493150100, 1493150160, 1493150220, 1493150280, 1493150340, 1493150400, 1493150460, 1493150520))
)
fit.xts <- auto.arima(qxts$y, seasonal = TRUE)
forecast_length <- 10
fore.xts <- forecast(fit.xts, h=forecast_length)
return(fore.xts)
})
#Auto.Arima Seasonal data
output$arima.seasonalData <- renderTable({
if (input$seasonal == "Seasonal"){
seasonaldata()
}
})
}
shinyApp(ui, server)
I would like the dashboard to wait for me to hit the predict button before displaying the chart and the forecasted data in the boxes based on selected radio-button option.
input$seasonal is a reactive; when it changes value it triggers an update in any reactive or observers that uses it (which also includes renderXXX). Your actionbutton input$click is also a reactive, and follows the exact same rules. Its value is a simple integer counting how many times the button has been clicked in that session.
Use isolate(input$seasonal) if you want to use the value of the radio buttons, but not have changes trigger an update.
To then have a reactive, observer or renderXXX to be triggered by a click on the action button, simply include input$click somewhere, i.e.:
output$plot <- renderPlot({
if (isolate(input$seasonal) == 'foo') {
data = "foo"
} else {
data = "bar"
}
input$click # <--- this is where we react to the action button
plot(data)
})
I have two datasets below, ItemListFiltered and ItemListAll. ItemListFiltered is just a subset of ItemListAll. I'm having trouble confirming the results of the apriori function.
ItemListAll returns a rule...
{V2=C15OMB3CH, V3=C15KRLMRO} => {V1=C15MORL} with a confidence of 37%
ItemListFiltered returns a rule...
{V1=C15KRLMRO, V3=C15OMB3CH} => {V2=C15MORL} with a confidence of 80%
I was expecting the confidence calculation to be...
Confidence = Count of Orders from X&Y / Count of Orders from X
I assume X is the left hand side of (=>) and Y is the right hand side of (=>)
C15OMB3CH, C15KRLMRO, and C15MORL are on 4 orders together, where as C15OMB3CH, and C15KRLMRO are on 11 orders. 4/11 = 37%
How can these two return a difference in confidence?
R 3.5.1 code
library(Matrix)
library(arules)
library(plyr)
library(dplyr)
library(gridExtra)
dataset <- read.csv("ItemListAll.csv", header=TRUE, sep=",")
#dataset <- read.csv("ItemListFiltered.csv", header=TRUE, sep=",")
df_itemList <- ddply(dataset,c("SALESID"),function(df1)paste(df1$ITEMID));
df_itemList$V1 <- as.factor(df_itemList$V1);
df_itemList$V2 <- as.factor(df_itemList$V2);
df_itemList$V3 <- as.factor(df_itemList$V3);
basket_rules <- apriori(select(df_itemList,c("V1","V2","V3")), parameter = list(sup=0.001,conf=0.1,target="rules", minlen=3));
df_basket <- as(basket_rules,"data.frame");
df_basket$support <- ceiling(df_basket$support * 100);
df_basket$confidence<- ceiling(df_basket$confidence * 100);
df_basket$lift<- round(df_basket$lift, digits = 2);
df_basket <- df_basket[rev(order(df_basket$support)),];
View(df_basket);
ItemListFiltered.csv
SALESID ITEMID
SO0973259-1 C15ABPRRO
SO0972779-1 C15LGS5M
SO0971629-1 C15MORL
SO0972734-1 C15MORL
SO0973561-3 C15MORL
SO0973561-4 C15MORL
SO0974783-1 C15PECHRO
SO0974783-2 C15PECHRO
SO0973561-1 C15RPDERO
SO0973561-2 C15RPDERO
SO0974466-1 C15RPDERO
SO0973259-1 C15OMB3CH
SO0972779-1 C15OMB3CH
SO0971629-1 C15OMB3CH
SO0972734-1 C15OMB3CH
SO0973561-3 C15OMB3CH
SO0973561-4 C15OMB3CH
SO0974783-1 C15OMB3CH
SO0974783-2 C15OMB3CH
SO0973561-1 C15OMB3CH
SO0973561-2 C15OMB3CH
SO0974466-1 C15OMB3CH
SO0973259-1 C15KRLMRO
SO0972779-1 C15KRLMRO
SO0971629-1 C15KRLMRO
SO0972734-1 C15KRLMRO
SO0973561-3 C15KRLMRO
SO0973561-4 C15KRLMRO
SO0974783-1 C15KRLMRO
SO0974783-2 C15KRLMRO
SO0973561-1 C15KRLMRO
SO0973561-2 C15KRLMRO
SO0974466-1 C15KRLMRO
ItemListAll.csv
SALESID ITEMID
SO0971629-1 C15KRLMRO
SO0971629-1 C15MORL
SO0971629-1 C15OMB3CH
SO0971841-2 C15ARBSIM
SO0971841-2 C15LADCRO
SO0971841-2 C15MORL
SO0972404-1 C15ABPRRO
SO0972404-1 C15CAOBG
SO0972404-1 C15PV20
SO0972564-3 C15ABPRRO
SO0972564-3 C15CAINM550
SO0972564-3 C15DLMT
SO0972710-1 C15DLMT
SO0972710-1 C15LGS5M
SO0972710-1 C15OMB3CH
SO0972734-1 C15KRLMRO
SO0972734-1 C15MORL
SO0972734-1 C15OMB3CH
SO0972744-1 C15CAINM550
SO0972744-1 C15KRLMRO
SO0972744-1 C15LGS5M
SO0972779-1 C15KRLMRO
SO0972779-1 C15LGS5M
SO0972779-1 C15OMB3CH
SO0972882-1 C15CAOBG
SO0972882-1 C15LADCRO
SO0972882-1 C15RPDERO
SO0972914-2 C15KRLMRO
SO0972914-2 C15LADCRO
SO0972914-2 C15LGS5M
SO0972948-1 C15ARBSIM
SO0972948-1 C15EPCNWCB
SO0972948-1 C15LGS5M
SO0972981-1 C15CAOBG
SO0972981-1 C15LADCRO
SO0972981-1 C15PECHRO
SO0972992-1 C15CAOBG
SO0972992-1 C15OMB3CH
SO0972992-1 C15PECHRO
SO0972995-2 C15CAOBG
SO0972995-2 C15MORL
SO0972995-2 C15OMB3CH
SO0973002-1 C15KRLMRO
SO0973002-1 C15LADCRO
SO0973002-1 C15MORL
SO0973010-1 C15ABPRRO
SO0973010-1 C15CAOBG
SO0973010-1 C15LADCRO
SO0973065-3 C15ABPRRO
SO0973065-3 C15CAOBG
SO0973065-3 C15LADCRO
SO0973127-1 C15ARBSIM
SO0973127-1 C15LADCRO
SO0973127-1 C15PECHRO
SO0973195-1 C15CAOBG
SO0973195-1 C15LADCRO
SO0973195-1 C15LGS5M
SO0973207-1 C15ARBSIM
SO0973207-1 C15OMB3CH
SO0973207-1 C15RPDERO
SO0973208-1 C15ABPRRO
SO0973208-1 C15NUBCA
SO0973208-1 C15OMB3CH
SO0973209-1 C15LADCRO
SO0973209-1 C15NUBCA
SO0973209-1 C15RPDERO
SO0973223-1 C15CAOBG
SO0973223-1 C15LGS5M
SO0973223-1 C15OMB3CH
SO0973251-1 C15ARBSIM
SO0973251-1 C15PV20
SO0973251-1 C15RPDERO
SO0973259-1 C15ABPRRO
SO0973259-1 C15KRLMRO
SO0973259-1 C15OMB3CH
SO0973266-1 C15CAOBG
SO0973266-1 C15MORL
SO0973266-1 C15PV20
SO0973277-1 C15EPCNWCB
SO0973277-1 C15MORL
SO0973277-1 C15NUBCA
SO0973307-2 C15CAOBG
SO0973307-2 C15OMB3CH
SO0973307-2 C15RPDERO
SO0973311-1 C15ABPRRO
SO0973311-1 C15CAOBG
SO0973311-1 C15OMB3CH
SO0973414-1 C15CAOBG
SO0973414-1 C15MORL
SO0973414-1 C15OMB3CH
SO0973457-1 C15ABPRRO
SO0973457-1 C15KRLMRO
SO0973457-1 C15LADCRO
SO0973460-1 C15DLMT
SO0973460-1 C15MORL
SO0973460-1 C15OMB3CH
SO0973466-1 C15ABPRRO
SO0973466-1 C15KRLMRO
SO0973466-1 C15LADCRO
SO0973516-1 C15DLMT
SO0973516-1 C15MORL
SO0973516-1 C15OMB3CH
SO0973541-1 C15CAOBG
SO0973541-1 C15PV20
SO0973541-1 C15RPDERO
SO0973547-1 C15CAOBG
SO0973547-1 C15MORL
SO0973547-1 C15OMB3CH
SO0973555-1 C15KRLMRO
SO0973555-1 C15LADCRO
SO0973555-1 C15MORL
SO0973561-1 C15KRLMRO
SO0973561-1 C15OMB3CH
SO0973561-1 C15RPDERO
SO0973561-2 C15KRLMRO
SO0973561-2 C15OMB3CH
SO0973561-2 C15RPDERO
SO0973561-3 C15KRLMRO
SO0973561-3 C15MORL
SO0973561-3 C15OMB3CH
SO0973561-4 C15KRLMRO
SO0973561-4 C15MORL
SO0973561-4 C15OMB3CH
SO0973739-1 C15ABPRRO
SO0973739-1 C15CAOBG
SO0973739-1 C15PV20
SO0974027-2 C15CAINM550
SO0974027-2 C15CAOBG
SO0974027-2 C15LGS5M
SO0974027-3 C15CAINM550
SO0974027-3 C15CAOBG
SO0974027-3 C15LGS5M
SO0974093-1 C15CAOBG
SO0974093-1 C15LGS5M
SO0974093-1 C15OMB3CH
SO0974104-1 C15CAOBG
SO0974104-1 C15LADCRO
SO0974104-1 C15MORL
SO0974114-1 C15KRLMRO
SO0974114-1 C15LADCRO
SO0974114-1 C15LGS5M
SO0974233-3 C15NUBCA
SO0974233-3 C15OMB3CH
SO0974233-3 C15PECHRO
SO0974437-3 C15DLMT
SO0974437-3 C15LADCRO
SO0974437-3 C15LGS5M
SO0974466-1 C15KRLMRO
SO0974466-1 C15OMB3CH
SO0974466-1 C15RPDERO
SO0974471-1 C15DLMT
SO0974471-1 C15MORL
SO0974471-1 C15OMB3CH
SO0974538-8 C15DLMT
SO0974538-8 C15MORL
SO0974538-8 C15OMB3CH
SO0974539-1 C15CAOBG
SO0974539-1 C15EPCNWCB
SO0974539-1 C15PECHRO
SO0974558-1 C15ABPRRO
SO0974558-1 C15CAOBG
SO0974558-1 C15PV20
SO0974581-1 C15MORL
SO0974581-1 C15NUBCA
SO0974581-1 C15PV20
SO0974649-1 C15CAOBG
SO0974649-1 C15LADCRO
SO0974649-1 C15PECHRO
SO0974663-1 C15CAINM550
SO0974663-1 C15KRLMRO
SO0974663-1 C15RPDERO
SO0974702-1 C15CAOBG
SO0974702-1 C15PECHRO
SO0974702-1 C15PV20
SO0974715-1 C15DLMT
SO0974715-1 C15OMB3CH
SO0974715-1 C15RPDERO
SO0974753-1 C15DLMT
SO0974753-1 C15OMB3CH
SO0974753-1 C15RPDERO
SO0974783-1 C15KRLMRO
SO0974783-1 C15OMB3CH
SO0974783-1 C15PECHRO
SO0974783-2 C15KRLMRO
SO0974783-2 C15OMB3CH
SO0974783-2 C15PECHRO
SO0974923-2 C15DLMT
SO0974923-2 C15LGS5M
SO0974923-2 C15OMB3CH
Here is updated code I used as well with references to here and here
#Main
library(Matrix)
library(arules)
library(plyr)
library(dplyr)
library(gridExtra)
#itemList <- read.csv("ItemListAll.csv", header=TRUE, sep=",")
itemList <- read.csv("ItemListFiltered.csv", header=TRUE, sep=",")
# Converting to a Matrix ####
itemList$const = TRUE
# Remove duplicates
dim(itemList)
orders <- unique(itemList)
dim(itemList)
# Need to reshape the matrix
itemList_max_prep <- reshape(data = itemList,
idvar = "SALESID",
timevar = "ITEMID",
direction = "wide")
# Drop the SALESID
itemList_matrix <- as.matrix(itemList_max_prep[,-1])
# Clean up the missing values to be FALSE
itemList_matrix[is.na(itemList_matrix)] <- FALSE
# Clean up names
colnames(itemList_matrix) <- gsub(x=colnames(itemList_matrix),
pattern="const\\.", replacement="")
itemList_trans <- as(itemList_matrix,"transactions")
inspect(itemList_trans)
basket_rules <- apriori(itemList_trans, parameter = list(sup=0.01,conf=0.5,target="rules", minlen=3));
df_basket <- as(basket_rules,"data.frame")
df_basket$support <- ceiling(df_basket$support * 100)
df_basket$confidence<- ceiling(df_basket$confidence * 100)
df_basket$lift<- round(df_basket$lift, digits = 2)
df_basket <- df_basket[rev(order(df_basket$support)),];
write.xlsx(df_basket, "Basket_Analysis.xlsx", col.names=TRUE, row.names=TRUE)
#grid.table(head(df_basket));
I have this zoo series.
> str(ALLXbm)
‘zoo’ series from 1998-01-01 to 2017-12-01
Data: num [1:240, 1:19] 0.00289 0.00992 -0.00424 0.07958 -0.035 ...
- attr(*, "dimnames")=List of 2
..$ : NULL
..$ : chr [1:19] "DE" "DB" "DM" "IE" ...
Index: Date[1:240], format: "1998-01-01" "1998-02-01" "1998-03-01" "1998-04-01" "1998-05-01" "1998-06-01" "1998-07-01" ...
The followings are the dput() results. The dependent variable is
> dput(ALLX$DE)
structure(c(" 0.0028923077", " 0.0099160000", "-0.0042440000",
" 0.0795800000", "-0.0350041667", "-0.2480500000", " 0.0042880000",
"-0.1620888889", " 0.0637518519", " 0.0928178571", "-0.0974379310",
"-0.1012612903", "-0.0141133333", " 0.0373500000", " 0.0420323529",
" 0.0406323529", "-0.0479314286", " 0.0393514286", " 0.0134594595",
"-0.0042162162", " 0.0282142857", "-0.0467717949", " 0.0271487179",
" 0.0484756757", " 0.0388105263", "-0.0232324324", "-0.1534487179",
"-0.0553538462", "-0.0127976190", " 0.0248613636", "-0.0069093023",
" 0.0016234043", "-0.0004617021", "-0.0747617021", " 0.0362914894",
" 0.0773816327", "-0.0038936170", " 0.0876829787", "-0.0198040000",
"-0.0078372549", " 0.0343978723", " 0.0350833333", "-0.0382480000",
" 0.0824346939", "-0.0280615385", "-0.0312228070", "-0.1020928571",
"-0.1011178571", " 0.0460186441", " 0.0407403509", "-0.0200607143",
" 0.1692310345", " 0.1254927273", "-0.0012947368", " 0.0410393443",
"-0.0060650000", " 0.1175457627", "-0.0219952381", " 0.0538709677",
" 0.1302564516", " 0.0177306452", " 0.0465440678", "-0.0318777778",
"-0.0027257576", "-0.0808714286", " 0.0753223881", " 0.0421338235",
" 0.0192617647", " 0.0645835821", " 0.0178743243", " 0.0917106667",
"-0.0504720000", "-0.0605917808", " 0.0670115942", " 0.0841644737",
"-0.0862937500", " 0.0583458333", " 0.0341283951", " 0.0363376623",
"-0.0072297619", " 0.0697188235", " 0.0708900000", " 0.0964263736",
" 0.0265102273", "-0.0310302326", " 0.0419810000", "-0.0719272727",
" 0.0243521739", "-0.0945943925", " 0.0339018868", " 0.0440350000",
" 0.0503681818", " 0.0629381356", "-0.0510975000", " 0.0249760684",
" 0.0827214876", " 0.0892224806", "-0.0168313043", " 0.0291585366",
" 0.0273647541", "-0.0975674242", "-0.0746062016", " 0.0238595238",
"-0.0278552000", "-0.0925223077", " 0.0659656934", " 0.0456731343",
" 0.0466132450", "-0.0135568493", " 0.0414873016", " 0.0340628205",
" 0.0267903846", "-0.0494906667", " 0.0180333333", " 0.0387062500",
"-0.0283124183", " 0.0306898810", " 0.0886619048", "-0.0198162921",
"-0.0154431818", "-0.0838285714", " 0.0658326531", "-0.0614076503",
" 0.0424984211", "-0.0170115000", " 0.0205202970", " 0.0596765258",
"-0.0538638298", "-0.0453000000", "-0.0493330317", " 0.0545885714",
" 0.1193183168", "-0.0335423810", " 0.0599138393", " 0.0207066964",
" 0.0554773756", " 0.0606587156", " 0.0448098039", "-0.0354082569",
" 0.0057316038", " 0.0019272727", " 0.0214976303", "-0.0108309179",
" 0.0166442211", " 0.0028019512", "-0.0548807910", " 0.0297556962",
"-0.0161016854", " 0.0100739726", " 0.0538203704", " 0.0304168724",
" 0.0246617647", " 0.0294436975", "-0.0192055300", "-0.0221397906",
" 0.0429194215", "-0.1356025974", " 0.0155836364", " 0.0421339806",
" 0.0193581152", "-0.0316510526", " 0.0171187192", " 0.0299633484",
" 0.0136308017", "-0.0803899225", "-0.0462198413", "-0.0337478469",
"-0.0044057018", " 0.0297085185", " 0.0273300000", "-0.0814074689",
" 0.0194720183", "-0.0640064885", " 0.0194782979", " 0.0047418868",
"-0.0393941634", " 0.0119372093", "-0.0092268116", "-0.0221629956",
" 0.0803489130", "-0.0976939502", "-0.0391824701", "-0.0672646154",
"-0.0124784247", "-0.1038085627", "-0.0105592476", "-0.0279538710",
" 0.0196077181", " 0.0338205788", "-0.0189602996", "-0.0860247273",
"-0.0478760714", "-0.0510344720", " 0.0275029586", " 0.0428622291",
" 0.0166286232", "-0.0288036101", "-0.0180476923", " 0.0393909483",
"-0.0641453901", "-0.0841982143", " 0.0047675900", "-0.0240106312",
"-0.0700313793", " 0.0733418994", "-0.0533325905", "-0.0331509589",
" 0.0489809399", "-0.0736352078", " 0.0240526646", "-0.0939243553",
"-0.0042405995", "-0.0387828042", "-0.0470409756", "-0.0931115294",
"-0.0910497738", " 0.0013111111", " 0.0348456633", " 0.0705459082",
" 0.0758425000", "-0.1478202222", " 0.0717040179", " 0.0569382114",
"-0.0945091623", " 0.1044561713", " 0.0293398396", "-0.1272769802",
" 0.0325193750", " 0.0277184066", "-0.0341782369", "-0.0250598160",
" 0.0189187500", "-0.0074535385", "-0.0359813953", " 0.0160150115",
" 0.0476877551", "-0.1018837398", "-0.0427113689", " 0.0191961637",
" 0.1059326403"), index = structure(c(883612800, 886291200, 888710400,
891388800, 893980800, 896659200, 899251200, 901929600, 904608000,
907200000, 909878400, 912470400, 915148800, 917827200, 920246400,
922924800, 925516800, 928195200, 930787200, 933465600, 936144000,
938736000, 941414400, 944006400, 946684800, 949363200, 951868800,
954547200, 957139200, 959817600, 962409600, 965088000, 967766400,
970358400, 973036800, 975628800, 978307200, 980985600, 983404800,
986083200, 988675200, 991353600, 993945600, 996624000, 999302400,
1001894400, 1004572800, 1007164800, 1009843200, 1012521600, 1014940800,
1017619200, 1020211200, 1022889600, 1025481600, 1028160000, 1030838400,
1033430400, 1036108800, 1038700800, 1041379200, 1044057600, 1046476800,
1049155200, 1051747200, 1054425600, 1057017600, 1059696000, 1062374400,
1064966400, 1067644800, 1070236800, 1072915200, 1075593600, 1078099200,
1080777600, 1083369600, 1086048000, 1088640000, 1091318400, 1093996800,
1096588800, 1099267200, 1101859200, 1104537600, 1107216000, 1109635200,
1112313600, 1114905600, 1117584000, 1120176000, 1122854400, 1125532800,
1128124800, 1130803200, 1133395200, 1136073600, 1138752000, 1141171200,
1143849600, 1146441600, 1149120000, 1151712000, 1154390400, 1157068800,
1159660800, 1162339200, 1164931200, 1167609600, 1170288000, 1172707200,
1175385600, 1177977600, 1180656000, 1183248000, 1185926400, 1188604800,
1191196800, 1193875200, 1196467200, 1199145600, 1201824000, 1204329600,
1207008000, 1209600000, 1212278400, 1214870400, 1217548800, 1220227200,
1222819200, 1225497600, 1228089600, 1230768000, 1233446400, 1235865600,
1238544000, 1241136000, 1243814400, 1246406400, 1249084800, 1251763200,
1254355200, 1257033600, 1259625600, 1262304000, 1264982400, 1267401600,
1270080000, 1272672000, 1275350400, 1277942400, 1280620800, 1283299200,
1285891200, 1288569600, 1291161600, 1293840000, 1296518400, 1298937600,
1301616000, 1304208000, 1306886400, 1309478400, 1312156800, 1314835200,
1317427200, 1320105600, 1322697600, 1325376000, 1328054400, 1330560000,
1333238400, 1335830400, 1338508800, 1341100800, 1343779200, 1346457600,
1349049600, 1351728000, 1354320000, 1356998400, 1359676800, 1362096000,
1364774400, 1367366400, 1370044800, 1372636800, 1375315200, 1377993600,
1380585600, 1383264000, 1385856000, 1388534400, 1391212800, 1393632000,
1396310400, 1398902400, 1401580800, 1404172800, 1406851200, 1409529600,
1412121600, 1414800000, 1417392000, 1420070400, 1422748800, 1425168000,
1427846400, 1430438400, 1433116800, 1435708800, 1438387200, 1441065600,
1443657600, 1446336000, 1448928000, 1451606400, 1454284800, 1456790400,
1459468800, 1462060800, 1464739200, 1467331200, 1470009600, 1472688000,
1475280000, 1477958400, 1480550400, 1483228800, 1485907200, 1488326400,
1491004800, 1493596800, 1496275200, 1498867200, 1501545600, 1504224000,
1506816000, 1509494400, 1512086400), tzone = "UTC", tclass = "Date"), class = c("xts",
"zoo"), .indexCLASS = "Date", tclass = "Date", .indexTZ = "UTC", tzone = "UTC", .Dim = c(240L,
1L), .Dimnames = list(NULL, "DE"))
and the independent variable dput() is
> dput(ALLX$SMB_L)
structure(c("-0.0465320170", "-0.0697526439", " 0.0273518415",
"-0.0740531793", " 0.0316742119", " 0.0152537452", "-0.0264616225",
" 0.0288440266", "-0.1158052743", "-0.1629207669", "-0.0321036909",
"-0.0155047450", "-0.0826901559", "-0.0958867805", "-0.1071221681",
"-0.0448685814", "-0.0246406184", "-0.1212194655", "-0.0753136163",
" 0.0410887618", "-0.0472807197", "-0.0821359482", "-0.0538328997",
"-0.1193170755", " 0.0451348940", "-0.0666430865", " 0.1030261973",
"-0.0512679603", "-0.0858501881", "-0.1100767741", " 0.2831559983",
" 0.0049661923", "-0.0466122487", "-0.0566777104", "-0.0268275721",
"-0.0592022576", " 0.0095647556", "-0.0369526728", "-0.0246719322",
" 0.0051186559", "-0.0462739541", "-0.1027257463", " 0.0041844398",
"-0.0443327227", " 0.0513661492", "-0.0672995604", "-0.0509956017",
"-0.0377391126", " 0.1208954628", " 0.0656631973", "-0.0109091553",
"-0.0319298318", " 0.0762361370", " 0.0082987739", " 0.0744911252",
" 0.0268991873", " 0.0106712270", "-0.0162705327", "-0.0030339799",
" 0.0437060161", "-0.0484248103", "-0.0131028206", "-0.0551732558",
"-0.0827144280", " 0.0412271313", "-0.0625801905", "-0.0531687861",
"-0.0194387563", " 0.0651252536", " 0.0185203178", " 0.0212495415",
"-0.0058409877", " 0.0005737324", " 0.0033661115", " 0.0166057656",
" 0.0273215778", "-0.0160116096", "-0.0516442850", " 0.0622634418",
" 0.0059385646", "-0.1501443652", " 0.0901923236", " 0.0033375512",
"-0.0753535560", " 0.0112127728", " 0.0174845326", " 0.0523481171",
"-0.0687008257", " 0.0364307873", " 0.0136230636", " 0.0058142113",
" 0.0554519027", " 0.0519022715", "-0.0570389442", "-0.0170328466",
" 0.0497557753", "-0.0022714909", " 0.0319902902", " 0.0567367912",
"-0.0807706219", " 0.0858961234", " 0.0222949488", " 0.0242892126",
" 0.0215915873", " 0.0537398523", "-0.0632469276", " 0.0111250970",
" 0.0532290206", "-0.0146115567", "-0.0497947429", "-0.0261614454",
"-0.0231347005", " 0.1003858901", " 0.0445487378", "-0.0150134676",
"-0.0457789040", "-0.0034933110", "-0.0668239816", "-0.0159024243",
"-0.0077140317", " 3.6823096028", "-0.0901045024", "-0.0143629036",
"-0.0098836289", "-0.0311347099", "-0.0151321886", "-0.0233462522",
"-0.0458991158", "-0.0282613742", "-0.0779876682", " 0.0117855889",
"-0.0893392649", "-0.0328418515", "-0.0090601777", "-0.0322611560",
"-0.0641631336", "-0.0595636188", " 0.0206379342", "-0.0891543562",
" 0.0717742257", "-0.0461537130", "-0.0954444539", "-0.0543154861",
"-0.0714790350", " 0.0474770475", "-0.0037168520", "-0.0077110744",
"-0.0355079545", " 0.0354765095", "-0.0152485350", "-0.0037694163",
" 0.0088693973", "-0.0013576659", "-0.0158975649", "-0.0569949742",
"-0.0047187899", " 0.0359832832", " 0.0048156392", "-0.0213823507",
"-0.0253481502", "-0.0305909946", "-0.0459944589", " 0.0697664274",
"-0.0053565569", " 0.0119311291", "-0.0068067904", "-0.0046527312",
"-0.0317481424", " 0.0322190631", " 0.0095508469", " 0.0009492415",
" 0.0109837401", " 0.0129883296", "-0.0230918819", " 0.0291994624",
" 0.0053826778", " 0.0116937939", " 0.7079728243", " 0.0103552944",
" 0.0016701093", " 0.0283116307", " 0.0367686262", " 0.0352952954",
" 0.0251452293", " 0.6522285618", " 0.0005676648", "-0.0355068478",
" 0.0383934967", "-0.0054156326", "-0.0288061771", " 0.0076047198",
"-0.0261335344", " 0.0587687305", " 0.0635220884", "-0.0035395762",
" 0.0638110053", "-0.0043276134", "-0.0183498085", "-0.0306437649",
"-0.0153833213", " 0.0165571651", "-0.0014514427", "-0.0111594992",
" 0.0061678879", "-0.0898811021", "-0.0106517096", " 0.0582380903",
" 0.0471135599", " 0.0946734825", " 0.0518059714", " 0.0075574758",
" 0.0078389969", " 0.0998944849", "-0.0259425031", "-0.0031939886",
"-0.0078991041", "-0.0583307843", "-0.0664592297", "-0.0428550515",
"-0.0282970576", " 0.0085063220", " 0.0027516516", "-0.0413816309",
"-0.0080472825", " 0.0001809490", "-0.0719654720", " 0.0250187336",
" 0.0017591069", "-0.0352040828", " 0.0032731177", "-0.0182234806",
"-0.0276918415", " 0.0246782744", " 0.0390556980", "-0.0240913343",
"-0.0207629243", " 0.0289632858", "-0.0398579730", " 0.0048926744",
"-0.0397687581"), index = structure(c(883612800, 886291200, 888710400,
891388800, 893980800, 896659200, 899251200, 901929600, 904608000,
907200000, 909878400, 912470400, 915148800, 917827200, 920246400,
922924800, 925516800, 928195200, 930787200, 933465600, 936144000,
938736000, 941414400, 944006400, 946684800, 949363200, 951868800,
954547200, 957139200, 959817600, 962409600, 965088000, 967766400,
970358400, 973036800, 975628800, 978307200, 980985600, 983404800,
986083200, 988675200, 991353600, 993945600, 996624000, 999302400,
1001894400, 1004572800, 1007164800, 1009843200, 1012521600, 1014940800,
1017619200, 1020211200, 1022889600, 1025481600, 1028160000, 1030838400,
1033430400, 1036108800, 1038700800, 1041379200, 1044057600, 1046476800,
1049155200, 1051747200, 1054425600, 1057017600, 1059696000, 1062374400,
1064966400, 1067644800, 1070236800, 1072915200, 1075593600, 1078099200,
1080777600, 1083369600, 1086048000, 1088640000, 1091318400, 1093996800,
1096588800, 1099267200, 1101859200, 1104537600, 1107216000, 1109635200,
1112313600, 1114905600, 1117584000, 1120176000, 1122854400, 1125532800,
1128124800, 1130803200, 1133395200, 1136073600, 1138752000, 1141171200,
1143849600, 1146441600, 1149120000, 1151712000, 1154390400, 1157068800,
1159660800, 1162339200, 1164931200, 1167609600, 1170288000, 1172707200,
1175385600, 1177977600, 1180656000, 1183248000, 1185926400, 1188604800,
1191196800, 1193875200, 1196467200, 1199145600, 1201824000, 1204329600,
1207008000, 1209600000, 1212278400, 1214870400, 1217548800, 1220227200,
1222819200, 1225497600, 1228089600, 1230768000, 1233446400, 1235865600,
1238544000, 1241136000, 1243814400, 1246406400, 1249084800, 1251763200,
1254355200, 1257033600, 1259625600, 1262304000, 1264982400, 1267401600,
1270080000, 1272672000, 1275350400, 1277942400, 1280620800, 1283299200,
1285891200, 1288569600, 1291161600, 1293840000, 1296518400, 1298937600,
1301616000, 1304208000, 1306886400, 1309478400, 1312156800, 1314835200,
1317427200, 1320105600, 1322697600, 1325376000, 1328054400, 1330560000,
1333238400, 1335830400, 1338508800, 1341100800, 1343779200, 1346457600,
1349049600, 1351728000, 1354320000, 1356998400, 1359676800, 1362096000,
1364774400, 1367366400, 1370044800, 1372636800, 1375315200, 1377993600,
1380585600, 1383264000, 1385856000, 1388534400, 1391212800, 1393632000,
1396310400, 1398902400, 1401580800, 1404172800, 1406851200, 1409529600,
1412121600, 1414800000, 1417392000, 1420070400, 1422748800, 1425168000,
1427846400, 1430438400, 1433116800, 1435708800, 1438387200, 1441065600,
1443657600, 1446336000, 1448928000, 1451606400, 1454284800, 1456790400,
1459468800, 1462060800, 1464739200, 1467331200, 1470009600, 1472688000,
1475280000, 1477958400, 1480550400, 1483228800, 1485907200, 1488326400,
1491004800, 1493596800, 1496275200, 1498867200, 1501545600, 1504224000,
1506816000, 1509494400, 1512086400), tzone = "UTC", tclass = "Date"), class = c("xts",
"zoo"), .indexCLASS = "Date", tclass = "Date", .indexTZ = "UTC", tzone = "UTC", .Dim = c(240L,
1L), .Dimnames = list(NULL, "SMB_L"))
I wanted a gmm rolling regression and i set up the code below.
dogmm <- function(x) {gmm(ALLX$DE~ ALLX$SMB_L, x = ALLXbm)}
kk=rollapply(dogmm, ALLXbm,width=24, by.column = FALSE)
Error message is
Error in zoo(rval, index(x)[i]) :
“x” : attempt to define invalid zoo object
I really can't figure out even if i spent much time. when i tried the gmm code alone, it does work
gmm(ALLXbm$DE~ ALLXbm$SMB_L, x = ALLXbm). I'm thinking that, the problem is in the rollapply code. I'm sorry this's my first attempt to GMM. Please some help?
I have a simple problem of trying to produce an HTML file with rmarkdown. I want to include a table which has one column which has a long list of genes.
All I want to do is set the column width of each column, both columns will require wrapping.
I have gone with pander and it seems to be working when called in R, but when I create the HTML only the first row is included in the table.
pander(df, split.cells = c(30,70), split.table = Inf)
Below is the HTML output any ideas what is happening?
dput(df)
structure(c("Cell Cycle, Mitotic", "Cell Cycle", "Adaptive Immune System",
"Mitotic G1-G1/S phases", "Cytokine Signaling in Immune system",
"G1/S Transition", "DNA Replication", "G2/M Checkpoints", "Activation of ATR in response to replication stress",
"S Phase", "CDKN2A, SKP2, CCNE1, NCAPD2, LPIN1, MCM10, TFDP1, NUP155, MCM3, E2F3, OPTN, RFC4, NUP160, CENPN, TFDP2, CCNA1, NUP153, PSMB10, PRKCA, CKAP5, KIF2C, PSME4, GINS4, FOXM1, SEC13, RB1, MCM4, MCM5, CCNA2, PSMB9, RFC3, NCAPG2, CDC7, MCM7, NUP205, CENPA, PRKCB, BIRC5, SMC2, BUB1, CDCA8, PRIM2, NDC80, CDKN2C, PLK1, PPP1CB, CETN2, CEP57, NUP50, PSME1, CCNE2, PSMB8, NCAPH2, VRK1, CDC20, PCM1, MCM2, E2F5, TUBA4A, PSMB2, CDT1, SMC4, ORC1, CUL1, AURKB, PSMD2, CDC25B, KIF18A, UBE2D1, RFC2, TUBG2, SMC1A, PSMB4, NCAPD3, BTRC, CDK6, DNA2, NUP93, POLD1, SPDL1, RBL1, CENPI, PCNA, MCM6, NUP62, POLE2, SKA1, CLASP2, ARPP19, CKS1B, AKAP9, PSMC2, WAPL, BORA, DBF4, CENPF, PLK4, NEK2, SKP1, ORC6, NUMA1, PAFAH1B1, TYMS, KIF23, ZWILCH, FBXO5, FEN1, PSMF1, RRM2",
"CDKN2A, SFN, SKP2, CCNE1, NCAPD2, LPIN1, MCM10, TFDP1, NUP155, CHEK1, MCM3, E2F3, OPTN, RFC4, NUP160, CENPN, TFDP2, CCNA1, STAG3, NUP153, PSMB10, PRKCA, CKAP5, KIF2C, PSME4, GINS4, FOXM1, SEC13, RB1, MCM4, MCM5, CCNA2, PSMB9, RFC3, NCAPG2, CDC7, MCM7, NUP205, CENPA, HSPA2, PRKCB, BIRC5, SMC2, BUB1, CDCA8, PRIM2, NDC80, CDKN2C, PLK1, PPP1CB, CETN2, CEP57, NUP50, PSME1, CCNE2, PSMB8, MLH1, NCAPH2, VRK1, CDC20, PCM1, MCM2, E2F5, TUBA4A, MRE11A, PSMB2, CDT1, SMC4, ORC1, CUL1, AURKB, PSMD2, CHEK2, CDC25B, KIF18A, UBE2D1, RFC2, TUBG2, SMC1A, HJURP, PSMB4, NCAPD3, HUS1, BTRC, BLM, YWHAB, CDK6, DNA2, NUP93, POLD1, SPDL1, RBL1, CENPI, YWHAQ, PCNA, MCM6, RAD1, NUP62, BRCA2, POLE2, SKA1, CLASP2, ARPP19, CKS1B, AKAP9, PSMC2, WAPL, SUN2, BORA, DBF4, CENPF, PLK4, NEK2, SKP1, ORC6, NUMA1, PAFAH1B1, TYMS, KIF23, ZWILCH, FBXO5, FEN1, MIS18A, PSMF1, RRM2",
"CTSV, HLA-DOB, PLCG2, NCK1, DYNC1I1, ERBB4, TAP1, LYN, TAP2, ICAM1, RAC1, CYBB, ICOS, CD19, NFKBIE, PSMB10, SOS1, YES1, SELL, KIF2C, AGO1, ITPR1, ITK, PSME4, CD247, SEC13, INPP5D, CD79A, PIK3CD, RIPK2, HLA-F, PSMB9, CD8B, CD96, CTSS, LCK, NCF4, FYB, RASGRP2, PRKCB, ZAP70, PTPRC, TRAT1, LILRB2, PVR, ITPR3, CD8A, CD40, PRKCQ, CD3D, RASGRP3, LCP2, KLRD1, ITGB5, CRTAM, MALT1, CD79B, PSME1, BTK, PSMB8, VAV1, PDGFRA, PSMB2, PAK3, NRAS, CUL1, HLA-DMB, AP1M2, PSMD2, KIF18A, SYK, VCAM1, CD28, NCF2, CD86, CD3E, CD3G, EGFR, FGFR3, PSMB4, LILRB1, BTRC, YWHAB, CD4, KIF3C, PAK2, SEC24B, IKBKB, NPEPPS, RAPGEF3, ICAM2, ICOSLG, ITGB7, DCTN4, PSMC2, KIF3A, RAPGEF4, CSK, CTLA4, FYN, SKP1, ERBB3, KIF23, NFKBIB, ITGB2, FCGR1B, PSMF1",
"CDKN2A, SKP2, CCNE1, MCM10, TFDP1, MCM3, E2F3, TFDP2, CCNA1, PSMB10, RB1, MCM4, MCM5, PSMB9, CDC7, MCM7, PRIM2, CDKN2C, PSME1, CCNE2, PSMB8, MCM2, E2F5, PSMB2, CDT1, ORC1, CUL1, PSMD2, PSMB4, CDK6, RBL1, PCNA, MCM6, POLE2, CKS1B, PSMC2, DBF4, SKP1, ORC6, TYMS, FBXO5, PSMF1, RRM2",
"IRAK1, IL1R2, IL6ST, NUP155, IFNG, IRF1, CSF2RB, LYN, ICAM1, IRF8, MAP2K4, NUP160, ADAM17, SOCS2, IRF4, SOCS1, NUP153, SOS1, IL2RG, PTPN2, INPP5D, PIK3CD, IL1RAP, RIPK2, GBP1, HLA-F, TRIM38, PML, NUP205, IFNGR1, CIITA, MYD88, HERC5, MID1, ISG20, IFNAR2, IL7, NUP50, PSMB8, CASP1, CSF2RA, VAV1, CAMK2B, GBP2, KPNA4, NRAS, STAT1, CUL1, SYK, PRLR, IL2RB, VCAM1, IL7R, USP18, CISH, MX2, IL2RA, IL6R, BTRC, YWHAB, LGALS9, NUP93, PTAFR, KPNA3, MAPK1, IKBKB, RSAD2, NUP62, OASL, EIF4E2, SAMHD1, TRIM8, KPNA5, FYN, SKP1, FCGR1B",
"SKP2, CCNE1, MCM10, TFDP1, MCM3, CCNA1, PSMB10, RB1, MCM4, MCM5, PSMB9, CDC7, MCM7, PRIM2, PSME1, CCNE2, PSMB8, MCM2, PSMB2, CDT1, ORC1, CUL1, PSMD2, PSMB4, PCNA, MCM6, POLE2, CKS1B, PSMC2, DBF4, SKP1, ORC6, TYMS, FBXO5, PSMF1, RRM2",
"MCM10, MCM3, RFC4, PSMB10, PSME4, GINS4, RB1, MCM4, MCM5, PSMB9, RFC3, CDC7, MCM7, PRIM2, PSME1, PSMB8, MCM2, PSMB2, CDT1, ORC1, PSMD2, RFC2, PSMB4, DNA2, POLD1, PCNA, MCM6, POLE2, PSMC2, DBF4, ORC6, FEN1, PSMF1",
"SFN, MCM10, CHEK1, MCM3, RFC4, MCM4, MCM5, RFC3, CDC7, MCM7, MCM2, ORC1, CHEK2, RFC2, HUS1, YWHAB, YWHAQ, MCM6, RAD1, DBF4, ORC6",
"MCM10, CHEK1, MCM3, RFC4, MCM4, MCM5, RFC3, CDC7, MCM7, MCM2, ORC1, RFC2, HUS1, MCM6, RAD1, DBF4, ORC6",
"SKP2, MCM3, RFC4, PSMB10, PSME4, GINS4, RB1, MCM4, MCM5, PSMB9, RFC3, MCM7, PRIM2, PSME1, PSMB8, MCM2, PSMB2, CDT1, ORC1, CUL1, PSMD2, RFC2, SMC1A, PSMB4, DNA2, POLD1, PCNA, MCM6, POLE2, CKS1B, PSMC2, WAPL, SKP1, ORC6, FEN1, PSMF1"
), .Dim = c(10L, 2L), .Dimnames = list(c("1", "2", "3", "4",
"5", "6", "7", "8", "9", "10"), c("Pathway", "GeneID")))
view of R console showing splitting working.