plotUnitsECDF Generic plot function to plot the ECDF

Plot ECDF curve

Usage

plotUnitsECDF.qraLm(x, ...)

Arguments

x

qraLm object see Lot2LotGen()

...

Optional plot parameters passed to the function

Author

Vasco Cadavez

Examples

prod <- Lot2LotGen(
                   nLots = 1000,
                   sizeLot = 1000,
                   unitSize = 500,
                   betaAlpha = 0.5112,
                   betaBeta = 9.959,
                   C0MeanLog = 1.023,
                   C0SdLog = 0.3267,
                   propVarInter = 0.7
                   )
prod1 <- Lot2LotGen(
                    nLots = 1000,
                    sizeLot = 1000,
                    unitSize = 500,
                    betaAlpha = 0.9112,
                    betaBeta = 2.959,
                    C0MeanLog = 2.023,
                    C0SdLog = 0.5267,
                    propVarInter = 0.7
                    )
str(prod1)
#> List of 10
#>  $ Lot2LotGenParameters:List of 9
#>   ..$ nLots       : num 1000
#>   ..$ sizeLot     : num 1000
#>   ..$ unitSize    : num 500
#>   ..$ betaAlpha   : num 0.911
#>   ..$ betaBeta    : num 2.96
#>   ..$ C0MeanLog   : num 2.02
#>   ..$ C0SdLog     : num 0.527
#>   ..$ propVarInter: num 0.7
#>   ..$ Poisson     : logi FALSE
#>  $ lotMeans            : num [1:1000] 10.3 23.4 57.2 99.4 14.3 ...
#>  $ unitsCounts         : num [1:1000000] 0 0 0 0 0 0 0 0 0 0 ...
#>  $ N                   : num [1:1000, 1:1000] 0 0 0 0 0 0 0 0 0 0 ...
#>  $ ProbUnitPos         : num [1:1000] 1 1 1 1 1 ...
#>  $ P                   : num 0.995
#>  $ betaGen             : num [1:1000] 0.1896 0.2501 0.4956 0.3845 0.0302 ...
#>  $ nLots               : num 1000
#>  $ sizeLot             : num 1000
#>  $ unitSize            : num 500
#>  - attr(*, "class")= chr "qraLm"
plotUnitsECDF.qraLm(prod1)


p1 <- plotUnitsECDF.qraLm(prod)
p2 <- plotUnitsECDF.qraLm(prod1)
par(mfrow = c(1, 2))
p1

p2

par(mfrow = c(1, 1)) # reset this parameter