Comparing the Sage code base to the UCLA-based package, we see a difference in the selection of the 'soft threshold' (beta). The former picked a value of 4.0 and the latter 4.5. Both packages save the details of the optimization process:
This is the output from the UCLA algorithm:
Power |
SFT.R.sq |
slope |
truncated.R.sq |
mean.k. |
median.k. |
max.k. |
|---|---|---|---|---|---|---|
1 |
0.736501290349614 |
1.13578366310746 |
0.994431751751329 |
699.572313034948 |
699.648966320864 |
1116.02918233716 |
1.5 |
0.251519543919548 |
0.297533023504088 |
0.97937139813633 |
439.899479746776 |
423.612496822198 |
840.47949387766 |
2 |
0.119730110387487 |
-0.163692811558184 |
0.91713459200463 |
295.588422458769 |
268.94730085273 |
667.945223302198 |
2.5 |
0.622896623436384 |
-0.491086993626836 |
0.900740857562872 |
208.537727229792 |
177.79413914117 |
549.155319758776 |
3 |
0.788460807596207 |
-0.68266797986082 |
0.903728452209731 |
152.800323593864 |
120.243286725551 |
463.29082261878 |
3.5 |
0.863229344305799 |
-0.85737992139182 |
0.889320106067623 |
115.448359595823 |
83.7476291500883 |
399.156919651325 |
4 |
0.894834165325891 |
-0.951860705585135 |
0.894194460202712 |
89.495687627513 |
59.4342800786908 |
350.910943247122 |
4.5 |
0.906475253699557 |
-1.04323304332767 |
0.88699005785024 |
70.9224191593446 |
42.2148824089728 |
312.992004152579 |
5 |
0.917140610429364 |
-1.08867663568657 |
0.893779632076047 |
57.2968574266249 |
30.4851790719122 |
282.652362917002 |
5.5 |
0.902963257559623 |
-1.12874378525066 |
0.876283735962727 |
47.0867813585159 |
22.5831735990022 |
257.887940247821 |
6 |
0.893657667984861 |
-1.138641996747 |
0.872082886273499 |
39.2934112572874 |
16.7925249162364 |
237.362271809525 |
6.5 |
0.877096891199445 |
-1.16563222664169 |
0.860765166816465 |
33.2467401837027 |
12.5431252833217 |
220.116544465302 |
7 |
0.867836927064159 |
-1.17261710771242 |
0.86332246600206 |
28.4861111749774 |
9.75194877632007 |
205.448215921732 |
7.5 |
0.868548947706417 |
-1.15942033025233 |
0.888456324814093 |
24.6879610646154 |
7.43724266918352 |
192.83304720061 |
8 |
0.85993341256465 |
-1.14979252940141 |
0.889730024616774 |
21.6207148555113 |
5.72770680062926 |
181.87355758458 |
8.5 |
0.839526774373017 |
-1.1454863855111 |
0.889555649590024 |
19.1158525842939 |
4.42968785517107 |
172.264089318976 |
9 |
0.828450334938823 |
-1.12976930526696 |
0.896114026911798 |
17.0489014282755 |
3.48574207332117 |
163.766610421789 |
9.5 |
0.824533609394867 |
-1.10327784158547 |
0.912422650093755 |
15.3266812705319 |
2.74605181254734 |
156.193631947154 |
10 |
0.805943284944799 |
-1.08884406876998 |
0.914807175372563 |
13.8785832450612 |
2.18328519701562 |
149.395941841956 |
10.5 |
0.800039898486579 |
-1.0763656575682 |
0.920931944553393 |
12.6505039786172 |
1.74424089765862 |
143.266795277355 |
11 |
0.788808310243576 |
-1.05656008289991 |
0.933254569845571 |
11.6005616530767 |
1.40321696471758 |
137.745622540377 |
11.5 |
0.779205608927366 |
-1.03924741920651 |
0.935482608107987 |
10.696027985041 |
1.14073738236037 |
132.707798864861 |
12 |
0.779321654012786 |
-1.02460870779617 |
0.947336705575843 |
9.91110282257499 |
0.925552317403238 |
128.118730788473 |
Add Comment