
this function prints an array of the bias, variance, risk constants R and CR for three estimators (parametric, non parametric and De Vielder) with different side effect probability for p
Table(file, nb = 10, mod)
the file in wich the simulation will be done
the number of simulation for each estimators
first mod makes simulation for an estimator in different cases (arg2) second mod makes simulation for a case with different estimators third mod makes simulation for a case with the parametric and non parametric estimatorse
Null.
# NOT RUN {
#Table("d.rda",100,3)
# plot the following result
#[1] "arg2Exp"
# nonpar par
#bias 4.407856e-02 -3.503516e-06
#var 6.001352e-05 1.972828e-05
#R 9.405885e-02 4.997679e-02
#CR 8.492920e+00 1.455619e+01
#[1] "arg2Cst"
# nonpar par
#bias 1.575410e-02 4.573187e-04
#var 1.714712e-05 5.675577e-06
#R 3.242554e-02 1.712875e-02
#CR 2.197198e+01 4.066553e+01
#[1] "arg2Comp"
# nonpar par
#bias 3.247060e-02 -5.558817e-05
#var 4.621880e-05 1.006148e-05
#R 6.880894e-02 3.628275e-02
#CR 1.101930e+01 1.967524e+01
#[1] "arg2Gamma"
# nonpar par
#bias 1.832645e-02 -7.552530e-04
#var 2.267637e-05 7.363311e-06
#R 3.778277e-02 1.870107e-02
#CR 1.911443e+01 3.773611e+01
#[1] "arg2Lnorm"
# nonpar par
#bias 4.556463e-02 1.017082e-03
#var 7.013414e-05 3.104221e-05
#R 1.005464e-01 5.599889e-02
#CR 8.101552e+00 1.318976e+01
#[1] "arg2Unif"
# nonpar par
#bias 2.514426e-02 2.598740e-04
#var 3.878331e-05 1.073282e-05
#R 5.236710e-02 2.748271e-02
#CR 1.408411e+01 2.580822e+01
#[1] "arg2Weib"
# nonpar par
#bias 4.588473e-03 1.540856e-04
#var 5.405285e-06 1.589171e-06
#R 9.367470e-03 4.933083e-03
#CR 7.709814e+01 1.468196e+02
#
# }
Run the code above in your browser using DataLab