Tukey = Function[x,
   l = Length[x];
   xx = Sort[x];
   xx[[Ceiling[l/2]]] +
    0.185*(xx[[Ceiling[0.95*l]]] +
       xx[[Ceiling[0.05*l]]] -
       2*xx[[Ceiling[0.5*l]]])
   ];

n = 50;
p1 = 0;
p2 = 1;
eeps = 5;
eps = 0.1;

dov = 5; (* Доверительная вероятность в процентах *)
dov = dov/2;

v1 = List[];
v1 = Join[RandomVariate[NormalDistribution[p1, p2], n]];

v2 = List[];
v2 = Join[v2, RandomVariate[NormalDistribution[p1, p2], n - n*eps]];
v2 = Join[v2, RandomVariate[NormalDistribution[p1 + eeps, p2], n*eps]];
v2 = RandomSample[v2];

v3 = List[];
v3 = Join[v3, RandomVariate[NormalDistribution[p1, p2], n - n*eps]];
v3 = Join[v3, RandomVariate[NormalDistribution[p1, p2 + eeps], n*eps]];
v3 = RandomSample[v3];

bv11 = List[];
bv12 = List[];
bv21 = List[];
bv22 = List[];
bv31 = List[];
bv32 = List[];

For[i = 1, i <= n, i++,
 bv11 = Append[bv11, Tukey[RandomChoice[v1, Length[v1]]]];
 bv12 = Append[bv12, Median[RandomChoice[v1, Length[v1]]]];

 bv21 = Append[bv21, Tukey[RandomChoice[v2, Length[v2]]]];
 bv22 = Append[bv22, Median[RandomChoice[v2, Length[v2]]]];

 bv31 = Append[bv31, Tukey[RandomChoice[v3, Length[v3]]]];
 bv32 = Append[bv32, Median[RandomChoice[v3, Length[v3]]]];
 ]

text = "";

Print["Оценка Тьюки - ", Mean[bv11], " - ", Variance[bv11],
  " - {", Sort[bv11][[Ceiling[(0 + (dov/100))*Length[bv11]]]], "; ",
  Sort[bv11][[Ceiling[(1 - (dov/100))*Length[bv11]]]], "}"];
Print["Выборочная медиана - ", Mean[bv12], " - ", Variance[bv12],
  " - {", Sort[bv12][[Ceiling[(0 + (dov/100))*Length[bv12]]]], "; ",
  Sort[bv12][[Ceiling[(1 - (dov/100))*Length[bv12]]]], "}"];
Print[];
Print["Оценка Тьюки(ассиметричный выброс) - ", Mean[bv21], " - ",
  Variance[bv21],
  " - {", Sort[bv21][[Ceiling[(0 + (dov/100))*Length[bv21]]]], "; ",
  Sort[bv21][[Ceiling[(1 - (dov/100))*Length[bv21]]]], "}"];
Print["Выборочная медиана(ассиметричный выброс) - ", Mean[bv22],
  " - ", Variance[bv22],
  " - {", Sort[bv22][[Ceiling[(0 + (dov/100))*Length[bv22]]]], "; ",
  Sort[bv22][[Ceiling[(1 - (dov/100))*Length[bv22]]]], "}"];
Print[]
Print["Оценка Тьюки(симметричный выброс) - ", Mean[bv31], " - ",
  Variance[bv31],
  " - {", Sort[bv31][[Ceiling[(0 + (dov/100))*Length[bv31]]]], "; ",
  Sort[bv31][[Ceiling[(1 - (dov/100))*Length[bv31]]]], "}"];
Print["Выборочная медиана(симметричный выброс) - ", Mean[bv32], " - ",
   Variance[bv32],
  " - {", Sort[bv32][[Ceiling[(0 + (dov/100))*Length[bv32]]]], "; ",
  Sort[bv32][[Ceiling[(1 - (dov/100))*Length[bv32]]]], "}"];

\[ScriptCapitalD] = SmoothKernelDistribution[v1];
Print[Table[
   Plot[f[\[ScriptCapitalD], x], {x, p1 - 8, p1 + 8},
    PlotLabel -> f], {f, {PDF, CDF}}]];
\[ScriptCapitalD] = SmoothKernelDistribution[v2];
Print[Table[
   Plot[f[\[ScriptCapitalD], x], {x, p1 - 8, p1 + 8},
    PlotLabel -> f], {f, {PDF, CDF}}]];
\[ScriptCapitalD] = SmoothKernelDistribution[v3];
Print[Table[
   Plot[f[\[ScriptCapitalD], x], {x, p1 - 8, p1 + 8},
    PlotLabel -> f], {f, {PDF, CDF}}]];