Nothing
library(RTDE)
#####
# (1) simulation
n <- 100
x <- dataRTDE(simu.nb=n, simu.marg="ufrechet", simu.cop="indep")
print(x)
summary(x)
plot(x, xlab="x", ylab="y")
plot(x, xlab="x", ylab="y", which=2)
n <- 100
y <- dataRTDE(simu.nb=n, simu.marg="ufrechet", simu.cop="indep", contamin.eps=1/100,
contamin.method="max+", contamin.marg="ufrechet", contamin.cop="indep")
print(y)
summary(y)
plot(y, xlab="x", ylab="y")
legend("right", pch=21, col=c("black", "red"), leg=c("orig", "cont"))
plot(y, xlab="x", ylab="y", which=2)
legend("right", pch=21, col=c("black", "red"), leg=c("orig", "cont"))
#####
# (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)
obs <- dataRTDE(cbind(x1, x2))
obs
summary(obs)
plot(obs)
obs <- dataRTDE(cbind(x1, x2), contamin.eps=2/100,
contamin.method="max+", contamin.marg="ufrechet", contamin.cop="indep")
obs
summary(obs)
plot(obs)
plot(obs, which=2)
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.