Daily Means of NOx (mono-nitrogen oxides) in air
ambientNOxCH.RdThis dataset contains daily means (from midnight to midnight) of NOx, i.e., mono-nitrogen oxides, in [ppb] at 13 sites in central Switzerland and Aarau for the year 2004.
data(ambientNOxCH, package="robustbase")Format
A data frame with 366 observations on the following 14 variables.
datedate of day, of class
"Date".adSite is located north of Altdorf 100 meters east of motorway A2, on an open field at the beginning of a more than 2000m deep valley (690.175, 193.55; 438; inLuft)
baSite is located in the centre of the little town of Baden in a residential area. Baden has 34'000 inhabitants and is situated on the swiss plateau (666.075, 257.972; 377; inLuft).
efSite is located 6 km south of altdorf and 800 m north of the village of Erstfeld. The motorway A2 passes 5 m west of the measuring site. Over 8 million vehicles have passed Erstfeld in 2004 where 13% of the counts were attributed to trucks (691.43, 187.69; 457; MFM-U).
laSite is located on a wooded hill in a rural area called Laegern, about 190 m above Baden, which is about 5 km away (669.8, 259; 690; NABEL).
luSite is located in the center of town of Lucerne, which has 57'000 inhabitants (666.19, 211.975; 460; inLuft).
reSite is located 1 km west of Reiden on the Swiss plateau. The motorway A2 passes 5 m west of the measuring site (639.56, 232.11; 462; MFM-U).
riSite is located at Rigi Seebodenalp, 649 m above the lake of Lucerne on an alp with half a dozen small houses (677.9, 213.5; 1030; NABEL).
seSite is located in Sedel next to town of Lucerne 35m above and 250m south of motorway A14 from Zug to Lucerne on a low hill with free 360° panorama (665.5, 213.41; 484; inLuft).
siSite is located at the border of a small industrial area in Sisseln, 300 m east of a main road (640.725, 266.25; 305; inLuft).
stSite is located at the south east border of Stans with 7'000 inhabitants (670.85, 201.025; 438; inLuft).
suSite is located in the center of Suhr (8700 inhabitants), 10 m from the main road (648.49, 246.985; 403; inLuft).
szSite is located in Schwyz (14'200 inhabitants) near a shopping center (691.92, 208.03; 470; inLuft).
zgSite is located in the centre of Zug with 22'000 inhabitants, 24 m from the main road (681.625, 224.625; 420; inLuft).
Details
The 13 sites are part of one of the three air quality monitoring networks:
inLuft (regional authorities of central Switzerland and canton Aargau)
NABEL (Swiss federal network)
MFM-U (Monitoring flankierende Massnahmen Umwelt), special Swiss
federal network along transit motorways A2 and A13 from Germany to
Italy through Switzerland
The information within the brackets means: Swiss
coordinates km east, km north; m above sea level; network
When the measuring sites are exposed to the same atmospheric condition
and when there is no singular emission event at any site,
log(mean(NOx) of a specific day at each site) is a linear
function of log(yearly.mean(NOx) at the corresponding site). The
offset and the slope of the straight
line reflects the atmospheric conditions at this specific day. During
winter time, often an inversion prevents the emissions from being
diluted vertically, so that there evolve two separate atmospheric
compartements: One below the inversion boundary with polluted air and one
above with relatively clean air. In our example below, Rigi Seebodenalp
is above the inversion boundary between December 10th and 12th.
Source
http://www.in-luft.ch/
http://www.empa.ch/plugin/template/empa/*/6794
http://www.bafu.admin.ch/umweltbeobachtung/02272/02280
See also
another NOx dataset, NOxEmissions.
Examples
data(ambientNOxCH)
str (ambientNOxCH)
#> 'data.frame': 366 obs. of 14 variables:
#> $ date: Date, format: "2004-01-01" "2004-01-02" ...
#> $ ad : num 12 17 12.3 13.5 47 ...
#> $ ba : num 14.66 18.71 8.75 19.35 57.47 ...
#> $ ef : num 17.3 40.4 35.1 26.3 53.3 ...
#> $ la : num 7.73 8.01 6.83 11.84 10.74 ...
#> $ lu : num 22.7 21.2 26.8 26.5 58.6 ...
#> $ re : num 20.1 45.2 8.4 37.4 110.4 ...
#> $ ri : num 6.32 8.73 10.31 7 1.4 ...
#> $ se : num 11.3 15.4 11 19.4 NA ...
#> $ si : num 13.46 15.76 9.82 13.02 37.27 ...
#> $ st : num 13.4 19.3 16.5 15 27.8 ...
#> $ su : num 27 35.8 34.9 34.8 78.7 ...
#> $ sz : num 16.2 21.2 16.4 18.1 42.7 ...
#> $ zg : num 25.5 34.5 29.1 26.2 76.9 ...
yearly <- log(colMeans(ambientNOxCH[,-1], na.rm=TRUE))
xlim <- range(yearly)
lNOx <- log(ambientNOxCH[, -1])
days <- ambientNOxCH[, "date"]
## Subset of 9 days starting at April 4:
idays <- seq(which(ambientNOxCH$date=="2004-12-04"), length=9)
ylim <- range(lNOx[idays,],na.rm=TRUE)
op <- par(mfrow=c(3,3),mar=rep(1,4), oma = c(0,0,2,0))
for (id in idays) {
daily <- unlist(lNOx[id,])
plot(NA, xlim=xlim,ylim=ylim, ann=FALSE, type = "n")
abline(0:1, col="light gray")
abline(lmrob(daily~yearly, na.action=na.exclude),
col="red", lwd=2)
text(yearly, daily, names(yearly), col="blue")
mtext(days[id], side=1, line=-1.2, cex=.75, adj=.98)
}
#> Warning: find_scale() did not converge in 'maxit.scale' (= 200) iterations with tol=1e-10, last rel.diff=0
#> Warning: find_scale() did not converge in 'maxit.scale' (= 200) iterations with tol=1e-10, last rel.diff=0
#> Warning: find_scale() did not converge in 'maxit.scale' (= 200) iterations with tol=1e-10, last rel.diff=0
#> Warning: find_scale() did not converge in 'maxit.scale' (= 200) iterations with tol=1e-10, last rel.diff=0
#> Warning: find_scale() did not converge in 'maxit.scale' (= 200) iterations with tol=1e-10, last rel.diff=0
#> Warning: find_scale() did not converge in 'maxit.scale' (= 200) iterations with tol=1e-10, last rel.diff=0
#> Warning: find_scale() did not converge in 'maxit.scale' (= 200) iterations with tol=1e-10, last rel.diff=0
mtext("Daily ~ Yearly log( NOx mean values ) at 13 Swiss locations",
outer=TRUE)
par(op)
## do all 366 regressions: Least Squares and Robust:
LS <- lapply(1:nrow(ambientNOxCH), function(id)
lm(unlist(lNOx[id,]) ~ yearly,
na.action = na.exclude))
R <- lapply(1:nrow(ambientNOxCH),
function(id) lmrob(unlist(lNOx[id,]) ~ yearly,
na.action = na.exclude))
#> Warning: M-step did NOT converge. Returning unconverged SM-estimate
#> Warning: find_scale() did not converge in 'maxit.scale' (= 200) iterations with tol=1e-10, last rel.diff=0
#> Warning: find_scale() did not converge in 'maxit.scale' (= 200) iterations with tol=1e-10, last rel.diff=0
#> Warning: find_scale() did not converge in 'maxit.scale' (= 200) iterations with tol=1e-10, last rel.diff=0
#> Warning: find_scale() did not converge in 'maxit.scale' (= 200) iterations with tol=1e-10, last rel.diff=0
#> Warning: find_scale() did not converge in 'maxit.scale' (= 200) iterations with tol=1e-10, last rel.diff=0
#> Warning: find_scale() did not converge in 'maxit.scale' (= 200) iterations with tol=1e-10, last rel.diff=0
#> Warning: find_scale() did not converge in 'maxit.scale' (= 200) iterations with tol=1e-10, last rel.diff=0
#> Warning: find_scale() did not converge in 'maxit.scale' (= 200) iterations with tol=1e-10, last rel.diff=0
#> Warning: M-step did NOT converge. Returning unconverged SM-estimate
#> Warning: M-step did NOT converge. Returning unconverged SM-estimate
#> Warning: find_scale() did not converge in 'maxit.scale' (= 200) iterations with tol=1e-10, last rel.diff=0
#> Warning: find_scale() did not converge in 'maxit.scale' (= 200) iterations with tol=1e-10, last rel.diff=0
#> Warning: find_scale() did not converge in 'maxit.scale' (= 200) iterations with tol=1e-10, last rel.diff=0
#> Warning: find_scale() did not converge in 'maxit.scale' (= 200) iterations with tol=1e-10, last rel.diff=0
#> Warning: find_scale() did not converge in 'maxit.scale' (= 200) iterations with tol=1e-10, last rel.diff=0
#> Warning: find_scale() did not converge in 'maxit.scale' (= 200) iterations with tol=1e-10, last rel.diff=0
#> Warning: find_scale() did not converge in 'maxit.scale' (= 200) iterations with tol=1e-10, last rel.diff=0
#> Warning: find_scale() did not converge in 'maxit.scale' (= 200) iterations with tol=1e-10, last rel.diff=0
#> Warning: find_scale() did not converge in 'maxit.scale' (= 200) iterations with tol=1e-10, last rel.diff=0
#> Warning: find_scale() did not converge in 'maxit.scale' (= 200) iterations with tol=1e-10, last rel.diff=0
#> Warning: M-step did NOT converge. Returning unconverged SM-estimate
#> Warning: find_scale() did not converge in 'maxit.scale' (= 200) iterations with tol=1e-10, last rel.diff=0
#> Warning: find_scale() did not converge in 'maxit.scale' (= 200) iterations with tol=1e-10, last rel.diff=0
#> Warning: find_scale() did not converge in 'maxit.scale' (= 200) iterations with tol=1e-10, last rel.diff=0
#> Warning: find_scale() did not converge in 'maxit.scale' (= 200) iterations with tol=1e-10, last rel.diff=0
#> Warning: find_scale() did not converge in 'maxit.scale' (= 200) iterations with tol=1e-10, last rel.diff=0
#> Warning: find_scale() did not converge in 'maxit.scale' (= 200) iterations with tol=1e-10, last rel.diff=0
#> Warning: find_scale() did not converge in 'maxit.scale' (= 200) iterations with tol=1e-10, last rel.diff=0
#> Warning: find_scale() did not converge in 'maxit.scale' (= 200) iterations with tol=1e-10, last rel.diff=0
#> Warning: find_scale() did not converge in 'maxit.scale' (= 200) iterations with tol=1e-10, last rel.diff=0
#> Warning: find_scale() did not converge in 'maxit.scale' (= 200) iterations with tol=1e-10, last rel.diff=0
#> Warning: find_scale() did not converge in 'maxit.scale' (= 200) iterations with tol=1e-10, last rel.diff=0
#> Warning: find_scale() did not converge in 'maxit.scale' (= 200) iterations with tol=1e-10, last rel.diff=0
#> Warning: find_scale() did not converge in 'maxit.scale' (= 200) iterations with tol=1e-10, last rel.diff=0
#> Warning: find_scale() did not converge in 'maxit.scale' (= 200) iterations with tol=1e-10, last rel.diff=0
#> Warning: find_scale() did not converge in 'maxit.scale' (= 200) iterations with tol=1e-10, last rel.diff=0
#> Warning: find_scale() did not converge in 'maxit.scale' (= 200) iterations with tol=1e-10, last rel.diff=0
#> Warning: find_scale() did not converge in 'maxit.scale' (= 200) iterations with tol=1e-10, last rel.diff=0
## currently 4 warnings about non-convergence;
## which ones?
days[notOk <- ! sapply(R, `[[`, "converged") ]
#> [1] "2004-01-10" "2004-05-12" "2004-05-16" "2004-11-16"
## "2004-01-10" "2004-05-12" "2004-05-16" "2004-11-16"
## first problematic case:
daily <- unlist(lNOx[which(notOk)[1],])
plot(daily ~ yearly,
main = paste("lmrob() non-convergent:",days[notOk[1]]))
rr <- lmrob(daily ~ yearly, na.action = na.exclude,
control = lmrob.control(trace=3, max.it = 100))
#> lmrob_S(n = 13, nRes = 500): fast_s() [non-large n]:
#> fast_s(*, s_y=2.68976, n=13, p=2, ipsi=1, ..) before INIT_WLS():
#> Subsampling 500 times to find candidate betas:
#> Sample[ 0]: idc = 1 7
#> b^[] = -13.1628 5.179
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -13.1628 5.179
#> after refine_(*, conv=F):
#> beta_ref : -0.962999 1.2121
#> with ||beta_ref - beta_cand|| = 12.8285446507, --> sc = 2.75529751980209
#> find_scale(*, ini.scale = 2.7552975198, tol=1e-10):
#> it | new scale
#> 0 | 0.6887294089
#> 1 | 0.5514682195
#> 2 | 0.4873381196
#> 3 | 0.4498934328
#> 4 | 0.4268842462
#> 5 | 0.4125538979
#> 6 | 0.4036216805
#> 7 | 0.3980495439
#> 8 | 0.3945701010
#> 9 | 0.3923956398
#> 10 | 0.3910359122
#> 11 | 0.3901853054
#> 12 | 0.3896530471
#> 13 | 0.3893199342
#> 14 | 0.3891114330
#> 15 | 0.3889809193
#> 16 | 0.3888992191
#> 17 | 0.3888480742
#> 18 | 0.3888160566
#> 19 | 0.3887960128
#> 20 | 0.3887834648
#> 21 | 0.3887756093
#> 22 | 0.3887706915
#> 23 | 0.3887676129
#> 24 | 0.3887656855
#> 25 | 0.3887644789
#> 26 | 0.3887637235
#> 27 | 0.3887632506
#> 28 | 0.3887629546
#> 29 | 0.3887627692
#> 30 | 0.3887626532
#> 31 | 0.3887625806
#> 32 | 0.3887625351
#> 33 | 0.3887625066
#> 34 | 0.3887624888
#> 35 | 0.3887624776
#> 36 | 0.3887624707
#> 37 | 0.3887624663
#> 38 | 0.3887624635
#> 39 | 0.3887624618
#> 40 | 0.3887624608
#> 41 | 0.3887624601
#> 42 | 0.3887624597
#> 43 | 0.3887624594
#> 44 | 0.3887624592
#> 45 | 0.3887624591
#> 46 | 0.3887624591
#> 47 | 0.3887624590
#> 48 | 0.3887624590
#> found new candidate with scale 0.3887625 in 48 iter (worst sc=2.6898e+20)
#> Sample[ 1]: idc = 4 1
#> b^[] = 2.25191 0.271184
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 2.25191 0.271184
#> after refine_(*, conv=F):
#> beta_ref : 0.0105456 0.920859
#> with ||beta_ref - beta_cand|| = 2.33362366183, --> sc = 0.813737418280548
#> find_scale(*, ini.scale =0.81373741828, tol=1e-10):
#> it | new scale
#> 0 | 0.5999288225
#> 1 | 0.5069437881
#> 2 | 0.4606430377
#> 3 | 0.4367117751
#> 4 | 0.4239104087
#> 5 | 0.4168979926
#> 6 | 0.4129998849
#> 7 | 0.4108141207
#> 8 | 0.4095823474
#> 9 | 0.4088861910
#> 10 | 0.4084921012
#> 11 | 0.4082688009
#> 12 | 0.4081422069
#> 13 | 0.4080704161
#> 14 | 0.4080296970
#> 15 | 0.4080065993
#> 16 | 0.4079934965
#> 17 | 0.4079860633
#> 18 | 0.4079818465
#> 19 | 0.4079794542
#> 20 | 0.4079780970
#> 21 | 0.4079773271
#> 22 | 0.4079768903
#> 23 | 0.4079766424
#> 24 | 0.4079765019
#> 25 | 0.4079764221
#> 26 | 0.4079763769
#> 27 | 0.4079763512
#> 28 | 0.4079763366
#> 29 | 0.4079763284
#> 30 | 0.4079763237
#> 31 | 0.4079763210
#> 32 | 0.4079763195
#> 33 | 0.4079763187
#> 34 | 0.4079763182
#> 35 | 0.4079763179
#> 36 | 0.4079763177
#> 37 | 0.4079763176
#> 38 | 0.4079763176
#> 39 | 0.4079763176
#> found new candidate with scale 0.4079763 in 39 iter (worst sc=0.40798)
#> Sample[ 2]: idc = 2 6
#> b^[] = -1.78547 1.24971
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -1.78547 1.24971
#> after refine_(*, conv=F):
#> beta_ref : -1.44982 1.30483
#> with ||beta_ref - beta_cand|| = 0.340142551919, --> sc = 1.10255399571346
#> Sample[ 3]: idc = 6 9
#> b^[] = -2.77944 1.94279
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -2.77944 1.94279
#> after refine_(*, conv=F):
#> beta_ref : -1.85666 1.53242
#> with ||beta_ref - beta_cand|| = 1.00990893711, --> sc = 0.885734481714681
#> find_scale(*, ini.scale =0.88573448171, tol=1e-10):
#> it | new scale
#> 0 | 0.6425722266
#> 1 | 0.5357202428
#> 2 | 0.4804169373
#> 3 | 0.4503098092
#> 4 | 0.4331560654
#> 5 | 0.4230703440
#> 6 | 0.4170211694
#> 7 | 0.4133482914
#> 8 | 0.4111014368
#> 9 | 0.4097206062
#> 10 | 0.4088695999
#> 11 | 0.4083442109
#> 12 | 0.4080195010
#> 13 | 0.4078186852
#> 14 | 0.4076944402
#> 15 | 0.4076175503
#> 16 | 0.4075699589
#> 17 | 0.4075404991
#> 18 | 0.4075222620
#> 19 | 0.4075109718
#> 20 | 0.4075039822
#> 21 | 0.4074996550
#> 22 | 0.4074969759
#> 23 | 0.4074953173
#> 24 | 0.4074942905
#> 25 | 0.4074936547
#> 26 | 0.4074932611
#> 27 | 0.4074930175
#> 28 | 0.4074928666
#> 29 | 0.4074927732
#> 30 | 0.4074927154
#> 31 | 0.4074926796
#> 32 | 0.4074926574
#> 33 | 0.4074926437
#> 34 | 0.4074926352
#> 35 | 0.4074926299
#> 36 | 0.4074926267
#> 37 | 0.4074926247
#> 38 | 0.4074926234
#> 39 | 0.4074926226
#> 40 | 0.4074926222
#> 41 | 0.4074926219
#> 42 | 0.4074926217
#> 43 | 0.4074926216
#> 44 | 0.4074926215
#> 45 | 0.4074926214
#> 46 | 0.4074926214
#> found new candidate with scale 0.4074926 in 46 iter (worst sc=0.40749)
#> Sample[ 4]: idc = 2 12
#> b^[] = 6.21824 -0.743034
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 6.21824 -0.743034
#> after refine_(*, conv=F):
#> beta_ref : -0.199361 0.979666
#> with ||beta_ref - beta_cand|| = 6.64479808942, --> sc = 1.90760566523583
#> find_scale(*, ini.scale = 1.9076056652, tol=1e-10):
#> it | new scale
#> 0 | 0.7573874141
#> 1 | 0.5728868666
#> 2 | 0.4889139924
#> 3 | 0.4457215402
#> 4 | 0.4224382771
#> 5 | 0.4093795666
#> 6 | 0.4018555477
#> 7 | 0.3974464540
#> 8 | 0.3948358419
#> 9 | 0.3932803961
#> 10 | 0.3923501302
#> 11 | 0.3917924997
#> 12 | 0.3914577810
#> 13 | 0.3912567002
#> 14 | 0.3911358419
#> 15 | 0.3910631792
#> 16 | 0.3910194849
#> 17 | 0.3909932074
#> 18 | 0.3909774033
#> 19 | 0.3909678977
#> 20 | 0.3909621804
#> 21 | 0.3909587416
#> 22 | 0.3909566732
#> 23 | 0.3909554290
#> 24 | 0.3909546807
#> 25 | 0.3909542306
#> 26 | 0.3909539599
#> 27 | 0.3909537970
#> 28 | 0.3909536991
#> 29 | 0.3909536401
#> 30 | 0.3909536047
#> 31 | 0.3909535834
#> 32 | 0.3909535706
#> 33 | 0.3909535629
#> 34 | 0.3909535582
#> 35 | 0.3909535554
#> 36 | 0.3909535537
#> 37 | 0.3909535527
#> 38 | 0.3909535521
#> 39 | 0.3909535518
#> 40 | 0.3909535515
#> 41 | 0.3909535514
#> 42 | 0.3909535513
#> 43 | 0.3909535513
#> 44 | 0.3909535513
#> found new candidate with scale 0.3909536 in 44 iter (worst sc=0.39095)
#> Sample[ 5]: idc = 12 1
#> b^[] = -0.484845 1.14253
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -0.484845 1.14253
#> after refine_(*, conv=F):
#> beta_ref : -0.587028 1.14935
#> with ||beta_ref - beta_cand|| = 0.102411328665, --> sc = 0.312564686828127
#> find_scale(*, ini.scale =0.31256468683, tol=1e-10):
#> it | new scale
#> 0 | 0.3033762364
#> 1 | 0.2969165179
#> 2 | 0.2923364690
#> 3 | 0.2890682994
#> 4 | 0.2867251665
#> 5 | 0.2850393926
#> 6 | 0.2838234747
#> 7 | 0.2829448327
#> 8 | 0.2823090578
#> 9 | 0.2818485692
#> 10 | 0.2815148029
#> 11 | 0.2812727614
#> 12 | 0.2810971715
#> 13 | 0.2809697545
#> 14 | 0.2808772760
#> 15 | 0.2808101460
#> 16 | 0.2807614114
#> 17 | 0.2807260287
#> 18 | 0.2807003385
#> 19 | 0.2806816849
#> 20 | 0.2806681402
#> 21 | 0.2806583049
#> 22 | 0.2806511631
#> 23 | 0.2806459770
#> 24 | 0.2806422111
#> 25 | 0.2806394765
#> 26 | 0.2806374906
#> 27 | 0.2806360486
#> 28 | 0.2806350014
#> 29 | 0.2806342410
#> 30 | 0.2806336888
#> 31 | 0.2806332878
#> 32 | 0.2806329967
#> 33 | 0.2806327852
#> 34 | 0.2806326317
#> 35 | 0.2806325202
#> 36 | 0.2806324392
#> 37 | 0.2806323804
#> 38 | 0.2806323377
#> 39 | 0.2806323067
#> 40 | 0.2806322842
#> 41 | 0.2806322678
#> 42 | 0.2806322559
#> 43 | 0.2806322473
#> 44 | 0.2806322411
#> 45 | 0.2806322365
#> 46 | 0.2806322332
#> 47 | 0.2806322308
#> 48 | 0.2806322291
#> 49 | 0.2806322278
#> 50 | 0.2806322269
#> 51 | 0.2806322262
#> 52 | 0.2806322257
#> 53 | 0.2806322254
#> 54 | 0.2806322251
#> 55 | 0.2806322249
#> 56 | 0.2806322248
#> 57 | 0.2806322247
#> 58 | 0.2806322246
#> 59 | 0.2806322246
#> 60 | 0.2806322246
#> 61 | 0.2806322245
#> found new candidate with scale 0.2806322 in 61 iter (worst sc=0.38876)
#> Sample[ 6]: idc = 1 12
#> b^[] = -0.484845 1.14253
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -0.484845 1.14253
#> after refine_(*, conv=F):
#> beta_ref : -0.587028 1.14935
#> with ||beta_ref - beta_cand|| = 0.102411328665, --> sc = 0.312564686828127
#> find_scale(*, ini.scale =0.31256468683, tol=1e-10):
#> it | new scale
#> 0 | 0.3033762364
#> 1 | 0.2969165179
#> 2 | 0.2923364690
#> 3 | 0.2890682994
#> 4 | 0.2867251665
#> 5 | 0.2850393926
#> 6 | 0.2838234747
#> 7 | 0.2829448327
#> 8 | 0.2823090578
#> 9 | 0.2818485692
#> 10 | 0.2815148029
#> 11 | 0.2812727614
#> 12 | 0.2810971715
#> 13 | 0.2809697545
#> 14 | 0.2808772760
#> 15 | 0.2808101460
#> 16 | 0.2807614114
#> 17 | 0.2807260287
#> 18 | 0.2807003385
#> 19 | 0.2806816849
#> 20 | 0.2806681402
#> 21 | 0.2806583049
#> 22 | 0.2806511631
#> 23 | 0.2806459770
#> 24 | 0.2806422111
#> 25 | 0.2806394765
#> 26 | 0.2806374906
#> 27 | 0.2806360486
#> 28 | 0.2806350014
#> 29 | 0.2806342410
#> 30 | 0.2806336888
#> 31 | 0.2806332878
#> 32 | 0.2806329967
#> 33 | 0.2806327852
#> 34 | 0.2806326317
#> 35 | 0.2806325202
#> 36 | 0.2806324392
#> 37 | 0.2806323804
#> 38 | 0.2806323377
#> 39 | 0.2806323067
#> 40 | 0.2806322842
#> 41 | 0.2806322678
#> 42 | 0.2806322559
#> 43 | 0.2806322473
#> 44 | 0.2806322411
#> 45 | 0.2806322365
#> 46 | 0.2806322332
#> 47 | 0.2806322308
#> 48 | 0.2806322291
#> 49 | 0.2806322278
#> 50 | 0.2806322269
#> 51 | 0.2806322262
#> 52 | 0.2806322257
#> 53 | 0.2806322254
#> 54 | 0.2806322251
#> 55 | 0.2806322249
#> 56 | 0.2806322248
#> 57 | 0.2806322247
#> 58 | 0.2806322246
#> 59 | 0.2806322246
#> 60 | 0.2806322246
#> 61 | 0.2806322245
#> found new candidate with scale 0.2806322 in 61 iter (worst sc=0.28063)
#> Sample[ 7]: idc = 1 11
#> b^[] = -8.08546 3.56245
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -8.08546 3.56245
#> after refine_(*, conv=F):
#> beta_ref : -1.05124 1.24888
#> with ||beta_ref - beta_cand|| = 7.40491897223, --> sc = 1.76350297340318
#> Sample[ 8]: idc = 4 8
#> b^[] = 0.327489 0.842506
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 0.327489 0.842506
#> after refine_(*, conv=F):
#> beta_ref : -0.488056 1.09604
#> with ||beta_ref - beta_cand|| = 0.854043995939, --> sc = 0.406506642125284
#> Sample[ 9]: idc = 7 1
#> b^[] = -13.1628 5.179
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -13.1628 5.179
#> after refine_(*, conv=F):
#> beta_ref : -0.962999 1.2121
#> with ||beta_ref - beta_cand|| = 12.8285446507, --> sc = 2.75529751980209
#> Sample[ 10]: idc = 4 5
#> b^[] = -0.83181 1.18668
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -0.83181 1.18668
#> after refine_(*, conv=F):
#> beta_ref : -0.707317 1.16638
#> with ||beta_ref - beta_cand|| = 0.126137655448, --> sc = 0.307085307609751
#> Sample[ 11]: idc = 6 2
#> b^[] = -1.78547 1.24971
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -1.78547 1.24971
#> after refine_(*, conv=F):
#> beta_ref : -1.44982 1.30483
#> with ||beta_ref - beta_cand|| = 0.340142551919, --> sc = 1.10255399571346
#> Sample[ 12]: idc = 10 3
#> b^[] = -0.579864 1.17238
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -0.579864 1.17238
#> after refine_(*, conv=F):
#> beta_ref : -0.596172 1.15246
#> with ||beta_ref - beta_cand|| = 0.0257405973865, --> sc = 0.30900223531099
#> Sample[ 13]: idc = 5 11
#> b^[] = -1.19798 1.27918
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -1.19798 1.27918
#> after refine_(*, conv=F):
#> beta_ref : -0.898455 1.21902
#> with ||beta_ref - beta_cand|| = 0.305504698816, --> sc = 0.367535728188709
#> Sample[ 14]: idc = 4 7
#> b^[] = -3.13442 1.87028
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -3.13442 1.87028
#> after refine_(*, conv=F):
#> beta_ref : -2.01174 1.5504
#> with ||beta_ref - beta_cand|| = 1.16736173632, --> sc = 0.618257718400311
#> Sample[ 15]: idc = 10 1
#> b^[] = -0.57229 1.17037
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -0.57229 1.17037
#> after refine_(*, conv=F):
#> beta_ref : -0.595111 1.15223
#> with ||beta_ref - beta_cand|| = 0.0291487502429, --> sc = 0.308992048351462
#> Sample[ 16]: idc = 0 1
#> b^[] = 24.2493 -6.73246
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 24.2493 -6.73246
#> after refine_(*, conv=F):
#> beta_ref : -0.505791 1.05334
#> with ||beta_ref - beta_cand|| = 25.9506189658, --> sc = 5.22853378335815
#> Sample[ 17]: idc = 9 7
#> b^[] = 0.815075 0.5672
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 0.815075 0.5672
#> after refine_(*, conv=F):
#> beta_ref : 0.224676 0.802773
#> with ||beta_ref - beta_cand|| = 0.635661711031, --> sc = 0.679685535172079
#> Sample[ 18]: idc = 4 10
#> b^[] = -2.48714 1.67811
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -2.48714 1.67811
#> after refine_(*, conv=F):
#> beta_ref : -2.13634 1.58637
#> with ||beta_ref - beta_cand|| = 0.362596815943, --> sc = 0.436315397353933
#> Sample[ 19]: idc = 5 0
#> b^[] = -6.1579 2.53215
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -6.1579 2.53215
#> after refine_(*, conv=F):
#> beta_ref : -0.763906 1.12366
#> with ||beta_ref - beta_cand|| = 5.57485454363, --> sc = 1.64015945823134
#> Sample[ 20]: idc = 1 0
#> b^[] = 24.2493 -6.73246
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 24.2493 -6.73246
#> after refine_(*, conv=F):
#> beta_ref : -0.505791 1.05334
#> with ||beta_ref - beta_cand|| = 25.9506189658, --> sc = 5.22853378335815
#> Sample[ 21]: idc = 10 0
#> b^[] = -9.17664 3.45192
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -9.17664 3.45192
#> after refine_(*, conv=F):
#> beta_ref : -0.840132 1.1533
#> with ||beta_ref - beta_cand|| = 8.64759632726, --> sc = 2.16009257867085
#> Sample[ 22]: idc = 5 4
#> b^[] = -0.83181 1.18668
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -0.83181 1.18668
#> after refine_(*, conv=F):
#> beta_ref : -0.707317 1.16638
#> with ||beta_ref - beta_cand|| = 0.126137655448, --> sc = 0.307085307609751
#> Sample[ 23]: idc = 6 9
#> b^[] = -2.77944 1.94279
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -2.77944 1.94279
#> after refine_(*, conv=F):
#> beta_ref : -1.85666 1.53242
#> with ||beta_ref - beta_cand|| = 1.00990893711, --> sc = 0.885734481714681
#> Sample[ 24]: idc = 2 11
#> b^[] = 0.931408 0.573269
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 0.931408 0.573269
#> after refine_(*, conv=F):
#> beta_ref : 0.290339 0.794023
#> with ||beta_ref - beta_cand|| = 0.678012796036, --> sc = 0.523721624955058
#> Sample[ 25]: idc = 1 0
#> b^[] = 24.2493 -6.73246
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 24.2493 -6.73246
#> after refine_(*, conv=F):
#> beta_ref : -0.505791 1.05334
#> with ||beta_ref - beta_cand|| = 25.9506189658, --> sc = 5.22853378335815
#> Sample[ 26]: idc = 9 2
#> b^[] = 0.553075 0.667465
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 0.553075 0.667465
#> after refine_(*, conv=F):
#> beta_ref : 0.207049 0.814328
#> with ||beta_ref - beta_cand|| = 0.375903533209, --> sc = 0.577742866999082
#> Sample[ 27]: idc = 10 2
#> b^[] = 13.188 -2.47835
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 13.188 -2.47835
#> after refine_(*, conv=F):
#> beta_ref : -0.333051 1.01284
#> with ||beta_ref - beta_cand|| = 13.9645138921, --> sc = 4.27637683340342
#> Sample[ 28]: idc = 2 1
#> b^[] = 2.63647 0.148748
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 2.63647 0.148748
#> after refine_(*, conv=F):
#> beta_ref : 0.060579 0.901992
#> with ||beta_ref - beta_cand|| = 2.68375998949, --> sc = 0.857121045415407
#> Sample[ 29]: idc = 3 4
#> b^[] = -0.284513 1.0242
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -0.284513 1.0242
#> after refine_(*, conv=F):
#> beta_ref : -0.565247 1.12529
#> with ||beta_ref - beta_cand|| = 0.298379314564, --> sc = 0.280155905387692
#> Sample[ 30]: idc = 7 8
#> b^[] = 8.05358 -1.82104
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 8.05358 -1.82104
#> after refine_(*, conv=F):
#> beta_ref : 0.205065 0.812545
#> with ||beta_ref - beta_cand|| = 8.2785869972, --> sc = 2.06063624179738
#> Sample[ 31]: idc = 3 4
#> b^[] = -0.284513 1.0242
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -0.284513 1.0242
#> after refine_(*, conv=F):
#> beta_ref : -0.565247 1.12529
#> with ||beta_ref - beta_cand|| = 0.298379314564, --> sc = 0.280155905387692
#> Sample[ 32]: idc = 11 2
#> b^[] = 0.931408 0.573269
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 0.931408 0.573269
#> after refine_(*, conv=F):
#> beta_ref : 0.290339 0.794023
#> with ||beta_ref - beta_cand|| = 0.678012796036, --> sc = 0.523721624955058
#> Sample[ 33]: idc = 3 8
#> b^[] = -0.471158 1.11784
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -0.471158 1.11784
#> after refine_(*, conv=F):
#> beta_ref : -0.548769 1.1376
#> with ||beta_ref - beta_cand|| = 0.08008750246, --> sc = 0.213513497916964
#> Sample[ 34]: idc = 1 12
#> b^[] = -0.484845 1.14253
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -0.484845 1.14253
#> after refine_(*, conv=F):
#> beta_ref : -0.587028 1.14935
#> with ||beta_ref - beta_cand|| = 0.102411328665, --> sc = 0.312564686828127
#> find_scale(*, ini.scale =0.31256468683, tol=1e-10):
#> it | new scale
#> 0 | 0.3033762364
#> 1 | 0.2969165179
#> 2 | 0.2923364690
#> 3 | 0.2890682994
#> 4 | 0.2867251665
#> 5 | 0.2850393926
#> 6 | 0.2838234747
#> 7 | 0.2829448327
#> 8 | 0.2823090578
#> 9 | 0.2818485692
#> 10 | 0.2815148029
#> 11 | 0.2812727614
#> 12 | 0.2810971715
#> 13 | 0.2809697545
#> 14 | 0.2808772760
#> 15 | 0.2808101460
#> 16 | 0.2807614114
#> 17 | 0.2807260287
#> 18 | 0.2807003385
#> 19 | 0.2806816849
#> 20 | 0.2806681402
#> 21 | 0.2806583049
#> 22 | 0.2806511631
#> 23 | 0.2806459770
#> 24 | 0.2806422111
#> 25 | 0.2806394765
#> 26 | 0.2806374906
#> 27 | 0.2806360486
#> 28 | 0.2806350014
#> 29 | 0.2806342410
#> 30 | 0.2806336888
#> 31 | 0.2806332878
#> 32 | 0.2806329967
#> 33 | 0.2806327852
#> 34 | 0.2806326317
#> 35 | 0.2806325202
#> 36 | 0.2806324392
#> 37 | 0.2806323804
#> 38 | 0.2806323377
#> 39 | 0.2806323067
#> 40 | 0.2806322842
#> 41 | 0.2806322678
#> 42 | 0.2806322559
#> 43 | 0.2806322473
#> 44 | 0.2806322411
#> 45 | 0.2806322365
#> 46 | 0.2806322332
#> 47 | 0.2806322308
#> 48 | 0.2806322291
#> 49 | 0.2806322278
#> 50 | 0.2806322269
#> 51 | 0.2806322262
#> 52 | 0.2806322257
#> 53 | 0.2806322254
#> 54 | 0.2806322251
#> 55 | 0.2806322249
#> 56 | 0.2806322248
#> 57 | 0.2806322247
#> 58 | 0.2806322246
#> 59 | 0.2806322246
#> 60 | 0.2806322246
#> 61 | 0.2806322245
#> found new candidate with scale 0.2806322 in 61 iter (worst sc=0.28063)
#> Sample[ 35]: idc = 6 0
#> b^[] = -1.6587 1.16132
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -1.6587 1.16132
#> after refine_(*, conv=F):
#> beta_ref : -1.38655 1.28782
#> with ||beta_ref - beta_cand|| = 0.300121734322, --> sc = 1.32303815737056
#> Sample[ 36]: idc = 0 4
#> b^[] = -36.3623 11.735
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -36.3623 11.735
#> after refine_(*, conv=F):
#> beta_ref : -1.03377 1.21864
#> with ||beta_ref - beta_cand|| = 36.8605343085, --> sc = 6.75365155222264
#> Sample[ 37]: idc = 5 12
#> b^[] = 1.03191 0.715867
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 1.03191 0.715867
#> after refine_(*, conv=F):
#> beta_ref : -0.376319 1.0606
#> with ||beta_ref - beta_cand|| = 1.4498103741, --> sc = 0.787749674878659
#> Sample[ 38]: idc = 1 4
#> b^[] = 2.25191 0.271184
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 2.25191 0.271184
#> after refine_(*, conv=F):
#> beta_ref : 0.0105456 0.920859
#> with ||beta_ref - beta_cand|| = 2.33362366183, --> sc = 0.813737418280548
#> Sample[ 39]: idc = 0 3
#> b^[] = 1.14468 0.30717
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 1.14468 0.30717
#> after refine_(*, conv=F):
#> beta_ref : -0.307508 0.964066
#> with ||beta_ref - beta_cand|| = 1.59385397212, --> sc = 1.34390327405485
#> Sample[ 40]: idc = 3 4
#> b^[] = -0.284513 1.0242
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -0.284513 1.0242
#> after refine_(*, conv=F):
#> beta_ref : -0.565247 1.12529
#> with ||beta_ref - beta_cand|| = 0.298379314564, --> sc = 0.280155905387692
#> Sample[ 41]: idc = 6 2
#> b^[] = -1.78547 1.24971
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -1.78547 1.24971
#> after refine_(*, conv=F):
#> beta_ref : -1.44982 1.30483
#> with ||beta_ref - beta_cand|| = 0.340142551919, --> sc = 1.10255399571346
#> Sample[ 42]: idc = 11 6
#> b^[] = -2.39849 1.67716
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -2.39849 1.67716
#> after refine_(*, conv=F):
#> beta_ref : -2.18647 1.6102
#> with ||beta_ref - beta_cand|| = 0.222349315351, --> sc = 0.365943372441649
#> Sample[ 43]: idc = 1 9
#> b^[] = -1.69558 1.52801
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -1.69558 1.52801
#> after refine_(*, conv=F):
#> beta_ref : -1.38063 1.38331
#> with ||beta_ref - beta_cand|| = 0.346599646578, --> sc = 0.467390770067497
#> Sample[ 44]: idc = 5 4
#> b^[] = -0.83181 1.18668
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -0.83181 1.18668
#> after refine_(*, conv=F):
#> beta_ref : -0.707317 1.16638
#> with ||beta_ref - beta_cand|| = 0.126137655448, --> sc = 0.307085307609751
#> Sample[ 45]: idc = 11 7
#> b^[] = 29.1679 -8.78739
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 29.1679 -8.78739
#> after refine_(*, conv=F):
#> beta_ref : -0.244111 0.966143
#> with ||beta_ref - beta_cand|| = 30.9870463767, --> sc = 6.48063687575041
#> Sample[ 46]: idc = 0 3
#> b^[] = 1.14468 0.30717
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 1.14468 0.30717
#> after refine_(*, conv=F):
#> beta_ref : -0.307508 0.964066
#> with ||beta_ref - beta_cand|| = 1.59385397212, --> sc = 1.34390327405485
#> Sample[ 47]: idc = 1 7
#> b^[] = -13.1628 5.179
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -13.1628 5.179
#> after refine_(*, conv=F):
#> beta_ref : -0.962999 1.2121
#> with ||beta_ref - beta_cand|| = 12.8285446507, --> sc = 2.75529751980209
#> Sample[ 48]: idc = 0 10
#> b^[] = -9.17664 3.45192
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -9.17664 3.45192
#> after refine_(*, conv=F):
#> beta_ref : -0.840132 1.1533
#> with ||beta_ref - beta_cand|| = 8.64759632726, --> sc = 2.16009257867086
#> Sample[ 49]: idc = 11 5
#> b^[] = -1.19798 1.27918
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -1.19798 1.27918
#> after refine_(*, conv=F):
#> beta_ref : -0.898455 1.21902
#> with ||beta_ref - beta_cand|| = 0.305504698816, --> sc = 0.367535728188709
#> Sample[ 50]: idc = 1 6
#> b^[] = -2.59549 1.81452
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -2.59549 1.81452
#> after refine_(*, conv=F):
#> beta_ref : -2.14459 1.61896
#> with ||beta_ref - beta_cand|| = 0.491482052855, --> sc = 0.553607074103745
#> Sample[ 51]: idc = 2 8
#> b^[] = 1.56876 0.414582
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 1.56876 0.414582
#> after refine_(*, conv=F):
#> beta_ref : 0.260798 0.824069
#> with ||beta_ref - beta_cand|| = 1.37056338354, --> sc = 0.594013641561037
#> Sample[ 52]: idc = 7 11
#> b^[] = 29.1679 -8.78739
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 29.1679 -8.78739
#> after refine_(*, conv=F):
#> beta_ref : -0.244111 0.966143
#> with ||beta_ref - beta_cand|| = 30.9870463767, --> sc = 6.48063687575061
#> Sample[ 53]: idc = 3 12
#> b^[] = -0.565704 1.16527
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -0.565704 1.16527
#> after refine_(*, conv=F):
#> beta_ref : -0.589765 1.1512
#> with ||beta_ref - beta_cand|| = 0.0278742526664, --> sc = 0.289862591554884
#> Sample[ 54]: idc = 11 12
#> b^[] = -2.47152 1.70137
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -2.47152 1.70137
#> after refine_(*, conv=F):
#> beta_ref : -2.20041 1.61542
#> with ||beta_ref - beta_cand|| = 0.284411253958, --> sc = 0.366625751619072
#> Sample[ 55]: idc = 0 12
#> b^[] = -14.974 5.21828
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -14.974 5.21828
#> after refine_(*, conv=F):
#> beta_ref : -1.04426 1.21824
#> with ||beta_ref - beta_cand|| = 14.4926713752, --> sc = 2.94185072560702
#> Sample[ 56]: idc = 12 10
#> b^[] = -0.773295 1.22367
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -0.773295 1.22367
#> after refine_(*, conv=F):
#> beta_ref : -0.62546 1.16042
#> with ||beta_ref - beta_cand|| = 0.160794576384, --> sc = 0.287874239981522
#> find_scale(*, ini.scale =0.28787423998, tol=1e-10):
#> it | new scale
#> 0 | 0.2855822164
#> 1 | 0.2839452738
#> 2 | 0.2827733665
#> 3 | 0.2819329210
#> 4 | 0.2813294280
#> 5 | 0.2808956885
#> 6 | 0.2805837492
#> 7 | 0.2803593010
#> 8 | 0.2801977500
#> 9 | 0.2800814420
#> 10 | 0.2799976917
#> 11 | 0.2799373777
#> 12 | 0.2798939377
#> 13 | 0.2798626488
#> 14 | 0.2798401110
#> 15 | 0.2798238762
#> 16 | 0.2798121813
#> 17 | 0.2798037568
#> 18 | 0.2797976879
#> 19 | 0.2797933160
#> 20 | 0.2797901665
#> 21 | 0.2797878976
#> 22 | 0.2797862632
#> 23 | 0.2797850857
#> 24 | 0.2797842374
#> 25 | 0.2797836264
#> 26 | 0.2797831862
#> 27 | 0.2797828690
#> 28 | 0.2797826406
#> 29 | 0.2797824760
#> 30 | 0.2797823574
#> 31 | 0.2797822720
#> 32 | 0.2797822105
#> 33 | 0.2797821661
#> 34 | 0.2797821342
#> 35 | 0.2797821112
#> 36 | 0.2797820946
#> 37 | 0.2797820827
#> 38 | 0.2797820741
#> 39 | 0.2797820679
#> 40 | 0.2797820634
#> 41 | 0.2797820602
#> 42 | 0.2797820579
#> 43 | 0.2797820562
#> 44 | 0.2797820550
#> 45 | 0.2797820541
#> 46 | 0.2797820535
#> 47 | 0.2797820531
#> 48 | 0.2797820528
#> 49 | 0.2797820525
#> 50 | 0.2797820523
#> 51 | 0.2797820522
#> 52 | 0.2797820521
#> 53 | 0.2797820521
#> 54 | 0.2797820520
#> 55 | 0.2797820520
#> 56 | 0.2797820520
#> found new candidate with scale 0.2797821 in 56 iter (worst sc=0.28063)
#> Sample[ 57]: idc = 0 12
#> b^[] = -14.974 5.21828
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -14.974 5.21828
#> after refine_(*, conv=F):
#> beta_ref : -1.04426 1.21824
#> with ||beta_ref - beta_cand|| = 14.4926713752, --> sc = 2.94185072560702
#> Sample[ 58]: idc = 8 6
#> b^[] = -2.69667 1.88507
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -2.69667 1.88507
#> after refine_(*, conv=F):
#> beta_ref : -1.94237 1.55803
#> with ||beta_ref - beta_cand|| = 0.822144059656, --> sc = 0.744918088188394
#> Sample[ 59]: idc = 8 6
#> b^[] = -2.69667 1.88507
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -2.69667 1.88507
#> after refine_(*, conv=F):
#> beta_ref : -1.94237 1.55803
#> with ||beta_ref - beta_cand|| = 0.822144059656, --> sc = 0.744918088188394
#> Sample[ 60]: idc = 11 1
#> b^[] = -8.08546 3.56245
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -8.08546 3.56245
#> after refine_(*, conv=F):
#> beta_ref : -1.05124 1.24888
#> with ||beta_ref - beta_cand|| = 7.40491897223, --> sc = 1.76350297340317
#> Sample[ 61]: idc = 3 4
#> b^[] = -0.284513 1.0242
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -0.284513 1.0242
#> after refine_(*, conv=F):
#> beta_ref : -0.565247 1.12529
#> with ||beta_ref - beta_cand|| = 0.298379314564, --> sc = 0.280155905387692
#> Sample[ 62]: idc = 2 0
#> b^[] = -2.67885 1.47214
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -2.67885 1.47214
#> after refine_(*, conv=F):
#> beta_ref : -1.4879 1.31816
#> with ||beta_ref - beta_cand|| = 1.20086135032, --> sc = 1.37247806701654
#> Sample[ 63]: idc = 7 3
#> b^[] = 0.263938 0.74904
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 0.263938 0.74904
#> after refine_(*, conv=F):
#> beta_ref : 0.145859 0.83463
#> with ||beta_ref - beta_cand|| = 0.145837154087, --> sc = 0.597892842089404
#> Sample[ 64]: idc = 2 4
#> b^[] = 2.8091 0.105765
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 2.8091 0.105765
#> after refine_(*, conv=F):
#> beta_ref : 0.030844 0.911664
#> with ||beta_ref - beta_cand|| = 2.89278369939, --> sc = 0.908806025778877
#> Sample[ 65]: idc = 3 1
#> b^[] = -0.582063 1.17348
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -0.582063 1.17348
#> after refine_(*, conv=F):
#> beta_ref : -0.597198 1.15265
#> with ||beta_ref - beta_cand|| = 0.0257497340318, --> sc = 0.312109152133832
#> Sample[ 66]: idc = 9 0
#> b^[] = 2.86111 -0.215798
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 2.86111 -0.215798
#> after refine_(*, conv=F):
#> beta_ref : 0.171545 0.81517
#> with ||beta_ref - beta_cand|| = 2.88038890444, --> sc = 1.45332402557657
#> Sample[ 67]: idc = 0 3
#> b^[] = 1.14468 0.30717
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 1.14468 0.30717
#> after refine_(*, conv=F):
#> beta_ref : -0.307508 0.964066
#> with ||beta_ref - beta_cand|| = 1.59385397212, --> sc = 1.34390327405485
#> Sample[ 68]: idc = 2 6
#> b^[] = -1.78547 1.24971
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -1.78547 1.24971
#> after refine_(*, conv=F):
#> beta_ref : -1.44982 1.30483
#> with ||beta_ref - beta_cand|| = 0.340142551919, --> sc = 1.10255399571346
#> Sample[ 69]: idc = 7 10
#> b^[] = -2.81746 1.7657
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -2.81746 1.7657
#> after refine_(*, conv=F):
#> beta_ref : -2.10053 1.5736
#> with ||beta_ref - beta_cand|| = 0.742227256343, --> sc = 0.492521121454336
#> Sample[ 70]: idc = 8 9
#> b^[] = -2.01045 1.6485
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -2.01045 1.6485
#> after refine_(*, conv=F):
#> beta_ref : -1.74249 1.49589
#> with ||beta_ref - beta_cand|| = 0.308377559274, --> sc = 0.579998989395441
#> Sample[ 71]: idc = 4 10
#> b^[] = -2.48714 1.67811
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -2.48714 1.67811
#> after refine_(*, conv=F):
#> beta_ref : -2.13634 1.58637
#> with ||beta_ref - beta_cand|| = 0.362596815943, --> sc = 0.436315397353933
#> Sample[ 72]: idc = 5 2
#> b^[] = 47.0534 -10.9101
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 47.0534 -10.9101
#> after refine_(*, conv=F):
#> beta_ref : -0.397375 1.02753
#> with ||beta_ref - beta_cand|| = 48.929345941, --> sc = 15.1481999432314
#> Sample[ 73]: idc = 3 2
#> b^[] = 0.301869 0.73001
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 0.301869 0.73001
#> after refine_(*, conv=F):
#> beta_ref : 0.14998 0.831077
#> with ||beta_ref - beta_cand|| = 0.182440710486, --> sc = 0.631676782286114
#> Sample[ 74]: idc = 6 8
#> b^[] = -2.69667 1.88507
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -2.69667 1.88507
#> after refine_(*, conv=F):
#> beta_ref : -1.94237 1.55803
#> with ||beta_ref - beta_cand|| = 0.822144059656, --> sc = 0.744918088188393
#> Sample[ 75]: idc = 6 7
#> b^[] = -2.26322 1.58284
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -2.26322 1.58284
#> after refine_(*, conv=F):
#> beta_ref : -2.18436 1.58611
#> with ||beta_ref - beta_cand|| = 0.078929299088, --> sc = 0.361556853894617
#> Sample[ 76]: idc = 10 12
#> b^[] = -0.773295 1.22367
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -0.773295 1.22367
#> after refine_(*, conv=F):
#> beta_ref : -0.62546 1.16042
#> with ||beta_ref - beta_cand|| = 0.160794576384, --> sc = 0.287874239981522
#> find_scale(*, ini.scale =0.28787423998, tol=1e-10):
#> it | new scale
#> 0 | 0.2855822164
#> 1 | 0.2839452738
#> 2 | 0.2827733665
#> 3 | 0.2819329210
#> 4 | 0.2813294280
#> 5 | 0.2808956885
#> 6 | 0.2805837492
#> 7 | 0.2803593010
#> 8 | 0.2801977500
#> 9 | 0.2800814420
#> 10 | 0.2799976917
#> 11 | 0.2799373777
#> 12 | 0.2798939377
#> 13 | 0.2798626488
#> 14 | 0.2798401110
#> 15 | 0.2798238762
#> 16 | 0.2798121813
#> 17 | 0.2798037568
#> 18 | 0.2797976879
#> 19 | 0.2797933160
#> 20 | 0.2797901665
#> 21 | 0.2797878976
#> 22 | 0.2797862632
#> 23 | 0.2797850857
#> 24 | 0.2797842374
#> 25 | 0.2797836264
#> 26 | 0.2797831862
#> 27 | 0.2797828690
#> 28 | 0.2797826406
#> 29 | 0.2797824760
#> 30 | 0.2797823574
#> 31 | 0.2797822720
#> 32 | 0.2797822105
#> 33 | 0.2797821661
#> 34 | 0.2797821342
#> 35 | 0.2797821112
#> 36 | 0.2797820946
#> 37 | 0.2797820827
#> 38 | 0.2797820741
#> 39 | 0.2797820679
#> 40 | 0.2797820634
#> 41 | 0.2797820602
#> 42 | 0.2797820579
#> 43 | 0.2797820562
#> 44 | 0.2797820550
#> 45 | 0.2797820541
#> 46 | 0.2797820535
#> 47 | 0.2797820531
#> 48 | 0.2797820528
#> 49 | 0.2797820525
#> 50 | 0.2797820523
#> 51 | 0.2797820522
#> 52 | 0.2797820521
#> 53 | 0.2797820521
#> 54 | 0.2797820520
#> 55 | 0.2797820520
#> 56 | 0.2797820520
#> found new candidate with scale 0.2797821 in 56 iter (worst sc=0.27978)
#> Sample[ 77]: idc = 9 12
#> b^[] = -1.25272 1.35853
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -1.25272 1.35853
#> after refine_(*, conv=F):
#> beta_ref : -0.819674 1.21313
#> with ||beta_ref - beta_cand|| = 0.456802788525, --> sc = 0.301200329328142
#> Sample[ 78]: idc = 9 6
#> b^[] = -2.77944 1.94279
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -2.77944 1.94279
#> after refine_(*, conv=F):
#> beta_ref : -1.85666 1.53242
#> with ||beta_ref - beta_cand|| = 1.00990893711, --> sc = 0.885734481714681
#> Sample[ 79]: idc = 7 0
#> b^[] = 7.13585 -1.51825
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 7.13585 -1.51825
#> after refine_(*, conv=F):
#> beta_ref : 0.224305 0.805575
#> with ||beta_ref - beta_cand|| = 7.29174469856, --> sc = 1.88768873167888
#> Sample[ 80]: idc = 11 7
#> b^[] = 29.1679 -8.78739
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 29.1679 -8.78739
#> after refine_(*, conv=F):
#> beta_ref : -0.244111 0.966143
#> with ||beta_ref - beta_cand|| = 30.9870463767, --> sc = 6.48063687575041
#> Sample[ 81]: idc = 5 8
#> b^[] = -0.229469 1.03452
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -0.229469 1.03452
#> after refine_(*, conv=F):
#> beta_ref : -0.556193 1.13403
#> with ||beta_ref - beta_cand|| = 0.341541122565, --> sc = 0.270706285978508
#> find_scale(*, ini.scale =0.27070628598, tol=1e-10):
#> it | new scale
#> 0 | 0.2724486979
#> 1 | 0.2736969547
#> 2 | 0.2745894824
#> 3 | 0.2752267857
#> 4 | 0.2756814065
#> 5 | 0.2760054871
#> 6 | 0.2762363979
#> 7 | 0.2764008668
#> 8 | 0.2765179827
#> 9 | 0.2766013645
#> 10 | 0.2766607216
#> 11 | 0.2767029724
#> 12 | 0.2767330449
#> 13 | 0.2767544485
#> 14 | 0.2767696815
#> 15 | 0.2767805227
#> 16 | 0.2767882382
#> 17 | 0.2767937290
#> 18 | 0.2767976367
#> 19 | 0.2768004176
#> 20 | 0.2768023967
#> 21 | 0.2768038051
#> 22 | 0.2768048074
#> 23 | 0.2768055207
#> 24 | 0.2768060283
#> 25 | 0.2768063896
#> 26 | 0.2768066467
#> 27 | 0.2768068296
#> 28 | 0.2768069598
#> 29 | 0.2768070525
#> 30 | 0.2768071184
#> 31 | 0.2768071653
#> 32 | 0.2768071987
#> 33 | 0.2768072225
#> 34 | 0.2768072394
#> 35 | 0.2768072515
#> 36 | 0.2768072600
#> 37 | 0.2768072661
#> 38 | 0.2768072705
#> 39 | 0.2768072735
#> 40 | 0.2768072757
#> 41 | 0.2768072773
#> 42 | 0.2768072784
#> 43 | 0.2768072792
#> 44 | 0.2768072798
#> 45 | 0.2768072802
#> 46 | 0.2768072805
#> 47 | 0.2768072807
#> 48 | 0.2768072808
#> 49 | 0.2768072809
#> 50 | 0.2768072810
#> 51 | 0.2768072810
#> 52 | 0.2768072811
#> 53 | 0.2768072811
#> found new candidate with scale 0.2768073 in 53 iter (worst sc=0.27978)
#> Sample[ 82]: idc = 12 4
#> b^[] = -4.26167 2.20494
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -4.26167 2.20494
#> after refine_(*, conv=F):
#> beta_ref : -1.65432 1.43182
#> with ||beta_ref - beta_cand|| = 2.71955618857, --> sc = 0.91889370623179
#> Sample[ 83]: idc = 5 12
#> b^[] = 1.03191 0.715867
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 1.03191 0.715867
#> after refine_(*, conv=F):
#> beta_ref : -0.376319 1.0606
#> with ||beta_ref - beta_cand|| = 1.4498103741, --> sc = 0.787749674878659
#> Sample[ 84]: idc = 3 6
#> b^[] = -4.48261 3.13038
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -4.48261 3.13038
#> after refine_(*, conv=F):
#> beta_ref : -1.5156 1.39333
#> with ||beta_ref - beta_cand|| = 3.43809929004, --> sc = 3.95137685252739
#> Sample[ 85]: idc = 9 0
#> b^[] = 2.86111 -0.215798
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 2.86111 -0.215798
#> after refine_(*, conv=F):
#> beta_ref : 0.171545 0.81517
#> with ||beta_ref - beta_cand|| = 2.88038890444, --> sc = 1.45332402557657
#> Sample[ 86]: idc = 1 2
#> b^[] = 2.63647 0.148748
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 2.63647 0.148748
#> after refine_(*, conv=F):
#> beta_ref : 0.060579 0.901992
#> with ||beta_ref - beta_cand|| = 2.68375998949, --> sc = 0.857121045415408
#> Sample[ 87]: idc = 0 4
#> b^[] = -36.3623 11.735
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -36.3623 11.735
#> after refine_(*, conv=F):
#> beta_ref : -1.03377 1.21864
#> with ||beta_ref - beta_cand|| = 36.8605343085, --> sc = 6.75365155222264
#> Sample[ 88]: idc = 0 9
#> b^[] = 2.86111 -0.215798
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 2.86111 -0.215798
#> after refine_(*, conv=F):
#> beta_ref : 0.171545 0.81517
#> with ||beta_ref - beta_cand|| = 2.88038890444, --> sc = 1.45332402557657
#> Sample[ 89]: idc = 4 9
#> b^[] = -0.706312 1.14942
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -0.706312 1.14942
#> after refine_(*, conv=F):
#> beta_ref : -0.662089 1.15324
#> with ||beta_ref - beta_cand|| = 0.0443876719904, --> sc = 0.302818894470935
#> Sample[ 90]: idc = 0 8
#> b^[] = 7.47514 -1.62162
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 7.47514 -1.62162
#> after refine_(*, conv=F):
#> beta_ref : 0.245606 0.798807
#> with ||beta_ref - beta_cand|| = 7.62395703783, --> sc = 1.90882061471928
#> Sample[ 91]: idc = 12 8
#> b^[] = -0.799579 1.23106
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -0.799579 1.23106
#> after refine_(*, conv=F):
#> beta_ref : -0.631179 1.16203
#> with ||beta_ref - beta_cand|| = 0.181998057096, --> sc = 0.286399245856023
#> find_scale(*, ini.scale =0.28639924586, tol=1e-10):
#> it | new scale
#> 0 | 0.2844989854
#> 1 | 0.2831418046
#> 2 | 0.2821705552
#> 3 | 0.2814744887
#> 4 | 0.2809751185
#> 5 | 0.2806165930
#> 6 | 0.2803590489
#> 7 | 0.2801739720
#> 8 | 0.2800409345
#> 9 | 0.2799452850
#> 10 | 0.2798765061
#> 11 | 0.2798270441
#> 12 | 0.2797914710
#> 13 | 0.2797658854
#> 14 | 0.2797474826
#> 15 | 0.2797342457
#> 16 | 0.2797247244
#> 17 | 0.2797178756
#> 18 | 0.2797129491
#> 19 | 0.2797094054
#> 20 | 0.2797068563
#> 21 | 0.2797050227
#> 22 | 0.2797037036
#> 23 | 0.2797027548
#> 24 | 0.2797020723
#> 25 | 0.2797015814
#> 26 | 0.2797012282
#> 27 | 0.2797009742
#> 28 | 0.2797007914
#> 29 | 0.2797006600
#> 30 | 0.2797005654
#> 31 | 0.2797004974
#> 32 | 0.2797004485
#> 33 | 0.2797004133
#> 34 | 0.2797003880
#> 35 | 0.2797003698
#> 36 | 0.2797003567
#> 37 | 0.2797003472
#> 38 | 0.2797003404
#> 39 | 0.2797003356
#> 40 | 0.2797003321
#> 41 | 0.2797003295
#> 42 | 0.2797003277
#> 43 | 0.2797003264
#> 44 | 0.2797003255
#> 45 | 0.2797003248
#> 46 | 0.2797003243
#> 47 | 0.2797003240
#> 48 | 0.2797003237
#> 49 | 0.2797003235
#> 50 | 0.2797003234
#> 51 | 0.2797003233
#> 52 | 0.2797003232
#> 53 | 0.2797003232
#> 54 | 0.2797003232
#> 55 | 0.2797003231
#> found new candidate with scale 0.2797003 in 55 iter (worst sc=0.2797)
#> Sample[ 92]: idc = 6 10
#> b^[] = -2.3463 1.64077
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -2.3463 1.64077
#> after refine_(*, conv=F):
#> beta_ref : -2.13209 1.58647
#> with ||beta_ref - beta_cand|| = 0.220987541723, --> sc = 0.417145691249323
#> Sample[ 93]: idc = 7 1
#> b^[] = -13.1628 5.179
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -13.1628 5.179
#> after refine_(*, conv=F):
#> beta_ref : -0.962999 1.2121
#> with ||beta_ref - beta_cand|| = 12.8285446507, --> sc = 2.75529751980209
#> Sample[ 94]: idc = 10 7
#> b^[] = -2.81746 1.7657
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -2.81746 1.7657
#> after refine_(*, conv=F):
#> beta_ref : -2.10053 1.5736
#> with ||beta_ref - beta_cand|| = 0.742227256343, --> sc = 0.492521121454337
#> Sample[ 95]: idc = 11 9
#> b^[] = -0.0569989 0.900935
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -0.0569989 0.900935
#> after refine_(*, conv=F):
#> beta_ref : -0.177777 0.971152
#> with ||beta_ref - beta_cand|| = 0.139706152895, --> sc = 0.486048896285347
#> Sample[ 96]: idc = 4 6
#> b^[] = -2.33514 1.63299
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -2.33514 1.63299
#> after refine_(*, conv=F):
#> beta_ref : -2.14415 1.58792
#> with ||beta_ref - beta_cand|| = 0.196236405624, --> sc = 0.399338780038278
#> Sample[ 97]: idc = 4 9
#> b^[] = -0.706312 1.14942
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -0.706312 1.14942
#> after refine_(*, conv=F):
#> beta_ref : -0.662089 1.15324
#> with ||beta_ref - beta_cand|| = 0.0443876719904, --> sc = 0.302818894470935
#> Sample[ 98]: idc = 9 11
#> b^[] = -0.0569989 0.900935
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -0.0569989 0.900935
#> after refine_(*, conv=F):
#> beta_ref : -0.177777 0.971152
#> with ||beta_ref - beta_cand|| = 0.139706152895, --> sc = 0.486048896285346
#> Sample[ 99]: idc = 4 2
#> b^[] = 2.8091 0.105765
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 2.8091 0.105765
#> after refine_(*, conv=F):
#> beta_ref : 0.030844 0.911664
#> with ||beta_ref - beta_cand|| = 2.89278369939, --> sc = 0.908806025778876
#> Sample[100]: idc = 12 11
#> b^[] = -2.47152 1.70137
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -2.47152 1.70137
#> after refine_(*, conv=F):
#> beta_ref : -2.20041 1.61542
#> with ||beta_ref - beta_cand|| = 0.284411253958, --> sc = 0.366625751619073
#> Sample[101]: idc = 6 11
#> b^[] = -2.39849 1.67716
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -2.39849 1.67716
#> after refine_(*, conv=F):
#> beta_ref : -2.18647 1.6102
#> with ||beta_ref - beta_cand|| = 0.222349315351, --> sc = 0.36594337244165
#> Sample[102]: idc = 5 9
#> b^[] = -0.749016 1.16576
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -0.749016 1.16576
#> after refine_(*, conv=F):
#> beta_ref : -0.670506 1.15677
#> with ||beta_ref - beta_cand|| = 0.0790237477718, --> sc = 0.28633536578424
#> Sample[103]: idc = 6 10
#> b^[] = -2.3463 1.64077
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -2.3463 1.64077
#> after refine_(*, conv=F):
#> beta_ref : -2.13209 1.58647
#> with ||beta_ref - beta_cand|| = 0.220987541723, --> sc = 0.417145691249323
#> Sample[104]: idc = 3 7
#> b^[] = 0.263938 0.74904
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 0.263938 0.74904
#> after refine_(*, conv=F):
#> beta_ref : 0.145859 0.83463
#> with ||beta_ref - beta_cand|| = 0.145837154087, --> sc = 0.597892842089404
#> Sample[105]: idc = 3 11
#> b^[] = -0.00341923 0.883173
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -0.00341923 0.883173
#> after refine_(*, conv=F):
#> beta_ref : -0.142669 0.957455
#> with ||beta_ref - beta_cand|| = 0.157824084971, --> sc = 0.473157356422567
#> Sample[106]: idc = 6 5
#> b^[] = -2.18551 1.52865
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -2.18551 1.52865
#> after refine_(*, conv=F):
#> beta_ref : -2.06748 1.54965
#> with ||beta_ref - beta_cand|| = 0.119878917294, --> sc = 0.539116852816471
#> Sample[107]: idc = 0 7
#> b^[] = 7.13585 -1.51825
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 7.13585 -1.51825
#> after refine_(*, conv=F):
#> beta_ref : 0.224305 0.805575
#> with ||beta_ref - beta_cand|| = 7.29174469856, --> sc = 1.88768873167888
#> Sample[108]: idc = 2 1
#> b^[] = 2.63647 0.148748
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 2.63647 0.148748
#> after refine_(*, conv=F):
#> beta_ref : 0.060579 0.901992
#> with ||beta_ref - beta_cand|| = 2.68375998949, --> sc = 0.857121045415407
#> Sample[109]: idc = 7 12
#> b^[] = -3.49557 1.98944
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -3.49557 1.98944
#> after refine_(*, conv=F):
#> beta_ref : -1.92479 1.5228
#> with ||beta_ref - beta_cand|| = 1.63862904714, --> sc = 0.741363565741794
#> Sample[110]: idc = 3 2
#> b^[] = 0.301869 0.73001
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 0.301869 0.73001
#> after refine_(*, conv=F):
#> beta_ref : 0.14998 0.831077
#> with ||beta_ref - beta_cand|| = 0.182440710486, --> sc = 0.631676782286114
#> Sample[111]: idc = 9 8
#> b^[] = -2.01045 1.6485
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -2.01045 1.6485
#> after refine_(*, conv=F):
#> beta_ref : -1.74249 1.49589
#> with ||beta_ref - beta_cand|| = 0.308377559274, --> sc = 0.579998989395441
#> Sample[112]: idc = 12 9
#> b^[] = -1.25272 1.35853
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -1.25272 1.35853
#> after refine_(*, conv=F):
#> beta_ref : -0.819674 1.21313
#> with ||beta_ref - beta_cand|| = 0.456802788525, --> sc = 0.301200329328143
#> Sample[113]: idc = 12 2
#> b^[] = 6.21824 -0.743034
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 6.21824 -0.743034
#> after refine_(*, conv=F):
#> beta_ref : -0.199361 0.979666
#> with ||beta_ref - beta_cand|| = 6.64479808942, --> sc = 1.90760566523583
#> Sample[114]: idc = 8 2
#> b^[] = 1.56876 0.414582
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 1.56876 0.414582
#> after refine_(*, conv=F):
#> beta_ref : 0.260798 0.824069
#> with ||beta_ref - beta_cand|| = 1.37056338354, --> sc = 0.594013641561037
#> Sample[115]: idc = 1 6
#> b^[] = -2.59549 1.81452
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -2.59549 1.81452
#> after refine_(*, conv=F):
#> beta_ref : -2.14459 1.61896
#> with ||beta_ref - beta_cand|| = 0.491482052855, --> sc = 0.553607074103745
#> Sample[116]: idc = 8 6
#> b^[] = -2.69667 1.88507
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -2.69667 1.88507
#> after refine_(*, conv=F):
#> beta_ref : -1.94237 1.55803
#> with ||beta_ref - beta_cand|| = 0.822144059656, --> sc = 0.744918088188394
#> Sample[117]: idc = 6 9
#> b^[] = -2.77944 1.94279
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -2.77944 1.94279
#> after refine_(*, conv=F):
#> beta_ref : -1.85666 1.53242
#> with ||beta_ref - beta_cand|| = 1.00990893711, --> sc = 0.885734481714681
#> Sample[118]: idc = 6 5
#> b^[] = -2.18551 1.52865
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -2.18551 1.52865
#> after refine_(*, conv=F):
#> beta_ref : -2.06748 1.54965
#> with ||beta_ref - beta_cand|| = 0.119878917294, --> sc = 0.539116852816471
#> Sample[119]: idc = 0 11
#> b^[] = 8.42881 -1.91219
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 8.42881 -1.91219
#> after refine_(*, conv=F):
#> beta_ref : 0.302657 0.780686
#> with ||beta_ref - beta_cand|| = 8.56072293339, --> sc = 1.9669857749306
#> Sample[120]: idc = 11 5
#> b^[] = -1.19798 1.27918
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -1.19798 1.27918
#> after refine_(*, conv=F):
#> beta_ref : -0.898455 1.21902
#> with ||beta_ref - beta_cand|| = 0.305504698816, --> sc = 0.367535728188709
#> Sample[121]: idc = 4 1
#> b^[] = 2.25191 0.271184
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 2.25191 0.271184
#> after refine_(*, conv=F):
#> beta_ref : 0.0105456 0.920859
#> with ||beta_ref - beta_cand|| = 2.33362366183, --> sc = 0.813737418280548
#> Sample[122]: idc = 0 11
#> b^[] = 8.42881 -1.91219
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 8.42881 -1.91219
#> after refine_(*, conv=F):
#> beta_ref : 0.302657 0.780686
#> with ||beta_ref - beta_cand|| = 8.56072293339, --> sc = 1.9669857749306
#> Sample[123]: idc = 4 7
#> b^[] = -3.13442 1.87028
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -3.13442 1.87028
#> after refine_(*, conv=F):
#> beta_ref : -2.01174 1.5504
#> with ||beta_ref - beta_cand|| = 1.16736173632, --> sc = 0.618257718400311
#> Sample[124]: idc = 7 6
#> b^[] = -2.26322 1.58284
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -2.26322 1.58284
#> after refine_(*, conv=F):
#> beta_ref : -2.18436 1.58611
#> with ||beta_ref - beta_cand|| = 0.078929299088, --> sc = 0.361556853894617
#> Sample[125]: idc = 10 11
#> b^[] = -2.05854 1.56446
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -2.05854 1.56446
#> after refine_(*, conv=F):
#> beta_ref : -2.07016 1.57129
#> with ||beta_ref - beta_cand|| = 0.0134780521353, --> sc = 0.37327449031415
#> Sample[126]: idc = 1 3
#> b^[] = -0.582063 1.17348
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -0.582063 1.17348
#> after refine_(*, conv=F):
#> beta_ref : -0.597198 1.15265
#> with ||beta_ref - beta_cand|| = 0.0257497340318, --> sc = 0.312109152133831
#> Sample[127]: idc = 4 12
#> b^[] = -4.26167 2.20494
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -4.26167 2.20494
#> after refine_(*, conv=F):
#> beta_ref : -1.65432 1.43182
#> with ||beta_ref - beta_cand|| = 2.71955618857, --> sc = 0.918893706231788
#> Sample[128]: idc = 5 6
#> b^[] = -2.18551 1.52865
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -2.18551 1.52865
#> after refine_(*, conv=F):
#> beta_ref : -2.06748 1.54965
#> with ||beta_ref - beta_cand|| = 0.119878917294, --> sc = 0.539116852816472
#> Sample[129]: idc = 10 0
#> b^[] = -9.17664 3.45192
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -9.17664 3.45192
#> after refine_(*, conv=F):
#> beta_ref : -0.840132 1.1533
#> with ||beta_ref - beta_cand|| = 8.64759632726, --> sc = 2.16009257867085
#> Sample[130]: idc = 7 9
#> b^[] = 0.815075 0.5672
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 0.815075 0.5672
#> after refine_(*, conv=F):
#> beta_ref : 0.224676 0.802773
#> with ||beta_ref - beta_cand|| = 0.635661711031, --> sc = 0.67968553517208
#> Sample[131]: idc = 2 9
#> b^[] = 0.553075 0.667465
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 0.553075 0.667465
#> after refine_(*, conv=F):
#> beta_ref : 0.207049 0.814328
#> with ||beta_ref - beta_cand|| = 0.375903533209, --> sc = 0.577742866999083
#> Sample[132]: idc = 10 9
#> b^[] = -1.1869 1.33334
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -1.1869 1.33334
#> after refine_(*, conv=F):
#> beta_ref : -0.809481 1.2096
#> with ||beta_ref - beta_cand|| = 0.39718852225, --> sc = 0.284340124502365
#> Sample[133]: idc = 5 4
#> b^[] = -0.83181 1.18668
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -0.83181 1.18668
#> after refine_(*, conv=F):
#> beta_ref : -0.707317 1.16638
#> with ||beta_ref - beta_cand|| = 0.126137655448, --> sc = 0.307085307609751
#> Sample[134]: idc = 9 5
#> b^[] = -0.749016 1.16576
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -0.749016 1.16576
#> after refine_(*, conv=F):
#> beta_ref : -0.670506 1.15677
#> with ||beta_ref - beta_cand|| = 0.0790237477718, --> sc = 0.28633536578424
#> Sample[135]: idc = 10 7
#> b^[] = -2.81746 1.7657
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -2.81746 1.7657
#> after refine_(*, conv=F):
#> beta_ref : -2.10053 1.5736
#> with ||beta_ref - beta_cand|| = 0.742227256343, --> sc = 0.492521121454337
#> Sample[136]: idc = 10 12
#> b^[] = -0.773295 1.22367
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -0.773295 1.22367
#> after refine_(*, conv=F):
#> beta_ref : -0.62546 1.16042
#> with ||beta_ref - beta_cand|| = 0.160794576384, --> sc = 0.287874239981522
#> Sample[137]: idc = 2 7
#> b^[] = 0.382345 0.709973
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 0.382345 0.709973
#> after refine_(*, conv=F):
#> beta_ref : 0.167816 0.825709
#> with ||beta_ref - beta_cand|| = 0.243757115723, --> sc = 0.613778594777883
#> Sample[138]: idc = 5 9
#> b^[] = -0.749016 1.16576
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -0.749016 1.16576
#> after refine_(*, conv=F):
#> beta_ref : -0.670506 1.15677
#> with ||beta_ref - beta_cand|| = 0.0790237477718, --> sc = 0.28633536578424
#> Sample[139]: idc = 3 5
#> b^[] = -0.381765 1.07299
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -0.381765 1.07299
#> after refine_(*, conv=F):
#> beta_ref : -0.547363 1.12925
#> with ||beta_ref - beta_cand|| = 0.17489423735, --> sc = 0.235441204593026
#> find_scale(*, ini.scale =0.23544120459, tol=1e-10):
#> it | new scale
#> 0 | 0.2468443051
#> 1 | 0.2552730356
#> 2 | 0.2614232664
#> 3 | 0.2658666606
#> 4 | 0.2690540256
#> 5 | 0.2713288789
#> 6 | 0.2729467122
#> 7 | 0.2740944264
#> 8 | 0.2749072117
#> 9 | 0.2754821033
#> 10 | 0.2758883802
#> 11 | 0.2761753230
#> 12 | 0.2763778968
#> 13 | 0.2765208655
#> 14 | 0.2766217459
#> 15 | 0.2766929176
#> 16 | 0.2767431245
#> 17 | 0.2767785394
#> 18 | 0.2768035191
#> 19 | 0.2768211377
#> 20 | 0.2768335640
#> 21 | 0.2768423282
#> 22 | 0.2768485093
#> 23 | 0.2768528687
#> 24 | 0.2768559432
#> 25 | 0.2768581116
#> 26 | 0.2768596408
#> 27 | 0.2768607194
#> 28 | 0.2768614800
#> 29 | 0.2768620164
#> 30 | 0.2768623948
#> 31 | 0.2768626616
#> 32 | 0.2768628498
#> 33 | 0.2768629825
#> 34 | 0.2768630761
#> 35 | 0.2768631421
#> 36 | 0.2768631886
#> 37 | 0.2768632215
#> 38 | 0.2768632446
#> 39 | 0.2768632609
#> 40 | 0.2768632725
#> 41 | 0.2768632806
#> 42 | 0.2768632863
#> 43 | 0.2768632903
#> 44 | 0.2768632932
#> 45 | 0.2768632952
#> 46 | 0.2768632966
#> 47 | 0.2768632976
#> 48 | 0.2768632983
#> 49 | 0.2768632988
#> 50 | 0.2768632992
#> 51 | 0.2768632994
#> 52 | 0.2768632996
#> 53 | 0.2768632997
#> 54 | 0.2768632998
#> 55 | 0.2768632999
#> 56 | 0.2768632999
#> 57 | 0.2768632999
#> 58 | 0.2768633000
#> found new candidate with scale 0.2768633 in 58 iter (worst sc=0.27686)
#> Sample[140]: idc = 7 2
#> b^[] = 0.382345 0.709973
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 0.382345 0.709973
#> after refine_(*, conv=F):
#> beta_ref : 0.167816 0.825709
#> with ||beta_ref - beta_cand|| = 0.243757115723, --> sc = 0.613778594777883
#> Sample[141]: idc = 5 9
#> b^[] = -0.749016 1.16576
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -0.749016 1.16576
#> after refine_(*, conv=F):
#> beta_ref : -0.670506 1.15677
#> with ||beta_ref - beta_cand|| = 0.0790237477718, --> sc = 0.28633536578424
#> Sample[142]: idc = 1 8
#> b^[] = -1.24237 1.38371
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -1.24237 1.38371
#> after refine_(*, conv=F):
#> beta_ref : -0.710666 1.18096
#> with ||beta_ref - beta_cand|| = 0.569048775686, --> sc = 0.404424556196255
#> Sample[143]: idc = 12 5
#> b^[] = 1.03191 0.715867
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 1.03191 0.715867
#> after refine_(*, conv=F):
#> beta_ref : -0.376319 1.0606
#> with ||beta_ref - beta_cand|| = 1.4498103741, --> sc = 0.787749674878658
#> Sample[144]: idc = 1 6
#> b^[] = -2.59549 1.81452
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -2.59549 1.81452
#> after refine_(*, conv=F):
#> beta_ref : -2.14459 1.61896
#> with ||beta_ref - beta_cand|| = 0.491482052855, --> sc = 0.553607074103745
#> Sample[145]: idc = 5 10
#> b^[] = 3.35407 0.129244
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 3.35407 0.129244
#> after refine_(*, conv=F):
#> beta_ref : -0.22268 0.997273
#> with ||beta_ref - beta_cand|| = 3.68057059515, --> sc = 1.49558843265863
#> Sample[146]: idc = 3 2
#> b^[] = 0.301869 0.73001
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 0.301869 0.73001
#> after refine_(*, conv=F):
#> beta_ref : 0.14998 0.831077
#> with ||beta_ref - beta_cand|| = 0.182440710486, --> sc = 0.631676782286114
#> Sample[147]: idc = 2 8
#> b^[] = 1.56876 0.414582
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 1.56876 0.414582
#> after refine_(*, conv=F):
#> beta_ref : 0.260798 0.824069
#> with ||beta_ref - beta_cand|| = 1.37056338354, --> sc = 0.594013641561037
#> Sample[148]: idc = 11 2
#> b^[] = 0.931408 0.573269
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 0.931408 0.573269
#> after refine_(*, conv=F):
#> beta_ref : 0.290339 0.794023
#> with ||beta_ref - beta_cand|| = 0.678012796036, --> sc = 0.523721624955058
#> Sample[149]: idc = 7 2
#> b^[] = 0.382345 0.709973
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 0.382345 0.709973
#> after refine_(*, conv=F):
#> beta_ref : 0.167816 0.825709
#> with ||beta_ref - beta_cand|| = 0.243757115723, --> sc = 0.613778594777883
#> Sample[150]: idc = 3 9
#> b^[] = 0.0196235 0.871612
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 0.0196235 0.871612
#> after refine_(*, conv=F):
#> beta_ref : -0.0873117 0.935189
#> with ||beta_ref - beta_cand|| = 0.124406911988, --> sc = 0.456445784317417
#> Sample[151]: idc = 11 3
#> b^[] = -0.00341923 0.883173
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -0.00341923 0.883173
#> after refine_(*, conv=F):
#> beta_ref : -0.142669 0.957455
#> with ||beta_ref - beta_cand|| = 0.157824084971, --> sc = 0.473157356422566
#> Sample[152]: idc = 8 11
#> b^[] = 5.54256 -0.955372
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 5.54256 -0.955372
#> after refine_(*, conv=F):
#> beta_ref : 0.410243 0.743998
#> with ||beta_ref - beta_cand|| = 5.4063409827, --> sc = 1.40713960042917
#> Sample[153]: idc = 7 2
#> b^[] = 0.382345 0.709973
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 0.382345 0.709973
#> after refine_(*, conv=F):
#> beta_ref : 0.167816 0.825709
#> with ||beta_ref - beta_cand|| = 0.243757115723, --> sc = 0.613778594777883
#> Sample[154]: idc = 6 0
#> b^[] = -1.6587 1.16132
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -1.6587 1.16132
#> after refine_(*, conv=F):
#> beta_ref : -1.38655 1.28782
#> with ||beta_ref - beta_cand|| = 0.300121734322, --> sc = 1.32303815737056
#> Sample[155]: idc = 2 7
#> b^[] = 0.382345 0.709973
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 0.382345 0.709973
#> after refine_(*, conv=F):
#> beta_ref : 0.167816 0.825709
#> with ||beta_ref - beta_cand|| = 0.243757115723, --> sc = 0.613778594777883
#> Sample[156]: idc = 5 12
#> b^[] = 1.03191 0.715867
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 1.03191 0.715867
#> after refine_(*, conv=F):
#> beta_ref : -0.376319 1.0606
#> with ||beta_ref - beta_cand|| = 1.4498103741, --> sc = 0.787749674878659
#> Sample[157]: idc = 9 12
#> b^[] = -1.25272 1.35853
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -1.25272 1.35853
#> after refine_(*, conv=F):
#> beta_ref : -0.819674 1.21313
#> with ||beta_ref - beta_cand|| = 0.456802788525, --> sc = 0.301200329328142
#> Sample[158]: idc = 7 3
#> b^[] = 0.263938 0.74904
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 0.263938 0.74904
#> after refine_(*, conv=F):
#> beta_ref : 0.145859 0.83463
#> with ||beta_ref - beta_cand|| = 0.145837154087, --> sc = 0.597892842089404
#> Sample[159]: idc = 1 8
#> b^[] = -1.24237 1.38371
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -1.24237 1.38371
#> after refine_(*, conv=F):
#> beta_ref : -0.710666 1.18096
#> with ||beta_ref - beta_cand|| = 0.569048775686, --> sc = 0.404424556196255
#> Sample[160]: idc = 2 12
#> b^[] = 6.21824 -0.743034
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 6.21824 -0.743034
#> after refine_(*, conv=F):
#> beta_ref : -0.199361 0.979666
#> with ||beta_ref - beta_cand|| = 6.64479808942, --> sc = 1.90760566523583
#> Sample[161]: idc = 7 1
#> b^[] = -13.1628 5.179
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -13.1628 5.179
#> after refine_(*, conv=F):
#> beta_ref : -0.962999 1.2121
#> with ||beta_ref - beta_cand|| = 12.8285446507, --> sc = 2.75529751980209
#> Sample[162]: idc = 4 3
#> b^[] = -0.284513 1.0242
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -0.284513 1.0242
#> after refine_(*, conv=F):
#> beta_ref : -0.565247 1.12529
#> with ||beta_ref - beta_cand|| = 0.298379314564, --> sc = 0.280155905387692
#> Sample[163]: idc = 2 11
#> b^[] = 0.931408 0.573269
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 0.931408 0.573269
#> after refine_(*, conv=F):
#> beta_ref : 0.290339 0.794023
#> with ||beta_ref - beta_cand|| = 0.678012796036, --> sc = 0.523721624955058
#> Sample[164]: idc = 1 5
#> b^[] = 0.17655 0.931947
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 0.17655 0.931947
#> after refine_(*, conv=F):
#> beta_ref : -0.62749 1.15422
#> with ||beta_ref - beta_cand|| = 0.834196899796, --> sc = 0.429182774060505
#> find_scale(*, ini.scale =0.42918277406, tol=1e-10):
#> it | new scale
#> 0 | 0.3783009807
#> 1 | 0.3447919807
#> 2 | 0.3226261856
#> 3 | 0.3078259959
#> 4 | 0.2978417806
#> 5 | 0.2910448963
#> 6 | 0.2863841378
#> 7 | 0.2831706086
#> 8 | 0.2809460201
#> 9 | 0.2794015839
#> 10 | 0.2783271441
#> 11 | 0.2775785869
#> 12 | 0.2770565367
#> 13 | 0.2766921926
#> 14 | 0.2764377848
#> 15 | 0.2762600787
#> 16 | 0.2761359186
#> 17 | 0.2760491552
#> 18 | 0.2759885173
#> 19 | 0.2759461346
#> 20 | 0.2759165096
#> 21 | 0.2758958012
#> 22 | 0.2758813253
#> 23 | 0.2758712059
#> 24 | 0.2758641317
#> 25 | 0.2758591865
#> 26 | 0.2758557293
#> 27 | 0.2758533125
#> 28 | 0.2758516230
#> 29 | 0.2758504419
#> 30 | 0.2758496162
#> 31 | 0.2758490390
#> 32 | 0.2758486354
#> 33 | 0.2758483533
#> 34 | 0.2758481561
#> 35 | 0.2758480183
#> 36 | 0.2758479219
#> 37 | 0.2758478545
#> 38 | 0.2758478074
#> 39 | 0.2758477745
#> 40 | 0.2758477514
#> 41 | 0.2758477353
#> 42 | 0.2758477241
#> 43 | 0.2758477162
#> 44 | 0.2758477107
#> 45 | 0.2758477069
#> 46 | 0.2758477042
#> 47 | 0.2758477023
#> 48 | 0.2758477010
#> 49 | 0.2758477001
#> 50 | 0.2758476994
#> 51 | 0.2758476990
#> 52 | 0.2758476987
#> 53 | 0.2758476985
#> 54 | 0.2758476983
#> 55 | 0.2758476982
#> 56 | 0.2758476981
#> 57 | 0.2758476981
#> 58 | 0.2758476980
#> 59 | 0.2758476980
#> found new candidate with scale 0.2758477 in 59 iter (worst sc=0.27681)
#> Sample[165]: idc = 8 7
#> b^[] = 8.05358 -1.82104
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 8.05358 -1.82104
#> after refine_(*, conv=F):
#> beta_ref : 0.205065 0.812545
#> with ||beta_ref - beta_cand|| = 8.2785869972, --> sc = 2.06063624179737
#> Sample[166]: idc = 2 6
#> b^[] = -1.78547 1.24971
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -1.78547 1.24971
#> after refine_(*, conv=F):
#> beta_ref : -1.44982 1.30483
#> with ||beta_ref - beta_cand|| = 0.340142551919, --> sc = 1.10255399571346
#> Sample[167]: idc = 5 9
#> b^[] = -0.749016 1.16576
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -0.749016 1.16576
#> after refine_(*, conv=F):
#> beta_ref : -0.670506 1.15677
#> with ||beta_ref - beta_cand|| = 0.0790237477718, --> sc = 0.28633536578424
#> Sample[168]: idc = 3 8
#> b^[] = -0.471158 1.11784
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -0.471158 1.11784
#> after refine_(*, conv=F):
#> beta_ref : -0.548769 1.1376
#> with ||beta_ref - beta_cand|| = 0.08008750246, --> sc = 0.213513497916964
#> Sample[169]: idc = 8 10
#> b^[] = -0.79425 1.22922
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -0.79425 1.22922
#> after refine_(*, conv=F):
#> beta_ref : -0.629979 1.16172
#> with ||beta_ref - beta_cand|| = 0.177599485206, --> sc = 0.285059041656877
#> Sample[170]: idc = 3 2
#> b^[] = 0.301869 0.73001
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 0.301869 0.73001
#> after refine_(*, conv=F):
#> beta_ref : 0.14998 0.831077
#> with ||beta_ref - beta_cand|| = 0.182440710486, --> sc = 0.631676782286114
#> Sample[171]: idc = 2 1
#> b^[] = 2.63647 0.148748
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 2.63647 0.148748
#> after refine_(*, conv=F):
#> beta_ref : 0.060579 0.901992
#> with ||beta_ref - beta_cand|| = 2.68375998949, --> sc = 0.857121045415407
#> Sample[172]: idc = 2 9
#> b^[] = 0.553075 0.667465
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 0.553075 0.667465
#> after refine_(*, conv=F):
#> beta_ref : 0.207049 0.814328
#> with ||beta_ref - beta_cand|| = 0.375903533209, --> sc = 0.577742866999083
#> Sample[173]: idc = 2 12
#> b^[] = 6.21824 -0.743034
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 6.21824 -0.743034
#> after refine_(*, conv=F):
#> beta_ref : -0.199361 0.979666
#> with ||beta_ref - beta_cand|| = 6.64479808942, --> sc = 1.90760566523583
#> Sample[174]: idc = 0 1
#> b^[] = 24.2493 -6.73246
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 24.2493 -6.73246
#> after refine_(*, conv=F):
#> beta_ref : -0.505791 1.05334
#> with ||beta_ref - beta_cand|| = 25.9506189658, --> sc = 5.22853378335815
#> Sample[175]: idc = 6 0
#> b^[] = -1.6587 1.16132
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -1.6587 1.16132
#> after refine_(*, conv=F):
#> beta_ref : -1.38655 1.28782
#> with ||beta_ref - beta_cand|| = 0.300121734322, --> sc = 1.32303815737056
#> Sample[176]: idc = 0 7
#> b^[] = 7.13585 -1.51825
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 7.13585 -1.51825
#> after refine_(*, conv=F):
#> beta_ref : 0.224305 0.805575
#> with ||beta_ref - beta_cand|| = 7.29174469856, --> sc = 1.88768873167888
#> Sample[177]: idc = 7 3
#> b^[] = 0.263938 0.74904
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 0.263938 0.74904
#> after refine_(*, conv=F):
#> beta_ref : 0.145859 0.83463
#> with ||beta_ref - beta_cand|| = 0.145837154087, --> sc = 0.597892842089404
#> Sample[178]: idc = 9 4
#> b^[] = -0.706312 1.14942
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -0.706312 1.14942
#> after refine_(*, conv=F):
#> beta_ref : -0.662089 1.15324
#> with ||beta_ref - beta_cand|| = 0.0443876719904, --> sc = 0.302818894470934
#> Sample[179]: idc = 9 8
#> b^[] = -2.01045 1.6485
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -2.01045 1.6485
#> after refine_(*, conv=F):
#> beta_ref : -1.74249 1.49589
#> with ||beta_ref - beta_cand|| = 0.308377559274, --> sc = 0.579998989395441
#> Sample[180]: idc = 11 1
#> b^[] = -8.08546 3.56245
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -8.08546 3.56245
#> after refine_(*, conv=F):
#> beta_ref : -1.05124 1.24888
#> with ||beta_ref - beta_cand|| = 7.40491897223, --> sc = 1.76350297340317
#> Sample[181]: idc = 5 4
#> b^[] = -0.83181 1.18668
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -0.83181 1.18668
#> after refine_(*, conv=F):
#> beta_ref : -0.707317 1.16638
#> with ||beta_ref - beta_cand|| = 0.126137655448, --> sc = 0.307085307609751
#> Sample[182]: idc = 8 10
#> b^[] = -0.79425 1.22922
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -0.79425 1.22922
#> after refine_(*, conv=F):
#> beta_ref : -0.629979 1.16172
#> with ||beta_ref - beta_cand|| = 0.177599485206, --> sc = 0.285059041656877
#> Sample[183]: idc = 1 9
#> b^[] = -1.69558 1.52801
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -1.69558 1.52801
#> after refine_(*, conv=F):
#> beta_ref : -1.38063 1.38331
#> with ||beta_ref - beta_cand|| = 0.346599646578, --> sc = 0.467390770067497
#> Sample[184]: idc = 7 12
#> b^[] = -3.49557 1.98944
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -3.49557 1.98944
#> after refine_(*, conv=F):
#> beta_ref : -1.92479 1.5228
#> with ||beta_ref - beta_cand|| = 1.63862904714, --> sc = 0.741363565741794
#> Sample[185]: idc = 12 4
#> b^[] = -4.26167 2.20494
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -4.26167 2.20494
#> after refine_(*, conv=F):
#> beta_ref : -1.65432 1.43182
#> with ||beta_ref - beta_cand|| = 2.71955618857, --> sc = 0.91889370623179
#> Sample[186]: idc = 2 9
#> b^[] = 0.553075 0.667465
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 0.553075 0.667465
#> after refine_(*, conv=F):
#> beta_ref : 0.207049 0.814328
#> with ||beta_ref - beta_cand|| = 0.375903533209, --> sc = 0.577742866999083
#> Sample[187]: idc = 5 6
#> b^[] = -2.18551 1.52865
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -2.18551 1.52865
#> after refine_(*, conv=F):
#> beta_ref : -2.06748 1.54965
#> with ||beta_ref - beta_cand|| = 0.119878917294, --> sc = 0.539116852816472
#> Sample[188]: idc = 11 9
#> b^[] = -0.0569989 0.900935
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -0.0569989 0.900935
#> after refine_(*, conv=F):
#> beta_ref : -0.177777 0.971152
#> with ||beta_ref - beta_cand|| = 0.139706152895, --> sc = 0.486048896285347
#> Sample[189]: idc = 7 12
#> b^[] = -3.49557 1.98944
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -3.49557 1.98944
#> after refine_(*, conv=F):
#> beta_ref : -1.92479 1.5228
#> with ||beta_ref - beta_cand|| = 1.63862904714, --> sc = 0.741363565741794
#> Sample[190]: idc = 11 12
#> b^[] = -2.47152 1.70137
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -2.47152 1.70137
#> after refine_(*, conv=F):
#> beta_ref : -2.20041 1.61542
#> with ||beta_ref - beta_cand|| = 0.284411253958, --> sc = 0.366625751619072
#> Sample[191]: idc = 5 8
#> b^[] = -0.229469 1.03452
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -0.229469 1.03452
#> after refine_(*, conv=F):
#> beta_ref : -0.556193 1.13403
#> with ||beta_ref - beta_cand|| = 0.341541122565, --> sc = 0.270706285978508
#> Sample[192]: idc = 4 6
#> b^[] = -2.33514 1.63299
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -2.33514 1.63299
#> after refine_(*, conv=F):
#> beta_ref : -2.14415 1.58792
#> with ||beta_ref - beta_cand|| = 0.196236405624, --> sc = 0.399338780038278
#> Sample[193]: idc = 3 9
#> b^[] = 0.0196235 0.871612
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 0.0196235 0.871612
#> after refine_(*, conv=F):
#> beta_ref : -0.0873117 0.935189
#> with ||beta_ref - beta_cand|| = 0.124406911988, --> sc = 0.456445784317417
#> Sample[194]: idc = 1 7
#> b^[] = -13.1628 5.179
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -13.1628 5.179
#> after refine_(*, conv=F):
#> beta_ref : -0.962999 1.2121
#> with ||beta_ref - beta_cand|| = 12.8285446507, --> sc = 2.75529751980209
#> Sample[195]: idc = 9 8
#> b^[] = -2.01045 1.6485
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -2.01045 1.6485
#> after refine_(*, conv=F):
#> beta_ref : -1.74249 1.49589
#> with ||beta_ref - beta_cand|| = 0.308377559274, --> sc = 0.579998989395441
#> Sample[196]: idc = 11 4
#> b^[] = -1.66598 1.43433
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -1.66598 1.43433
#> after refine_(*, conv=F):
#> beta_ref : -1.67229 1.44956
#> with ||beta_ref - beta_cand|| = 0.0164864177799, --> sc = 0.348265804600926
#> Sample[197]: idc = 12 2
#> b^[] = 6.21824 -0.743034
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 6.21824 -0.743034
#> after refine_(*, conv=F):
#> beta_ref : -0.199361 0.979666
#> with ||beta_ref - beta_cand|| = 6.64479808942, --> sc = 1.90760566523583
#> Sample[198]: idc = 2 11
#> b^[] = 0.931408 0.573269
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 0.931408 0.573269
#> after refine_(*, conv=F):
#> beta_ref : 0.290339 0.794023
#> with ||beta_ref - beta_cand|| = 0.678012796036, --> sc = 0.523721624955058
#> Sample[199]: idc = 3 5
#> b^[] = -0.381765 1.07299
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -0.381765 1.07299
#> after refine_(*, conv=F):
#> beta_ref : -0.547363 1.12925
#> with ||beta_ref - beta_cand|| = 0.17489423735, --> sc = 0.235441204593026
#> Sample[200]: idc = 11 6
#> b^[] = -2.39849 1.67716
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -2.39849 1.67716
#> after refine_(*, conv=F):
#> beta_ref : -2.18647 1.6102
#> with ||beta_ref - beta_cand|| = 0.222349315351, --> sc = 0.365943372441649
#> Sample[201]: idc = 6 12
#> b^[] = -2.40731 1.68331
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -2.40731 1.68331
#> after refine_(*, conv=F):
#> beta_ref : -2.20087 1.61624
#> with ||beta_ref - beta_cand|| = 0.2170640588, --> sc = 0.35543328734486
#> Sample[202]: idc = 0 3
#> b^[] = 1.14468 0.30717
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 1.14468 0.30717
#> after refine_(*, conv=F):
#> beta_ref : -0.307508 0.964066
#> with ||beta_ref - beta_cand|| = 1.59385397212, --> sc = 1.34390327405485
#> Sample[203]: idc = 6 10
#> b^[] = -2.3463 1.64077
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -2.3463 1.64077
#> after refine_(*, conv=F):
#> beta_ref : -2.13209 1.58647
#> with ||beta_ref - beta_cand|| = 0.220987541723, --> sc = 0.417145691249323
#> Sample[204]: idc = 9 10
#> b^[] = -1.1869 1.33334
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -1.1869 1.33334
#> after refine_(*, conv=F):
#> beta_ref : -0.809481 1.2096
#> with ||beta_ref - beta_cand|| = 0.39718852225, --> sc = 0.284340124502366
#> Sample[205]: idc = 12 10
#> b^[] = -0.773295 1.22367
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -0.773295 1.22367
#> after refine_(*, conv=F):
#> beta_ref : -0.62546 1.16042
#> with ||beta_ref - beta_cand|| = 0.160794576384, --> sc = 0.287874239981522
#> Sample[206]: idc = 10 6
#> b^[] = -2.3463 1.64077
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -2.3463 1.64077
#> after refine_(*, conv=F):
#> beta_ref : -2.13209 1.58647
#> with ||beta_ref - beta_cand|| = 0.220987541723, --> sc = 0.417145691249323
#> Sample[207]: idc = 7 9
#> b^[] = 0.815075 0.5672
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 0.815075 0.5672
#> after refine_(*, conv=F):
#> beta_ref : 0.224676 0.802773
#> with ||beta_ref - beta_cand|| = 0.635661711031, --> sc = 0.67968553517208
#> Sample[208]: idc = 9 7
#> b^[] = 0.815075 0.5672
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 0.815075 0.5672
#> after refine_(*, conv=F):
#> beta_ref : 0.224676 0.802773
#> with ||beta_ref - beta_cand|| = 0.635661711031, --> sc = 0.679685535172079
#> Sample[209]: idc = 10 11
#> b^[] = -2.05854 1.56446
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -2.05854 1.56446
#> after refine_(*, conv=F):
#> beta_ref : -2.07016 1.57129
#> with ||beta_ref - beta_cand|| = 0.0134780521353, --> sc = 0.37327449031415
#> Sample[210]: idc = 10 8
#> b^[] = -0.79425 1.22922
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -0.79425 1.22922
#> after refine_(*, conv=F):
#> beta_ref : -0.629979 1.16172
#> with ||beta_ref - beta_cand|| = 0.177599485206, --> sc = 0.285059041656878
#> Sample[211]: idc = 12 4
#> b^[] = -4.26167 2.20494
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -4.26167 2.20494
#> after refine_(*, conv=F):
#> beta_ref : -1.65432 1.43182
#> with ||beta_ref - beta_cand|| = 2.71955618857, --> sc = 0.91889370623179
#> Sample[212]: idc = 1 4
#> b^[] = 2.25191 0.271184
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 2.25191 0.271184
#> after refine_(*, conv=F):
#> beta_ref : 0.0105456 0.920859
#> with ||beta_ref - beta_cand|| = 2.33362366183, --> sc = 0.813737418280548
#> Sample[213]: idc = 3 5
#> b^[] = -0.381765 1.07299
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -0.381765 1.07299
#> after refine_(*, conv=F):
#> beta_ref : -0.547363 1.12925
#> with ||beta_ref - beta_cand|| = 0.17489423735, --> sc = 0.235441204593026
#> Sample[214]: idc = 0 8
#> b^[] = 7.47514 -1.62162
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 7.47514 -1.62162
#> after refine_(*, conv=F):
#> beta_ref : 0.245606 0.798807
#> with ||beta_ref - beta_cand|| = 7.62395703783, --> sc = 1.90882061471928
#> Sample[215]: idc = 8 2
#> b^[] = 1.56876 0.414582
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 1.56876 0.414582
#> after refine_(*, conv=F):
#> beta_ref : 0.260798 0.824069
#> with ||beta_ref - beta_cand|| = 1.37056338354, --> sc = 0.594013641561037
#> Sample[216]: idc = 6 7
#> b^[] = -2.26322 1.58284
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -2.26322 1.58284
#> after refine_(*, conv=F):
#> beta_ref : -2.18436 1.58611
#> with ||beta_ref - beta_cand|| = 0.078929299088, --> sc = 0.361556853894617
#> Sample[217]: idc = 0 2
#> b^[] = -2.67885 1.47214
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -2.67885 1.47214
#> after refine_(*, conv=F):
#> beta_ref : -1.4879 1.31816
#> with ||beta_ref - beta_cand|| = 1.20086135032, --> sc = 1.37247806701654
#> Sample[218]: idc = 3 8
#> b^[] = -0.471158 1.11784
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -0.471158 1.11784
#> after refine_(*, conv=F):
#> beta_ref : -0.548769 1.1376
#> with ||beta_ref - beta_cand|| = 0.08008750246, --> sc = 0.213513497916964
#> Sample[219]: idc = 3 7
#> b^[] = 0.263938 0.74904
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 0.263938 0.74904
#> after refine_(*, conv=F):
#> beta_ref : 0.145859 0.83463
#> with ||beta_ref - beta_cand|| = 0.145837154087, --> sc = 0.597892842089404
#> Sample[220]: idc = 7 4
#> b^[] = -3.13442 1.87028
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -3.13442 1.87028
#> after refine_(*, conv=F):
#> beta_ref : -2.01174 1.5504
#> with ||beta_ref - beta_cand|| = 1.16736173632, --> sc = 0.618257718400314
#> Sample[221]: idc = 11 5
#> b^[] = -1.19798 1.27918
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -1.19798 1.27918
#> after refine_(*, conv=F):
#> beta_ref : -0.898455 1.21902
#> with ||beta_ref - beta_cand|| = 0.305504698816, --> sc = 0.367535728188709
#> Sample[222]: idc = 9 7
#> b^[] = 0.815075 0.5672
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 0.815075 0.5672
#> after refine_(*, conv=F):
#> beta_ref : 0.224676 0.802773
#> with ||beta_ref - beta_cand|| = 0.635661711031, --> sc = 0.679685535172079
#> Sample[223]: idc = 8 7
#> b^[] = 8.05358 -1.82104
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 8.05358 -1.82104
#> after refine_(*, conv=F):
#> beta_ref : 0.205065 0.812545
#> with ||beta_ref - beta_cand|| = 8.2785869972, --> sc = 2.06063624179737
#> Sample[224]: idc = 4 5
#> b^[] = -0.83181 1.18668
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -0.83181 1.18668
#> after refine_(*, conv=F):
#> beta_ref : -0.707317 1.16638
#> with ||beta_ref - beta_cand|| = 0.126137655448, --> sc = 0.307085307609751
#> Sample[225]: idc = 12 8
#> b^[] = -0.799579 1.23106
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -0.799579 1.23106
#> after refine_(*, conv=F):
#> beta_ref : -0.631179 1.16203
#> with ||beta_ref - beta_cand|| = 0.181998057096, --> sc = 0.286399245856023
#> Sample[226]: idc = 6 2
#> b^[] = -1.78547 1.24971
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -1.78547 1.24971
#> after refine_(*, conv=F):
#> beta_ref : -1.44982 1.30483
#> with ||beta_ref - beta_cand|| = 0.340142551919, --> sc = 1.10255399571346
#> Sample[227]: idc = 9 1
#> b^[] = -1.69558 1.52801
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -1.69558 1.52801
#> after refine_(*, conv=F):
#> beta_ref : -1.38063 1.38331
#> with ||beta_ref - beta_cand|| = 0.346599646578, --> sc = 0.467390770067498
#> Sample[228]: idc = 8 4
#> b^[] = 0.327489 0.842506
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 0.327489 0.842506
#> after refine_(*, conv=F):
#> beta_ref : -0.488056 1.09604
#> with ||beta_ref - beta_cand|| = 0.854043995939, --> sc = 0.406506642125283
#> Sample[229]: idc = 7 11
#> b^[] = 29.1679 -8.78739
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 29.1679 -8.78739
#> after refine_(*, conv=F):
#> beta_ref : -0.244111 0.966143
#> with ||beta_ref - beta_cand|| = 30.9870463767, --> sc = 6.48063687575061
#> Sample[230]: idc = 9 0
#> b^[] = 2.86111 -0.215798
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 2.86111 -0.215798
#> after refine_(*, conv=F):
#> beta_ref : 0.171545 0.81517
#> with ||beta_ref - beta_cand|| = 2.88038890444, --> sc = 1.45332402557657
#> Sample[231]: idc = 4 2
#> b^[] = 2.8091 0.105765
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 2.8091 0.105765
#> after refine_(*, conv=F):
#> beta_ref : 0.030844 0.911664
#> with ||beta_ref - beta_cand|| = 2.89278369939, --> sc = 0.908806025778876
#> Sample[232]: idc = 0 7
#> b^[] = 7.13585 -1.51825
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 7.13585 -1.51825
#> after refine_(*, conv=F):
#> beta_ref : 0.224305 0.805575
#> with ||beta_ref - beta_cand|| = 7.29174469856, --> sc = 1.88768873167888
#> Sample[233]: idc = 1 4
#> b^[] = 2.25191 0.271184
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 2.25191 0.271184
#> after refine_(*, conv=F):
#> beta_ref : 0.0105456 0.920859
#> with ||beta_ref - beta_cand|| = 2.33362366183, --> sc = 0.813737418280548
#> Sample[234]: idc = 3 1
#> b^[] = -0.582063 1.17348
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -0.582063 1.17348
#> after refine_(*, conv=F):
#> beta_ref : -0.597198 1.15265
#> with ||beta_ref - beta_cand|| = 0.0257497340318, --> sc = 0.312109152133832
#> Sample[235]: idc = 7 4
#> b^[] = -3.13442 1.87028
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -3.13442 1.87028
#> after refine_(*, conv=F):
#> beta_ref : -2.01174 1.5504
#> with ||beta_ref - beta_cand|| = 1.16736173632, --> sc = 0.618257718400314
#> Sample[236]: idc = 10 11
#> b^[] = -2.05854 1.56446
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -2.05854 1.56446
#> after refine_(*, conv=F):
#> beta_ref : -2.07016 1.57129
#> with ||beta_ref - beta_cand|| = 0.0134780521353, --> sc = 0.37327449031415
#> Sample[237]: idc = 7 11
#> b^[] = 29.1679 -8.78739
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 29.1679 -8.78739
#> after refine_(*, conv=F):
#> beta_ref : -0.244111 0.966143
#> with ||beta_ref - beta_cand|| = 30.9870463767, --> sc = 6.48063687575061
#> Sample[238]: idc = 7 9
#> b^[] = 0.815075 0.5672
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 0.815075 0.5672
#> after refine_(*, conv=F):
#> beta_ref : 0.224676 0.802773
#> with ||beta_ref - beta_cand|| = 0.635661711031, --> sc = 0.67968553517208
#> Sample[239]: idc = 2 9
#> b^[] = 0.553075 0.667465
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 0.553075 0.667465
#> after refine_(*, conv=F):
#> beta_ref : 0.207049 0.814328
#> with ||beta_ref - beta_cand|| = 0.375903533209, --> sc = 0.577742866999083
#> Sample[240]: idc = 8 5
#> b^[] = -0.229469 1.03452
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -0.229469 1.03452
#> after refine_(*, conv=F):
#> beta_ref : -0.556193 1.13403
#> with ||beta_ref - beta_cand|| = 0.341541122565, --> sc = 0.270706285978508
#> Sample[241]: idc = 10 9
#> b^[] = -1.1869 1.33334
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -1.1869 1.33334
#> after refine_(*, conv=F):
#> beta_ref : -0.809481 1.2096
#> with ||beta_ref - beta_cand|| = 0.39718852225, --> sc = 0.284340124502365
#> Sample[242]: idc = 7 11
#> b^[] = 29.1679 -8.78739
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 29.1679 -8.78739
#> after refine_(*, conv=F):
#> beta_ref : -0.244111 0.966143
#> with ||beta_ref - beta_cand|| = 30.9870463767, --> sc = 6.48063687575061
#> Sample[243]: idc = 11 8
#> b^[] = 5.54256 -0.955372
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 5.54256 -0.955372
#> after refine_(*, conv=F):
#> beta_ref : 0.410243 0.743998
#> with ||beta_ref - beta_cand|| = 5.4063409827, --> sc = 1.40713960042918
#> Sample[244]: idc = 0 5
#> b^[] = -6.1579 2.53215
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -6.1579 2.53215
#> after refine_(*, conv=F):
#> beta_ref : -0.763906 1.12366
#> with ||beta_ref - beta_cand|| = 5.57485454363, --> sc = 1.64015945823134
#> Sample[245]: idc = 5 12
#> b^[] = 1.03191 0.715867
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 1.03191 0.715867
#> after refine_(*, conv=F):
#> beta_ref : -0.376319 1.0606
#> with ||beta_ref - beta_cand|| = 1.4498103741, --> sc = 0.787749674878659
#> Sample[246]: idc = 10 6
#> b^[] = -2.3463 1.64077
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -2.3463 1.64077
#> after refine_(*, conv=F):
#> beta_ref : -2.13209 1.58647
#> with ||beta_ref - beta_cand|| = 0.220987541723, --> sc = 0.417145691249323
#> Sample[247]: idc = 5 7
#> b^[] = -1.81626 1.43537
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -1.81626 1.43537
#> after refine_(*, conv=F):
#> beta_ref : -2.02342 1.54672
#> with ||beta_ref - beta_cand|| = 0.235186747649, --> sc = 0.472769193580754
#> Sample[248]: idc = 5 12
#> b^[] = 1.03191 0.715867
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 1.03191 0.715867
#> after refine_(*, conv=F):
#> beta_ref : -0.376319 1.0606
#> with ||beta_ref - beta_cand|| = 1.4498103741, --> sc = 0.787749674878659
#> Sample[249]: idc = 6 0
#> b^[] = -1.6587 1.16132
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -1.6587 1.16132
#> after refine_(*, conv=F):
#> beta_ref : -1.38655 1.28782
#> with ||beta_ref - beta_cand|| = 0.300121734322, --> sc = 1.32303815737056
#> Sample[250]: idc = 8 9
#> b^[] = -2.01045 1.6485
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -2.01045 1.6485
#> after refine_(*, conv=F):
#> beta_ref : -1.74249 1.49589
#> with ||beta_ref - beta_cand|| = 0.308377559274, --> sc = 0.579998989395441
#> Sample[251]: idc = 3 6
#> b^[] = -4.48261 3.13038
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -4.48261 3.13038
#> after refine_(*, conv=F):
#> beta_ref : -1.5156 1.39333
#> with ||beta_ref - beta_cand|| = 3.43809929004, --> sc = 3.95137685252739
#> Sample[252]: idc = 0 2
#> b^[] = -2.67885 1.47214
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -2.67885 1.47214
#> after refine_(*, conv=F):
#> beta_ref : -1.4879 1.31816
#> with ||beta_ref - beta_cand|| = 1.20086135032, --> sc = 1.37247806701654
#> Sample[253]: idc = 10 8
#> b^[] = -0.79425 1.22922
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -0.79425 1.22922
#> after refine_(*, conv=F):
#> beta_ref : -0.629979 1.16172
#> with ||beta_ref - beta_cand|| = 0.177599485206, --> sc = 0.285059041656878
#> Sample[254]: idc = 5 10
#> b^[] = 3.35407 0.129244
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 3.35407 0.129244
#> after refine_(*, conv=F):
#> beta_ref : -0.22268 0.997273
#> with ||beta_ref - beta_cand|| = 3.68057059515, --> sc = 1.49558843265863
#> Sample[255]: idc = 12 1
#> b^[] = -0.484845 1.14253
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -0.484845 1.14253
#> after refine_(*, conv=F):
#> beta_ref : -0.587028 1.14935
#> with ||beta_ref - beta_cand|| = 0.102411328665, --> sc = 0.312564686828127
#> Sample[256]: idc = 7 2
#> b^[] = 0.382345 0.709973
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 0.382345 0.709973
#> after refine_(*, conv=F):
#> beta_ref : 0.167816 0.825709
#> with ||beta_ref - beta_cand|| = 0.243757115723, --> sc = 0.613778594777883
#> Sample[257]: idc = 12 6
#> b^[] = -2.40731 1.68331
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -2.40731 1.68331
#> after refine_(*, conv=F):
#> beta_ref : -2.20087 1.61624
#> with ||beta_ref - beta_cand|| = 0.2170640588, --> sc = 0.355433287344859
#> Sample[258]: idc = 8 12
#> b^[] = -0.799579 1.23106
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -0.799579 1.23106
#> after refine_(*, conv=F):
#> beta_ref : -0.631179 1.16203
#> with ||beta_ref - beta_cand|| = 0.181998057096, --> sc = 0.286399245856023
#> Sample[259]: idc = 6 7
#> b^[] = -2.26322 1.58284
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -2.26322 1.58284
#> after refine_(*, conv=F):
#> beta_ref : -2.18436 1.58611
#> with ||beta_ref - beta_cand|| = 0.078929299088, --> sc = 0.361556853894617
#> Sample[260]: idc = 4 3
#> b^[] = -0.284513 1.0242
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -0.284513 1.0242
#> after refine_(*, conv=F):
#> beta_ref : -0.565247 1.12529
#> with ||beta_ref - beta_cand|| = 0.298379314564, --> sc = 0.280155905387692
#> Sample[261]: idc = 6 3
#> b^[] = -4.48261 3.13038
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -4.48261 3.13038
#> after refine_(*, conv=F):
#> beta_ref : -1.5156 1.39333
#> with ||beta_ref - beta_cand|| = 3.43809929004, --> sc = 3.95137685252739
#> Sample[262]: idc = 4 9
#> b^[] = -0.706312 1.14942
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -0.706312 1.14942
#> after refine_(*, conv=F):
#> beta_ref : -0.662089 1.15324
#> with ||beta_ref - beta_cand|| = 0.0443876719904, --> sc = 0.302818894470935
#> Sample[263]: idc = 11 2
#> b^[] = 0.931408 0.573269
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 0.931408 0.573269
#> after refine_(*, conv=F):
#> beta_ref : 0.290339 0.794023
#> with ||beta_ref - beta_cand|| = 0.678012796036, --> sc = 0.523721624955058
#> Sample[264]: idc = 3 9
#> b^[] = 0.0196235 0.871612
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 0.0196235 0.871612
#> after refine_(*, conv=F):
#> beta_ref : -0.0873117 0.935189
#> with ||beta_ref - beta_cand|| = 0.124406911988, --> sc = 0.456445784317417
#> Sample[265]: idc = 10 8
#> b^[] = -0.79425 1.22922
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -0.79425 1.22922
#> after refine_(*, conv=F):
#> beta_ref : -0.629979 1.16172
#> with ||beta_ref - beta_cand|| = 0.177599485206, --> sc = 0.285059041656878
#> Sample[266]: idc = 12 6
#> b^[] = -2.40731 1.68331
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -2.40731 1.68331
#> after refine_(*, conv=F):
#> beta_ref : -2.20087 1.61624
#> with ||beta_ref - beta_cand|| = 0.2170640588, --> sc = 0.355433287344859
#> Sample[267]: idc = 1 8
#> b^[] = -1.24237 1.38371
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -1.24237 1.38371
#> after refine_(*, conv=F):
#> beta_ref : -0.710666 1.18096
#> with ||beta_ref - beta_cand|| = 0.569048775686, --> sc = 0.404424556196255
#> Sample[268]: idc = 1 3
#> b^[] = -0.582063 1.17348
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -0.582063 1.17348
#> after refine_(*, conv=F):
#> beta_ref : -0.597198 1.15265
#> with ||beta_ref - beta_cand|| = 0.0257497340318, --> sc = 0.312109152133831
#> Sample[269]: idc = 1 3
#> b^[] = -0.582063 1.17348
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -0.582063 1.17348
#> after refine_(*, conv=F):
#> beta_ref : -0.597198 1.15265
#> with ||beta_ref - beta_cand|| = 0.0257497340318, --> sc = 0.312109152133831
#> Sample[270]: idc = 5 11
#> b^[] = -1.19798 1.27918
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -1.19798 1.27918
#> after refine_(*, conv=F):
#> beta_ref : -0.898455 1.21902
#> with ||beta_ref - beta_cand|| = 0.305504698816, --> sc = 0.367535728188709
#> Sample[271]: idc = 7 2
#> b^[] = 0.382345 0.709973
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 0.382345 0.709973
#> after refine_(*, conv=F):
#> beta_ref : 0.167816 0.825709
#> with ||beta_ref - beta_cand|| = 0.243757115723, --> sc = 0.613778594777883
#> Sample[272]: idc = 9 4
#> b^[] = -0.706312 1.14942
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -0.706312 1.14942
#> after refine_(*, conv=F):
#> beta_ref : -0.662089 1.15324
#> with ||beta_ref - beta_cand|| = 0.0443876719904, --> sc = 0.302818894470934
#> Sample[273]: idc = 3 8
#> b^[] = -0.471158 1.11784
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -0.471158 1.11784
#> after refine_(*, conv=F):
#> beta_ref : -0.548769 1.1376
#> with ||beta_ref - beta_cand|| = 0.08008750246, --> sc = 0.213513497916964
#> Sample[274]: idc = 10 0
#> b^[] = -9.17664 3.45192
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -9.17664 3.45192
#> after refine_(*, conv=F):
#> beta_ref : -0.840132 1.1533
#> with ||beta_ref - beta_cand|| = 8.64759632726, --> sc = 2.16009257867085
#> Sample[275]: idc = 2 10
#> b^[] = 13.188 -2.47835
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 13.188 -2.47835
#> after refine_(*, conv=F):
#> beta_ref : -0.333051 1.01284
#> with ||beta_ref - beta_cand|| = 13.9645138921, --> sc = 4.27637683340343
#> Sample[276]: idc = 6 1
#> b^[] = -2.59549 1.81452
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -2.59549 1.81452
#> after refine_(*, conv=F):
#> beta_ref : -2.14459 1.61896
#> with ||beta_ref - beta_cand|| = 0.491482052855, --> sc = 0.553607074103743
#> Sample[277]: idc = 4 7
#> b^[] = -3.13442 1.87028
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -3.13442 1.87028
#> after refine_(*, conv=F):
#> beta_ref : -2.01174 1.5504
#> with ||beta_ref - beta_cand|| = 1.16736173632, --> sc = 0.618257718400311
#> Sample[278]: idc = 11 1
#> b^[] = -8.08546 3.56245
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -8.08546 3.56245
#> after refine_(*, conv=F):
#> beta_ref : -1.05124 1.24888
#> with ||beta_ref - beta_cand|| = 7.40491897223, --> sc = 1.76350297340317
#> Sample[279]: idc = 2 10
#> b^[] = 13.188 -2.47835
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 13.188 -2.47835
#> after refine_(*, conv=F):
#> beta_ref : -0.333051 1.01284
#> with ||beta_ref - beta_cand|| = 13.9645138921, --> sc = 4.27637683340343
#> Sample[280]: idc = 7 0
#> b^[] = 7.13585 -1.51825
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 7.13585 -1.51825
#> after refine_(*, conv=F):
#> beta_ref : 0.224305 0.805575
#> with ||beta_ref - beta_cand|| = 7.29174469856, --> sc = 1.88768873167888
#> Sample[281]: idc = 6 1
#> b^[] = -2.59549 1.81452
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -2.59549 1.81452
#> after refine_(*, conv=F):
#> beta_ref : -2.14459 1.61896
#> with ||beta_ref - beta_cand|| = 0.491482052855, --> sc = 0.553607074103743
#> Sample[282]: idc = 9 12
#> b^[] = -1.25272 1.35853
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -1.25272 1.35853
#> after refine_(*, conv=F):
#> beta_ref : -0.819674 1.21313
#> with ||beta_ref - beta_cand|| = 0.456802788525, --> sc = 0.301200329328142
#> Sample[283]: idc = 12 7
#> b^[] = -3.49557 1.98944
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -3.49557 1.98944
#> after refine_(*, conv=F):
#> beta_ref : -1.92479 1.5228
#> with ||beta_ref - beta_cand|| = 1.63862904714, --> sc = 0.741363565741794
#> Sample[284]: idc = 9 0
#> b^[] = 2.86111 -0.215798
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 2.86111 -0.215798
#> after refine_(*, conv=F):
#> beta_ref : 0.171545 0.81517
#> with ||beta_ref - beta_cand|| = 2.88038890444, --> sc = 1.45332402557657
#> Sample[285]: idc = 6 5
#> b^[] = -2.18551 1.52865
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -2.18551 1.52865
#> after refine_(*, conv=F):
#> beta_ref : -2.06748 1.54965
#> with ||beta_ref - beta_cand|| = 0.119878917294, --> sc = 0.539116852816471
#> Sample[286]: idc = 7 11
#> b^[] = 29.1679 -8.78739
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 29.1679 -8.78739
#> after refine_(*, conv=F):
#> beta_ref : -0.244111 0.966143
#> with ||beta_ref - beta_cand|| = 30.9870463767, --> sc = 6.48063687575061
#> Sample[287]: idc = 10 8
#> b^[] = -0.79425 1.22922
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -0.79425 1.22922
#> after refine_(*, conv=F):
#> beta_ref : -0.629979 1.16172
#> with ||beta_ref - beta_cand|| = 0.177599485206, --> sc = 0.285059041656878
#> Sample[288]: idc = 7 1
#> b^[] = -13.1628 5.179
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -13.1628 5.179
#> after refine_(*, conv=F):
#> beta_ref : -0.962999 1.2121
#> with ||beta_ref - beta_cand|| = 12.8285446507, --> sc = 2.75529751980209
#> Sample[289]: idc = 6 8
#> b^[] = -2.69667 1.88507
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -2.69667 1.88507
#> after refine_(*, conv=F):
#> beta_ref : -1.94237 1.55803
#> with ||beta_ref - beta_cand|| = 0.822144059656, --> sc = 0.744918088188393
#> Sample[290]: idc = 5 3
#> b^[] = -0.381765 1.07299
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -0.381765 1.07299
#> after refine_(*, conv=F):
#> beta_ref : -0.547363 1.12925
#> with ||beta_ref - beta_cand|| = 0.17489423735, --> sc = 0.235441204593026
#> Sample[291]: idc = 12 8
#> b^[] = -0.799579 1.23106
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -0.799579 1.23106
#> after refine_(*, conv=F):
#> beta_ref : -0.631179 1.16203
#> with ||beta_ref - beta_cand|| = 0.181998057096, --> sc = 0.286399245856023
#> Sample[292]: idc = 8 6
#> b^[] = -2.69667 1.88507
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -2.69667 1.88507
#> after refine_(*, conv=F):
#> beta_ref : -1.94237 1.55803
#> with ||beta_ref - beta_cand|| = 0.822144059656, --> sc = 0.744918088188394
#> Sample[293]: idc = 8 2
#> b^[] = 1.56876 0.414582
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 1.56876 0.414582
#> after refine_(*, conv=F):
#> beta_ref : 0.260798 0.824069
#> with ||beta_ref - beta_cand|| = 1.37056338354, --> sc = 0.594013641561037
#> Sample[294]: idc = 8 2
#> b^[] = 1.56876 0.414582
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 1.56876 0.414582
#> after refine_(*, conv=F):
#> beta_ref : 0.260798 0.824069
#> with ||beta_ref - beta_cand|| = 1.37056338354, --> sc = 0.594013641561037
#> Sample[295]: idc = 11 12
#> b^[] = -2.47152 1.70137
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -2.47152 1.70137
#> after refine_(*, conv=F):
#> beta_ref : -2.20041 1.61542
#> with ||beta_ref - beta_cand|| = 0.284411253958, --> sc = 0.366625751619072
#> Sample[296]: idc = 10 12
#> b^[] = -0.773295 1.22367
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -0.773295 1.22367
#> after refine_(*, conv=F):
#> beta_ref : -0.62546 1.16042
#> with ||beta_ref - beta_cand|| = 0.160794576384, --> sc = 0.287874239981522
#> Sample[297]: idc = 11 7
#> b^[] = 29.1679 -8.78739
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 29.1679 -8.78739
#> after refine_(*, conv=F):
#> beta_ref : -0.244111 0.966143
#> with ||beta_ref - beta_cand|| = 30.9870463767, --> sc = 6.48063687575041
#> Sample[298]: idc = 11 10
#> b^[] = -2.05854 1.56446
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -2.05854 1.56446
#> after refine_(*, conv=F):
#> beta_ref : -2.07016 1.57129
#> with ||beta_ref - beta_cand|| = 0.0134780521353, --> sc = 0.37327449031415
#> Sample[299]: idc = 6 7
#> b^[] = -2.26322 1.58284
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -2.26322 1.58284
#> after refine_(*, conv=F):
#> beta_ref : -2.18436 1.58611
#> with ||beta_ref - beta_cand|| = 0.078929299088, --> sc = 0.361556853894617
#> Sample[300]: idc = 3 10
#> b^[] = -0.579864 1.17238
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -0.579864 1.17238
#> after refine_(*, conv=F):
#> beta_ref : -0.596172 1.15246
#> with ||beta_ref - beta_cand|| = 0.0257405973865, --> sc = 0.309002235310989
#> Sample[301]: idc = 10 6
#> b^[] = -2.3463 1.64077
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -2.3463 1.64077
#> after refine_(*, conv=F):
#> beta_ref : -2.13209 1.58647
#> with ||beta_ref - beta_cand|| = 0.220987541723, --> sc = 0.417145691249323
#> Sample[302]: idc = 8 0
#> b^[] = 7.47514 -1.62162
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 7.47514 -1.62162
#> after refine_(*, conv=F):
#> beta_ref : 0.245606 0.798807
#> with ||beta_ref - beta_cand|| = 7.62395703783, --> sc = 1.90882061471928
#> Sample[303]: idc = 1 4
#> b^[] = 2.25191 0.271184
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 2.25191 0.271184
#> after refine_(*, conv=F):
#> beta_ref : 0.0105456 0.920859
#> with ||beta_ref - beta_cand|| = 2.33362366183, --> sc = 0.813737418280548
#> Sample[304]: idc = 2 12
#> b^[] = 6.21824 -0.743034
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 6.21824 -0.743034
#> after refine_(*, conv=F):
#> beta_ref : -0.199361 0.979666
#> with ||beta_ref - beta_cand|| = 6.64479808942, --> sc = 1.90760566523583
#> Sample[305]: idc = 4 6
#> b^[] = -2.33514 1.63299
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -2.33514 1.63299
#> after refine_(*, conv=F):
#> beta_ref : -2.14415 1.58792
#> with ||beta_ref - beta_cand|| = 0.196236405624, --> sc = 0.399338780038278
#> Sample[306]: idc = 2 5
#> b^[] = 47.0534 -10.9101
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 47.0534 -10.9101
#> after refine_(*, conv=F):
#> beta_ref : -0.397375 1.02753
#> with ||beta_ref - beta_cand|| = 48.929345941, --> sc = 15.1481999432315
#> Sample[307]: idc = 7 0
#> b^[] = 7.13585 -1.51825
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 7.13585 -1.51825
#> after refine_(*, conv=F):
#> beta_ref : 0.224305 0.805575
#> with ||beta_ref - beta_cand|| = 7.29174469856, --> sc = 1.88768873167888
#> Sample[308]: idc = 2 11
#> b^[] = 0.931408 0.573269
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 0.931408 0.573269
#> after refine_(*, conv=F):
#> beta_ref : 0.290339 0.794023
#> with ||beta_ref - beta_cand|| = 0.678012796036, --> sc = 0.523721624955058
#> Sample[309]: idc = 7 3
#> b^[] = 0.263938 0.74904
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 0.263938 0.74904
#> after refine_(*, conv=F):
#> beta_ref : 0.145859 0.83463
#> with ||beta_ref - beta_cand|| = 0.145837154087, --> sc = 0.597892842089404
#> Sample[310]: idc = 4 11
#> b^[] = -1.66598 1.43433
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -1.66598 1.43433
#> after refine_(*, conv=F):
#> beta_ref : -1.67229 1.44956
#> with ||beta_ref - beta_cand|| = 0.0164864177799, --> sc = 0.348265804600927
#> Sample[311]: idc = 2 9
#> b^[] = 0.553075 0.667465
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 0.553075 0.667465
#> after refine_(*, conv=F):
#> beta_ref : 0.207049 0.814328
#> with ||beta_ref - beta_cand|| = 0.375903533209, --> sc = 0.577742866999083
#> Sample[312]: idc = 7 0
#> b^[] = 7.13585 -1.51825
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 7.13585 -1.51825
#> after refine_(*, conv=F):
#> beta_ref : 0.224305 0.805575
#> with ||beta_ref - beta_cand|| = 7.29174469856, --> sc = 1.88768873167888
#> Sample[313]: idc = 4 5
#> b^[] = -0.83181 1.18668
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -0.83181 1.18668
#> after refine_(*, conv=F):
#> beta_ref : -0.707317 1.16638
#> with ||beta_ref - beta_cand|| = 0.126137655448, --> sc = 0.307085307609751
#> Sample[314]: idc = 5 9
#> b^[] = -0.749016 1.16576
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -0.749016 1.16576
#> after refine_(*, conv=F):
#> beta_ref : -0.670506 1.15677
#> with ||beta_ref - beta_cand|| = 0.0790237477718, --> sc = 0.28633536578424
#> Sample[315]: idc = 11 6
#> b^[] = -2.39849 1.67716
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -2.39849 1.67716
#> after refine_(*, conv=F):
#> beta_ref : -2.18647 1.6102
#> with ||beta_ref - beta_cand|| = 0.222349315351, --> sc = 0.365943372441649
#> Sample[316]: idc = 0 2
#> b^[] = -2.67885 1.47214
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -2.67885 1.47214
#> after refine_(*, conv=F):
#> beta_ref : -1.4879 1.31816
#> with ||beta_ref - beta_cand|| = 1.20086135032, --> sc = 1.37247806701654
#> Sample[317]: idc = 8 6
#> b^[] = -2.69667 1.88507
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -2.69667 1.88507
#> after refine_(*, conv=F):
#> beta_ref : -1.94237 1.55803
#> with ||beta_ref - beta_cand|| = 0.822144059656, --> sc = 0.744918088188394
#> Sample[318]: idc = 11 12
#> b^[] = -2.47152 1.70137
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -2.47152 1.70137
#> after refine_(*, conv=F):
#> beta_ref : -2.20041 1.61542
#> with ||beta_ref - beta_cand|| = 0.284411253958, --> sc = 0.366625751619072
#> Sample[319]: idc = 6 4
#> b^[] = -2.33514 1.63299
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -2.33514 1.63299
#> after refine_(*, conv=F):
#> beta_ref : -2.14415 1.58792
#> with ||beta_ref - beta_cand|| = 0.196236405624, --> sc = 0.399338780038278
#> Sample[320]: idc = 9 11
#> b^[] = -0.0569989 0.900935
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -0.0569989 0.900935
#> after refine_(*, conv=F):
#> beta_ref : -0.177777 0.971152
#> with ||beta_ref - beta_cand|| = 0.139706152895, --> sc = 0.486048896285346
#> Sample[321]: idc = 6 7
#> b^[] = -2.26322 1.58284
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -2.26322 1.58284
#> after refine_(*, conv=F):
#> beta_ref : -2.18436 1.58611
#> with ||beta_ref - beta_cand|| = 0.078929299088, --> sc = 0.361556853894617
#> Sample[322]: idc = 9 6
#> b^[] = -2.77944 1.94279
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -2.77944 1.94279
#> after refine_(*, conv=F):
#> beta_ref : -1.85666 1.53242
#> with ||beta_ref - beta_cand|| = 1.00990893711, --> sc = 0.885734481714681
#> Sample[323]: idc = 11 5
#> b^[] = -1.19798 1.27918
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -1.19798 1.27918
#> after refine_(*, conv=F):
#> beta_ref : -0.898455 1.21902
#> with ||beta_ref - beta_cand|| = 0.305504698816, --> sc = 0.367535728188709
#> Sample[324]: idc = 8 9
#> b^[] = -2.01045 1.6485
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -2.01045 1.6485
#> after refine_(*, conv=F):
#> beta_ref : -1.74249 1.49589
#> with ||beta_ref - beta_cand|| = 0.308377559274, --> sc = 0.579998989395441
#> Sample[325]: idc = 5 7
#> b^[] = -1.81626 1.43537
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -1.81626 1.43537
#> after refine_(*, conv=F):
#> beta_ref : -2.02342 1.54672
#> with ||beta_ref - beta_cand|| = 0.235186747649, --> sc = 0.472769193580754
#> Sample[326]: idc = 8 0
#> b^[] = 7.47514 -1.62162
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 7.47514 -1.62162
#> after refine_(*, conv=F):
#> beta_ref : 0.245606 0.798807
#> with ||beta_ref - beta_cand|| = 7.62395703783, --> sc = 1.90882061471928
#> Sample[327]: idc = 2 6
#> b^[] = -1.78547 1.24971
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -1.78547 1.24971
#> after refine_(*, conv=F):
#> beta_ref : -1.44982 1.30483
#> with ||beta_ref - beta_cand|| = 0.340142551919, --> sc = 1.10255399571346
#> Sample[328]: idc = 5 12
#> b^[] = 1.03191 0.715867
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 1.03191 0.715867
#> after refine_(*, conv=F):
#> beta_ref : -0.376319 1.0606
#> with ||beta_ref - beta_cand|| = 1.4498103741, --> sc = 0.787749674878659
#> Sample[329]: idc = 5 9
#> b^[] = -0.749016 1.16576
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -0.749016 1.16576
#> after refine_(*, conv=F):
#> beta_ref : -0.670506 1.15677
#> with ||beta_ref - beta_cand|| = 0.0790237477718, --> sc = 0.28633536578424
#> Sample[330]: idc = 5 4
#> b^[] = -0.83181 1.18668
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -0.83181 1.18668
#> after refine_(*, conv=F):
#> beta_ref : -0.707317 1.16638
#> with ||beta_ref - beta_cand|| = 0.126137655448, --> sc = 0.307085307609751
#> Sample[331]: idc = 2 8
#> b^[] = 1.56876 0.414582
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 1.56876 0.414582
#> after refine_(*, conv=F):
#> beta_ref : 0.260798 0.824069
#> with ||beta_ref - beta_cand|| = 1.37056338354, --> sc = 0.594013641561037
#> Sample[332]: idc = 1 9
#> b^[] = -1.69558 1.52801
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -1.69558 1.52801
#> after refine_(*, conv=F):
#> beta_ref : -1.38063 1.38331
#> with ||beta_ref - beta_cand|| = 0.346599646578, --> sc = 0.467390770067497
#> Sample[333]: idc = 9 10
#> b^[] = -1.1869 1.33334
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -1.1869 1.33334
#> after refine_(*, conv=F):
#> beta_ref : -0.809481 1.2096
#> with ||beta_ref - beta_cand|| = 0.39718852225, --> sc = 0.284340124502366
#> Sample[334]: idc = 11 8
#> b^[] = 5.54256 -0.955372
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 5.54256 -0.955372
#> after refine_(*, conv=F):
#> beta_ref : 0.410243 0.743998
#> with ||beta_ref - beta_cand|| = 5.4063409827, --> sc = 1.40713960042918
#> Sample[335]: idc = 9 3
#> b^[] = 0.0196235 0.871612
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 0.0196235 0.871612
#> after refine_(*, conv=F):
#> beta_ref : -0.0873117 0.935189
#> with ||beta_ref - beta_cand|| = 0.124406911988, --> sc = 0.456445784317417
#> Sample[336]: idc = 10 12
#> b^[] = -0.773295 1.22367
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -0.773295 1.22367
#> after refine_(*, conv=F):
#> beta_ref : -0.62546 1.16042
#> with ||beta_ref - beta_cand|| = 0.160794576384, --> sc = 0.287874239981522
#> Sample[337]: idc = 11 12
#> b^[] = -2.47152 1.70137
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -2.47152 1.70137
#> after refine_(*, conv=F):
#> beta_ref : -2.20041 1.61542
#> with ||beta_ref - beta_cand|| = 0.284411253958, --> sc = 0.366625751619072
#> Sample[338]: idc = 9 6
#> b^[] = -2.77944 1.94279
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -2.77944 1.94279
#> after refine_(*, conv=F):
#> beta_ref : -1.85666 1.53242
#> with ||beta_ref - beta_cand|| = 1.00990893711, --> sc = 0.885734481714681
#> Sample[339]: idc = 0 7
#> b^[] = 7.13585 -1.51825
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 7.13585 -1.51825
#> after refine_(*, conv=F):
#> beta_ref : 0.224305 0.805575
#> with ||beta_ref - beta_cand|| = 7.29174469856, --> sc = 1.88768873167888
#> Sample[340]: idc = 11 1
#> b^[] = -8.08546 3.56245
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -8.08546 3.56245
#> after refine_(*, conv=F):
#> beta_ref : -1.05124 1.24888
#> with ||beta_ref - beta_cand|| = 7.40491897223, --> sc = 1.76350297340317
#> Sample[341]: idc = 10 6
#> b^[] = -2.3463 1.64077
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -2.3463 1.64077
#> after refine_(*, conv=F):
#> beta_ref : -2.13209 1.58647
#> with ||beta_ref - beta_cand|| = 0.220987541723, --> sc = 0.417145691249323
#> Sample[342]: idc = 0 3
#> b^[] = 1.14468 0.30717
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 1.14468 0.30717
#> after refine_(*, conv=F):
#> beta_ref : -0.307508 0.964066
#> with ||beta_ref - beta_cand|| = 1.59385397212, --> sc = 1.34390327405485
#> Sample[343]: idc = 4 1
#> b^[] = 2.25191 0.271184
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 2.25191 0.271184
#> after refine_(*, conv=F):
#> beta_ref : 0.0105456 0.920859
#> with ||beta_ref - beta_cand|| = 2.33362366183, --> sc = 0.813737418280548
#> Sample[344]: idc = 9 7
#> b^[] = 0.815075 0.5672
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 0.815075 0.5672
#> after refine_(*, conv=F):
#> beta_ref : 0.224676 0.802773
#> with ||beta_ref - beta_cand|| = 0.635661711031, --> sc = 0.679685535172079
#> Sample[345]: idc = 4 12
#> b^[] = -4.26167 2.20494
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -4.26167 2.20494
#> after refine_(*, conv=F):
#> beta_ref : -1.65432 1.43182
#> with ||beta_ref - beta_cand|| = 2.71955618857, --> sc = 0.918893706231788
#> Sample[346]: idc = 10 11
#> b^[] = -2.05854 1.56446
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -2.05854 1.56446
#> after refine_(*, conv=F):
#> beta_ref : -2.07016 1.57129
#> with ||beta_ref - beta_cand|| = 0.0134780521353, --> sc = 0.37327449031415
#> Sample[347]: idc = 6 5
#> b^[] = -2.18551 1.52865
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -2.18551 1.52865
#> after refine_(*, conv=F):
#> beta_ref : -2.06748 1.54965
#> with ||beta_ref - beta_cand|| = 0.119878917294, --> sc = 0.539116852816471
#> Sample[348]: idc = 5 3
#> b^[] = -0.381765 1.07299
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -0.381765 1.07299
#> after refine_(*, conv=F):
#> beta_ref : -0.547363 1.12925
#> with ||beta_ref - beta_cand|| = 0.17489423735, --> sc = 0.235441204593026
#> Sample[349]: idc = 7 3
#> b^[] = 0.263938 0.74904
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 0.263938 0.74904
#> after refine_(*, conv=F):
#> beta_ref : 0.145859 0.83463
#> with ||beta_ref - beta_cand|| = 0.145837154087, --> sc = 0.597892842089404
#> Sample[350]: idc = 10 5
#> b^[] = 3.35407 0.129244
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 3.35407 0.129244
#> after refine_(*, conv=F):
#> beta_ref : -0.22268 0.997273
#> with ||beta_ref - beta_cand|| = 3.68057059515, --> sc = 1.49558843265863
#> Sample[351]: idc = 9 0
#> b^[] = 2.86111 -0.215798
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 2.86111 -0.215798
#> after refine_(*, conv=F):
#> beta_ref : 0.171545 0.81517
#> with ||beta_ref - beta_cand|| = 2.88038890444, --> sc = 1.45332402557657
#> Sample[352]: idc = 11 6
#> b^[] = -2.39849 1.67716
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -2.39849 1.67716
#> after refine_(*, conv=F):
#> beta_ref : -2.18647 1.6102
#> with ||beta_ref - beta_cand|| = 0.222349315351, --> sc = 0.365943372441649
#> Sample[353]: idc = 4 10
#> b^[] = -2.48714 1.67811
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -2.48714 1.67811
#> after refine_(*, conv=F):
#> beta_ref : -2.13634 1.58637
#> with ||beta_ref - beta_cand|| = 0.362596815943, --> sc = 0.436315397353933
#> Sample[354]: idc = 1 2
#> b^[] = 2.63647 0.148748
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 2.63647 0.148748
#> after refine_(*, conv=F):
#> beta_ref : 0.060579 0.901992
#> with ||beta_ref - beta_cand|| = 2.68375998949, --> sc = 0.857121045415408
#> Sample[355]: idc = 6 11
#> b^[] = -2.39849 1.67716
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -2.39849 1.67716
#> after refine_(*, conv=F):
#> beta_ref : -2.18647 1.6102
#> with ||beta_ref - beta_cand|| = 0.222349315351, --> sc = 0.36594337244165
#> Sample[356]: idc = 3 10
#> b^[] = -0.579864 1.17238
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -0.579864 1.17238
#> after refine_(*, conv=F):
#> beta_ref : -0.596172 1.15246
#> with ||beta_ref - beta_cand|| = 0.0257405973865, --> sc = 0.309002235310989
#> Sample[357]: idc = 2 6
#> b^[] = -1.78547 1.24971
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -1.78547 1.24971
#> after refine_(*, conv=F):
#> beta_ref : -1.44982 1.30483
#> with ||beta_ref - beta_cand|| = 0.340142551919, --> sc = 1.10255399571346
#> Sample[358]: idc = 7 10
#> b^[] = -2.81746 1.7657
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -2.81746 1.7657
#> after refine_(*, conv=F):
#> beta_ref : -2.10053 1.5736
#> with ||beta_ref - beta_cand|| = 0.742227256343, --> sc = 0.492521121454336
#> Sample[359]: idc = 6 4
#> b^[] = -2.33514 1.63299
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -2.33514 1.63299
#> after refine_(*, conv=F):
#> beta_ref : -2.14415 1.58792
#> with ||beta_ref - beta_cand|| = 0.196236405624, --> sc = 0.399338780038278
#> Sample[360]: idc = 10 5
#> b^[] = 3.35407 0.129244
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 3.35407 0.129244
#> after refine_(*, conv=F):
#> beta_ref : -0.22268 0.997273
#> with ||beta_ref - beta_cand|| = 3.68057059515, --> sc = 1.49558843265863
#> Sample[361]: idc = 2 10
#> b^[] = 13.188 -2.47835
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 13.188 -2.47835
#> after refine_(*, conv=F):
#> beta_ref : -0.333051 1.01284
#> with ||beta_ref - beta_cand|| = 13.9645138921, --> sc = 4.27637683340343
#> Sample[362]: idc = 9 6
#> b^[] = -2.77944 1.94279
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -2.77944 1.94279
#> after refine_(*, conv=F):
#> beta_ref : -1.85666 1.53242
#> with ||beta_ref - beta_cand|| = 1.00990893711, --> sc = 0.885734481714681
#> Sample[363]: idc = 0 2
#> b^[] = -2.67885 1.47214
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -2.67885 1.47214
#> after refine_(*, conv=F):
#> beta_ref : -1.4879 1.31816
#> with ||beta_ref - beta_cand|| = 1.20086135032, --> sc = 1.37247806701654
#> Sample[364]: idc = 8 3
#> b^[] = -0.471158 1.11784
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -0.471158 1.11784
#> after refine_(*, conv=F):
#> beta_ref : -0.548769 1.1376
#> with ||beta_ref - beta_cand|| = 0.08008750246, --> sc = 0.213513497916963
#> Sample[365]: idc = 11 6
#> b^[] = -2.39849 1.67716
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -2.39849 1.67716
#> after refine_(*, conv=F):
#> beta_ref : -2.18647 1.6102
#> with ||beta_ref - beta_cand|| = 0.222349315351, --> sc = 0.365943372441649
#> Sample[366]: idc = 4 3
#> b^[] = -0.284513 1.0242
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -0.284513 1.0242
#> after refine_(*, conv=F):
#> beta_ref : -0.565247 1.12529
#> with ||beta_ref - beta_cand|| = 0.298379314564, --> sc = 0.280155905387692
#> Sample[367]: idc = 9 6
#> b^[] = -2.77944 1.94279
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -2.77944 1.94279
#> after refine_(*, conv=F):
#> beta_ref : -1.85666 1.53242
#> with ||beta_ref - beta_cand|| = 1.00990893711, --> sc = 0.885734481714681
#> Sample[368]: idc = 0 7
#> b^[] = 7.13585 -1.51825
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 7.13585 -1.51825
#> after refine_(*, conv=F):
#> beta_ref : 0.224305 0.805575
#> with ||beta_ref - beta_cand|| = 7.29174469856, --> sc = 1.88768873167888
#> Sample[369]: idc = 9 2
#> b^[] = 0.553075 0.667465
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 0.553075 0.667465
#> after refine_(*, conv=F):
#> beta_ref : 0.207049 0.814328
#> with ||beta_ref - beta_cand|| = 0.375903533209, --> sc = 0.577742866999082
#> Sample[370]: idc = 10 8
#> b^[] = -0.79425 1.22922
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -0.79425 1.22922
#> after refine_(*, conv=F):
#> beta_ref : -0.629979 1.16172
#> with ||beta_ref - beta_cand|| = 0.177599485206, --> sc = 0.285059041656878
#> Sample[371]: idc = 11 6
#> b^[] = -2.39849 1.67716
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -2.39849 1.67716
#> after refine_(*, conv=F):
#> beta_ref : -2.18647 1.6102
#> with ||beta_ref - beta_cand|| = 0.222349315351, --> sc = 0.365943372441649
#> Sample[372]: idc = 0 3
#> b^[] = 1.14468 0.30717
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 1.14468 0.30717
#> after refine_(*, conv=F):
#> beta_ref : -0.307508 0.964066
#> with ||beta_ref - beta_cand|| = 1.59385397212, --> sc = 1.34390327405485
#> Sample[373]: idc = 6 8
#> b^[] = -2.69667 1.88507
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -2.69667 1.88507
#> after refine_(*, conv=F):
#> beta_ref : -1.94237 1.55803
#> with ||beta_ref - beta_cand|| = 0.822144059656, --> sc = 0.744918088188393
#> Sample[374]: idc = 12 11
#> b^[] = -2.47152 1.70137
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -2.47152 1.70137
#> after refine_(*, conv=F):
#> beta_ref : -2.20041 1.61542
#> with ||beta_ref - beta_cand|| = 0.284411253958, --> sc = 0.366625751619073
#> Sample[375]: idc = 11 6
#> b^[] = -2.39849 1.67716
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -2.39849 1.67716
#> after refine_(*, conv=F):
#> beta_ref : -2.18647 1.6102
#> with ||beta_ref - beta_cand|| = 0.222349315351, --> sc = 0.365943372441649
#> Sample[376]: idc = 1 9
#> b^[] = -1.69558 1.52801
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -1.69558 1.52801
#> after refine_(*, conv=F):
#> beta_ref : -1.38063 1.38331
#> with ||beta_ref - beta_cand|| = 0.346599646578, --> sc = 0.467390770067497
#> Sample[377]: idc = 0 7
#> b^[] = 7.13585 -1.51825
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 7.13585 -1.51825
#> after refine_(*, conv=F):
#> beta_ref : 0.224305 0.805575
#> with ||beta_ref - beta_cand|| = 7.29174469856, --> sc = 1.88768873167888
#> Sample[378]: idc = 5 10
#> b^[] = 3.35407 0.129244
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 3.35407 0.129244
#> after refine_(*, conv=F):
#> beta_ref : -0.22268 0.997273
#> with ||beta_ref - beta_cand|| = 3.68057059515, --> sc = 1.49558843265863
#> Sample[379]: idc = 10 3
#> b^[] = -0.579864 1.17238
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -0.579864 1.17238
#> after refine_(*, conv=F):
#> beta_ref : -0.596172 1.15246
#> with ||beta_ref - beta_cand|| = 0.0257405973865, --> sc = 0.30900223531099
#> Sample[380]: idc = 12 0
#> b^[] = -14.974 5.21828
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -14.974 5.21828
#> after refine_(*, conv=F):
#> beta_ref : -1.04426 1.21824
#> with ||beta_ref - beta_cand|| = 14.4926713752, --> sc = 2.94185072560702
#> Sample[381]: idc = 3 5
#> b^[] = -0.381765 1.07299
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -0.381765 1.07299
#> after refine_(*, conv=F):
#> beta_ref : -0.547363 1.12925
#> with ||beta_ref - beta_cand|| = 0.17489423735, --> sc = 0.235441204593026
#> Sample[382]: idc = 5 8
#> b^[] = -0.229469 1.03452
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -0.229469 1.03452
#> after refine_(*, conv=F):
#> beta_ref : -0.556193 1.13403
#> with ||beta_ref - beta_cand|| = 0.341541122565, --> sc = 0.270706285978508
#> Sample[383]: idc = 5 12
#> b^[] = 1.03191 0.715867
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 1.03191 0.715867
#> after refine_(*, conv=F):
#> beta_ref : -0.376319 1.0606
#> with ||beta_ref - beta_cand|| = 1.4498103741, --> sc = 0.787749674878659
#> Sample[384]: idc = 3 4
#> b^[] = -0.284513 1.0242
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -0.284513 1.0242
#> after refine_(*, conv=F):
#> beta_ref : -0.565247 1.12529
#> with ||beta_ref - beta_cand|| = 0.298379314564, --> sc = 0.280155905387692
#> Sample[385]: idc = 12 0
#> b^[] = -14.974 5.21828
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -14.974 5.21828
#> after refine_(*, conv=F):
#> beta_ref : -1.04426 1.21824
#> with ||beta_ref - beta_cand|| = 14.4926713752, --> sc = 2.94185072560702
#> Sample[386]: idc = 7 9
#> b^[] = 0.815075 0.5672
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 0.815075 0.5672
#> after refine_(*, conv=F):
#> beta_ref : 0.224676 0.802773
#> with ||beta_ref - beta_cand|| = 0.635661711031, --> sc = 0.67968553517208
#> Sample[387]: idc = 12 1
#> b^[] = -0.484845 1.14253
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -0.484845 1.14253
#> after refine_(*, conv=F):
#> beta_ref : -0.587028 1.14935
#> with ||beta_ref - beta_cand|| = 0.102411328665, --> sc = 0.312564686828127
#> Sample[388]: idc = 1 9
#> b^[] = -1.69558 1.52801
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -1.69558 1.52801
#> after refine_(*, conv=F):
#> beta_ref : -1.38063 1.38331
#> with ||beta_ref - beta_cand|| = 0.346599646578, --> sc = 0.467390770067497
#> Sample[389]: idc = 0 5
#> b^[] = -6.1579 2.53215
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -6.1579 2.53215
#> after refine_(*, conv=F):
#> beta_ref : -0.763906 1.12366
#> with ||beta_ref - beta_cand|| = 5.57485454363, --> sc = 1.64015945823134
#> Sample[390]: idc = 11 3
#> b^[] = -0.00341923 0.883173
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -0.00341923 0.883173
#> after refine_(*, conv=F):
#> beta_ref : -0.142669 0.957455
#> with ||beta_ref - beta_cand|| = 0.157824084971, --> sc = 0.473157356422566
#> Sample[391]: idc = 7 10
#> b^[] = -2.81746 1.7657
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -2.81746 1.7657
#> after refine_(*, conv=F):
#> beta_ref : -2.10053 1.5736
#> with ||beta_ref - beta_cand|| = 0.742227256343, --> sc = 0.492521121454336
#> Sample[392]: idc = 5 3
#> b^[] = -0.381765 1.07299
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -0.381765 1.07299
#> after refine_(*, conv=F):
#> beta_ref : -0.547363 1.12925
#> with ||beta_ref - beta_cand|| = 0.17489423735, --> sc = 0.235441204593026
#> Sample[393]: idc = 6 9
#> b^[] = -2.77944 1.94279
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -2.77944 1.94279
#> after refine_(*, conv=F):
#> beta_ref : -1.85666 1.53242
#> with ||beta_ref - beta_cand|| = 1.00990893711, --> sc = 0.885734481714681
#> Sample[394]: idc = 7 4
#> b^[] = -3.13442 1.87028
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -3.13442 1.87028
#> after refine_(*, conv=F):
#> beta_ref : -2.01174 1.5504
#> with ||beta_ref - beta_cand|| = 1.16736173632, --> sc = 0.618257718400314
#> Sample[395]: idc = 12 5
#> b^[] = 1.03191 0.715867
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 1.03191 0.715867
#> after refine_(*, conv=F):
#> beta_ref : -0.376319 1.0606
#> with ||beta_ref - beta_cand|| = 1.4498103741, --> sc = 0.787749674878658
#> Sample[396]: idc = 11 9
#> b^[] = -0.0569989 0.900935
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -0.0569989 0.900935
#> after refine_(*, conv=F):
#> beta_ref : -0.177777 0.971152
#> with ||beta_ref - beta_cand|| = 0.139706152895, --> sc = 0.486048896285347
#> Sample[397]: idc = 7 1
#> b^[] = -13.1628 5.179
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -13.1628 5.179
#> after refine_(*, conv=F):
#> beta_ref : -0.962999 1.2121
#> with ||beta_ref - beta_cand|| = 12.8285446507, --> sc = 2.75529751980209
#> Sample[398]: idc = 10 3
#> b^[] = -0.579864 1.17238
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -0.579864 1.17238
#> after refine_(*, conv=F):
#> beta_ref : -0.596172 1.15246
#> with ||beta_ref - beta_cand|| = 0.0257405973865, --> sc = 0.30900223531099
#> Sample[399]: idc = 5 2
#> b^[] = 47.0534 -10.9101
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 47.0534 -10.9101
#> after refine_(*, conv=F):
#> beta_ref : -0.397375 1.02753
#> with ||beta_ref - beta_cand|| = 48.929345941, --> sc = 15.1481999432314
#> Sample[400]: idc = 0 2
#> b^[] = -2.67885 1.47214
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -2.67885 1.47214
#> after refine_(*, conv=F):
#> beta_ref : -1.4879 1.31816
#> with ||beta_ref - beta_cand|| = 1.20086135032, --> sc = 1.37247806701654
#> Sample[401]: idc = 7 1
#> b^[] = -13.1628 5.179
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -13.1628 5.179
#> after refine_(*, conv=F):
#> beta_ref : -0.962999 1.2121
#> with ||beta_ref - beta_cand|| = 12.8285446507, --> sc = 2.75529751980209
#> Sample[402]: idc = 6 9
#> b^[] = -2.77944 1.94279
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -2.77944 1.94279
#> after refine_(*, conv=F):
#> beta_ref : -1.85666 1.53242
#> with ||beta_ref - beta_cand|| = 1.00990893711, --> sc = 0.885734481714681
#> Sample[403]: idc = 3 6
#> b^[] = -4.48261 3.13038
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -4.48261 3.13038
#> after refine_(*, conv=F):
#> beta_ref : -1.5156 1.39333
#> with ||beta_ref - beta_cand|| = 3.43809929004, --> sc = 3.95137685252739
#> Sample[404]: idc = 6 12
#> b^[] = -2.40731 1.68331
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -2.40731 1.68331
#> after refine_(*, conv=F):
#> beta_ref : -2.20087 1.61624
#> with ||beta_ref - beta_cand|| = 0.2170640588, --> sc = 0.35543328734486
#> Sample[405]: idc = 8 6
#> b^[] = -2.69667 1.88507
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -2.69667 1.88507
#> after refine_(*, conv=F):
#> beta_ref : -1.94237 1.55803
#> with ||beta_ref - beta_cand|| = 0.822144059656, --> sc = 0.744918088188394
#> Sample[406]: idc = 3 9
#> b^[] = 0.0196235 0.871612
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 0.0196235 0.871612
#> after refine_(*, conv=F):
#> beta_ref : -0.0873117 0.935189
#> with ||beta_ref - beta_cand|| = 0.124406911988, --> sc = 0.456445784317417
#> Sample[407]: idc = 11 0
#> b^[] = 8.42881 -1.91219
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 8.42881 -1.91219
#> after refine_(*, conv=F):
#> beta_ref : 0.302657 0.780686
#> with ||beta_ref - beta_cand|| = 8.56072293339, --> sc = 1.9669857749306
#> Sample[408]: idc = 8 12
#> b^[] = -0.799579 1.23106
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -0.799579 1.23106
#> after refine_(*, conv=F):
#> beta_ref : -0.631179 1.16203
#> with ||beta_ref - beta_cand|| = 0.181998057096, --> sc = 0.286399245856023
#> Sample[409]: idc = 4 1
#> b^[] = 2.25191 0.271184
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 2.25191 0.271184
#> after refine_(*, conv=F):
#> beta_ref : 0.0105456 0.920859
#> with ||beta_ref - beta_cand|| = 2.33362366183, --> sc = 0.813737418280548
#> Sample[410]: idc = 7 10
#> b^[] = -2.81746 1.7657
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -2.81746 1.7657
#> after refine_(*, conv=F):
#> beta_ref : -2.10053 1.5736
#> with ||beta_ref - beta_cand|| = 0.742227256343, --> sc = 0.492521121454336
#> Sample[411]: idc = 3 10
#> b^[] = -0.579864 1.17238
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -0.579864 1.17238
#> after refine_(*, conv=F):
#> beta_ref : -0.596172 1.15246
#> with ||beta_ref - beta_cand|| = 0.0257405973865, --> sc = 0.309002235310989
#> Sample[412]: idc = 8 12
#> b^[] = -0.799579 1.23106
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -0.799579 1.23106
#> after refine_(*, conv=F):
#> beta_ref : -0.631179 1.16203
#> with ||beta_ref - beta_cand|| = 0.181998057096, --> sc = 0.286399245856023
#> Sample[413]: idc = 3 5
#> b^[] = -0.381765 1.07299
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -0.381765 1.07299
#> after refine_(*, conv=F):
#> beta_ref : -0.547363 1.12925
#> with ||beta_ref - beta_cand|| = 0.17489423735, --> sc = 0.235441204593026
#> Sample[414]: idc = 10 5
#> b^[] = 3.35407 0.129244
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 3.35407 0.129244
#> after refine_(*, conv=F):
#> beta_ref : -0.22268 0.997273
#> with ||beta_ref - beta_cand|| = 3.68057059515, --> sc = 1.49558843265863
#> Sample[415]: idc = 1 2
#> b^[] = 2.63647 0.148748
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 2.63647 0.148748
#> after refine_(*, conv=F):
#> beta_ref : 0.060579 0.901992
#> with ||beta_ref - beta_cand|| = 2.68375998949, --> sc = 0.857121045415408
#> Sample[416]: idc = 0 1
#> b^[] = 24.2493 -6.73246
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 24.2493 -6.73246
#> after refine_(*, conv=F):
#> beta_ref : -0.505791 1.05334
#> with ||beta_ref - beta_cand|| = 25.9506189658, --> sc = 5.22853378335815
#> Sample[417]: idc = 0 5
#> b^[] = -6.1579 2.53215
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -6.1579 2.53215
#> after refine_(*, conv=F):
#> beta_ref : -0.763906 1.12366
#> with ||beta_ref - beta_cand|| = 5.57485454363, --> sc = 1.64015945823134
#> Sample[418]: idc = 1 8
#> b^[] = -1.24237 1.38371
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -1.24237 1.38371
#> after refine_(*, conv=F):
#> beta_ref : -0.710666 1.18096
#> with ||beta_ref - beta_cand|| = 0.569048775686, --> sc = 0.404424556196255
#> Sample[419]: idc = 9 12
#> b^[] = -1.25272 1.35853
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -1.25272 1.35853
#> after refine_(*, conv=F):
#> beta_ref : -0.819674 1.21313
#> with ||beta_ref - beta_cand|| = 0.456802788525, --> sc = 0.301200329328142
#> Sample[420]: idc = 1 7
#> b^[] = -13.1628 5.179
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -13.1628 5.179
#> after refine_(*, conv=F):
#> beta_ref : -0.962999 1.2121
#> with ||beta_ref - beta_cand|| = 12.8285446507, --> sc = 2.75529751980209
#> Sample[421]: idc = 3 12
#> b^[] = -0.565704 1.16527
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -0.565704 1.16527
#> after refine_(*, conv=F):
#> beta_ref : -0.589765 1.1512
#> with ||beta_ref - beta_cand|| = 0.0278742526664, --> sc = 0.289862591554884
#> Sample[422]: idc = 8 6
#> b^[] = -2.69667 1.88507
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -2.69667 1.88507
#> after refine_(*, conv=F):
#> beta_ref : -1.94237 1.55803
#> with ||beta_ref - beta_cand|| = 0.822144059656, --> sc = 0.744918088188394
#> Sample[423]: idc = 7 2
#> b^[] = 0.382345 0.709973
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 0.382345 0.709973
#> after refine_(*, conv=F):
#> beta_ref : 0.167816 0.825709
#> with ||beta_ref - beta_cand|| = 0.243757115723, --> sc = 0.613778594777883
#> Sample[424]: idc = 10 3
#> b^[] = -0.579864 1.17238
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -0.579864 1.17238
#> after refine_(*, conv=F):
#> beta_ref : -0.596172 1.15246
#> with ||beta_ref - beta_cand|| = 0.0257405973865, --> sc = 0.30900223531099
#> Sample[425]: idc = 5 1
#> b^[] = 0.17655 0.931947
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 0.17655 0.931947
#> after refine_(*, conv=F):
#> beta_ref : -0.62749 1.15422
#> with ||beta_ref - beta_cand|| = 0.834196899796, --> sc = 0.429182774060504
#> find_scale(*, ini.scale =0.42918277406, tol=1e-10):
#> it | new scale
#> 0 | 0.3783009807
#> 1 | 0.3447919807
#> 2 | 0.3226261856
#> 3 | 0.3078259959
#> 4 | 0.2978417806
#> 5 | 0.2910448963
#> 6 | 0.2863841378
#> 7 | 0.2831706086
#> 8 | 0.2809460201
#> 9 | 0.2794015839
#> 10 | 0.2783271441
#> 11 | 0.2775785869
#> 12 | 0.2770565367
#> 13 | 0.2766921926
#> 14 | 0.2764377848
#> 15 | 0.2762600787
#> 16 | 0.2761359186
#> 17 | 0.2760491552
#> 18 | 0.2759885173
#> 19 | 0.2759461346
#> 20 | 0.2759165096
#> 21 | 0.2758958012
#> 22 | 0.2758813253
#> 23 | 0.2758712059
#> 24 | 0.2758641317
#> 25 | 0.2758591865
#> 26 | 0.2758557293
#> 27 | 0.2758533125
#> 28 | 0.2758516230
#> 29 | 0.2758504419
#> 30 | 0.2758496162
#> 31 | 0.2758490390
#> 32 | 0.2758486354
#> 33 | 0.2758483533
#> 34 | 0.2758481561
#> 35 | 0.2758480183
#> 36 | 0.2758479219
#> 37 | 0.2758478545
#> 38 | 0.2758478074
#> 39 | 0.2758477745
#> 40 | 0.2758477514
#> 41 | 0.2758477353
#> 42 | 0.2758477241
#> 43 | 0.2758477162
#> 44 | 0.2758477107
#> 45 | 0.2758477069
#> 46 | 0.2758477042
#> 47 | 0.2758477023
#> 48 | 0.2758477010
#> 49 | 0.2758477001
#> 50 | 0.2758476994
#> 51 | 0.2758476990
#> 52 | 0.2758476987
#> 53 | 0.2758476985
#> 54 | 0.2758476983
#> 55 | 0.2758476982
#> 56 | 0.2758476981
#> 57 | 0.2758476981
#> 58 | 0.2758476980
#> 59 | 0.2758476980
#> found new candidate with scale 0.2758477 in 59 iter (worst sc=0.27585)
#> Sample[426]: idc = 1 8
#> b^[] = -1.24237 1.38371
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -1.24237 1.38371
#> after refine_(*, conv=F):
#> beta_ref : -0.710666 1.18096
#> with ||beta_ref - beta_cand|| = 0.569048775686, --> sc = 0.404424556196255
#> Sample[427]: idc = 4 12
#> b^[] = -4.26167 2.20494
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -4.26167 2.20494
#> after refine_(*, conv=F):
#> beta_ref : -1.65432 1.43182
#> with ||beta_ref - beta_cand|| = 2.71955618857, --> sc = 0.918893706231788
#> Sample[428]: idc = 4 2
#> b^[] = 2.8091 0.105765
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 2.8091 0.105765
#> after refine_(*, conv=F):
#> beta_ref : 0.030844 0.911664
#> with ||beta_ref - beta_cand|| = 2.89278369939, --> sc = 0.908806025778876
#> Sample[429]: idc = 5 2
#> b^[] = 47.0534 -10.9101
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 47.0534 -10.9101
#> after refine_(*, conv=F):
#> beta_ref : -0.397375 1.02753
#> with ||beta_ref - beta_cand|| = 48.929345941, --> sc = 15.1481999432314
#> Sample[430]: idc = 0 2
#> b^[] = -2.67885 1.47214
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -2.67885 1.47214
#> after refine_(*, conv=F):
#> beta_ref : -1.4879 1.31816
#> with ||beta_ref - beta_cand|| = 1.20086135032, --> sc = 1.37247806701654
#> Sample[431]: idc = 8 6
#> b^[] = -2.69667 1.88507
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -2.69667 1.88507
#> after refine_(*, conv=F):
#> beta_ref : -1.94237 1.55803
#> with ||beta_ref - beta_cand|| = 0.822144059656, --> sc = 0.744918088188394
#> Sample[432]: idc = 0 7
#> b^[] = 7.13585 -1.51825
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 7.13585 -1.51825
#> after refine_(*, conv=F):
#> beta_ref : 0.224305 0.805575
#> with ||beta_ref - beta_cand|| = 7.29174469856, --> sc = 1.88768873167888
#> Sample[433]: idc = 9 3
#> b^[] = 0.0196235 0.871612
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 0.0196235 0.871612
#> after refine_(*, conv=F):
#> beta_ref : -0.0873117 0.935189
#> with ||beta_ref - beta_cand|| = 0.124406911988, --> sc = 0.456445784317417
#> Sample[434]: idc = 2 10
#> b^[] = 13.188 -2.47835
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 13.188 -2.47835
#> after refine_(*, conv=F):
#> beta_ref : -0.333051 1.01284
#> with ||beta_ref - beta_cand|| = 13.9645138921, --> sc = 4.27637683340343
#> Sample[435]: idc = 12 6
#> b^[] = -2.40731 1.68331
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -2.40731 1.68331
#> after refine_(*, conv=F):
#> beta_ref : -2.20087 1.61624
#> with ||beta_ref - beta_cand|| = 0.2170640588, --> sc = 0.355433287344859
#> Sample[436]: idc = 6 7
#> b^[] = -2.26322 1.58284
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -2.26322 1.58284
#> after refine_(*, conv=F):
#> beta_ref : -2.18436 1.58611
#> with ||beta_ref - beta_cand|| = 0.078929299088, --> sc = 0.361556853894617
#> Sample[437]: idc = 0 5
#> b^[] = -6.1579 2.53215
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -6.1579 2.53215
#> after refine_(*, conv=F):
#> beta_ref : -0.763906 1.12366
#> with ||beta_ref - beta_cand|| = 5.57485454363, --> sc = 1.64015945823134
#> Sample[438]: idc = 10 6
#> b^[] = -2.3463 1.64077
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -2.3463 1.64077
#> after refine_(*, conv=F):
#> beta_ref : -2.13209 1.58647
#> with ||beta_ref - beta_cand|| = 0.220987541723, --> sc = 0.417145691249323
#> Sample[439]: idc = 7 3
#> b^[] = 0.263938 0.74904
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 0.263938 0.74904
#> after refine_(*, conv=F):
#> beta_ref : 0.145859 0.83463
#> with ||beta_ref - beta_cand|| = 0.145837154087, --> sc = 0.597892842089404
#> Sample[440]: idc = 12 6
#> b^[] = -2.40731 1.68331
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -2.40731 1.68331
#> after refine_(*, conv=F):
#> beta_ref : -2.20087 1.61624
#> with ||beta_ref - beta_cand|| = 0.2170640588, --> sc = 0.355433287344859
#> Sample[441]: idc = 8 12
#> b^[] = -0.799579 1.23106
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -0.799579 1.23106
#> after refine_(*, conv=F):
#> beta_ref : -0.631179 1.16203
#> with ||beta_ref - beta_cand|| = 0.181998057096, --> sc = 0.286399245856023
#> Sample[442]: idc = 9 4
#> b^[] = -0.706312 1.14942
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -0.706312 1.14942
#> after refine_(*, conv=F):
#> beta_ref : -0.662089 1.15324
#> with ||beta_ref - beta_cand|| = 0.0443876719904, --> sc = 0.302818894470934
#> Sample[443]: idc = 5 2
#> b^[] = 47.0534 -10.9101
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 47.0534 -10.9101
#> after refine_(*, conv=F):
#> beta_ref : -0.397375 1.02753
#> with ||beta_ref - beta_cand|| = 48.929345941, --> sc = 15.1481999432314
#> Sample[444]: idc = 5 3
#> b^[] = -0.381765 1.07299
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -0.381765 1.07299
#> after refine_(*, conv=F):
#> beta_ref : -0.547363 1.12925
#> with ||beta_ref - beta_cand|| = 0.17489423735, --> sc = 0.235441204593026
#> Sample[445]: idc = 5 1
#> b^[] = 0.17655 0.931947
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 0.17655 0.931947
#> after refine_(*, conv=F):
#> beta_ref : -0.62749 1.15422
#> with ||beta_ref - beta_cand|| = 0.834196899796, --> sc = 0.429182774060504
#> find_scale(*, ini.scale =0.42918277406, tol=1e-10):
#> it | new scale
#> 0 | 0.3783009807
#> 1 | 0.3447919807
#> 2 | 0.3226261856
#> 3 | 0.3078259959
#> 4 | 0.2978417806
#> 5 | 0.2910448963
#> 6 | 0.2863841378
#> 7 | 0.2831706086
#> 8 | 0.2809460201
#> 9 | 0.2794015839
#> 10 | 0.2783271441
#> 11 | 0.2775785869
#> 12 | 0.2770565367
#> 13 | 0.2766921926
#> 14 | 0.2764377848
#> 15 | 0.2762600787
#> 16 | 0.2761359186
#> 17 | 0.2760491552
#> 18 | 0.2759885173
#> 19 | 0.2759461346
#> 20 | 0.2759165096
#> 21 | 0.2758958012
#> 22 | 0.2758813253
#> 23 | 0.2758712059
#> 24 | 0.2758641317
#> 25 | 0.2758591865
#> 26 | 0.2758557293
#> 27 | 0.2758533125
#> 28 | 0.2758516230
#> 29 | 0.2758504419
#> 30 | 0.2758496162
#> 31 | 0.2758490390
#> 32 | 0.2758486354
#> 33 | 0.2758483533
#> 34 | 0.2758481561
#> 35 | 0.2758480183
#> 36 | 0.2758479219
#> 37 | 0.2758478545
#> 38 | 0.2758478074
#> 39 | 0.2758477745
#> 40 | 0.2758477514
#> 41 | 0.2758477353
#> 42 | 0.2758477241
#> 43 | 0.2758477162
#> 44 | 0.2758477107
#> 45 | 0.2758477069
#> 46 | 0.2758477042
#> 47 | 0.2758477023
#> 48 | 0.2758477010
#> 49 | 0.2758477001
#> 50 | 0.2758476994
#> 51 | 0.2758476990
#> 52 | 0.2758476987
#> 53 | 0.2758476985
#> 54 | 0.2758476983
#> 55 | 0.2758476982
#> 56 | 0.2758476981
#> 57 | 0.2758476981
#> 58 | 0.2758476980
#> 59 | 0.2758476980
#> found new candidate with scale 0.2758477 in 59 iter (worst sc=0.27585)
#> Sample[446]: idc = 5 2
#> b^[] = 47.0534 -10.9101
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 47.0534 -10.9101
#> after refine_(*, conv=F):
#> beta_ref : -0.397375 1.02753
#> with ||beta_ref - beta_cand|| = 48.929345941, --> sc = 15.1481999432314
#> Sample[447]: idc = 10 2
#> b^[] = 13.188 -2.47835
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 13.188 -2.47835
#> after refine_(*, conv=F):
#> beta_ref : -0.333051 1.01284
#> with ||beta_ref - beta_cand|| = 13.9645138921, --> sc = 4.27637683340342
#> Sample[448]: idc = 1 11
#> b^[] = -8.08546 3.56245
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -8.08546 3.56245
#> after refine_(*, conv=F):
#> beta_ref : -1.05124 1.24888
#> with ||beta_ref - beta_cand|| = 7.40491897223, --> sc = 1.76350297340318
#> Sample[449]: idc = 7 2
#> b^[] = 0.382345 0.709973
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 0.382345 0.709973
#> after refine_(*, conv=F):
#> beta_ref : 0.167816 0.825709
#> with ||beta_ref - beta_cand|| = 0.243757115723, --> sc = 0.613778594777883
#> Sample[450]: idc = 7 6
#> b^[] = -2.26322 1.58284
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -2.26322 1.58284
#> after refine_(*, conv=F):
#> beta_ref : -2.18436 1.58611
#> with ||beta_ref - beta_cand|| = 0.078929299088, --> sc = 0.361556853894617
#> Sample[451]: idc = 0 5
#> b^[] = -6.1579 2.53215
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -6.1579 2.53215
#> after refine_(*, conv=F):
#> beta_ref : -0.763906 1.12366
#> with ||beta_ref - beta_cand|| = 5.57485454363, --> sc = 1.64015945823134
#> Sample[452]: idc = 5 1
#> b^[] = 0.17655 0.931947
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 0.17655 0.931947
#> after refine_(*, conv=F):
#> beta_ref : -0.62749 1.15422
#> with ||beta_ref - beta_cand|| = 0.834196899796, --> sc = 0.429182774060504
#> find_scale(*, ini.scale =0.42918277406, tol=1e-10):
#> it | new scale
#> 0 | 0.3783009807
#> 1 | 0.3447919807
#> 2 | 0.3226261856
#> 3 | 0.3078259959
#> 4 | 0.2978417806
#> 5 | 0.2910448963
#> 6 | 0.2863841378
#> 7 | 0.2831706086
#> 8 | 0.2809460201
#> 9 | 0.2794015839
#> 10 | 0.2783271441
#> 11 | 0.2775785869
#> 12 | 0.2770565367
#> 13 | 0.2766921926
#> 14 | 0.2764377848
#> 15 | 0.2762600787
#> 16 | 0.2761359186
#> 17 | 0.2760491552
#> 18 | 0.2759885173
#> 19 | 0.2759461346
#> 20 | 0.2759165096
#> 21 | 0.2758958012
#> 22 | 0.2758813253
#> 23 | 0.2758712059
#> 24 | 0.2758641317
#> 25 | 0.2758591865
#> 26 | 0.2758557293
#> 27 | 0.2758533125
#> 28 | 0.2758516230
#> 29 | 0.2758504419
#> 30 | 0.2758496162
#> 31 | 0.2758490390
#> 32 | 0.2758486354
#> 33 | 0.2758483533
#> 34 | 0.2758481561
#> 35 | 0.2758480183
#> 36 | 0.2758479219
#> 37 | 0.2758478545
#> 38 | 0.2758478074
#> 39 | 0.2758477745
#> 40 | 0.2758477514
#> 41 | 0.2758477353
#> 42 | 0.2758477241
#> 43 | 0.2758477162
#> 44 | 0.2758477107
#> 45 | 0.2758477069
#> 46 | 0.2758477042
#> 47 | 0.2758477023
#> 48 | 0.2758477010
#> 49 | 0.2758477001
#> 50 | 0.2758476994
#> 51 | 0.2758476990
#> 52 | 0.2758476987
#> 53 | 0.2758476985
#> 54 | 0.2758476983
#> 55 | 0.2758476982
#> 56 | 0.2758476981
#> 57 | 0.2758476981
#> 58 | 0.2758476980
#> 59 | 0.2758476980
#> found new candidate with scale 0.2758477 in 59 iter (worst sc=0.27585)
#> Sample[453]: idc = 3 4
#> b^[] = -0.284513 1.0242
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -0.284513 1.0242
#> after refine_(*, conv=F):
#> beta_ref : -0.565247 1.12529
#> with ||beta_ref - beta_cand|| = 0.298379314564, --> sc = 0.280155905387692
#> Sample[454]: idc = 3 4
#> b^[] = -0.284513 1.0242
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -0.284513 1.0242
#> after refine_(*, conv=F):
#> beta_ref : -0.565247 1.12529
#> with ||beta_ref - beta_cand|| = 0.298379314564, --> sc = 0.280155905387692
#> Sample[455]: idc = 4 12
#> b^[] = -4.26167 2.20494
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -4.26167 2.20494
#> after refine_(*, conv=F):
#> beta_ref : -1.65432 1.43182
#> with ||beta_ref - beta_cand|| = 2.71955618857, --> sc = 0.918893706231788
#> Sample[456]: idc = 11 7
#> b^[] = 29.1679 -8.78739
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 29.1679 -8.78739
#> after refine_(*, conv=F):
#> beta_ref : -0.244111 0.966143
#> with ||beta_ref - beta_cand|| = 30.9870463767, --> sc = 6.48063687575041
#> Sample[457]: idc = 0 12
#> b^[] = -14.974 5.21828
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -14.974 5.21828
#> after refine_(*, conv=F):
#> beta_ref : -1.04426 1.21824
#> with ||beta_ref - beta_cand|| = 14.4926713752, --> sc = 2.94185072560702
#> Sample[458]: idc = 3 7
#> b^[] = 0.263938 0.74904
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 0.263938 0.74904
#> after refine_(*, conv=F):
#> beta_ref : 0.145859 0.83463
#> with ||beta_ref - beta_cand|| = 0.145837154087, --> sc = 0.597892842089404
#> Sample[459]: idc = 7 9
#> b^[] = 0.815075 0.5672
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 0.815075 0.5672
#> after refine_(*, conv=F):
#> beta_ref : 0.224676 0.802773
#> with ||beta_ref - beta_cand|| = 0.635661711031, --> sc = 0.67968553517208
#> Sample[460]: idc = 12 3
#> b^[] = -0.565704 1.16527
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -0.565704 1.16527
#> after refine_(*, conv=F):
#> beta_ref : -0.589765 1.1512
#> with ||beta_ref - beta_cand|| = 0.0278742526664, --> sc = 0.289862591554884
#> Sample[461]: idc = 11 4
#> b^[] = -1.66598 1.43433
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -1.66598 1.43433
#> after refine_(*, conv=F):
#> beta_ref : -1.67229 1.44956
#> with ||beta_ref - beta_cand|| = 0.0164864177799, --> sc = 0.348265804600926
#> Sample[462]: idc = 3 12
#> b^[] = -0.565704 1.16527
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -0.565704 1.16527
#> after refine_(*, conv=F):
#> beta_ref : -0.589765 1.1512
#> with ||beta_ref - beta_cand|| = 0.0278742526664, --> sc = 0.289862591554884
#> Sample[463]: idc = 5 0
#> b^[] = -6.1579 2.53215
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -6.1579 2.53215
#> after refine_(*, conv=F):
#> beta_ref : -0.763906 1.12366
#> with ||beta_ref - beta_cand|| = 5.57485454363, --> sc = 1.64015945823134
#> Sample[464]: idc = 0 2
#> b^[] = -2.67885 1.47214
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -2.67885 1.47214
#> after refine_(*, conv=F):
#> beta_ref : -1.4879 1.31816
#> with ||beta_ref - beta_cand|| = 1.20086135032, --> sc = 1.37247806701654
#> Sample[465]: idc = 1 5
#> b^[] = 0.17655 0.931947
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 0.17655 0.931947
#> after refine_(*, conv=F):
#> beta_ref : -0.62749 1.15422
#> with ||beta_ref - beta_cand|| = 0.834196899796, --> sc = 0.429182774060505
#> find_scale(*, ini.scale =0.42918277406, tol=1e-10):
#> it | new scale
#> 0 | 0.3783009807
#> 1 | 0.3447919807
#> 2 | 0.3226261856
#> 3 | 0.3078259959
#> 4 | 0.2978417806
#> 5 | 0.2910448963
#> 6 | 0.2863841378
#> 7 | 0.2831706086
#> 8 | 0.2809460201
#> 9 | 0.2794015839
#> 10 | 0.2783271441
#> 11 | 0.2775785869
#> 12 | 0.2770565367
#> 13 | 0.2766921926
#> 14 | 0.2764377848
#> 15 | 0.2762600787
#> 16 | 0.2761359186
#> 17 | 0.2760491552
#> 18 | 0.2759885173
#> 19 | 0.2759461346
#> 20 | 0.2759165096
#> 21 | 0.2758958012
#> 22 | 0.2758813253
#> 23 | 0.2758712059
#> 24 | 0.2758641317
#> 25 | 0.2758591865
#> 26 | 0.2758557293
#> 27 | 0.2758533125
#> 28 | 0.2758516230
#> 29 | 0.2758504419
#> 30 | 0.2758496162
#> 31 | 0.2758490390
#> 32 | 0.2758486354
#> 33 | 0.2758483533
#> 34 | 0.2758481561
#> 35 | 0.2758480183
#> 36 | 0.2758479219
#> 37 | 0.2758478545
#> 38 | 0.2758478074
#> 39 | 0.2758477745
#> 40 | 0.2758477514
#> 41 | 0.2758477353
#> 42 | 0.2758477241
#> 43 | 0.2758477162
#> 44 | 0.2758477107
#> 45 | 0.2758477069
#> 46 | 0.2758477042
#> 47 | 0.2758477023
#> 48 | 0.2758477010
#> 49 | 0.2758477001
#> 50 | 0.2758476994
#> 51 | 0.2758476990
#> 52 | 0.2758476987
#> 53 | 0.2758476985
#> 54 | 0.2758476983
#> 55 | 0.2758476982
#> 56 | 0.2758476981
#> 57 | 0.2758476981
#> 58 | 0.2758476980
#> 59 | 0.2758476980
#> found new candidate with scale 0.2758477 in 59 iter (worst sc=0.27585)
#> Sample[466]: idc = 1 9
#> b^[] = -1.69558 1.52801
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -1.69558 1.52801
#> after refine_(*, conv=F):
#> beta_ref : -1.38063 1.38331
#> with ||beta_ref - beta_cand|| = 0.346599646578, --> sc = 0.467390770067497
#> Sample[467]: idc = 2 8
#> b^[] = 1.56876 0.414582
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 1.56876 0.414582
#> after refine_(*, conv=F):
#> beta_ref : 0.260798 0.824069
#> with ||beta_ref - beta_cand|| = 1.37056338354, --> sc = 0.594013641561037
#> Sample[468]: idc = 3 9
#> b^[] = 0.0196235 0.871612
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 0.0196235 0.871612
#> after refine_(*, conv=F):
#> beta_ref : -0.0873117 0.935189
#> with ||beta_ref - beta_cand|| = 0.124406911988, --> sc = 0.456445784317417
#> Sample[469]: idc = 4 3
#> b^[] = -0.284513 1.0242
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -0.284513 1.0242
#> after refine_(*, conv=F):
#> beta_ref : -0.565247 1.12529
#> with ||beta_ref - beta_cand|| = 0.298379314564, --> sc = 0.280155905387692
#> Sample[470]: idc = 0 11
#> b^[] = 8.42881 -1.91219
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 8.42881 -1.91219
#> after refine_(*, conv=F):
#> beta_ref : 0.302657 0.780686
#> with ||beta_ref - beta_cand|| = 8.56072293339, --> sc = 1.9669857749306
#> Sample[471]: idc = 3 4
#> b^[] = -0.284513 1.0242
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -0.284513 1.0242
#> after refine_(*, conv=F):
#> beta_ref : -0.565247 1.12529
#> with ||beta_ref - beta_cand|| = 0.298379314564, --> sc = 0.280155905387692
#> Sample[472]: idc = 5 0
#> b^[] = -6.1579 2.53215
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -6.1579 2.53215
#> after refine_(*, conv=F):
#> beta_ref : -0.763906 1.12366
#> with ||beta_ref - beta_cand|| = 5.57485454363, --> sc = 1.64015945823134
#> Sample[473]: idc = 2 8
#> b^[] = 1.56876 0.414582
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 1.56876 0.414582
#> after refine_(*, conv=F):
#> beta_ref : 0.260798 0.824069
#> with ||beta_ref - beta_cand|| = 1.37056338354, --> sc = 0.594013641561037
#> Sample[474]: idc = 5 2
#> b^[] = 47.0534 -10.9101
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 47.0534 -10.9101
#> after refine_(*, conv=F):
#> beta_ref : -0.397375 1.02753
#> with ||beta_ref - beta_cand|| = 48.929345941, --> sc = 15.1481999432314
#> Sample[475]: idc = 4 0
#> b^[] = -36.3623 11.735
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -36.3623 11.735
#> after refine_(*, conv=F):
#> beta_ref : -1.03377 1.21864
#> with ||beta_ref - beta_cand|| = 36.8605343085, --> sc = 6.75365155222259
#> Sample[476]: idc = 5 8
#> b^[] = -0.229469 1.03452
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -0.229469 1.03452
#> after refine_(*, conv=F):
#> beta_ref : -0.556193 1.13403
#> with ||beta_ref - beta_cand|| = 0.341541122565, --> sc = 0.270706285978508
#> Sample[477]: idc = 7 0
#> b^[] = 7.13585 -1.51825
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 7.13585 -1.51825
#> after refine_(*, conv=F):
#> beta_ref : 0.224305 0.805575
#> with ||beta_ref - beta_cand|| = 7.29174469856, --> sc = 1.88768873167888
#> Sample[478]: idc = 2 4
#> b^[] = 2.8091 0.105765
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 2.8091 0.105765
#> after refine_(*, conv=F):
#> beta_ref : 0.030844 0.911664
#> with ||beta_ref - beta_cand|| = 2.89278369939, --> sc = 0.908806025778877
#> Sample[479]: idc = 11 2
#> b^[] = 0.931408 0.573269
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 0.931408 0.573269
#> after refine_(*, conv=F):
#> beta_ref : 0.290339 0.794023
#> with ||beta_ref - beta_cand|| = 0.678012796036, --> sc = 0.523721624955058
#> Sample[480]: idc = 11 5
#> b^[] = -1.19798 1.27918
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -1.19798 1.27918
#> after refine_(*, conv=F):
#> beta_ref : -0.898455 1.21902
#> with ||beta_ref - beta_cand|| = 0.305504698816, --> sc = 0.367535728188709
#> Sample[481]: idc = 1 6
#> b^[] = -2.59549 1.81452
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -2.59549 1.81452
#> after refine_(*, conv=F):
#> beta_ref : -2.14459 1.61896
#> with ||beta_ref - beta_cand|| = 0.491482052855, --> sc = 0.553607074103745
#> Sample[482]: idc = 11 5
#> b^[] = -1.19798 1.27918
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -1.19798 1.27918
#> after refine_(*, conv=F):
#> beta_ref : -0.898455 1.21902
#> with ||beta_ref - beta_cand|| = 0.305504698816, --> sc = 0.367535728188709
#> Sample[483]: idc = 11 8
#> b^[] = 5.54256 -0.955372
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 5.54256 -0.955372
#> after refine_(*, conv=F):
#> beta_ref : 0.410243 0.743998
#> with ||beta_ref - beta_cand|| = 5.4063409827, --> sc = 1.40713960042918
#> Sample[484]: idc = 0 10
#> b^[] = -9.17664 3.45192
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -9.17664 3.45192
#> after refine_(*, conv=F):
#> beta_ref : -0.840132 1.1533
#> with ||beta_ref - beta_cand|| = 8.64759632726, --> sc = 2.16009257867086
#> Sample[485]: idc = 12 3
#> b^[] = -0.565704 1.16527
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -0.565704 1.16527
#> after refine_(*, conv=F):
#> beta_ref : -0.589765 1.1512
#> with ||beta_ref - beta_cand|| = 0.0278742526664, --> sc = 0.289862591554884
#> Sample[486]: idc = 9 1
#> b^[] = -1.69558 1.52801
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -1.69558 1.52801
#> after refine_(*, conv=F):
#> beta_ref : -1.38063 1.38331
#> with ||beta_ref - beta_cand|| = 0.346599646578, --> sc = 0.467390770067498
#> Sample[487]: idc = 9 0
#> b^[] = 2.86111 -0.215798
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 2.86111 -0.215798
#> after refine_(*, conv=F):
#> beta_ref : 0.171545 0.81517
#> with ||beta_ref - beta_cand|| = 2.88038890444, --> sc = 1.45332402557657
#> Sample[488]: idc = 7 11
#> b^[] = 29.1679 -8.78739
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 29.1679 -8.78739
#> after refine_(*, conv=F):
#> beta_ref : -0.244111 0.966143
#> with ||beta_ref - beta_cand|| = 30.9870463767, --> sc = 6.48063687575061
#> Sample[489]: idc = 11 3
#> b^[] = -0.00341923 0.883173
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -0.00341923 0.883173
#> after refine_(*, conv=F):
#> beta_ref : -0.142669 0.957455
#> with ||beta_ref - beta_cand|| = 0.157824084971, --> sc = 0.473157356422566
#> Sample[490]: idc = 10 4
#> b^[] = -2.48714 1.67811
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -2.48714 1.67811
#> after refine_(*, conv=F):
#> beta_ref : -2.13634 1.58637
#> with ||beta_ref - beta_cand|| = 0.362596815943, --> sc = 0.436315397353933
#> Sample[491]: idc = 11 4
#> b^[] = -1.66598 1.43433
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -1.66598 1.43433
#> after refine_(*, conv=F):
#> beta_ref : -1.67229 1.44956
#> with ||beta_ref - beta_cand|| = 0.0164864177799, --> sc = 0.348265804600926
#> Sample[492]: idc = 2 7
#> b^[] = 0.382345 0.709973
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 0.382345 0.709973
#> after refine_(*, conv=F):
#> beta_ref : 0.167816 0.825709
#> with ||beta_ref - beta_cand|| = 0.243757115723, --> sc = 0.613778594777883
#> Sample[493]: idc = 0 4
#> b^[] = -36.3623 11.735
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -36.3623 11.735
#> after refine_(*, conv=F):
#> beta_ref : -1.03377 1.21864
#> with ||beta_ref - beta_cand|| = 36.8605343085, --> sc = 6.75365155222264
#> Sample[494]: idc = 3 7
#> b^[] = 0.263938 0.74904
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 0.263938 0.74904
#> after refine_(*, conv=F):
#> beta_ref : 0.145859 0.83463
#> with ||beta_ref - beta_cand|| = 0.145837154087, --> sc = 0.597892842089404
#> Sample[495]: idc = 3 10
#> b^[] = -0.579864 1.17238
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -0.579864 1.17238
#> after refine_(*, conv=F):
#> beta_ref : -0.596172 1.15246
#> with ||beta_ref - beta_cand|| = 0.0257405973865, --> sc = 0.309002235310989
#> Sample[496]: idc = 3 8
#> b^[] = -0.471158 1.11784
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -0.471158 1.11784
#> after refine_(*, conv=F):
#> beta_ref : -0.548769 1.1376
#> with ||beta_ref - beta_cand|| = 0.08008750246, --> sc = 0.213513497916964
#> Sample[497]: idc = 5 11
#> b^[] = -1.19798 1.27918
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -1.19798 1.27918
#> after refine_(*, conv=F):
#> beta_ref : -0.898455 1.21902
#> with ||beta_ref - beta_cand|| = 0.305504698816, --> sc = 0.367535728188709
#> Sample[498]: idc = 3 5
#> b^[] = -0.381765 1.07299
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= -0.381765 1.07299
#> after refine_(*, conv=F):
#> beta_ref : -0.547363 1.12925
#> with ||beta_ref - beta_cand|| = 0.17489423735, --> sc = 0.235441204593026
#> Sample[499]: idc = 8 0
#> b^[] = 7.47514 -1.62162
#> refine_fast_s(s0=-1, convChk=FALSE): beta_cand= 7.47514 -1.62162
#> after refine_(*, conv=F):
#> beta_ref : 0.245606 0.798807
#> with ||beta_ref - beta_cand|| = 7.62395703783, --> sc = 1.90882061471928
#> Now refine() to convergence for 2 very best ones:
#> refine_fast_s(s0=0.275848, convChk=TRUE): beta_cand= -0.62749 1.15422
#> refine_() converged after 17 iterations
#> Best[0]: convergence (17 iter.): -> improved scale to 0.275726265760473; bbeta= -0.60771 1.14791
#>
#> refine_fast_s(s0=0.275848, convChk=TRUE): beta_cand= -0.62749 1.15422
#> refine_() converged after 17 iterations
#> Best[1]: convergence (17 iter.)
#> lmrob.S(): scale = 0.275726; coeff.=
#> [1] -0.6077097 1.1479105
#> init converged (remaining method = "M") -> coef=
#> (Intercept) yearly
#> -0.6077097 1.1479105
#> lmrob_MM(): rwls():
#> it 1: L(b1) = 0.217004380265 ||b0 - b1||_1 = 0.0750855
#> it 2: L(b1) = 0.215546067920 ||b0 - b1||_1 = 0.0769669
#> it 3: L(b1) = 0.215065302148 ||b0 - b1||_1 = 0.0773737
#> it 4: L(b1) = 0.214717389542 ||b0 - b1||_1 = 0.0726872
#> it 5: L(b1) = 0.214438109382 ||b0 - b1||_1 = 0.0658683
#> it 6: L(b1) = 0.214220966235 ||b0 - b1||_1 = 0.0580807
#> it 7: L(b1) = 0.214059745231 ||b0 - b1||_1 = 0.0499815
#> it 8: L(b1) = 0.213945290245 ||b0 - b1||_1 = 0.0420661
#> it 9: L(b1) = 0.213867250585 ||b0 - b1||_1 = 0.0347089
#> it 10: L(b1) = 0.213815865828 ||b0 - b1||_1 = 0.028151
#> it 11: L(b1) = 0.213783008597 ||b0 - b1||_1 = 0.0225045
#> it 12: L(b1) = 0.213762496979 ||b0 - b1||_1 = 0.0177781
#> it 13: L(b1) = 0.213749937328 ||b0 - b1||_1 = 0.0139104
#> it 14: L(b1) = 0.213742363795 ||b0 - b1||_1 = 0.0108015
#> it 15: L(b1) = 0.213737851503 ||b0 - b1||_1 = 0.00833734
#> it 16: L(b1) = 0.213735188119 ||b0 - b1||_1 = 0.00640538
#> it 17: L(b1) = 0.213733627376 ||b0 - b1||_1 = 0.00490338
#> it 18: L(b1) = 0.213732717848 ||b0 - b1||_1 = 0.00374317
#> it 19: L(b1) = 0.213732190070 ||b0 - b1||_1 = 0.00285142
#> it 20: L(b1) = 0.213731884805 ||b0 - b1||_1 = 0.00216858
#> it 21: L(b1) = 0.213731708678 ||b0 - b1||_1 = 0.00164722
#> it 22: L(b1) = 0.213731607251 ||b0 - b1||_1 = 0.00125003
#> it 23: L(b1) = 0.213731548924 ||b0 - b1||_1 = 0.000947929
#> it 24: L(b1) = 0.213731515420 ||b0 - b1||_1 = 0.000718449
#> it 25: L(b1) = 0.213731496189 ||b0 - b1||_1 = 0.000544298
#> it 26: L(b1) = 0.213731485159 ||b0 - b1||_1 = 0.000412232
#> it 27: L(b1) = 0.213731478835 ||b0 - b1||_1 = 0.000312136
#> it 28: L(b1) = 0.213731475210 ||b0 - b1||_1 = 0.000236303
#> it 29: L(b1) = 0.213731473134 ||b0 - b1||_1 = 0.000178869
#> it 30: L(b1) = 0.213731471944 ||b0 - b1||_1 = 0.00013538
#> it 31: L(b1) = 0.213731471263 ||b0 - b1||_1 = 0.000102457
#> it 32: L(b1) = 0.213731470872 ||b0 - b1||_1 = 7.7536e-05
#> it 33: L(b1) = 0.213731470649 ||b0 - b1||_1 = 5.86739e-05
#> it 34: L(b1) = 0.213731470521 ||b0 - b1||_1 = 4.43989e-05
#> it 35: L(b1) = 0.213731470448 ||b0 - b1||_1 = 3.35961e-05
#> it 36: L(b1) = 0.213731470406 ||b0 - b1||_1 = 2.54212e-05
#> it 37: L(b1) = 0.213731470382 ||b0 - b1||_1 = 1.92353e-05
#> it 38: L(b1) = 0.213731470368 ||b0 - b1||_1 = 1.45544e-05
#> it 39: L(b1) = 0.213731470360 ||b0 - b1||_1 = 1.10125e-05
#> it 40: L(b1) = 0.213731470356 ||b0 - b1||_1 = 8.33255e-06
#> it 41: L(b1) = 0.213731470353 ||b0 - b1||_1 = 6.30473e-06
#> it 42: L(b1) = 0.213731470352 ||b0 - b1||_1 = 4.77038e-06
#> it 43: L(b1) = 0.213731470351 ||b0 - b1||_1 = 3.60943e-06
#> it 44: L(b1) = 0.213731470350 ||b0 - b1||_1 = 2.73101e-06
#> it 45: L(b1) = 0.213731470350 ||b0 - b1||_1 = 2.06636e-06
#> it 46: L(b1) = 0.213731470350 ||b0 - b1||_1 = 1.56347e-06
#> it 47: L(b1) = 0.213731470350 ||b0 - b1||_1 = 1.18297e-06
#> it 48: L(b1) = 0.213731470350 ||b0 - b1||_1 = 8.95067e-07
#> it 49: L(b1) = 0.213731470350 ||b0 - b1||_1 = 6.77233e-07
#> it 50: L(b1) = 0.213731470350 ||b0 - b1||_1 = 5.12413e-07
#> it 51: L(b1) = 0.213731470350 ||b0 - b1||_1 = 3.87706e-07
#> it 52: L(b1) = 0.213731470350 ||b0 - b1||_1 = 2.93349e-07
#> it 53: L(b1) = 0.213731470350 ||b0 - b1||_1 = 2.21956e-07
#> rwls() used 53 it.; last ||b0 - b1||_1 = 2.21956e-07, L(b1) = 0.21373147035; convergence
#> lmrob..MM..fit(*, obj) --> updating .. qr(x * rweights) -> rank= 2, outlierStats()
#> step "M" -> new coef=
#> (Intercept) yearly
#> -1.149445 1.277453
##-> 53 iter.
## Look at all coefficients:
R.cf <- t(sapply(R, coef))
C.cf <- t(sapply(LS, coef))
plot(C.cf, xlim=range(C.cf[,1],R.cf[,1]),
ylim=range(C.cf[,2],R.cf[,2]))
mD1 <- rowMeans(abs(C.cf - R.cf))
lrg <- mD1 > quantile(mD1, 0.80)
arrows(C.cf[lrg,1], C.cf[lrg,2],
R.cf[lrg,1], R.cf[lrg,2], length=.1, col="light gray")
points(R.cf, col=2)
## All robustness weights
aW <- t(sapply(R, weights, type="robustness"))
colnames(aW) <- names(yearly)
summary(aW)
#> ad ba ef la
#> Min. :0.0000 Min. :0.0184 Min. :0.0000 Min. :0.0000
#> 1st Qu.:0.8183 1st Qu.:0.8741 1st Qu.:0.8546 1st Qu.:0.8661
#> Median :0.9315 Median :0.9441 Median :0.9555 Median :0.9529
#> Mean :0.8397 Mean :0.8877 Mean :0.8675 Mean :0.8727
#> 3rd Qu.:0.9820 3rd Qu.:0.9819 3rd Qu.:0.9905 3rd Qu.:0.9884
#> Max. :1.0000 Max. :1.0000 Max. :1.0000 Max. :1.0000
#> NA's :12 NA's :9 NA's :19 NA's :16
#> lu re ri se
#> Min. :0.5979 Min. :0.0000 Min. :0.0000 Min. :0.000092
#> 1st Qu.:0.9520 1st Qu.:0.8688 1st Qu.:0.8003 1st Qu.:0.923078
#> Median :0.9819 Median :0.9513 Median :0.9317 Median :0.966730
#> Mean :0.9631 Mean :0.8897 Mean :0.7835 Mean :0.933577
#> 3rd Qu.:0.9943 3rd Qu.:0.9911 3rd Qu.:0.9779 3rd Qu.:0.993047
#> Max. :1.0000 Max. :1.0000 Max. :1.0000 Max. :0.999999
#> NA's :7 NA's :18 NA's :14 NA's :18
#> si st su sz
#> Min. :0.0000 Min. :0.05688 Min. :0.001426 Min. :0.0000
#> 1st Qu.:0.8988 1st Qu.:0.92226 1st Qu.:0.884378 1st Qu.:0.9332
#> Median :0.9604 Median :0.97311 Median :0.954295 Median :0.9765
#> Mean :0.9114 Mean :0.93317 Mean :0.911515 Mean :0.9494
#> 3rd Qu.:0.9897 3rd Qu.:0.99242 3rd Qu.:0.985824 3rd Qu.:0.9947
#> Max. :1.0000 Max. :0.99997 Max. :0.999996 Max. :1.0000
#> NA's :10 NA's :6 NA's :7 NA's :20
#> zg
#> Min. :0.1744
#> 1st Qu.:0.9262
#> Median :0.9707
#> Mean :0.9389
#> 3rd Qu.:0.9905
#> Max. :1.0000
#> NA's :15
sort(colSums(aW < 0.05, na.rm = TRUE)) # how often "clear outlier":
#> lu st zg ba se sz su si re la ef ad ri
#> 0 0 0 1 1 1 2 3 4 10 14 17 48
# lu st zg ba se sz su si re la ef ad ri
# 0 0 0 1 1 1 2 3 4 10 14 17 48
lattice::levelplot(aW, asp=1/2, main="Robustness weights",
xlab= "day", ylab= "site")
