A numeric value indicating an offset to avoid the initial
iterations. This is usefull to see the likelihood evolution in the last
iterations where the changes in likelihood are of different order of
magnitud than those at the starting iterations.
Examples
# -- apply EMbC to the example path --
mybcp <- stbc(expth)
# -- inspect the likelihood evolution --
lkhp(mybcp)
# -- avoid the initial values --
lkhp(mybcp,10)
Results
R version 3.3.1 (2016-06-21) -- "Bug in Your Hair"
Copyright (C) 2016 The R Foundation for Statistical Computing
Platform: x86_64-pc-linux-gnu (64-bit)
R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.
R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.
Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.
> library(EMbC)
Loading required package: move
Loading required package: geosphere
Loading required package: sp
Loading required package: raster
Loading required package: rgdal
rgdal: version: 1.1-10, (SVN revision 622)
Geospatial Data Abstraction Library extensions to R successfully loaded
Loaded GDAL runtime: GDAL 1.11.3, released 2015/09/16
Path to GDAL shared files: /usr/share/gdal/1.11
Loaded PROJ.4 runtime: Rel. 4.9.2, 08 September 2015, [PJ_VERSION: 492]
Path to PROJ.4 shared files: (autodetected)
Linking to sp version: 1.2-3
> png(filename="/home/ddbj/snapshot/RGM3/R_CC/result/EMbC/lkhp.Rd_%03d_medium.png", width=480, height=480)
> ### Name: lkhp
> ### Title: Likelihood profile plots
> ### Aliases: lkhp lkhp,binClst-method lkhp,list-method
>
> ### ** Examples
>
> # -- apply EMbC to the example path --
> mybcp <- stbc(expth)
... computing starting delimiters:
... delimiter lL : 1.488106
... delimiter lH : 1.488106
... delimiter hL : 1.488106
... delimiter hH : 1.488106
... delimiter Ll : 1.982173
... delimiter Lh : 1.982173
... delimiter Hl : 0.8682456
... delimiter Hh : 0.8682456
[1] 0 -0.0000e+00 4 600
[1] 1 -4.2724e+00 4 447
[1] 2 -4.1181e+00 4 35
[1] 3 -3.6766e+00 4 34
[1] 4 -3.2744e+00 4 24
[1] 5 -3.1224e+00 4 21
[1] 6 -3.0675e+00 4 15
[1] 7 -3.0520e+00 4 7
[1] 8 -3.0456e+00 4 4
[1] 9 -3.0416e+00 4 0
[1] 10 -3.0395e+00 4 0
[1] 11 -3.0384e+00 4 0
[1] 12 -3.0379e+00 4 0
[1] 13 -3.0377e+00 4 0
[1] 14 -3.0376e+00 4 0
[1] ... Stable clustering
> # -- inspect the likelihood evolution --
> lkhp(mybcp)
> # -- avoid the initial values --
> lkhp(mybcp,10)
>
>
>
>
>
> dev.off()
null device
1
>