-
Notifications
You must be signed in to change notification settings - Fork 1
/
.Rhistory
102 lines (102 loc) · 5.68 KB
/
.Rhistory
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
library(hawkes)
help(gof)
devtools::document(hawkes_em_functions.R)
devtools::document(R/hawkes_em_functions.R)
devtools::document("R/hawkes_em_functions.R")
devtools::document("R")
library(hawkes)
gof
help(gof)
help(compensator)
help(pp)
library(hawkes)
help(pp)
library(hawkes)
help(pp)
stats::ts
devtools::document("R")
library(hawkes)
help(pp)
help(ts)
help(pp)
help(plot.pp)
help(print.pp)
help(plot.pp)
plot.pp <-function(pp_obj) {
m <- attr(pp_obj, "n_margins")
n <- attr(pp_obj, "n_events")
plot(pp_obj[[1]], 1:n[1], type = "n", main = "Number of Events by Time",
ylab = "Events", xlab = "Time", ylim = c(0, max(n)),
xlim = c(0, max(unlist(pp_obj))), las = 1)
axis(side = 4, las = 1)
par(mar = c(5,4,4,4) + 1)
#mtext("Events", side = 4, line = 3, cex = par("cex.lab"))
legend(x = "topleft", legend = attr(pp_obj, "names"),
col = c(3:(m+2)), lty = 1, lwd = 2, border = "n")
for (i in 1:m) {
points(pp_obj[[i]], 1:n[i], type = "l", col = i+2, lwd = 2)
}
# Under construction: Event times
# plot(pp[[1]], rep(1, length(pp[[1]])), type = "n", main = "Event Times", ylab = "", xlab = "Time", pch = 124, yaxt = "n", ylim = c(-.5, -(m+.5)))
#
# for (i in 1:m) {
# points(pp[[i]], rep(-i, length(pp[[i]])), pch = 124, col = i+2)
# }
}
devtools::document("R")
library(hawkes)
help(res)
pp_obj <- pp(n_margins = 2, end_time = 100)
pp_obj
attr(pp_obj, "n_margins")
""
pp_obj$n_margins
help(gof)
devtools::document("R")
library(hawkes)
devtools::load_all(".")
help(set_kernel)
help(set_parms)
devtools::document("R")
library(hawkes)
help(set_parms)
help(get_parms)
help(compute_I)
devtools::document("R")
library(hawkes)
help(compute_I)
help(intensity)
help(probs)
devtools::document("R")
library(hawkes)
devtools::document("R")
library(hawkes)
help(compute_Q)
help(Q-fun)
help(Q_fun)
devtools::document("R")
library(hawkes)
help(Q_fun)
help(compute_mu)
help(compute_alpha)
help(logL)
help(EM)
devtools::document("R")
library(hawkes)
tau1 = c(63.217,105.783,282.633,345.717,345.783,403.733,418.517,517.933,522.85,550.417,596.583,618.067,686.75,689.233,876.45,884.85,926.033,930.867,951.183,953.017,1050.983,1074.95,1074.967,1091.467,1191.7,1215.55,1215.583,1235.333,1239.2,1258.7,1258.933,1379.583,1380.35,1428.417,1428.483,1450.283,1523.433,1595.7,1601.567,1720.55,1740.767,1746.3,1764.117,1859.017,1884.583,1887.55,1930.95,2030.883,2098.183,2221.383,2224.95,2339.6,2459.533,2462.85,2557.033,2560.967,2584.383,2584.833,2628.45,2734.95,2752.35,2797.267,2891.167,3038.733,3083.15,3136,3225.933,3274.2,3303.433,3399.217,3421.317,3422.1,3449.183,3536.95,3763.467,3780.05,3780.167,3972.583,3978.4,4028.717,4051.383,4094.8,4137.95,4145.317,4218.217,4262.017,4262.4,4289.3,4713.2,4728.15,4742.217,4744.9,4770.983,4794.417,4932.8,4935.75,4964.383,5055.2,5057.283,5079.7,5206.367,5218.717,5225.95,5233.933,5536.983,5554.65,5586.45,6060.35,6085.55,6116.05,6119.85,6205.567,6233.067,6258.95,6594.567,6609.4,6621.483,6642.7,6647.183,6719.367,6733.367,6759.433,6763.1,6788.033,7069.65,7070.067,7217.233,7240.6,7269.867,7404.517,7515.05,7516.3,7604.783,7608.767,7619.883,7921.417,8104.65,8245.65,8245.667,8246.05,8248.15,8276.65,8296.85,8317.467,8320.4,8465.85,8485.667,8507.167,8583.817,8628.3,8652.1,8677.55,8680.25,8748.733,8781.383,8938.95,8939.367,9087.433,9110.45,9113.283,9114.3,9116.517,9131.817,9251.4,9257.25,9281.1,9420.333,9447.467,9494.083,9518.35,9521.583,9615.717,9639.317,9644.2,9781,9804.2,9804.95,9806.35,9852.517,9981.467,9985.383,9997.667,10089.167,10093.467,10095.583,10117.4,10120.483,10160.75,10163.35,10334.45,10357.967,10451.683,10453.817,10530.4,10769.317,10933.633,10941.417,11004.817,11127.25,11151.467,11152.583,11156.283,11174.433,11267.017,11271.067,11277.167,11317.667,11424.917,11442.283,11442.633,11488.183,11603.133,11654.233,11676.2,11850.017,11872.35,11878,11926.133,12089.25,12107.283,12264.967,12275.6,12307.35,12349.617,12351.183,12352.717)
tau2 = c(16.817,60.067,66.9,86.083,183.45,426.3,525.133,543.467,566.067,585.683,593.033,601.083,611.067,671.283,688.5,695.767,878,928.067,930.617,952.883,1000.2,1000.5,1030.3,1070,1073.45,1077.067,1121.767,1175.067,1175.083,1198.233,1198.683,1199.983,1221.083,1238.733,1256.967,1379.033,1437.6,1511.1,1523.017,1530.633,1583.517,1725.517,1741.2,1743.8,1745.2,1763.467,1769.517,1840.8,1841.8,1868.417,1883.867,1884.367,1884.6,1935.167,2045.417,2196.333,2205.733,2223.333,2345.233,2422.8,2465.767,2555.55,2632.15,2709.15,2709.267,2729.6,2735.833,2775.433,2782.05,2801.517,3062.85,3113.933,3143.3,3227.233,3274.483,3303.467,3399.25,3422.017,3448.067,3468.133,3763.933,3780.083,3978.517,4007.817,4048.783,4049.283,4137.967,4222.317,4262.217,4262.4,4290.067,4727.683,4739.633,4774.117,4798.05,4798.167,4932.017,4956.833,4985.867,5040.85,5060.25,5174.467,5197.85,5218.55,5219.283,5231.15,5241.333,5400.45,5489.833,5542.633,6065.583,6113.667,6167.6,6238.9,6590.233,6596.4,6622.267,6644.283,6720.233,6761.233,6792.117,7069.867,7080.267,7185.55,7218.617,7279.3,7420.733,7516.017,7516.683,7608.317,7608.35,7618.867,7620.35,7648.817,7930.683,8246,8265.733,8297.033,8302.317,8481.533,8493.867,8624.667,8626.817,8627.433,8630.867,8637.2,8652.733,8674.05,8674.067,8708.733,8710.117,8711.333,8727.15,8749.55,8781.833,8940.383,8948.133,9092.467,9113.25,9113.667,9114.733,9251.3,9251.817,9254.6,9281.133,9429.233,9474.733,9616.683,9640.233,9786.367,9804.767,9813.367,9856.033,9982.233,9999.733,10006.417,10154.5,10324.85,10334.667,10499.367,10532.083,10934.4,11175.517,11271.033,11301.233,11442.317,11538.55,11660.633,11872.367,11922.567,12065.817,12090.95,12273.667,12322.833,12350.45)
D = list(tau1, tau2)
names(D) = c("Employee", "Employer")
pp(D)
D
devtools:use_data(D)
devtools::use_data(D)
devtools::use_data()
helP(devtools::use_data)
help(devtools::use_data)
help(use_data)
email <- D
devtools::use_data(email)
library(hawkes)
help(git_install)
library(hawkes)