Nothing
library(RTDE)
#####
# (1) simulation
n <- 100
x <- RTDE(simu=list(nb=n, marg="ufrechet", cop="indep", replicate=1),
nbpoint=10:11, alpha=0, omega=1/2)
x
summary(x)
prob(x, 1:5)
y <- RTDE(simu=list(nb=n, marg="ufrechet", cop="indep", replicate=1),
contamin=list(eps=1/100, method="max+",
marg="ufrechet", cop="indep"),
nbpoint=10:15, alpha=0, omega=1/2)
y
summary(y)
prob(y, 1:5)
#####
# (2) part of the workers' compensation dataset
x1 <- c(
21.798086, 22.640528, 22.572010, 24.789710, 25.876764, 28.033613,
22.525887, 12.004031, 12.713178, 13.596610, 14.811727, 12.774073,
20.245789, 24.242468, 50.216515, 56.099793, 58.109747, 67.807105,
73.852437, 84.208474, 83.604216, 19.507341, 20.810822, 23.838122,
24.212193, 25.367578, 35.401344, 37.580989, 12.428727, 13.492474,
13.303600, 13.228920, 14.060473, 15.227696, 17.857683, 6.829432,
6.219655, 5.696577, 4.542884, 3.663628, 2.632086, 1.932376,
3.724269, 4.149906, 3.696872, 0.487892, 0.594629, 0.560057,
0.724571, 2.494260, 3.394857, 3.314038, 10.595191, 12.863955,
2.789484, 2.426480, 3.140245, 3.964679, 4.899610, 5.247896,
6.280799, 7.721616, 9.324806, 30.035148, 32.618878, 36.071278,
42.494763, 49.524793, 48.872610, 40.636397, 211.719199, 232.681695,
230.007534, 234.635378, 251.852677, 207.512194, 215.038202, 22.525133,
23.545828, 24.979337, 27.390359, 27.337832, 28.268401, 25.281128,
11.939588, 10.566079, 19.387104, 23.062327, 23.089131, 23.799825,
21.317594, 28.171269, 32.008767, 31.300733, 33.507831, 32.207902,
23.471988, 24.101833, 24.766193, 26.078216)
x2 <- c(
0.538707, 0.439184, 1.059775, 0.560013, 1.004997, 1.097314, 0.609833, 0.270222,
0.229566, 0.596850, 0.196539, 0.134248, 0.489312, 0.418218, 0.769208, 0.649707,
0.503919, 0.675466, 0.545745, 1.562266, 0.931762, 0.291125, 0.499927, 0.151084,
0.141910, 0.300373, 0.119761, 0.141300, 0.377662, 0.169574, 0.243585, 0.061215,
0.055272, 0.312816, 0.160196, 0.623029, 0.280707, 0.174422, 0.176666, 0.153907,
0.101874, 0.000000, 0.000000, 0.018082, 0.408802, 0.010500, 0.000000, 0.000000,
0.037500, 0.011234, 0.042146, 0.086599, 0.041500, 0.115592, 0.000000, 0.170623,
0.078436, 0.047501, 0.051499, 0.144621, 0.008998, 0.259875, 0.548369, 0.383658,
0.420764, 0.588094, 0.717523, 1.133547, 1.372065, 0.601273, 2.140846, 2.382040,
2.929906, 3.165114, 3.595126, 4.284873, 2.739517, 0.389158, 0.374650, 0.684493,
0.143933, 1.236701, 0.734178, 0.850133, 0.348156, 0.331977, 0.116060, 0.266805,
0.310604, 0.278189, 0.079680, 1.238292, 0.576512, 0.831404, 1.303705, 0.673907,
0.605122, 0.664457, 0.348918, 0.370878)
z <- RTDE(obs=cbind(x1, x2), nbpoint=10, alpha=0, omega=1/2:3)
z
summary(z)
prob(z, 1:5)
z1 <- RTDE(obs=cbind(x1, x2), contamin=list(eps=1/100, method="max+",
marg="ufrechet", cop="indep"), nbpoint=10:15, alpha=0, omega=1/2)
z1
summary(z1)
prob(z1, 1:5)
z2 <- RTDE(obs=cbind(x1, x2), contamin=list(eps=1/100, method="max+",
marg="ufrechet", cop="indep"), nbpoint=10:15, alpha=0:1, omega=1/2)
z2
summary(z2)
prob(z2, 1:5)
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.