mlica: Maximum likelihood implementation of independent component...
In mlica2: Independent Component Analysis using Maximum Likelihood
Description
Usage
Arguments
Value
Author(s)
References
Examples
Description
This function performs ICA using a maximum likelihood framework
and takes as arguments parameters to control the number of
algorithm runs and convergence criteria.
Usage
1
mlica(prNCP, nruns = 10, tol = 1e-04, maxit = 300, fail.th = 5, learn.mu = 1)
Arguments
prNCP
The output object from 'proposeNCP'.
nruns
The number of converged algorithm runs sought (function
returns the best solution according to the log-likelihood
value).
tol
Tolerance level for establishing convergence of run.
maxit
Maximum number of iterations to allow per run.
fail.th
A threshold on the number of consecutive runs that fail to
converge.
learn.mu
Learning parameter for fixed point algorithm (note that this
need not be changed since it has already been optimised).
Value
A list with following components:
A: Estimate of the mixing matrix.
B: Estimate of the inverse mixing matrix.
S: Estimate of the source matrix.
X: Normalised data matrix.
ncp: Number of independent components.
NC: Binary number specifying whether best run converged or
not.(=1 indicates convergence,=0 indicates no convergence).
LL: Log likelihood value of best run.
Author(s)
Andrew Teschendorff a.teschendorff@ucl.ac.uk
References
Hyvaerinen A., Karhunen J., and Oja E.: Independent Component
Analysis, John Wiley and Sons, New York, (2001).
Kreil D. and MacKay D. (2003): Reproducibility Assessment of
Independent Component Analysis of Expression Ratios from DNA
microarrays, Comparative and Functional Genomics *4* (3),300-317.
Liebermeister W. (2002): Linear Modes of gene expression determined
by independent component analysis, Bioinformatics *18*, no.1,
51-60.
Chiappetta P., Roubaud MC. and Torresani B.: Blind source separation
and the analysis of microarray data, J. Comput. Biol. 2004;
11(6):1090-109.
Examples
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
## Not run:
data(simMAdata);
dataX <- simMAdata[[1]];
prPCA <- PriorNormPCA(dataX);
prNCP <- proposeNCP(prPCA,0.1);
a.best.l <- list();
for( i in 1:5){
a.best.l[[i]] <- mlica(prNCP,nruns=5);
}
checkICA <- CheckStability(a.best.l,0.7);
sourceS <- simMAdata[[3]];
print(cor(a.best.l[[1]]$S,sourceS));
sModes <- SortModes(a.best.l[[1]],c.val=0.5);
## End(Not run)
## The function is currently defined as
function (prNCP, nruns = 10, tol = 1e-04, maxit = 300, fail.th = 5,
learn.mu = 1)
{
print("Entering mlica")
print("Performing preliminary run")
a <- mlicaMAIN(prNCP, tol = 1e-04, maxit = 10, mu = learn.mu)
ncp <- dim(a$S)[2]
max.logL <- a$LL
a.best <- a
print("Finished preliminary run")
print("Starting runs")
run.n <- 0
fail.count <- 0
v.logL <- vector()
v.NC <- vector()
while (run.n < nruns) {
a <- mlicaMAIN(prNCP, tol = 1e-04, maxit = maxit, mu = learn.mu)
v.logL <- c(v.logL, a$LL)
v.NC <- c(v.NC, a$NC)
if (a$NC == 0) {
fail.count <- 0
run.n <- run.n + 1
if (a$LL > max.logL) {
a.best <- a
}
}
else {
fail.count <- fail.count + 1
}
if (fail.count >= fail.th) {
print("Stopping: Five consecutive runs failed to converge!")
print("Consider either increasing the threshold for pca eigenvalues to perform ICA on a smaller subspace or increasing maxit")
stop
}
}
print("End of runs")
return(a.best)
}
Example output
[1] "Performing SVD"
Warning message:
In svd(X, LINPACK = FALSE) :
the LINPACK argument has been defunct since R 3.1.0
[1] "About to find ncp"
[1] "Entering mlica"
[1] "Performing preliminary run"
[1] "Entered MLica"
[1] "Entering iteration loop " "0"
[1] "Checkpt1"
[1] "AvDev=" "0.290219314280639"
[1] "iterated logL"
[1] -3590.578
[1] "Entering iteration loop " "1"
[1] "Checkpt1"
[1] "AvDev=" "0.105827278122874"
[1] "iterated logL"
[1] -3472.839
[1] "Entering iteration loop " "2"
[1] "Checkpt1"
[1] "AvDev=" "0.0471382644915915"
[1] "iterated logL"
[1] -3421.788
[1] "Entering iteration loop " "3"
[1] "Checkpt1"
[1] "AvDev=" "0.0264127579726401"
[1] "iterated logL"
[1] -3401.135
[1] "Entering iteration loop " "4"
[1] "Checkpt1"
[1] "AvDev=" "0.0142316603037475"
[1] "iterated logL"
[1] -3395.756
[1] "Entering iteration loop " "5"
[1] "Checkpt1"
[1] "AvDev=" "0.0048090862414591"
[1] "iterated logL"
[1] -3395.325
[1] "Entering iteration loop " "6"
[1] "Checkpt1"
[1] "AvDev=" "0.00127812851800342"
[1] "iterated logL"
[1] -3395.483
[1] "Entering iteration loop " "7"
[1] "Checkpt1"
[1] "AvDev=" "0.000311215663334201"
[1] "iterated logL"
[1] -3395.535
[1] "Entering iteration loop " "8"
[1] "Checkpt1"
[1] "AvDev=" "7.47703362656165e-05"
[1] "iterated logL"
[1] -3395.548
[1] "Entering iteration loop " "9"
[1] "Checkpt1"
[1] "AvDev=" "1.78819501199913e-05"
[1] "iterated logL"
[1] -3395.551
[1] "Finished preliminary run"
[1] "Starting runs"
[1] "Entered MLica"
[1] "Entering iteration loop " "0"
[1] "Checkpt1"
[1] "AvDev=" "0.400979712633759"
[1] "iterated logL"
[1] -4423.601
[1] "Entering iteration loop " "1"
[1] "Checkpt1"
[1] "AvDev=" "0.179739366038515"
[1] "iterated logL"
[1] -3473.989
[1] "Entering iteration loop " "2"
[1] "Checkpt1"
[1] "AvDev=" "0.0444486673678738"
[1] "iterated logL"
[1] -3398.851
[1] "Entering iteration loop " "3"
[1] "Checkpt1"
[1] "AvDev=" "0.00487559331142786"
[1] "iterated logL"
[1] -3395.584
[1] "Entering iteration loop " "4"
[1] "Checkpt1"
[1] "AvDev=" "0.000858328754907789"
[1] "iterated logL"
[1] -3395.607
[1] "Entering iteration loop " "5"
[1] "Checkpt1"
[1] "AvDev=" "0.000207039556067208"
[1] "iterated logL"
[1] -3395.554
[1] "Entering iteration loop " "6"
[1] "Checkpt1"
[1] "AvDev=" "4.13203110180985e-05"
[1] "iterated logL"
[1] -3395.554
[1] "Entering iteration loop " "7"
[1] "Checkpt1"
[1] "AvDev=" "1.05896434797988e-05"
[1] "iterated logL"
[1] -3395.552
[1] "Entered MLica"
[1] "Entering iteration loop " "0"
[1] "Checkpt1"
[1] "AvDev=" "0.463498076168755"
[1] "iterated logL"
[1] -3667.793
[1] "Entering iteration loop " "1"
[1] "Checkpt1"
[1] "AvDev=" "0.0764780178146065"
[1] "iterated logL"
[1] -3510.967
[1] "Entering iteration loop " "2"
[1] "Checkpt1"
[1] "AvDev=" "0.0842512230811305"
[1] "iterated logL"
[1] -3424.638
[1] "Entering iteration loop " "3"
[1] "Checkpt1"
[1] "AvDev=" "0.0322542264738013"
[1] "iterated logL"
[1] -3397.497
[1] "Entering iteration loop " "4"
[1] "Checkpt1"
[1] "AvDev=" "0.00461314705312814"
[1] "iterated logL"
[1] -3395.5
[1] "Entering iteration loop " "5"
[1] "Checkpt1"
[1] "AvDev=" "0.00065829292806568"
[1] "iterated logL"
[1] -3395.574
[1] "Entering iteration loop " "6"
[1] "Checkpt1"
[1] "AvDev=" "0.000113873449825682"
[1] "iterated logL"
[1] -3395.549
[1] "Entering iteration loop " "7"
[1] "Checkpt1"
[1] "AvDev=" "2.12979921296829e-05"
[1] "iterated logL"
[1] -3395.553
[1] "Entering iteration loop " "8"
[1] "Checkpt1"
[1] "AvDev=" "3.54673673779438e-06"
[1] "iterated logL"
[1] -3395.552
[1] "Entered MLica"
[1] "Entering iteration loop " "0"
[1] "Checkpt1"
[1] "AvDev=" "0.36629807731381"
[1] "iterated logL"
[1] -3992.956
[1] "Entering iteration loop " "1"
[1] "Checkpt1"
[1] "AvDev=" "0.100107945866601"
[1] "iterated logL"
[1] -3800.148
[1] "Entering iteration loop " "2"
[1] "Checkpt1"
[1] "AvDev=" "0.0747871283230001"
[1] "iterated logL"
[1] -3577.054
[1] "Entering iteration loop " "3"
[1] "Checkpt1"
[1] "AvDev=" "0.0508704848282896"
[1] "iterated logL"
[1] -3478.123
[1] "Entering iteration loop " "4"
[1] "Checkpt1"
[1] "AvDev=" "0.0516390652426347"
[1] "iterated logL"
[1] -3418.488
[1] "Entering iteration loop " "5"
[1] "Checkpt1"
[1] "AvDev=" "0.0270751636990737"
[1] "iterated logL"
[1] -3399.461
[1] "Entering iteration loop " "6"
[1] "Checkpt1"
[1] "AvDev=" "0.00507287047873504"
[1] "iterated logL"
[1] -3395.689
[1] "Entering iteration loop " "7"
[1] "Checkpt1"
[1] "AvDev=" "0.000827287456690507"
[1] "iterated logL"
[1] -3395.541
[1] "Entering iteration loop " "8"
[1] "Checkpt1"
[1] "AvDev=" "0.000156985449468821"
[1] "iterated logL"
[1] -3395.549
[1] "Entering iteration loop " "9"
[1] "Checkpt1"
[1] "AvDev=" "3.27897912097372e-05"
[1] "iterated logL"
[1] -3395.551
[1] "Entering iteration loop " "10"
[1] "Checkpt1"
[1] "AvDev=" "7.02164083348528e-06"
[1] "iterated logL"
[1] -3395.552
[1] "Entered MLica"
[1] "Entering iteration loop " "0"
[1] "Checkpt1"
[1] "AvDev=" "0.339996295702037"
[1] "iterated logL"
[1] -3734.831
[1] "Entering iteration loop " "1"
[1] "Checkpt1"
[1] "AvDev=" "0.115184530320527"
[1] "iterated logL"
[1] -3664.789
[1] "Entering iteration loop " "2"
[1] "Checkpt1"
[1] "AvDev=" "0.0761931714957583"
[1] "iterated logL"
[1] -3513.394
[1] "Entering iteration loop " "3"
[1] "Checkpt1"
[1] "AvDev=" "0.0395665341671519"
[1] "iterated logL"
[1] -3468.643
[1] "Entering iteration loop " "4"
[1] "Checkpt1"
[1] "AvDev=" "0.0329772639881574"
[1] "iterated logL"
[1] -3425.284
[1] "Entering iteration loop " "5"
[1] "Checkpt1"
[1] "AvDev=" "0.0131235875899943"
[1] "iterated logL"
[1] -3398.354
[1] "Entering iteration loop " "6"
[1] "Checkpt1"
[1] "AvDev=" "0.00317390696458019"
[1] "iterated logL"
[1] -3395.759
[1] "Entering iteration loop " "7"
[1] "Checkpt1"
[1] "AvDev=" "0.000725960885815795"
[1] "iterated logL"
[1] -3395.573
[1] "Entering iteration loop " "8"
[1] "Checkpt1"
[1] "AvDev=" "0.000163452040045647"
[1] "iterated logL"
[1] -3395.555
[1] "Entering iteration loop " "9"
[1] "Checkpt1"
[1] "AvDev=" "3.72448459157757e-05"
[1] "iterated logL"
[1] -3395.552
[1] "Entering iteration loop " "10"
[1] "Checkpt1"
[1] "AvDev=" "8.40711799153595e-06"
[1] "iterated logL"
[1] -3395.552
[1] "Entered MLica"
[1] "Entering iteration loop " "0"
[1] "Checkpt1"
[1] "AvDev=" "0.364821563279185"
[1] "iterated logL"
[1] -4328.342
[1] "Entering iteration loop " "1"
[1] "Checkpt1"
[1] "AvDev=" "0.115872175582269"
[1] "iterated logL"
[1] -3754.574
[1] "Entering iteration loop " "2"
[1] "Checkpt1"
[1] "AvDev=" "0.109689364039819"
[1] "iterated logL"
[1] -3431.945
[1] "Entering iteration loop " "3"
[1] "Checkpt1"
[1] "AvDev=" "0.0359576994516968"
[1] "iterated logL"
[1] -3395.668
[1] "Entering iteration loop " "4"
[1] "Checkpt1"
[1] "AvDev=" "0.00633261188333669"
[1] "iterated logL"
[1] -3395.861
[1] "Entering iteration loop " "5"
[1] "Checkpt1"
[1] "AvDev=" "0.00083586484581874"
[1] "iterated logL"
[1] -3395.485
[1] "Entering iteration loop " "6"
[1] "Checkpt1"
[1] "AvDev=" "0.000156792108413645"
[1] "iterated logL"
[1] -3395.561
[1] "Entering iteration loop " "7"
[1] "Checkpt1"
[1] "AvDev=" "2.85094680056924e-05"
[1] "iterated logL"
[1] -3395.549
[1] "Entering iteration loop " "8"
[1] "Checkpt1"
[1] "AvDev=" "5.68654754230662e-06"
[1] "iterated logL"
[1] -3395.552
[1] "End of runs"
[1] "Entering mlica"
[1] "Performing preliminary run"
[1] "Entered MLica"
[1] "Entering iteration loop " "0"
[1] "Checkpt1"
[1] "AvDev=" "0.386038204478516"
[1] "iterated logL"
[1] -4391.043
[1] "Entering iteration loop " "1"
[1] "Checkpt1"
[1] "AvDev=" "0.118722047951209"
[1] "iterated logL"
[1] -3787.804
[1] "Entering iteration loop " "2"
[1] "Checkpt1"
[1] "AvDev=" "0.1160421413945"
[1] "iterated logL"
[1] -3476.545
[1] "Entering iteration loop " "3"
[1] "Checkpt1"
[1] "AvDev=" "0.0606436361056195"
[1] "iterated logL"
[1] -3400.034
[1] "Entering iteration loop " "4"
[1] "Checkpt1"
[1] "AvDev=" "0.00988526357087198"
[1] "iterated logL"
[1] -3395.5
[1] "Entering iteration loop " "5"
[1] "Checkpt1"
[1] "AvDev=" "0.00112709094447862"
[1] "iterated logL"
[1] -3395.475
[1] "Entering iteration loop " "6"
[1] "Checkpt1"
[1] "AvDev=" "0.000248549371090569"
[1] "iterated logL"
[1] -3395.548
[1] "Entering iteration loop " "7"
[1] "Checkpt1"
[1] "AvDev=" "5.46497451578651e-05"
[1] "iterated logL"
[1] -3395.549
[1] "Entering iteration loop " "8"
[1] "Checkpt1"
[1] "AvDev=" "1.21802304175041e-05"
[1] "iterated logL"
[1] -3395.552
[1] "Finished preliminary run"
[1] "Starting runs"
[1] "Entered MLica"
[1] "Entering iteration loop " "0"
[1] "Checkpt1"
[1] "AvDev=" "0.356261911078639"
[1] "iterated logL"
[1] -3672.559
[1] "Entering iteration loop " "1"
[1] "Checkpt1"
[1] "AvDev=" "0.0883629996095928"
[1] "iterated logL"
[1] -3405.789
[1] "Entering iteration loop " "2"
[1] "Checkpt1"
[1] "AvDev=" "0.0115397758858259"
[1] "iterated logL"
[1] -3395.782
[1] "Entering iteration loop " "3"
[1] "Checkpt1"
[1] "AvDev=" "0.00159642887472422"
[1] "iterated logL"
[1] -3395.648
[1] "Entering iteration loop " "4"
[1] "Checkpt1"
[1] "AvDev=" "0.00036227500675863"
[1] "iterated logL"
[1] -3395.562
[1] "Entering iteration loop " "5"
[1] "Checkpt1"
[1] "AvDev=" "7.82784397283197e-05"
[1] "iterated logL"
[1] -3395.556
[1] "Entering iteration loop " "6"
[1] "Checkpt1"
[1] "AvDev=" "1.80340204448921e-05"
[1] "iterated logL"
[1] -3395.552
[1] "Entered MLica"
[1] "Entering iteration loop " "0"
[1] "Checkpt1"
[1] "AvDev=" "0.38765322127385"
[1] "iterated logL"
[1] -3673.776
[1] "Entering iteration loop " "1"
[1] "Checkpt1"
[1] "AvDev=" "0.11495659664447"
[1] "iterated logL"
[1] -3472.173
[1] "Entering iteration loop " "2"
[1] "Checkpt1"
[1] "AvDev=" "0.0386829672942188"
[1] "iterated logL"
[1] -3405.763
[1] "Entering iteration loop " "3"
[1] "Checkpt1"
[1] "AvDev=" "0.00826699766804716"
[1] "iterated logL"
[1] -3396.022
[1] "Entering iteration loop " "4"
[1] "Checkpt1"
[1] "AvDev=" "0.00153531227230042"
[1] "iterated logL"
[1] -3395.556
[1] "Entering iteration loop " "5"
[1] "Checkpt1"
[1] "AvDev=" "0.000352005254312435"
[1] "iterated logL"
[1] -3395.548
[1] "Entering iteration loop " "6"
[1] "Checkpt1"
[1] "AvDev=" "7.73710107153997e-05"
[1] "iterated logL"
[1] -3395.551
[1] "Entering iteration loop " "7"
[1] "Checkpt1"
[1] "AvDev=" "1.78830343199412e-05"
[1] "iterated logL"
[1] -3395.552
[1] "Entered MLica"
[1] "Entering iteration loop " "0"
[1] "Checkpt1"
[1] "AvDev=" "0.392189272097064"
[1] "iterated logL"
[1] -5432.134
[1] "Entering iteration loop " "1"
[1] "Checkpt1"
[1] "AvDev=" "0.1417022380835"
[1] "iterated logL"
[1] -4041.275
[1] "Entering iteration loop " "2"
[1] "Checkpt1"
[1] "AvDev=" "0.0796615817576913"
[1] "iterated logL"
[1] -3755.534
[1] "Entering iteration loop " "3"
[1] "Checkpt1"
[1] "AvDev=" "0.0693391715870864"
[1] "iterated logL"
[1] -3506.783
[1] "Entering iteration loop " "4"
[1] "Checkpt1"
[1] "AvDev=" "0.0486440891664531"
[1] "iterated logL"
[1] -3424.152
[1] "Entering iteration loop " "5"
[1] "Checkpt1"
[1] "AvDev=" "0.0221261664494838"
[1] "iterated logL"
[1] -3397.661
[1] "Entering iteration loop " "6"
[1] "Checkpt1"
[1] "AvDev=" "0.0028671746044226"
[1] "iterated logL"
[1] -3395.543
[1] "Entering iteration loop " "7"
[1] "Checkpt1"
[1] "AvDev=" "0.000574550590411959"
[1] "iterated logL"
[1] -3395.546
[1] "Entering iteration loop " "8"
[1] "Checkpt1"
[1] "AvDev=" "0.000125694040543018"
[1] "iterated logL"
[1] -3395.548
[1] "Entering iteration loop " "9"
[1] "Checkpt1"
[1] "AvDev=" "2.68782845113823e-05"
[1] "iterated logL"
[1] -3395.551
[1] "Entering iteration loop " "10"
[1] "Checkpt1"
[1] "AvDev=" "5.96964466208338e-06"
[1] "iterated logL"
[1] -3395.552
[1] "Entered MLica"
[1] "Entering iteration loop " "0"
[1] "Checkpt1"
[1] "AvDev=" "0.410912985876807"
[1] "iterated logL"
[1] -4043.816
[1] "Entering iteration loop " "1"
[1] "Checkpt1"
[1] "AvDev=" "0.107261300936882"
[1] "iterated logL"
[1] -3911.253
[1] "Entering iteration loop " "2"
[1] "Checkpt1"
[1] "AvDev=" "0.0600342739500788"
[1] "iterated logL"
[1] -3916.637
[1] "Entering iteration loop " "3"
[1] "Checkpt1"
[1] "AvDev=" "0.0441894098933856"
[1] "iterated logL"
[1] -3758.087
[1] "Entering iteration loop " "4"
[1] "Checkpt1"
[1] "AvDev=" "0.0518400801703871"
[1] "iterated logL"
[1] -3561.38
[1] "Entering iteration loop " "5"
[1] "Checkpt1"
[1] "AvDev=" "0.0501583042138534"
[1] "iterated logL"
[1] -3413.922
[1] "Entering iteration loop " "6"
[1] "Checkpt1"
[1] "AvDev=" "0.0218865273541719"
[1] "iterated logL"
[1] -3396.987
[1] "Entering iteration loop " "7"
[1] "Checkpt1"
[1] "AvDev=" "0.00430632415238061"
[1] "iterated logL"
[1] -3395.452
[1] "Entering iteration loop " "8"
[1] "Checkpt1"
[1] "AvDev=" "0.000742641245699041"
[1] "iterated logL"
[1] -3395.533
[1] "Entering iteration loop " "9"
[1] "Checkpt1"
[1] "AvDev=" "0.000150363951070377"
[1] "iterated logL"
[1] -3395.545
[1] "Entering iteration loop " "10"
[1] "Checkpt1"
[1] "AvDev=" "3.29787956548537e-05"
[1] "iterated logL"
[1] -3395.55
[1] "Entering iteration loop " "11"
[1] "Checkpt1"
[1] "AvDev=" "7.32844518741213e-06"
[1] "iterated logL"
[1] -3395.551
[1] "Entered MLica"
[1] "Entering iteration loop " "0"
[1] "Checkpt1"
[1] "AvDev=" "0.371274126606981"
[1] "iterated logL"
[1] -3648.194
[1] "Entering iteration loop " "1"
[1] "Checkpt1"
[1] "AvDev=" "0.0801941640029815"
[1] "iterated logL"
[1] -3513.596
[1] "Entering iteration loop " "2"
[1] "Checkpt1"
[1] "AvDev=" "0.0744793555550997"
[1] "iterated logL"
[1] -3448.159
[1] "Entering iteration loop " "3"
[1] "Checkpt1"
[1] "AvDev=" "0.0525100410417017"
[1] "iterated logL"
[1] -3414.35
[1] "Entering iteration loop " "4"
[1] "Checkpt1"
[1] "AvDev=" "0.0283072099719519"
[1] "iterated logL"
[1] -3400.12
[1] "Entering iteration loop " "5"
[1] "Checkpt1"
[1] "AvDev=" "0.0129825402674149"
[1] "iterated logL"
[1] -3395.6
[1] "Entering iteration loop " "6"
[1] "Checkpt1"
[1] "AvDev=" "0.00424652217030174"
[1] "iterated logL"
[1] -3395.338
[1] "Entering iteration loop " "7"
[1] "Checkpt1"
[1] "AvDev=" "0.00111200174326132"
[1] "iterated logL"
[1] -3395.493
[1] "Entering iteration loop " "8"
[1] "Checkpt1"
[1] "AvDev=" "0.000269207446763617"
[1] "iterated logL"
[1] -3395.537
[1] "Entering iteration loop " "9"
[1] "Checkpt1"
[1] "AvDev=" "6.47633209478762e-05"
[1] "iterated logL"
[1] -3395.548
[1] "Entering iteration loop " "10"
[1] "Checkpt1"
[1] "AvDev=" "1.5395072232993e-05"
[1] "iterated logL"
[1] -3395.551
[1] "End of runs"
[1] "Entering mlica"
[1] "Performing preliminary run"
[1] "Entered MLica"
[1] "Entering iteration loop " "0"
[1] "Checkpt1"
[1] "AvDev=" "0.296028470334294"
[1] "iterated logL"
[1] -3809.782
[1] "Entering iteration loop " "1"
[1] "Checkpt1"
[1] "AvDev=" "0.0894398874126622"
[1] "iterated logL"
[1] -3594.638
[1] "Entering iteration loop " "2"
[1] "Checkpt1"
[1] "AvDev=" "0.0591590009051705"
[1] "iterated logL"
[1] -3459.552
[1] "Entering iteration loop " "3"
[1] "Checkpt1"
[1] "AvDev=" "0.0200699441239609"
[1] "iterated logL"
[1] -3439.458
[1] "Entering iteration loop " "4"
[1] "Checkpt1"
[1] "AvDev=" "0.0242748564489362"
[1] "iterated logL"
[1] -3429.975
[1] "Entering iteration loop " "5"
[1] "Checkpt1"
[1] "AvDev=" "0.0308904719620776"
[1] "iterated logL"
[1] -3411.681
[1] "Entering iteration loop " "6"
[1] "Checkpt1"
[1] "AvDev=" "0.0218622093012664"
[1] "iterated logL"
[1] -3399.13
[1] "Entering iteration loop " "7"
[1] "Checkpt1"
[1] "AvDev=" "0.00865223324385435"
[1] "iterated logL"
[1] -3396.15
[1] "Entering iteration loop " "8"
[1] "Checkpt1"
[1] "AvDev=" "0.00230748595510655"
[1] "iterated logL"
[1] -3395.681
[1] "Entering iteration loop " "9"
[1] "Checkpt1"
[1] "AvDev=" "0.000523761338507808"
[1] "iterated logL"
[1] -3395.582
[1] "Finished preliminary run"
[1] "Starting runs"
[1] "Entered MLica"
[1] "Entering iteration loop " "0"
[1] "Checkpt1"
[1] "AvDev=" "0.280196120957882"
[1] "iterated logL"
[1] -4126.692
[1] "Entering iteration loop " "1"
[1] "Checkpt1"
[1] "AvDev=" "0.136967268241482"
[1] "iterated logL"
[1] -3528.859
[1] "Entering iteration loop " "2"
[1] "Checkpt1"
[1] "AvDev=" "0.063909758039693"
[1] "iterated logL"
[1] -3418.22
[1] "Entering iteration loop " "3"
[1] "Checkpt1"
[1] "AvDev=" "0.0136820341896695"
[1] "iterated logL"
[1] -3396.547
[1] "Entering iteration loop " "4"
[1] "Checkpt1"
[1] "AvDev=" "0.00229592432874156"
[1] "iterated logL"
[1] -3395.697
[1] "Entering iteration loop " "5"
[1] "Checkpt1"
[1] "AvDev=" "0.000501495866691987"
[1] "iterated logL"
[1] -3395.569
[1] "Entering iteration loop " "6"
[1] "Checkpt1"
[1] "AvDev=" "0.000114241962480851"
[1] "iterated logL"
[1] -3395.558
[1] "Entering iteration loop " "7"
[1] "Checkpt1"
[1] "AvDev=" "2.56100609335057e-05"
[1] "iterated logL"
[1] -3395.553
[1] "Entered MLica"
[1] "Entering iteration loop " "0"
[1] "Checkpt1"
[1] "AvDev=" "0.3591373454329"
[1] "iterated logL"
[1] -3649.387
[1] "Entering iteration loop " "1"
[1] "Checkpt1"
[1] "AvDev=" "0.0870331784836628"
[1] "iterated logL"
[1] -3474.703
[1] "Entering iteration loop " "2"
[1] "Checkpt1"
[1] "AvDev=" "0.0318251790046466"
[1] "iterated logL"
[1] -3442.092
[1] "Entering iteration loop " "3"
[1] "Checkpt1"
[1] "AvDev=" "0.026413729802557"
[1] "iterated logL"
[1] -3428.089
[1] "Entering iteration loop " "4"
[1] "Checkpt1"
[1] "AvDev=" "0.0305909886580975"
[1] "iterated logL"
[1] -3409.182
[1] "Entering iteration loop " "5"
[1] "Checkpt1"
[1] "AvDev=" "0.0206336363315325"
[1] "iterated logL"
[1] -3397.634
[1] "Entering iteration loop " "6"
[1] "Checkpt1"
[1] "AvDev=" "0.0089562560436709"
[1] "iterated logL"
[1] -3395.333
[1] "Entering iteration loop " "7"
[1] "Checkpt1"
[1] "AvDev=" "0.00269255449333509"
[1] "iterated logL"
[1] -3395.407
[1] "Entering iteration loop " "8"
[1] "Checkpt1"
[1] "AvDev=" "0.000679790417210284"
[1] "iterated logL"
[1] -3395.516
[1] "Entering iteration loop " "9"
[1] "Checkpt1"
[1] "AvDev=" "0.000163404609437547"
[1] "iterated logL"
[1] -3395.543
[1] "Entering iteration loop " "10"
[1] "Checkpt1"
[1] "AvDev=" "3.9450123678408e-05"
[1] "iterated logL"
[1] -3395.55
[1] "Entering iteration loop " "11"
[1] "Checkpt1"
[1] "AvDev=" "9.41764292642783e-06"
[1] "iterated logL"
[1] -3395.551
[1] "Entered MLica"
[1] "Entering iteration loop " "0"
[1] "Checkpt1"
[1] "AvDev=" "0.375329432885969"
[1] "iterated logL"
[1] -4102.084
[1] "Entering iteration loop " "1"
[1] "Checkpt1"
[1] "AvDev=" "0.0914017749959568"
[1] "iterated logL"
[1] -3852.777
[1] "Entering iteration loop " "2"
[1] "Checkpt1"
[1] "AvDev=" "0.0819096172436832"
[1] "iterated logL"
[1] -3614.669
[1] "Entering iteration loop " "3"
[1] "Checkpt1"
[1] "AvDev=" "0.0640237802367762"
[1] "iterated logL"
[1] -3468.583
[1] "Entering iteration loop " "4"
[1] "Checkpt1"
[1] "AvDev=" "0.0287827456740873"
[1] "iterated logL"
[1] -3410.449
[1] "Entering iteration loop " "5"
[1] "Checkpt1"
[1] "AvDev=" "0.00832633461415859"
[1] "iterated logL"
[1] -3396.645
[1] "Entering iteration loop " "6"
[1] "Checkpt1"
[1] "AvDev=" "0.00187009118538909"
[1] "iterated logL"
[1] -3395.638
[1] "Entering iteration loop " "7"
[1] "Checkpt1"
[1] "AvDev=" "0.000413676242951251"
[1] "iterated logL"
[1] -3395.563
[1] "Entering iteration loop " "8"
[1] "Checkpt1"
[1] "AvDev=" "9.46708401711417e-05"
[1] "iterated logL"
[1] -3395.554
[1] "Entering iteration loop " "9"
[1] "Checkpt1"
[1] "AvDev=" "2.13034932697826e-05"
[1] "iterated logL"
[1] -3395.552
[1] "Entering iteration loop " "10"
[1] "Checkpt1"
[1] "AvDev=" "4.80830754306893e-06"
[1] "iterated logL"
[1] -3395.552
[1] "Entered MLica"
[1] "Entering iteration loop " "0"
[1] "Checkpt1"
[1] "AvDev=" "0.414533771616944"
[1] "iterated logL"
[1] -3619.764
[1] "Entering iteration loop " "1"
[1] "Checkpt1"
[1] "AvDev=" "0.0998829443199397"
[1] "iterated logL"
[1] -3428.315
[1] "Entering iteration loop " "2"
[1] "Checkpt1"
[1] "AvDev=" "0.0371035179271871"
[1] "iterated logL"
[1] -3395.115
[1] "Entering iteration loop " "3"
[1] "Checkpt1"
[1] "AvDev=" "0.00585473854860562"
[1] "iterated logL"
[1] -3395.705
[1] "Entering iteration loop " "4"
[1] "Checkpt1"
[1] "AvDev=" "0.00102639880166599"
[1] "iterated logL"
[1] -3395.46
[1] "Entering iteration loop " "5"
[1] "Checkpt1"
[1] "AvDev=" "0.000233825810315878"
[1] "iterated logL"
[1] -3395.553
[1] "Entering iteration loop " "6"
[1] "Checkpt1"
[1] "AvDev=" "4.63066237271436e-05"
[1] "iterated logL"
[1] -3395.548
[1] "Entering iteration loop " "7"
[1] "Checkpt1"
[1] "AvDev=" "1.13079987055212e-05"
[1] "iterated logL"
[1] -3395.552
[1] "Entered MLica"
[1] "Entering iteration loop " "0"
[1] "Checkpt1"
[1] "AvDev=" "0.326979966868656"
[1] "iterated logL"
[1] -3635.347
[1] "Entering iteration loop " "1"
[1] "Checkpt1"
[1] "AvDev=" "0.0958728068461986"
[1] "iterated logL"
[1] -3460.853
[1] "Entering iteration loop " "2"
[1] "Checkpt1"
[1] "AvDev=" "0.0571714978175811"
[1] "iterated logL"
[1] -3405.342
[1] "Entering iteration loop " "3"
[1] "Checkpt1"
[1] "AvDev=" "0.011807503859014"
[1] "iterated logL"
[1] -3395.566
[1] "Entering iteration loop " "4"
[1] "Checkpt1"
[1] "AvDev=" "0.00163925954499247"
[1] "iterated logL"
[1] -3395.508
[1] "Entering iteration loop " "5"
[1] "Checkpt1"
[1] "AvDev=" "0.000385225711069099"
[1] "iterated logL"
[1] -3395.538
[1] "Entering iteration loop " "6"
[1] "Checkpt1"
[1] "AvDev=" "9.01734293590305e-05"
[1] "iterated logL"
[1] -3395.549
[1] "Entering iteration loop " "7"
[1] "Checkpt1"
[1] "AvDev=" "2.12590118020316e-05"
[1] "iterated logL"
[1] -3395.551
[1] "Entering iteration loop " "8"
[1] "Checkpt1"
[1] "AvDev=" "5.08094001730523e-06"
[1] "iterated logL"
[1] -3395.552
[1] "End of runs"
[1] "Entering mlica"
[1] "Performing preliminary run"
[1] "Entered MLica"
[1] "Entering iteration loop " "0"
[1] "Checkpt1"
[1] "AvDev=" "0.438000141407069"
[1] "iterated logL"
[1] -4276.765
[1] "Entering iteration loop " "1"
[1] "Checkpt1"
[1] "AvDev=" "0.154084115234504"
[1] "iterated logL"
[1] -3499.516
[1] "Entering iteration loop " "2"
[1] "Checkpt1"
[1] "AvDev=" "0.0596482103812549"
[1] "iterated logL"
[1] -3397.822
[1] "Entering iteration loop " "3"
[1] "Checkpt1"
[1] "AvDev=" "0.00981867306151161"
[1] "iterated logL"
[1] -3395.482
[1] "Entering iteration loop " "4"
[1] "Checkpt1"
[1] "AvDev=" "0.00217984430913429"
[1] "iterated logL"
[1] -3395.449
[1] "Entering iteration loop " "5"
[1] "Checkpt1"
[1] "AvDev=" "0.000498978167806657"
[1] "iterated logL"
[1] -3395.529
[1] "Entering iteration loop " "6"
[1] "Checkpt1"
[1] "AvDev=" "0.000111394576278191"
[1] "iterated logL"
[1] -3395.546
[1] "Entering iteration loop " "7"
[1] "Checkpt1"
[1] "AvDev=" "2.70277136988498e-05"
[1] "iterated logL"
[1] -3395.55
[1] "Entering iteration loop " "8"
[1] "Checkpt1"
[1] "AvDev=" "6.12100051378475e-06"
[1] "iterated logL"
[1] -3395.551
[1] "Finished preliminary run"
[1] "Starting runs"
[1] "Entered MLica"
[1] "Entering iteration loop " "0"
[1] "Checkpt1"
[1] "AvDev=" "0.311074578137404"
[1] "iterated logL"
[1] -4035.497
[1] "Entering iteration loop " "1"
[1] "Checkpt1"
[1] "AvDev=" "0.0958933207246391"
[1] "iterated logL"
[1] -3661.894
[1] "Entering iteration loop " "2"
[1] "Checkpt1"
[1] "AvDev=" "0.0637237256445566"
[1] "iterated logL"
[1] -3608.866
[1] "Entering iteration loop " "3"
[1] "Checkpt1"
[1] "AvDev=" "0.0664723119203546"
[1] "iterated logL"
[1] -3499.507
[1] "Entering iteration loop " "4"
[1] "Checkpt1"
[1] "AvDev=" "0.0524919873563722"
[1] "iterated logL"
[1] -3423.849
[1] "Entering iteration loop " "5"
[1] "Checkpt1"
[1] "AvDev=" "0.0266110652585197"
[1] "iterated logL"
[1] -3399.491
[1] "Entering iteration loop " "6"
[1] "Checkpt1"
[1] "AvDev=" "0.00948270616614102"
[1] "iterated logL"
[1] -3396.113
[1] "Entering iteration loop " "7"
[1] "Checkpt1"
[1] "AvDev=" "0.00236905376963106"
[1] "iterated logL"
[1] -3395.669
[1] "Entering iteration loop " "8"
[1] "Checkpt1"
[1] "AvDev=" "0.00053187290600404"
[1] "iterated logL"
[1] -3395.579
[1] "Entering iteration loop " "9"
[1] "Checkpt1"
[1] "AvDev=" "0.000127236444056749"
[1] "iterated logL"
[1] -3395.558
[1] "Entering iteration loop " "10"
[1] "Checkpt1"
[1] "AvDev=" "3.03095178792535e-05"
[1] "iterated logL"
[1] -3395.553
[1] "Entering iteration loop " "11"
[1] "Checkpt1"
[1] "AvDev=" "7.27387684576024e-06"
[1] "iterated logL"
[1] -3395.552
[1] "Entered MLica"
[1] "Entering iteration loop " "0"
[1] "Checkpt1"
[1] "AvDev=" "0.408283532394995"
[1] "iterated logL"
[1] -4144.323
[1] "Entering iteration loop " "1"
[1] "Checkpt1"
[1] "AvDev=" "0.131991485639095"
[1] "iterated logL"
[1] -3479.347
[1] "Entering iteration loop " "2"
[1] "Checkpt1"
[1] "AvDev=" "0.0481662412793996"
[1] "iterated logL"
[1] -3443.261
[1] "Entering iteration loop " "3"
[1] "Checkpt1"
[1] "AvDev=" "0.0410388259598368"
[1] "iterated logL"
[1] -3401.783
[1] "Entering iteration loop " "4"
[1] "Checkpt1"
[1] "AvDev=" "0.00960736728874848"
[1] "iterated logL"
[1] -3395.308
[1] "Entering iteration loop " "5"
[1] "Checkpt1"
[1] "AvDev=" "0.00159036836987607"
[1] "iterated logL"
[1] -3395.651
[1] "Entering iteration loop " "6"
[1] "Checkpt1"
[1] "AvDev=" "0.00030326441425986"
[1] "iterated logL"
[1] -3395.55
[1] "Entering iteration loop " "7"
[1] "Checkpt1"
[1] "AvDev=" "5.75657634615365e-05"
[1] "iterated logL"
[1] -3395.555
[1] "Entering iteration loop " "8"
[1] "Checkpt1"
[1] "AvDev=" "1.37327379236027e-05"
[1] "iterated logL"
[1] -3395.552
[1] "Entered MLica"
[1] "Entering iteration loop " "0"
[1] "Checkpt1"
[1] "AvDev=" "0.425815168841733"
[1] "iterated logL"
[1] -3667.934
[1] "Entering iteration loop " "1"
[1] "Checkpt1"
[1] "AvDev=" "0.0698379808222559"
[1] "iterated logL"
[1] -3524.008
[1] "Entering iteration loop " "2"
[1] "Checkpt1"
[1] "AvDev=" "0.0131519647708972"
[1] "iterated logL"
[1] -3509.167
[1] "Entering iteration loop " "3"
[1] "Checkpt1"
[1] "AvDev=" "0.00486706941295416"
[1] "iterated logL"
[1] -3507.437
[1] "Entering iteration loop " "4"
[1] "Checkpt1"
[1] "AvDev=" "0.00625192454533561"
[1] "iterated logL"
[1] -3509.155
[1] "Entering iteration loop " "5"
[1] "Checkpt1"
[1] "AvDev=" "0.0115161072696129"
[1] "iterated logL"
[1] -3510.774
[1] "Entering iteration loop " "6"
[1] "Checkpt1"
[1] "AvDev=" "0.0207314015437168"
[1] "iterated logL"
[1] -3506.719
[1] "Entering iteration loop " "7"
[1] "Checkpt1"
[1] "AvDev=" "0.0316288909278083"
[1] "iterated logL"
[1] -3479.767
[1] "Entering iteration loop " "8"
[1] "Checkpt1"
[1] "AvDev=" "0.0325973768372488"
[1] "iterated logL"
[1] -3425.318
[1] "Entering iteration loop " "9"
[1] "Checkpt1"
[1] "AvDev=" "0.018965538980085"
[1] "iterated logL"
[1] -3399.135
[1] "Entering iteration loop " "10"
[1] "Checkpt1"
[1] "AvDev=" "0.00534623296704367"
[1] "iterated logL"
[1] -3395.796
[1] "Entering iteration loop " "11"
[1] "Checkpt1"
[1] "AvDev=" "0.00113679433105254"
[1] "iterated logL"
[1] -3395.569
[1] "Entering iteration loop " "12"
[1] "Checkpt1"
[1] "AvDev=" "0.000209041273698903"
[1] "iterated logL"
[1] -3395.556
[1] "Entering iteration loop " "13"
[1] "Checkpt1"
[1] "AvDev=" "3.82955775854774e-05"
[1] "iterated logL"
[1] -3395.553
[1] "Entering iteration loop " "14"
[1] "Checkpt1"
[1] "AvDev=" "7.38487073760532e-06"
[1] "iterated logL"
[1] -3395.552
[1] "Entered MLica"
[1] "Entering iteration loop " "0"
[1] "Checkpt1"
[1] "AvDev=" "0.371359542049407"
[1] "iterated logL"
[1] -3876.2
[1] "Entering iteration loop " "1"
[1] "Checkpt1"
[1] "AvDev=" "0.142530946467836"
[1] "iterated logL"
[1] -3441.852
[1] "Entering iteration loop " "2"
[1] "Checkpt1"
[1] "AvDev=" "0.0334766990223169"
[1] "iterated logL"
[1] -3395.664
[1] "Entering iteration loop " "3"
[1] "Checkpt1"
[1] "AvDev=" "0.00382891203440433"
[1] "iterated logL"
[1] -3395.684
[1] "Entering iteration loop " "4"
[1] "Checkpt1"
[1] "AvDev=" "0.000594310499221397"
[1] "iterated logL"
[1] -3395.562
[1] "Entering iteration loop " "5"
[1] "Checkpt1"
[1] "AvDev=" "0.000116149447185206"
[1] "iterated logL"
[1] -3395.558
[1] "Entering iteration loop " "6"
[1] "Checkpt1"
[1] "AvDev=" "2.83853526208703e-05"
[1] "iterated logL"
[1] -3395.553
[1] "Entering iteration loop " "7"
[1] "Checkpt1"
[1] "AvDev=" "6.14471227920681e-06"
[1] "iterated logL"
[1] -3395.552
[1] "Entered MLica"
[1] "Entering iteration loop " "0"
[1] "Checkpt1"
[1] "AvDev=" "0.410470544121206"
[1] "iterated logL"
[1] -3761.94
[1] "Entering iteration loop " "1"
[1] "Checkpt1"
[1] "AvDev=" "0.0894245649461816"
[1] "iterated logL"
[1] -3558.636
[1] "Entering iteration loop " "2"
[1] "Checkpt1"
[1] "AvDev=" "0.0295823341135455"
[1] "iterated logL"
[1] -3497.191
[1] "Entering iteration loop " "3"
[1] "Checkpt1"
[1] "AvDev=" "0.0319827688603013"
[1] "iterated logL"
[1] -3458.716
[1] "Entering iteration loop " "4"
[1] "Checkpt1"
[1] "AvDev=" "0.0458463230582255"
[1] "iterated logL"
[1] -3414.967
[1] "Entering iteration loop " "5"
[1] "Checkpt1"
[1] "AvDev=" "0.0204951937399991"
[1] "iterated logL"
[1] -3398.147
[1] "Entering iteration loop " "6"
[1] "Checkpt1"
[1] "AvDev=" "0.00345410150464575"
[1] "iterated logL"
[1] -3395.399
[1] "Entering iteration loop " "7"
[1] "Checkpt1"
[1] "AvDev=" "0.000624447866132663"
[1] "iterated logL"
[1] -3395.583
[1] "Entering iteration loop " "8"
[1] "Checkpt1"
[1] "AvDev=" "0.000111555927307654"
[1] "iterated logL"
[1] -3395.548
[1] "Entering iteration loop " "9"
[1] "Checkpt1"
[1] "AvDev=" "2.68388826026896e-05"
[1] "iterated logL"
[1] -3395.552
[1] "Entering iteration loop " "10"
[1] "Checkpt1"
[1] "AvDev=" "5.42100137719388e-06"
[1] "iterated logL"
[1] -3395.552
[1] "End of runs"
[1] "Entering mlica"
[1] "Performing preliminary run"
[1] "Entered MLica"
[1] "Entering iteration loop " "0"
[1] "Checkpt1"
[1] "AvDev=" "0.366074296551384"
[1] "iterated logL"
[1] -4159.965
[1] "Entering iteration loop " "1"
[1] "Checkpt1"
[1] "AvDev=" "0.0990706074759644"
[1] "iterated logL"
[1] -3842.361
[1] "Entering iteration loop " "2"
[1] "Checkpt1"
[1] "AvDev=" "0.0925390792823253"
[1] "iterated logL"
[1] -3492.819
[1] "Entering iteration loop " "3"
[1] "Checkpt1"
[1] "AvDev=" "0.0526223128722005"
[1] "iterated logL"
[1] -3403.373
[1] "Entering iteration loop " "4"
[1] "Checkpt1"
[1] "AvDev=" "0.0181231017247815"
[1] "iterated logL"
[1] -3396.83
[1] "Entering iteration loop " "5"
[1] "Checkpt1"
[1] "AvDev=" "0.00445475815844939"
[1] "iterated logL"
[1] -3395.756
[1] "Entering iteration loop " "6"
[1] "Checkpt1"
[1] "AvDev=" "0.000901564622348875"
[1] "iterated logL"
[1] -3395.602
[1] "Entering iteration loop " "7"
[1] "Checkpt1"
[1] "AvDev=" "0.000213306212586903"
[1] "iterated logL"
[1] -3395.562
[1] "Entering iteration loop " "8"
[1] "Checkpt1"
[1] "AvDev=" "4.75160900748634e-05"
[1] "iterated logL"
[1] -3395.555
[1] "Entering iteration loop " "9"
[1] "Checkpt1"
[1] "AvDev=" "1.16181179872218e-05"
[1] "iterated logL"
[1] -3395.552
[1] "Finished preliminary run"
[1] "Starting runs"
[1] "Entered MLica"
[1] "Entering iteration loop " "0"
[1] "Checkpt1"
[1] "AvDev=" "0.438326263367479"
[1] "iterated logL"
[1] -4019.211
[1] "Entering iteration loop " "1"
[1] "Checkpt1"
[1] "AvDev=" "0.0999407970722891"
[1] "iterated logL"
[1] -3937.154
[1] "Entering iteration loop " "2"
[1] "Checkpt1"
[1] "AvDev=" "0.0857379663920665"
[1] "iterated logL"
[1] -3583.106
[1] "Entering iteration loop " "3"
[1] "Checkpt1"
[1] "AvDev=" "0.0749152696209581"
[1] "iterated logL"
[1] -3433.51
[1] "Entering iteration loop " "4"
[1] "Checkpt1"
[1] "AvDev=" "0.0259885513192216"
[1] "iterated logL"
[1] -3400.524
[1] "Entering iteration loop " "5"
[1] "Checkpt1"
[1] "AvDev=" "0.00445187199563686"
[1] "iterated logL"
[1] -3395.51
[1] "Entering iteration loop " "6"
[1] "Checkpt1"
[1] "AvDev=" "0.000831511651645693"
[1] "iterated logL"
[1] -3395.625
[1] "Entering iteration loop " "7"
[1] "Checkpt1"
[1] "AvDev=" "0.000192552944613526"
[1] "iterated logL"
[1] -3395.55
[1] "Entering iteration loop " "8"
[1] "Checkpt1"
[1] "AvDev=" "3.89059018286901e-05"
[1] "iterated logL"
[1] -3395.555
[1] "Entering iteration loop " "9"
[1] "Checkpt1"
[1] "AvDev=" "9.37276456141046e-06"
[1] "iterated logL"
[1] -3395.552
[1] "Entered MLica"
[1] "Entering iteration loop " "0"
[1] "Checkpt1"
[1] "AvDev=" "0.351009643350607"
[1] "iterated logL"
[1] -3735.173
[1] "Entering iteration loop " "1"
[1] "Checkpt1"
[1] "AvDev=" "0.109323958491958"
[1] "iterated logL"
[1] -3544.713
[1] "Entering iteration loop " "2"
[1] "Checkpt1"
[1] "AvDev=" "0.0707197937973701"
[1] "iterated logL"
[1] -3486.343
[1] "Entering iteration loop " "3"
[1] "Checkpt1"
[1] "AvDev=" "0.04351262542631"
[1] "iterated logL"
[1] -3414.097
[1] "Entering iteration loop " "4"
[1] "Checkpt1"
[1] "AvDev=" "0.0146397623893175"
[1] "iterated logL"
[1] -3396.041
[1] "Entering iteration loop " "5"
[1] "Checkpt1"
[1] "AvDev=" "0.00178476371045314"
[1] "iterated logL"
[1] -3395.552
[1] "Entering iteration loop " "6"
[1] "Checkpt1"
[1] "AvDev=" "0.000333464059947828"
[1] "iterated logL"
[1] -3395.562
[1] "Entering iteration loop " "7"
[1] "Checkpt1"
[1] "AvDev=" "7.88913069801145e-05"
[1] "iterated logL"
[1] -3395.554
[1] "Entering iteration loop " "8"
[1] "Checkpt1"
[1] "AvDev=" "1.78322360740409e-05"
[1] "iterated logL"
[1] -3395.552
[1] "Entered MLica"
[1] "Entering iteration loop " "0"
[1] "Checkpt1"
[1] "AvDev=" "0.49119825887111"
[1] "iterated logL"
[1] -4801.95
[1] "Entering iteration loop " "1"
[1] "Checkpt1"
[1] "AvDev=" "0.141448069534506"
[1] "iterated logL"
[1] -3756.975
[1] "Entering iteration loop " "2"
[1] "Checkpt1"
[1] "AvDev=" "0.113921181520134"
[1] "iterated logL"
[1] -3445.1
[1] "Entering iteration loop " "3"
[1] "Checkpt1"
[1] "AvDev=" "0.0451228348785501"
[1] "iterated logL"
[1] -3400.834
[1] "Entering iteration loop " "4"
[1] "Checkpt1"
[1] "AvDev=" "0.0120801507087052"
[1] "iterated logL"
[1] -3396.434
[1] "Entering iteration loop " "5"
[1] "Checkpt1"
[1] "AvDev=" "0.00329053928203262"
[1] "iterated logL"
[1] -3395.715
[1] "Entering iteration loop " "6"
[1] "Checkpt1"
[1] "AvDev=" "0.000666767357514407"
[1] "iterated logL"
[1] -3395.592
[1] "Entering iteration loop " "7"
[1] "Checkpt1"
[1] "AvDev=" "0.000162347881609018"
[1] "iterated logL"
[1] -3395.56
[1] "Entering iteration loop " "8"
[1] "Checkpt1"
[1] "AvDev=" "3.73592487872754e-05"
[1] "iterated logL"
[1] -3395.554
[1] "Entering iteration loop " "9"
[1] "Checkpt1"
[1] "AvDev=" "9.07690874152711e-06"
[1] "iterated logL"
[1] -3395.552
[1] "Entered MLica"
[1] "Entering iteration loop " "0"
[1] "Checkpt1"
[1] "AvDev=" "0.36746082558709"
[1] "iterated logL"
[1] -4123.382
[1] "Entering iteration loop " "1"
[1] "Checkpt1"
[1] "AvDev=" "0.13094735679392"
[1] "iterated logL"
[1] -3518.547
[1] "Entering iteration loop " "2"
[1] "Checkpt1"
[1] "AvDev=" "0.0521414764338787"
[1] "iterated logL"
[1] -3423.406
[1] "Entering iteration loop " "3"
[1] "Checkpt1"
[1] "AvDev=" "0.0150553718379829"
[1] "iterated logL"
[1] -3398.158
[1] "Entering iteration loop " "4"
[1] "Checkpt1"
[1] "AvDev=" "0.00365878582864727"
[1] "iterated logL"
[1] -3395.7
[1] "Entering iteration loop " "5"
[1] "Checkpt1"
[1] "AvDev=" "0.000857647373322417"
[1] "iterated logL"
[1] -3395.565
[1] "Entering iteration loop " "6"
[1] "Checkpt1"
[1] "AvDev=" "0.000194209693258142"
[1] "iterated logL"
[1] -3395.553
[1] "Entering iteration loop " "7"
[1] "Checkpt1"
[1] "AvDev=" "4.500676380397e-05"
[1] "iterated logL"
[1] -3395.552
[1] "Entering iteration loop " "8"
[1] "Checkpt1"
[1] "AvDev=" "1.02524134500583e-05"
[1] "iterated logL"
[1] -3395.552
[1] "Entered MLica"
[1] "Entering iteration loop " "0"
[1] "Checkpt1"
[1] "AvDev=" "0.381175722621573"
[1] "iterated logL"
[1] -4035.78
[1] "Entering iteration loop " "1"
[1] "Checkpt1"
[1] "AvDev=" "0.0828709227554443"
[1] "iterated logL"
[1] -3845.233
[1] "Entering iteration loop " "2"
[1] "Checkpt1"
[1] "AvDev=" "0.0935145589689355"
[1] "iterated logL"
[1] -3694.554
[1] "Entering iteration loop " "3"
[1] "Checkpt1"
[1] "AvDev=" "0.0844200429745166"
[1] "iterated logL"
[1] -3608.252
[1] "Entering iteration loop " "4"
[1] "Checkpt1"
[1] "AvDev=" "0.0663941351087329"
[1] "iterated logL"
[1] -3478.698
[1] "Entering iteration loop " "5"
[1] "Checkpt1"
[1] "AvDev=" "0.0446549747417372"
[1] "iterated logL"
[1] -3406.252
[1] "Entering iteration loop " "6"
[1] "Checkpt1"
[1] "AvDev=" "0.0140499598598769"
[1] "iterated logL"
[1] -3395.556
[1] "Entering iteration loop " "7"
[1] "Checkpt1"
[1] "AvDev=" "0.00392616751329586"
[1] "iterated logL"
[1] -3395.332
[1] "Entering iteration loop " "8"
[1] "Checkpt1"
[1] "AvDev=" "0.000990970829615874"
[1] "iterated logL"
[1] -3395.499
[1] "Entering iteration loop " "9"
[1] "Checkpt1"
[1] "AvDev=" "0.000237840071878959"
[1] "iterated logL"
[1] -3395.538
[1] "Entering iteration loop " "10"
[1] "Checkpt1"
[1] "AvDev=" "5.74841248261739e-05"
[1] "iterated logL"
[1] -3395.549
[1] "Entering iteration loop " "11"
[1] "Checkpt1"
[1] "AvDev=" "1.36830985332594e-05"
[1] "iterated logL"
[1] -3395.551
[1] "End of runs"
There were 50 or more warnings (use warnings() to see the first 50)
Warning message:
In solve.default(STS, LINPACK = FALSE) :
the LINPACK argument has been defunct since R 3.1.0
1 2 3 4 5
[1,] 1.707899e-02 0.006944718 -0.041919257 -0.996658231 -0.020187865
[2,] -1.244749e-02 0.045893055 -0.030323965 0.039829080 0.998275626
[3,] 7.363413e-05 0.032804710 0.998660069 -0.054682965 0.045595331
[4,] -4.938957e-02 -0.996255792 -0.000171045 -0.006949881 0.030553821
[5,] -9.985560e-01 0.065144815 0.001126950 0.045249163 -0.005074941
function (prNCP, nruns = 10, tol = 1e-04, maxit = 300, fail.th = 5,
learn.mu = 1)
{
print("Entering mlica")
print("Performing preliminary run")
a <- mlicaMAIN(prNCP, tol = 1e-04, maxit = 10, mu = learn.mu)
ncp <- dim(a$S)[2]
max.logL <- a$LL
a.best <- a
print("Finished preliminary run")
print("Starting runs")
run.n <- 0
fail.count <- 0
v.logL <- vector()
v.NC <- vector()
while (run.n < nruns) {
a <- mlicaMAIN(prNCP, tol = 1e-04, maxit = maxit, mu = learn.mu)
v.logL <- c(v.logL, a$LL)
v.NC <- c(v.NC, a$NC)
if (a$NC == 0) {
fail.count <- 0
run.n <- run.n + 1
if (a$LL > max.logL) {
a.best <- a
}
}
else {
fail.count <- fail.count + 1
}
if (fail.count >= fail.th) {
print("Stopping: Five consecutive runs failed to converge!")
print("Consider either increasing the threshold for pca eigenvalues to perform ICA on a smaller subspace or increasing maxit")
stop
}
}
print("End of runs")
return(a.best)
}
mlica2 documentation built on May 1, 2019, 10:45 p.m.
R Package Documentation
Browse R Packages
We want your feedback!
Note that we can't provide technical support on individual packages. You should contact the package authors for that.
Description Usage Arguments Value Author(s) References Examples
Description
This function performs ICA using a maximum likelihood framework and takes as arguments parameters to control the number of algorithm runs and convergence criteria.
Usage
1 | mlica(prNCP, nruns = 10, tol = 1e-04, maxit = 300, fail.th = 5, learn.mu = 1)
|
Arguments
prNCP |
The output object from 'proposeNCP'. |
nruns |
The number of converged algorithm runs sought (function returns the best solution according to the log-likelihood value). |
tol |
Tolerance level for establishing convergence of run. |
maxit |
Maximum number of iterations to allow per run. |
fail.th |
A threshold on the number of consecutive runs that fail to converge. |
learn.mu |
Learning parameter for fixed point algorithm (note that this need not be changed since it has already been optimised). |
Value
A list with following components:
A: Estimate of the mixing matrix.
B: Estimate of the inverse mixing matrix.
S: Estimate of the source matrix.
X: Normalised data matrix.
ncp: Number of independent components.
NC: Binary number specifying whether best run converged or not.(=1 indicates convergence,=0 indicates no convergence).
LL: Log likelihood value of best run.
Author(s)
Andrew Teschendorff a.teschendorff@ucl.ac.uk
References
Hyvaerinen A., Karhunen J., and Oja E.: Independent Component Analysis, John Wiley and Sons, New York, (2001).
Kreil D. and MacKay D. (2003): Reproducibility Assessment of Independent Component Analysis of Expression Ratios from DNA microarrays, Comparative and Functional Genomics *4* (3),300-317.
Liebermeister W. (2002): Linear Modes of gene expression determined by independent component analysis, Bioinformatics *18*, no.1, 51-60.
Chiappetta P., Roubaud MC. and Torresani B.: Blind source separation and the analysis of microarray data, J. Comput. Biol. 2004; 11(6):1090-109.
Examples
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 | ## Not run:
data(simMAdata);
dataX <- simMAdata[[1]];
prPCA <- PriorNormPCA(dataX);
prNCP <- proposeNCP(prPCA,0.1);
a.best.l <- list();
for( i in 1:5){
a.best.l[[i]] <- mlica(prNCP,nruns=5);
}
checkICA <- CheckStability(a.best.l,0.7);
sourceS <- simMAdata[[3]];
print(cor(a.best.l[[1]]$S,sourceS));
sModes <- SortModes(a.best.l[[1]],c.val=0.5);
## End(Not run)
## The function is currently defined as
function (prNCP, nruns = 10, tol = 1e-04, maxit = 300, fail.th = 5,
learn.mu = 1)
{
print("Entering mlica")
print("Performing preliminary run")
a <- mlicaMAIN(prNCP, tol = 1e-04, maxit = 10, mu = learn.mu)
ncp <- dim(a$S)[2]
max.logL <- a$LL
a.best <- a
print("Finished preliminary run")
print("Starting runs")
run.n <- 0
fail.count <- 0
v.logL <- vector()
v.NC <- vector()
while (run.n < nruns) {
a <- mlicaMAIN(prNCP, tol = 1e-04, maxit = maxit, mu = learn.mu)
v.logL <- c(v.logL, a$LL)
v.NC <- c(v.NC, a$NC)
if (a$NC == 0) {
fail.count <- 0
run.n <- run.n + 1
if (a$LL > max.logL) {
a.best <- a
}
}
else {
fail.count <- fail.count + 1
}
if (fail.count >= fail.th) {
print("Stopping: Five consecutive runs failed to converge!")
print("Consider either increasing the threshold for pca eigenvalues to perform ICA on a smaller subspace or increasing maxit")
stop
}
}
print("End of runs")
return(a.best)
}
|
Example output
[1] "Performing SVD"
Warning message:
In svd(X, LINPACK = FALSE) :
the LINPACK argument has been defunct since R 3.1.0
[1] "About to find ncp"
[1] "Entering mlica"
[1] "Performing preliminary run"
[1] "Entered MLica"
[1] "Entering iteration loop " "0"
[1] "Checkpt1"
[1] "AvDev=" "0.290219314280639"
[1] "iterated logL"
[1] -3590.578
[1] "Entering iteration loop " "1"
[1] "Checkpt1"
[1] "AvDev=" "0.105827278122874"
[1] "iterated logL"
[1] -3472.839
[1] "Entering iteration loop " "2"
[1] "Checkpt1"
[1] "AvDev=" "0.0471382644915915"
[1] "iterated logL"
[1] -3421.788
[1] "Entering iteration loop " "3"
[1] "Checkpt1"
[1] "AvDev=" "0.0264127579726401"
[1] "iterated logL"
[1] -3401.135
[1] "Entering iteration loop " "4"
[1] "Checkpt1"
[1] "AvDev=" "0.0142316603037475"
[1] "iterated logL"
[1] -3395.756
[1] "Entering iteration loop " "5"
[1] "Checkpt1"
[1] "AvDev=" "0.0048090862414591"
[1] "iterated logL"
[1] -3395.325
[1] "Entering iteration loop " "6"
[1] "Checkpt1"
[1] "AvDev=" "0.00127812851800342"
[1] "iterated logL"
[1] -3395.483
[1] "Entering iteration loop " "7"
[1] "Checkpt1"
[1] "AvDev=" "0.000311215663334201"
[1] "iterated logL"
[1] -3395.535
[1] "Entering iteration loop " "8"
[1] "Checkpt1"
[1] "AvDev=" "7.47703362656165e-05"
[1] "iterated logL"
[1] -3395.548
[1] "Entering iteration loop " "9"
[1] "Checkpt1"
[1] "AvDev=" "1.78819501199913e-05"
[1] "iterated logL"
[1] -3395.551
[1] "Finished preliminary run"
[1] "Starting runs"
[1] "Entered MLica"
[1] "Entering iteration loop " "0"
[1] "Checkpt1"
[1] "AvDev=" "0.400979712633759"
[1] "iterated logL"
[1] -4423.601
[1] "Entering iteration loop " "1"
[1] "Checkpt1"
[1] "AvDev=" "0.179739366038515"
[1] "iterated logL"
[1] -3473.989
[1] "Entering iteration loop " "2"
[1] "Checkpt1"
[1] "AvDev=" "0.0444486673678738"
[1] "iterated logL"
[1] -3398.851
[1] "Entering iteration loop " "3"
[1] "Checkpt1"
[1] "AvDev=" "0.00487559331142786"
[1] "iterated logL"
[1] -3395.584
[1] "Entering iteration loop " "4"
[1] "Checkpt1"
[1] "AvDev=" "0.000858328754907789"
[1] "iterated logL"
[1] -3395.607
[1] "Entering iteration loop " "5"
[1] "Checkpt1"
[1] "AvDev=" "0.000207039556067208"
[1] "iterated logL"
[1] -3395.554
[1] "Entering iteration loop " "6"
[1] "Checkpt1"
[1] "AvDev=" "4.13203110180985e-05"
[1] "iterated logL"
[1] -3395.554
[1] "Entering iteration loop " "7"
[1] "Checkpt1"
[1] "AvDev=" "1.05896434797988e-05"
[1] "iterated logL"
[1] -3395.552
[1] "Entered MLica"
[1] "Entering iteration loop " "0"
[1] "Checkpt1"
[1] "AvDev=" "0.463498076168755"
[1] "iterated logL"
[1] -3667.793
[1] "Entering iteration loop " "1"
[1] "Checkpt1"
[1] "AvDev=" "0.0764780178146065"
[1] "iterated logL"
[1] -3510.967
[1] "Entering iteration loop " "2"
[1] "Checkpt1"
[1] "AvDev=" "0.0842512230811305"
[1] "iterated logL"
[1] -3424.638
[1] "Entering iteration loop " "3"
[1] "Checkpt1"
[1] "AvDev=" "0.0322542264738013"
[1] "iterated logL"
[1] -3397.497
[1] "Entering iteration loop " "4"
[1] "Checkpt1"
[1] "AvDev=" "0.00461314705312814"
[1] "iterated logL"
[1] -3395.5
[1] "Entering iteration loop " "5"
[1] "Checkpt1"
[1] "AvDev=" "0.00065829292806568"
[1] "iterated logL"
[1] -3395.574
[1] "Entering iteration loop " "6"
[1] "Checkpt1"
[1] "AvDev=" "0.000113873449825682"
[1] "iterated logL"
[1] -3395.549
[1] "Entering iteration loop " "7"
[1] "Checkpt1"
[1] "AvDev=" "2.12979921296829e-05"
[1] "iterated logL"
[1] -3395.553
[1] "Entering iteration loop " "8"
[1] "Checkpt1"
[1] "AvDev=" "3.54673673779438e-06"
[1] "iterated logL"
[1] -3395.552
[1] "Entered MLica"
[1] "Entering iteration loop " "0"
[1] "Checkpt1"
[1] "AvDev=" "0.36629807731381"
[1] "iterated logL"
[1] -3992.956
[1] "Entering iteration loop " "1"
[1] "Checkpt1"
[1] "AvDev=" "0.100107945866601"
[1] "iterated logL"
[1] -3800.148
[1] "Entering iteration loop " "2"
[1] "Checkpt1"
[1] "AvDev=" "0.0747871283230001"
[1] "iterated logL"
[1] -3577.054
[1] "Entering iteration loop " "3"
[1] "Checkpt1"
[1] "AvDev=" "0.0508704848282896"
[1] "iterated logL"
[1] -3478.123
[1] "Entering iteration loop " "4"
[1] "Checkpt1"
[1] "AvDev=" "0.0516390652426347"
[1] "iterated logL"
[1] -3418.488
[1] "Entering iteration loop " "5"
[1] "Checkpt1"
[1] "AvDev=" "0.0270751636990737"
[1] "iterated logL"
[1] -3399.461
[1] "Entering iteration loop " "6"
[1] "Checkpt1"
[1] "AvDev=" "0.00507287047873504"
[1] "iterated logL"
[1] -3395.689
[1] "Entering iteration loop " "7"
[1] "Checkpt1"
[1] "AvDev=" "0.000827287456690507"
[1] "iterated logL"
[1] -3395.541
[1] "Entering iteration loop " "8"
[1] "Checkpt1"
[1] "AvDev=" "0.000156985449468821"
[1] "iterated logL"
[1] -3395.549
[1] "Entering iteration loop " "9"
[1] "Checkpt1"
[1] "AvDev=" "3.27897912097372e-05"
[1] "iterated logL"
[1] -3395.551
[1] "Entering iteration loop " "10"
[1] "Checkpt1"
[1] "AvDev=" "7.02164083348528e-06"
[1] "iterated logL"
[1] -3395.552
[1] "Entered MLica"
[1] "Entering iteration loop " "0"
[1] "Checkpt1"
[1] "AvDev=" "0.339996295702037"
[1] "iterated logL"
[1] -3734.831
[1] "Entering iteration loop " "1"
[1] "Checkpt1"
[1] "AvDev=" "0.115184530320527"
[1] "iterated logL"
[1] -3664.789
[1] "Entering iteration loop " "2"
[1] "Checkpt1"
[1] "AvDev=" "0.0761931714957583"
[1] "iterated logL"
[1] -3513.394
[1] "Entering iteration loop " "3"
[1] "Checkpt1"
[1] "AvDev=" "0.0395665341671519"
[1] "iterated logL"
[1] -3468.643
[1] "Entering iteration loop " "4"
[1] "Checkpt1"
[1] "AvDev=" "0.0329772639881574"
[1] "iterated logL"
[1] -3425.284
[1] "Entering iteration loop " "5"
[1] "Checkpt1"
[1] "AvDev=" "0.0131235875899943"
[1] "iterated logL"
[1] -3398.354
[1] "Entering iteration loop " "6"
[1] "Checkpt1"
[1] "AvDev=" "0.00317390696458019"
[1] "iterated logL"
[1] -3395.759
[1] "Entering iteration loop " "7"
[1] "Checkpt1"
[1] "AvDev=" "0.000725960885815795"
[1] "iterated logL"
[1] -3395.573
[1] "Entering iteration loop " "8"
[1] "Checkpt1"
[1] "AvDev=" "0.000163452040045647"
[1] "iterated logL"
[1] -3395.555
[1] "Entering iteration loop " "9"
[1] "Checkpt1"
[1] "AvDev=" "3.72448459157757e-05"
[1] "iterated logL"
[1] -3395.552
[1] "Entering iteration loop " "10"
[1] "Checkpt1"
[1] "AvDev=" "8.40711799153595e-06"
[1] "iterated logL"
[1] -3395.552
[1] "Entered MLica"
[1] "Entering iteration loop " "0"
[1] "Checkpt1"
[1] "AvDev=" "0.364821563279185"
[1] "iterated logL"
[1] -4328.342
[1] "Entering iteration loop " "1"
[1] "Checkpt1"
[1] "AvDev=" "0.115872175582269"
[1] "iterated logL"
[1] -3754.574
[1] "Entering iteration loop " "2"
[1] "Checkpt1"
[1] "AvDev=" "0.109689364039819"
[1] "iterated logL"
[1] -3431.945
[1] "Entering iteration loop " "3"
[1] "Checkpt1"
[1] "AvDev=" "0.0359576994516968"
[1] "iterated logL"
[1] -3395.668
[1] "Entering iteration loop " "4"
[1] "Checkpt1"
[1] "AvDev=" "0.00633261188333669"
[1] "iterated logL"
[1] -3395.861
[1] "Entering iteration loop " "5"
[1] "Checkpt1"
[1] "AvDev=" "0.00083586484581874"
[1] "iterated logL"
[1] -3395.485
[1] "Entering iteration loop " "6"
[1] "Checkpt1"
[1] "AvDev=" "0.000156792108413645"
[1] "iterated logL"
[1] -3395.561
[1] "Entering iteration loop " "7"
[1] "Checkpt1"
[1] "AvDev=" "2.85094680056924e-05"
[1] "iterated logL"
[1] -3395.549
[1] "Entering iteration loop " "8"
[1] "Checkpt1"
[1] "AvDev=" "5.68654754230662e-06"
[1] "iterated logL"
[1] -3395.552
[1] "End of runs"
[1] "Entering mlica"
[1] "Performing preliminary run"
[1] "Entered MLica"
[1] "Entering iteration loop " "0"
[1] "Checkpt1"
[1] "AvDev=" "0.386038204478516"
[1] "iterated logL"
[1] -4391.043
[1] "Entering iteration loop " "1"
[1] "Checkpt1"
[1] "AvDev=" "0.118722047951209"
[1] "iterated logL"
[1] -3787.804
[1] "Entering iteration loop " "2"
[1] "Checkpt1"
[1] "AvDev=" "0.1160421413945"
[1] "iterated logL"
[1] -3476.545
[1] "Entering iteration loop " "3"
[1] "Checkpt1"
[1] "AvDev=" "0.0606436361056195"
[1] "iterated logL"
[1] -3400.034
[1] "Entering iteration loop " "4"
[1] "Checkpt1"
[1] "AvDev=" "0.00988526357087198"
[1] "iterated logL"
[1] -3395.5
[1] "Entering iteration loop " "5"
[1] "Checkpt1"
[1] "AvDev=" "0.00112709094447862"
[1] "iterated logL"
[1] -3395.475
[1] "Entering iteration loop " "6"
[1] "Checkpt1"
[1] "AvDev=" "0.000248549371090569"
[1] "iterated logL"
[1] -3395.548
[1] "Entering iteration loop " "7"
[1] "Checkpt1"
[1] "AvDev=" "5.46497451578651e-05"
[1] "iterated logL"
[1] -3395.549
[1] "Entering iteration loop " "8"
[1] "Checkpt1"
[1] "AvDev=" "1.21802304175041e-05"
[1] "iterated logL"
[1] -3395.552
[1] "Finished preliminary run"
[1] "Starting runs"
[1] "Entered MLica"
[1] "Entering iteration loop " "0"
[1] "Checkpt1"
[1] "AvDev=" "0.356261911078639"
[1] "iterated logL"
[1] -3672.559
[1] "Entering iteration loop " "1"
[1] "Checkpt1"
[1] "AvDev=" "0.0883629996095928"
[1] "iterated logL"
[1] -3405.789
[1] "Entering iteration loop " "2"
[1] "Checkpt1"
[1] "AvDev=" "0.0115397758858259"
[1] "iterated logL"
[1] -3395.782
[1] "Entering iteration loop " "3"
[1] "Checkpt1"
[1] "AvDev=" "0.00159642887472422"
[1] "iterated logL"
[1] -3395.648
[1] "Entering iteration loop " "4"
[1] "Checkpt1"
[1] "AvDev=" "0.00036227500675863"
[1] "iterated logL"
[1] -3395.562
[1] "Entering iteration loop " "5"
[1] "Checkpt1"
[1] "AvDev=" "7.82784397283197e-05"
[1] "iterated logL"
[1] -3395.556
[1] "Entering iteration loop " "6"
[1] "Checkpt1"
[1] "AvDev=" "1.80340204448921e-05"
[1] "iterated logL"
[1] -3395.552
[1] "Entered MLica"
[1] "Entering iteration loop " "0"
[1] "Checkpt1"
[1] "AvDev=" "0.38765322127385"
[1] "iterated logL"
[1] -3673.776
[1] "Entering iteration loop " "1"
[1] "Checkpt1"
[1] "AvDev=" "0.11495659664447"
[1] "iterated logL"
[1] -3472.173
[1] "Entering iteration loop " "2"
[1] "Checkpt1"
[1] "AvDev=" "0.0386829672942188"
[1] "iterated logL"
[1] -3405.763
[1] "Entering iteration loop " "3"
[1] "Checkpt1"
[1] "AvDev=" "0.00826699766804716"
[1] "iterated logL"
[1] -3396.022
[1] "Entering iteration loop " "4"
[1] "Checkpt1"
[1] "AvDev=" "0.00153531227230042"
[1] "iterated logL"
[1] -3395.556
[1] "Entering iteration loop " "5"
[1] "Checkpt1"
[1] "AvDev=" "0.000352005254312435"
[1] "iterated logL"
[1] -3395.548
[1] "Entering iteration loop " "6"
[1] "Checkpt1"
[1] "AvDev=" "7.73710107153997e-05"
[1] "iterated logL"
[1] -3395.551
[1] "Entering iteration loop " "7"
[1] "Checkpt1"
[1] "AvDev=" "1.78830343199412e-05"
[1] "iterated logL"
[1] -3395.552
[1] "Entered MLica"
[1] "Entering iteration loop " "0"
[1] "Checkpt1"
[1] "AvDev=" "0.392189272097064"
[1] "iterated logL"
[1] -5432.134
[1] "Entering iteration loop " "1"
[1] "Checkpt1"
[1] "AvDev=" "0.1417022380835"
[1] "iterated logL"
[1] -4041.275
[1] "Entering iteration loop " "2"
[1] "Checkpt1"
[1] "AvDev=" "0.0796615817576913"
[1] "iterated logL"
[1] -3755.534
[1] "Entering iteration loop " "3"
[1] "Checkpt1"
[1] "AvDev=" "0.0693391715870864"
[1] "iterated logL"
[1] -3506.783
[1] "Entering iteration loop " "4"
[1] "Checkpt1"
[1] "AvDev=" "0.0486440891664531"
[1] "iterated logL"
[1] -3424.152
[1] "Entering iteration loop " "5"
[1] "Checkpt1"
[1] "AvDev=" "0.0221261664494838"
[1] "iterated logL"
[1] -3397.661
[1] "Entering iteration loop " "6"
[1] "Checkpt1"
[1] "AvDev=" "0.0028671746044226"
[1] "iterated logL"
[1] -3395.543
[1] "Entering iteration loop " "7"
[1] "Checkpt1"
[1] "AvDev=" "0.000574550590411959"
[1] "iterated logL"
[1] -3395.546
[1] "Entering iteration loop " "8"
[1] "Checkpt1"
[1] "AvDev=" "0.000125694040543018"
[1] "iterated logL"
[1] -3395.548
[1] "Entering iteration loop " "9"
[1] "Checkpt1"
[1] "AvDev=" "2.68782845113823e-05"
[1] "iterated logL"
[1] -3395.551
[1] "Entering iteration loop " "10"
[1] "Checkpt1"
[1] "AvDev=" "5.96964466208338e-06"
[1] "iterated logL"
[1] -3395.552
[1] "Entered MLica"
[1] "Entering iteration loop " "0"
[1] "Checkpt1"
[1] "AvDev=" "0.410912985876807"
[1] "iterated logL"
[1] -4043.816
[1] "Entering iteration loop " "1"
[1] "Checkpt1"
[1] "AvDev=" "0.107261300936882"
[1] "iterated logL"
[1] -3911.253
[1] "Entering iteration loop " "2"
[1] "Checkpt1"
[1] "AvDev=" "0.0600342739500788"
[1] "iterated logL"
[1] -3916.637
[1] "Entering iteration loop " "3"
[1] "Checkpt1"
[1] "AvDev=" "0.0441894098933856"
[1] "iterated logL"
[1] -3758.087
[1] "Entering iteration loop " "4"
[1] "Checkpt1"
[1] "AvDev=" "0.0518400801703871"
[1] "iterated logL"
[1] -3561.38
[1] "Entering iteration loop " "5"
[1] "Checkpt1"
[1] "AvDev=" "0.0501583042138534"
[1] "iterated logL"
[1] -3413.922
[1] "Entering iteration loop " "6"
[1] "Checkpt1"
[1] "AvDev=" "0.0218865273541719"
[1] "iterated logL"
[1] -3396.987
[1] "Entering iteration loop " "7"
[1] "Checkpt1"
[1] "AvDev=" "0.00430632415238061"
[1] "iterated logL"
[1] -3395.452
[1] "Entering iteration loop " "8"
[1] "Checkpt1"
[1] "AvDev=" "0.000742641245699041"
[1] "iterated logL"
[1] -3395.533
[1] "Entering iteration loop " "9"
[1] "Checkpt1"
[1] "AvDev=" "0.000150363951070377"
[1] "iterated logL"
[1] -3395.545
[1] "Entering iteration loop " "10"
[1] "Checkpt1"
[1] "AvDev=" "3.29787956548537e-05"
[1] "iterated logL"
[1] -3395.55
[1] "Entering iteration loop " "11"
[1] "Checkpt1"
[1] "AvDev=" "7.32844518741213e-06"
[1] "iterated logL"
[1] -3395.551
[1] "Entered MLica"
[1] "Entering iteration loop " "0"
[1] "Checkpt1"
[1] "AvDev=" "0.371274126606981"
[1] "iterated logL"
[1] -3648.194
[1] "Entering iteration loop " "1"
[1] "Checkpt1"
[1] "AvDev=" "0.0801941640029815"
[1] "iterated logL"
[1] -3513.596
[1] "Entering iteration loop " "2"
[1] "Checkpt1"
[1] "AvDev=" "0.0744793555550997"
[1] "iterated logL"
[1] -3448.159
[1] "Entering iteration loop " "3"
[1] "Checkpt1"
[1] "AvDev=" "0.0525100410417017"
[1] "iterated logL"
[1] -3414.35
[1] "Entering iteration loop " "4"
[1] "Checkpt1"
[1] "AvDev=" "0.0283072099719519"
[1] "iterated logL"
[1] -3400.12
[1] "Entering iteration loop " "5"
[1] "Checkpt1"
[1] "AvDev=" "0.0129825402674149"
[1] "iterated logL"
[1] -3395.6
[1] "Entering iteration loop " "6"
[1] "Checkpt1"
[1] "AvDev=" "0.00424652217030174"
[1] "iterated logL"
[1] -3395.338
[1] "Entering iteration loop " "7"
[1] "Checkpt1"
[1] "AvDev=" "0.00111200174326132"
[1] "iterated logL"
[1] -3395.493
[1] "Entering iteration loop " "8"
[1] "Checkpt1"
[1] "AvDev=" "0.000269207446763617"
[1] "iterated logL"
[1] -3395.537
[1] "Entering iteration loop " "9"
[1] "Checkpt1"
[1] "AvDev=" "6.47633209478762e-05"
[1] "iterated logL"
[1] -3395.548
[1] "Entering iteration loop " "10"
[1] "Checkpt1"
[1] "AvDev=" "1.5395072232993e-05"
[1] "iterated logL"
[1] -3395.551
[1] "End of runs"
[1] "Entering mlica"
[1] "Performing preliminary run"
[1] "Entered MLica"
[1] "Entering iteration loop " "0"
[1] "Checkpt1"
[1] "AvDev=" "0.296028470334294"
[1] "iterated logL"
[1] -3809.782
[1] "Entering iteration loop " "1"
[1] "Checkpt1"
[1] "AvDev=" "0.0894398874126622"
[1] "iterated logL"
[1] -3594.638
[1] "Entering iteration loop " "2"
[1] "Checkpt1"
[1] "AvDev=" "0.0591590009051705"
[1] "iterated logL"
[1] -3459.552
[1] "Entering iteration loop " "3"
[1] "Checkpt1"
[1] "AvDev=" "0.0200699441239609"
[1] "iterated logL"
[1] -3439.458
[1] "Entering iteration loop " "4"
[1] "Checkpt1"
[1] "AvDev=" "0.0242748564489362"
[1] "iterated logL"
[1] -3429.975
[1] "Entering iteration loop " "5"
[1] "Checkpt1"
[1] "AvDev=" "0.0308904719620776"
[1] "iterated logL"
[1] -3411.681
[1] "Entering iteration loop " "6"
[1] "Checkpt1"
[1] "AvDev=" "0.0218622093012664"
[1] "iterated logL"
[1] -3399.13
[1] "Entering iteration loop " "7"
[1] "Checkpt1"
[1] "AvDev=" "0.00865223324385435"
[1] "iterated logL"
[1] -3396.15
[1] "Entering iteration loop " "8"
[1] "Checkpt1"
[1] "AvDev=" "0.00230748595510655"
[1] "iterated logL"
[1] -3395.681
[1] "Entering iteration loop " "9"
[1] "Checkpt1"
[1] "AvDev=" "0.000523761338507808"
[1] "iterated logL"
[1] -3395.582
[1] "Finished preliminary run"
[1] "Starting runs"
[1] "Entered MLica"
[1] "Entering iteration loop " "0"
[1] "Checkpt1"
[1] "AvDev=" "0.280196120957882"
[1] "iterated logL"
[1] -4126.692
[1] "Entering iteration loop " "1"
[1] "Checkpt1"
[1] "AvDev=" "0.136967268241482"
[1] "iterated logL"
[1] -3528.859
[1] "Entering iteration loop " "2"
[1] "Checkpt1"
[1] "AvDev=" "0.063909758039693"
[1] "iterated logL"
[1] -3418.22
[1] "Entering iteration loop " "3"
[1] "Checkpt1"
[1] "AvDev=" "0.0136820341896695"
[1] "iterated logL"
[1] -3396.547
[1] "Entering iteration loop " "4"
[1] "Checkpt1"
[1] "AvDev=" "0.00229592432874156"
[1] "iterated logL"
[1] -3395.697
[1] "Entering iteration loop " "5"
[1] "Checkpt1"
[1] "AvDev=" "0.000501495866691987"
[1] "iterated logL"
[1] -3395.569
[1] "Entering iteration loop " "6"
[1] "Checkpt1"
[1] "AvDev=" "0.000114241962480851"
[1] "iterated logL"
[1] -3395.558
[1] "Entering iteration loop " "7"
[1] "Checkpt1"
[1] "AvDev=" "2.56100609335057e-05"
[1] "iterated logL"
[1] -3395.553
[1] "Entered MLica"
[1] "Entering iteration loop " "0"
[1] "Checkpt1"
[1] "AvDev=" "0.3591373454329"
[1] "iterated logL"
[1] -3649.387
[1] "Entering iteration loop " "1"
[1] "Checkpt1"
[1] "AvDev=" "0.0870331784836628"
[1] "iterated logL"
[1] -3474.703
[1] "Entering iteration loop " "2"
[1] "Checkpt1"
[1] "AvDev=" "0.0318251790046466"
[1] "iterated logL"
[1] -3442.092
[1] "Entering iteration loop " "3"
[1] "Checkpt1"
[1] "AvDev=" "0.026413729802557"
[1] "iterated logL"
[1] -3428.089
[1] "Entering iteration loop " "4"
[1] "Checkpt1"
[1] "AvDev=" "0.0305909886580975"
[1] "iterated logL"
[1] -3409.182
[1] "Entering iteration loop " "5"
[1] "Checkpt1"
[1] "AvDev=" "0.0206336363315325"
[1] "iterated logL"
[1] -3397.634
[1] "Entering iteration loop " "6"
[1] "Checkpt1"
[1] "AvDev=" "0.0089562560436709"
[1] "iterated logL"
[1] -3395.333
[1] "Entering iteration loop " "7"
[1] "Checkpt1"
[1] "AvDev=" "0.00269255449333509"
[1] "iterated logL"
[1] -3395.407
[1] "Entering iteration loop " "8"
[1] "Checkpt1"
[1] "AvDev=" "0.000679790417210284"
[1] "iterated logL"
[1] -3395.516
[1] "Entering iteration loop " "9"
[1] "Checkpt1"
[1] "AvDev=" "0.000163404609437547"
[1] "iterated logL"
[1] -3395.543
[1] "Entering iteration loop " "10"
[1] "Checkpt1"
[1] "AvDev=" "3.9450123678408e-05"
[1] "iterated logL"
[1] -3395.55
[1] "Entering iteration loop " "11"
[1] "Checkpt1"
[1] "AvDev=" "9.41764292642783e-06"
[1] "iterated logL"
[1] -3395.551
[1] "Entered MLica"
[1] "Entering iteration loop " "0"
[1] "Checkpt1"
[1] "AvDev=" "0.375329432885969"
[1] "iterated logL"
[1] -4102.084
[1] "Entering iteration loop " "1"
[1] "Checkpt1"
[1] "AvDev=" "0.0914017749959568"
[1] "iterated logL"
[1] -3852.777
[1] "Entering iteration loop " "2"
[1] "Checkpt1"
[1] "AvDev=" "0.0819096172436832"
[1] "iterated logL"
[1] -3614.669
[1] "Entering iteration loop " "3"
[1] "Checkpt1"
[1] "AvDev=" "0.0640237802367762"
[1] "iterated logL"
[1] -3468.583
[1] "Entering iteration loop " "4"
[1] "Checkpt1"
[1] "AvDev=" "0.0287827456740873"
[1] "iterated logL"
[1] -3410.449
[1] "Entering iteration loop " "5"
[1] "Checkpt1"
[1] "AvDev=" "0.00832633461415859"
[1] "iterated logL"
[1] -3396.645
[1] "Entering iteration loop " "6"
[1] "Checkpt1"
[1] "AvDev=" "0.00187009118538909"
[1] "iterated logL"
[1] -3395.638
[1] "Entering iteration loop " "7"
[1] "Checkpt1"
[1] "AvDev=" "0.000413676242951251"
[1] "iterated logL"
[1] -3395.563
[1] "Entering iteration loop " "8"
[1] "Checkpt1"
[1] "AvDev=" "9.46708401711417e-05"
[1] "iterated logL"
[1] -3395.554
[1] "Entering iteration loop " "9"
[1] "Checkpt1"
[1] "AvDev=" "2.13034932697826e-05"
[1] "iterated logL"
[1] -3395.552
[1] "Entering iteration loop " "10"
[1] "Checkpt1"
[1] "AvDev=" "4.80830754306893e-06"
[1] "iterated logL"
[1] -3395.552
[1] "Entered MLica"
[1] "Entering iteration loop " "0"
[1] "Checkpt1"
[1] "AvDev=" "0.414533771616944"
[1] "iterated logL"
[1] -3619.764
[1] "Entering iteration loop " "1"
[1] "Checkpt1"
[1] "AvDev=" "0.0998829443199397"
[1] "iterated logL"
[1] -3428.315
[1] "Entering iteration loop " "2"
[1] "Checkpt1"
[1] "AvDev=" "0.0371035179271871"
[1] "iterated logL"
[1] -3395.115
[1] "Entering iteration loop " "3"
[1] "Checkpt1"
[1] "AvDev=" "0.00585473854860562"
[1] "iterated logL"
[1] -3395.705
[1] "Entering iteration loop " "4"
[1] "Checkpt1"
[1] "AvDev=" "0.00102639880166599"
[1] "iterated logL"
[1] -3395.46
[1] "Entering iteration loop " "5"
[1] "Checkpt1"
[1] "AvDev=" "0.000233825810315878"
[1] "iterated logL"
[1] -3395.553
[1] "Entering iteration loop " "6"
[1] "Checkpt1"
[1] "AvDev=" "4.63066237271436e-05"
[1] "iterated logL"
[1] -3395.548
[1] "Entering iteration loop " "7"
[1] "Checkpt1"
[1] "AvDev=" "1.13079987055212e-05"
[1] "iterated logL"
[1] -3395.552
[1] "Entered MLica"
[1] "Entering iteration loop " "0"
[1] "Checkpt1"
[1] "AvDev=" "0.326979966868656"
[1] "iterated logL"
[1] -3635.347
[1] "Entering iteration loop " "1"
[1] "Checkpt1"
[1] "AvDev=" "0.0958728068461986"
[1] "iterated logL"
[1] -3460.853
[1] "Entering iteration loop " "2"
[1] "Checkpt1"
[1] "AvDev=" "0.0571714978175811"
[1] "iterated logL"
[1] -3405.342
[1] "Entering iteration loop " "3"
[1] "Checkpt1"
[1] "AvDev=" "0.011807503859014"
[1] "iterated logL"
[1] -3395.566
[1] "Entering iteration loop " "4"
[1] "Checkpt1"
[1] "AvDev=" "0.00163925954499247"
[1] "iterated logL"
[1] -3395.508
[1] "Entering iteration loop " "5"
[1] "Checkpt1"
[1] "AvDev=" "0.000385225711069099"
[1] "iterated logL"
[1] -3395.538
[1] "Entering iteration loop " "6"
[1] "Checkpt1"
[1] "AvDev=" "9.01734293590305e-05"
[1] "iterated logL"
[1] -3395.549
[1] "Entering iteration loop " "7"
[1] "Checkpt1"
[1] "AvDev=" "2.12590118020316e-05"
[1] "iterated logL"
[1] -3395.551
[1] "Entering iteration loop " "8"
[1] "Checkpt1"
[1] "AvDev=" "5.08094001730523e-06"
[1] "iterated logL"
[1] -3395.552
[1] "End of runs"
[1] "Entering mlica"
[1] "Performing preliminary run"
[1] "Entered MLica"
[1] "Entering iteration loop " "0"
[1] "Checkpt1"
[1] "AvDev=" "0.438000141407069"
[1] "iterated logL"
[1] -4276.765
[1] "Entering iteration loop " "1"
[1] "Checkpt1"
[1] "AvDev=" "0.154084115234504"
[1] "iterated logL"
[1] -3499.516
[1] "Entering iteration loop " "2"
[1] "Checkpt1"
[1] "AvDev=" "0.0596482103812549"
[1] "iterated logL"
[1] -3397.822
[1] "Entering iteration loop " "3"
[1] "Checkpt1"
[1] "AvDev=" "0.00981867306151161"
[1] "iterated logL"
[1] -3395.482
[1] "Entering iteration loop " "4"
[1] "Checkpt1"
[1] "AvDev=" "0.00217984430913429"
[1] "iterated logL"
[1] -3395.449
[1] "Entering iteration loop " "5"
[1] "Checkpt1"
[1] "AvDev=" "0.000498978167806657"
[1] "iterated logL"
[1] -3395.529
[1] "Entering iteration loop " "6"
[1] "Checkpt1"
[1] "AvDev=" "0.000111394576278191"
[1] "iterated logL"
[1] -3395.546
[1] "Entering iteration loop " "7"
[1] "Checkpt1"
[1] "AvDev=" "2.70277136988498e-05"
[1] "iterated logL"
[1] -3395.55
[1] "Entering iteration loop " "8"
[1] "Checkpt1"
[1] "AvDev=" "6.12100051378475e-06"
[1] "iterated logL"
[1] -3395.551
[1] "Finished preliminary run"
[1] "Starting runs"
[1] "Entered MLica"
[1] "Entering iteration loop " "0"
[1] "Checkpt1"
[1] "AvDev=" "0.311074578137404"
[1] "iterated logL"
[1] -4035.497
[1] "Entering iteration loop " "1"
[1] "Checkpt1"
[1] "AvDev=" "0.0958933207246391"
[1] "iterated logL"
[1] -3661.894
[1] "Entering iteration loop " "2"
[1] "Checkpt1"
[1] "AvDev=" "0.0637237256445566"
[1] "iterated logL"
[1] -3608.866
[1] "Entering iteration loop " "3"
[1] "Checkpt1"
[1] "AvDev=" "0.0664723119203546"
[1] "iterated logL"
[1] -3499.507
[1] "Entering iteration loop " "4"
[1] "Checkpt1"
[1] "AvDev=" "0.0524919873563722"
[1] "iterated logL"
[1] -3423.849
[1] "Entering iteration loop " "5"
[1] "Checkpt1"
[1] "AvDev=" "0.0266110652585197"
[1] "iterated logL"
[1] -3399.491
[1] "Entering iteration loop " "6"
[1] "Checkpt1"
[1] "AvDev=" "0.00948270616614102"
[1] "iterated logL"
[1] -3396.113
[1] "Entering iteration loop " "7"
[1] "Checkpt1"
[1] "AvDev=" "0.00236905376963106"
[1] "iterated logL"
[1] -3395.669
[1] "Entering iteration loop " "8"
[1] "Checkpt1"
[1] "AvDev=" "0.00053187290600404"
[1] "iterated logL"
[1] -3395.579
[1] "Entering iteration loop " "9"
[1] "Checkpt1"
[1] "AvDev=" "0.000127236444056749"
[1] "iterated logL"
[1] -3395.558
[1] "Entering iteration loop " "10"
[1] "Checkpt1"
[1] "AvDev=" "3.03095178792535e-05"
[1] "iterated logL"
[1] -3395.553
[1] "Entering iteration loop " "11"
[1] "Checkpt1"
[1] "AvDev=" "7.27387684576024e-06"
[1] "iterated logL"
[1] -3395.552
[1] "Entered MLica"
[1] "Entering iteration loop " "0"
[1] "Checkpt1"
[1] "AvDev=" "0.408283532394995"
[1] "iterated logL"
[1] -4144.323
[1] "Entering iteration loop " "1"
[1] "Checkpt1"
[1] "AvDev=" "0.131991485639095"
[1] "iterated logL"
[1] -3479.347
[1] "Entering iteration loop " "2"
[1] "Checkpt1"
[1] "AvDev=" "0.0481662412793996"
[1] "iterated logL"
[1] -3443.261
[1] "Entering iteration loop " "3"
[1] "Checkpt1"
[1] "AvDev=" "0.0410388259598368"
[1] "iterated logL"
[1] -3401.783
[1] "Entering iteration loop " "4"
[1] "Checkpt1"
[1] "AvDev=" "0.00960736728874848"
[1] "iterated logL"
[1] -3395.308
[1] "Entering iteration loop " "5"
[1] "Checkpt1"
[1] "AvDev=" "0.00159036836987607"
[1] "iterated logL"
[1] -3395.651
[1] "Entering iteration loop " "6"
[1] "Checkpt1"
[1] "AvDev=" "0.00030326441425986"
[1] "iterated logL"
[1] -3395.55
[1] "Entering iteration loop " "7"
[1] "Checkpt1"
[1] "AvDev=" "5.75657634615365e-05"
[1] "iterated logL"
[1] -3395.555
[1] "Entering iteration loop " "8"
[1] "Checkpt1"
[1] "AvDev=" "1.37327379236027e-05"
[1] "iterated logL"
[1] -3395.552
[1] "Entered MLica"
[1] "Entering iteration loop " "0"
[1] "Checkpt1"
[1] "AvDev=" "0.425815168841733"
[1] "iterated logL"
[1] -3667.934
[1] "Entering iteration loop " "1"
[1] "Checkpt1"
[1] "AvDev=" "0.0698379808222559"
[1] "iterated logL"
[1] -3524.008
[1] "Entering iteration loop " "2"
[1] "Checkpt1"
[1] "AvDev=" "0.0131519647708972"
[1] "iterated logL"
[1] -3509.167
[1] "Entering iteration loop " "3"
[1] "Checkpt1"
[1] "AvDev=" "0.00486706941295416"
[1] "iterated logL"
[1] -3507.437
[1] "Entering iteration loop " "4"
[1] "Checkpt1"
[1] "AvDev=" "0.00625192454533561"
[1] "iterated logL"
[1] -3509.155
[1] "Entering iteration loop " "5"
[1] "Checkpt1"
[1] "AvDev=" "0.0115161072696129"
[1] "iterated logL"
[1] -3510.774
[1] "Entering iteration loop " "6"
[1] "Checkpt1"
[1] "AvDev=" "0.0207314015437168"
[1] "iterated logL"
[1] -3506.719
[1] "Entering iteration loop " "7"
[1] "Checkpt1"
[1] "AvDev=" "0.0316288909278083"
[1] "iterated logL"
[1] -3479.767
[1] "Entering iteration loop " "8"
[1] "Checkpt1"
[1] "AvDev=" "0.0325973768372488"
[1] "iterated logL"
[1] -3425.318
[1] "Entering iteration loop " "9"
[1] "Checkpt1"
[1] "AvDev=" "0.018965538980085"
[1] "iterated logL"
[1] -3399.135
[1] "Entering iteration loop " "10"
[1] "Checkpt1"
[1] "AvDev=" "0.00534623296704367"
[1] "iterated logL"
[1] -3395.796
[1] "Entering iteration loop " "11"
[1] "Checkpt1"
[1] "AvDev=" "0.00113679433105254"
[1] "iterated logL"
[1] -3395.569
[1] "Entering iteration loop " "12"
[1] "Checkpt1"
[1] "AvDev=" "0.000209041273698903"
[1] "iterated logL"
[1] -3395.556
[1] "Entering iteration loop " "13"
[1] "Checkpt1"
[1] "AvDev=" "3.82955775854774e-05"
[1] "iterated logL"
[1] -3395.553
[1] "Entering iteration loop " "14"
[1] "Checkpt1"
[1] "AvDev=" "7.38487073760532e-06"
[1] "iterated logL"
[1] -3395.552
[1] "Entered MLica"
[1] "Entering iteration loop " "0"
[1] "Checkpt1"
[1] "AvDev=" "0.371359542049407"
[1] "iterated logL"
[1] -3876.2
[1] "Entering iteration loop " "1"
[1] "Checkpt1"
[1] "AvDev=" "0.142530946467836"
[1] "iterated logL"
[1] -3441.852
[1] "Entering iteration loop " "2"
[1] "Checkpt1"
[1] "AvDev=" "0.0334766990223169"
[1] "iterated logL"
[1] -3395.664
[1] "Entering iteration loop " "3"
[1] "Checkpt1"
[1] "AvDev=" "0.00382891203440433"
[1] "iterated logL"
[1] -3395.684
[1] "Entering iteration loop " "4"
[1] "Checkpt1"
[1] "AvDev=" "0.000594310499221397"
[1] "iterated logL"
[1] -3395.562
[1] "Entering iteration loop " "5"
[1] "Checkpt1"
[1] "AvDev=" "0.000116149447185206"
[1] "iterated logL"
[1] -3395.558
[1] "Entering iteration loop " "6"
[1] "Checkpt1"
[1] "AvDev=" "2.83853526208703e-05"
[1] "iterated logL"
[1] -3395.553
[1] "Entering iteration loop " "7"
[1] "Checkpt1"
[1] "AvDev=" "6.14471227920681e-06"
[1] "iterated logL"
[1] -3395.552
[1] "Entered MLica"
[1] "Entering iteration loop " "0"
[1] "Checkpt1"
[1] "AvDev=" "0.410470544121206"
[1] "iterated logL"
[1] -3761.94
[1] "Entering iteration loop " "1"
[1] "Checkpt1"
[1] "AvDev=" "0.0894245649461816"
[1] "iterated logL"
[1] -3558.636
[1] "Entering iteration loop " "2"
[1] "Checkpt1"
[1] "AvDev=" "0.0295823341135455"
[1] "iterated logL"
[1] -3497.191
[1] "Entering iteration loop " "3"
[1] "Checkpt1"
[1] "AvDev=" "0.0319827688603013"
[1] "iterated logL"
[1] -3458.716
[1] "Entering iteration loop " "4"
[1] "Checkpt1"
[1] "AvDev=" "0.0458463230582255"
[1] "iterated logL"
[1] -3414.967
[1] "Entering iteration loop " "5"
[1] "Checkpt1"
[1] "AvDev=" "0.0204951937399991"
[1] "iterated logL"
[1] -3398.147
[1] "Entering iteration loop " "6"
[1] "Checkpt1"
[1] "AvDev=" "0.00345410150464575"
[1] "iterated logL"
[1] -3395.399
[1] "Entering iteration loop " "7"
[1] "Checkpt1"
[1] "AvDev=" "0.000624447866132663"
[1] "iterated logL"
[1] -3395.583
[1] "Entering iteration loop " "8"
[1] "Checkpt1"
[1] "AvDev=" "0.000111555927307654"
[1] "iterated logL"
[1] -3395.548
[1] "Entering iteration loop " "9"
[1] "Checkpt1"
[1] "AvDev=" "2.68388826026896e-05"
[1] "iterated logL"
[1] -3395.552
[1] "Entering iteration loop " "10"
[1] "Checkpt1"
[1] "AvDev=" "5.42100137719388e-06"
[1] "iterated logL"
[1] -3395.552
[1] "End of runs"
[1] "Entering mlica"
[1] "Performing preliminary run"
[1] "Entered MLica"
[1] "Entering iteration loop " "0"
[1] "Checkpt1"
[1] "AvDev=" "0.366074296551384"
[1] "iterated logL"
[1] -4159.965
[1] "Entering iteration loop " "1"
[1] "Checkpt1"
[1] "AvDev=" "0.0990706074759644"
[1] "iterated logL"
[1] -3842.361
[1] "Entering iteration loop " "2"
[1] "Checkpt1"
[1] "AvDev=" "0.0925390792823253"
[1] "iterated logL"
[1] -3492.819
[1] "Entering iteration loop " "3"
[1] "Checkpt1"
[1] "AvDev=" "0.0526223128722005"
[1] "iterated logL"
[1] -3403.373
[1] "Entering iteration loop " "4"
[1] "Checkpt1"
[1] "AvDev=" "0.0181231017247815"
[1] "iterated logL"
[1] -3396.83
[1] "Entering iteration loop " "5"
[1] "Checkpt1"
[1] "AvDev=" "0.00445475815844939"
[1] "iterated logL"
[1] -3395.756
[1] "Entering iteration loop " "6"
[1] "Checkpt1"
[1] "AvDev=" "0.000901564622348875"
[1] "iterated logL"
[1] -3395.602
[1] "Entering iteration loop " "7"
[1] "Checkpt1"
[1] "AvDev=" "0.000213306212586903"
[1] "iterated logL"
[1] -3395.562
[1] "Entering iteration loop " "8"
[1] "Checkpt1"
[1] "AvDev=" "4.75160900748634e-05"
[1] "iterated logL"
[1] -3395.555
[1] "Entering iteration loop " "9"
[1] "Checkpt1"
[1] "AvDev=" "1.16181179872218e-05"
[1] "iterated logL"
[1] -3395.552
[1] "Finished preliminary run"
[1] "Starting runs"
[1] "Entered MLica"
[1] "Entering iteration loop " "0"
[1] "Checkpt1"
[1] "AvDev=" "0.438326263367479"
[1] "iterated logL"
[1] -4019.211
[1] "Entering iteration loop " "1"
[1] "Checkpt1"
[1] "AvDev=" "0.0999407970722891"
[1] "iterated logL"
[1] -3937.154
[1] "Entering iteration loop " "2"
[1] "Checkpt1"
[1] "AvDev=" "0.0857379663920665"
[1] "iterated logL"
[1] -3583.106
[1] "Entering iteration loop " "3"
[1] "Checkpt1"
[1] "AvDev=" "0.0749152696209581"
[1] "iterated logL"
[1] -3433.51
[1] "Entering iteration loop " "4"
[1] "Checkpt1"
[1] "AvDev=" "0.0259885513192216"
[1] "iterated logL"
[1] -3400.524
[1] "Entering iteration loop " "5"
[1] "Checkpt1"
[1] "AvDev=" "0.00445187199563686"
[1] "iterated logL"
[1] -3395.51
[1] "Entering iteration loop " "6"
[1] "Checkpt1"
[1] "AvDev=" "0.000831511651645693"
[1] "iterated logL"
[1] -3395.625
[1] "Entering iteration loop " "7"
[1] "Checkpt1"
[1] "AvDev=" "0.000192552944613526"
[1] "iterated logL"
[1] -3395.55
[1] "Entering iteration loop " "8"
[1] "Checkpt1"
[1] "AvDev=" "3.89059018286901e-05"
[1] "iterated logL"
[1] -3395.555
[1] "Entering iteration loop " "9"
[1] "Checkpt1"
[1] "AvDev=" "9.37276456141046e-06"
[1] "iterated logL"
[1] -3395.552
[1] "Entered MLica"
[1] "Entering iteration loop " "0"
[1] "Checkpt1"
[1] "AvDev=" "0.351009643350607"
[1] "iterated logL"
[1] -3735.173
[1] "Entering iteration loop " "1"
[1] "Checkpt1"
[1] "AvDev=" "0.109323958491958"
[1] "iterated logL"
[1] -3544.713
[1] "Entering iteration loop " "2"
[1] "Checkpt1"
[1] "AvDev=" "0.0707197937973701"
[1] "iterated logL"
[1] -3486.343
[1] "Entering iteration loop " "3"
[1] "Checkpt1"
[1] "AvDev=" "0.04351262542631"
[1] "iterated logL"
[1] -3414.097
[1] "Entering iteration loop " "4"
[1] "Checkpt1"
[1] "AvDev=" "0.0146397623893175"
[1] "iterated logL"
[1] -3396.041
[1] "Entering iteration loop " "5"
[1] "Checkpt1"
[1] "AvDev=" "0.00178476371045314"
[1] "iterated logL"
[1] -3395.552
[1] "Entering iteration loop " "6"
[1] "Checkpt1"
[1] "AvDev=" "0.000333464059947828"
[1] "iterated logL"
[1] -3395.562
[1] "Entering iteration loop " "7"
[1] "Checkpt1"
[1] "AvDev=" "7.88913069801145e-05"
[1] "iterated logL"
[1] -3395.554
[1] "Entering iteration loop " "8"
[1] "Checkpt1"
[1] "AvDev=" "1.78322360740409e-05"
[1] "iterated logL"
[1] -3395.552
[1] "Entered MLica"
[1] "Entering iteration loop " "0"
[1] "Checkpt1"
[1] "AvDev=" "0.49119825887111"
[1] "iterated logL"
[1] -4801.95
[1] "Entering iteration loop " "1"
[1] "Checkpt1"
[1] "AvDev=" "0.141448069534506"
[1] "iterated logL"
[1] -3756.975
[1] "Entering iteration loop " "2"
[1] "Checkpt1"
[1] "AvDev=" "0.113921181520134"
[1] "iterated logL"
[1] -3445.1
[1] "Entering iteration loop " "3"
[1] "Checkpt1"
[1] "AvDev=" "0.0451228348785501"
[1] "iterated logL"
[1] -3400.834
[1] "Entering iteration loop " "4"
[1] "Checkpt1"
[1] "AvDev=" "0.0120801507087052"
[1] "iterated logL"
[1] -3396.434
[1] "Entering iteration loop " "5"
[1] "Checkpt1"
[1] "AvDev=" "0.00329053928203262"
[1] "iterated logL"
[1] -3395.715
[1] "Entering iteration loop " "6"
[1] "Checkpt1"
[1] "AvDev=" "0.000666767357514407"
[1] "iterated logL"
[1] -3395.592
[1] "Entering iteration loop " "7"
[1] "Checkpt1"
[1] "AvDev=" "0.000162347881609018"
[1] "iterated logL"
[1] -3395.56
[1] "Entering iteration loop " "8"
[1] "Checkpt1"
[1] "AvDev=" "3.73592487872754e-05"
[1] "iterated logL"
[1] -3395.554
[1] "Entering iteration loop " "9"
[1] "Checkpt1"
[1] "AvDev=" "9.07690874152711e-06"
[1] "iterated logL"
[1] -3395.552
[1] "Entered MLica"
[1] "Entering iteration loop " "0"
[1] "Checkpt1"
[1] "AvDev=" "0.36746082558709"
[1] "iterated logL"
[1] -4123.382
[1] "Entering iteration loop " "1"
[1] "Checkpt1"
[1] "AvDev=" "0.13094735679392"
[1] "iterated logL"
[1] -3518.547
[1] "Entering iteration loop " "2"
[1] "Checkpt1"
[1] "AvDev=" "0.0521414764338787"
[1] "iterated logL"
[1] -3423.406
[1] "Entering iteration loop " "3"
[1] "Checkpt1"
[1] "AvDev=" "0.0150553718379829"
[1] "iterated logL"
[1] -3398.158
[1] "Entering iteration loop " "4"
[1] "Checkpt1"
[1] "AvDev=" "0.00365878582864727"
[1] "iterated logL"
[1] -3395.7
[1] "Entering iteration loop " "5"
[1] "Checkpt1"
[1] "AvDev=" "0.000857647373322417"
[1] "iterated logL"
[1] -3395.565
[1] "Entering iteration loop " "6"
[1] "Checkpt1"
[1] "AvDev=" "0.000194209693258142"
[1] "iterated logL"
[1] -3395.553
[1] "Entering iteration loop " "7"
[1] "Checkpt1"
[1] "AvDev=" "4.500676380397e-05"
[1] "iterated logL"
[1] -3395.552
[1] "Entering iteration loop " "8"
[1] "Checkpt1"
[1] "AvDev=" "1.02524134500583e-05"
[1] "iterated logL"
[1] -3395.552
[1] "Entered MLica"
[1] "Entering iteration loop " "0"
[1] "Checkpt1"
[1] "AvDev=" "0.381175722621573"
[1] "iterated logL"
[1] -4035.78
[1] "Entering iteration loop " "1"
[1] "Checkpt1"
[1] "AvDev=" "0.0828709227554443"
[1] "iterated logL"
[1] -3845.233
[1] "Entering iteration loop " "2"
[1] "Checkpt1"
[1] "AvDev=" "0.0935145589689355"
[1] "iterated logL"
[1] -3694.554
[1] "Entering iteration loop " "3"
[1] "Checkpt1"
[1] "AvDev=" "0.0844200429745166"
[1] "iterated logL"
[1] -3608.252
[1] "Entering iteration loop " "4"
[1] "Checkpt1"
[1] "AvDev=" "0.0663941351087329"
[1] "iterated logL"
[1] -3478.698
[1] "Entering iteration loop " "5"
[1] "Checkpt1"
[1] "AvDev=" "0.0446549747417372"
[1] "iterated logL"
[1] -3406.252
[1] "Entering iteration loop " "6"
[1] "Checkpt1"
[1] "AvDev=" "0.0140499598598769"
[1] "iterated logL"
[1] -3395.556
[1] "Entering iteration loop " "7"
[1] "Checkpt1"
[1] "AvDev=" "0.00392616751329586"
[1] "iterated logL"
[1] -3395.332
[1] "Entering iteration loop " "8"
[1] "Checkpt1"
[1] "AvDev=" "0.000990970829615874"
[1] "iterated logL"
[1] -3395.499
[1] "Entering iteration loop " "9"
[1] "Checkpt1"
[1] "AvDev=" "0.000237840071878959"
[1] "iterated logL"
[1] -3395.538
[1] "Entering iteration loop " "10"
[1] "Checkpt1"
[1] "AvDev=" "5.74841248261739e-05"
[1] "iterated logL"
[1] -3395.549
[1] "Entering iteration loop " "11"
[1] "Checkpt1"
[1] "AvDev=" "1.36830985332594e-05"
[1] "iterated logL"
[1] -3395.551
[1] "End of runs"
There were 50 or more warnings (use warnings() to see the first 50)
Warning message:
In solve.default(STS, LINPACK = FALSE) :
the LINPACK argument has been defunct since R 3.1.0
1 2 3 4 5
[1,] 1.707899e-02 0.006944718 -0.041919257 -0.996658231 -0.020187865
[2,] -1.244749e-02 0.045893055 -0.030323965 0.039829080 0.998275626
[3,] 7.363413e-05 0.032804710 0.998660069 -0.054682965 0.045595331
[4,] -4.938957e-02 -0.996255792 -0.000171045 -0.006949881 0.030553821
[5,] -9.985560e-01 0.065144815 0.001126950 0.045249163 -0.005074941
function (prNCP, nruns = 10, tol = 1e-04, maxit = 300, fail.th = 5,
learn.mu = 1)
{
print("Entering mlica")
print("Performing preliminary run")
a <- mlicaMAIN(prNCP, tol = 1e-04, maxit = 10, mu = learn.mu)
ncp <- dim(a$S)[2]
max.logL <- a$LL
a.best <- a
print("Finished preliminary run")
print("Starting runs")
run.n <- 0
fail.count <- 0
v.logL <- vector()
v.NC <- vector()
while (run.n < nruns) {
a <- mlicaMAIN(prNCP, tol = 1e-04, maxit = maxit, mu = learn.mu)
v.logL <- c(v.logL, a$LL)
v.NC <- c(v.NC, a$NC)
if (a$NC == 0) {
fail.count <- 0
run.n <- run.n + 1
if (a$LL > max.logL) {
a.best <- a
}
}
else {
fail.count <- fail.count + 1
}
if (fail.count >= fail.th) {
print("Stopping: Five consecutive runs failed to converge!")
print("Consider either increasing the threshold for pca eigenvalues to perform ICA on a smaller subspace or increasing maxit")
stop
}
}
print("End of runs")
return(a.best)
}
R Package Documentation
Browse R Packages
We want your feedback!
Note that we can't provide technical support on individual packages. You should contact the package authors for that.
Embedding an R snippet on your website
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.