R Under development (unstable) (2025-10-07 r88904 ucrt) -- "Unsuffered Consequences" Copyright (C) 2025 The R Foundation for Statistical Computing Platform: x86_64-w64-mingw32/x64 R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain conditions. Type 'license()' or 'licence()' for distribution details. R is a collaborative project with many contributors. Type 'contributors()' for more information and 'citation()' on how to cite R or R packages in publications. Type 'demo()' for some demos, 'help()' for on-line help, or 'help.start()' for an HTML browser interface to help. Type 'q()' to quit R. > # Clean environment and load libraries > rm(list = ls()) # Clear environment to avoid conflicts > > # Source files > library(MixStable) Attaching package: 'MixStable' The following object is masked from 'package:graphics': clip > > set.seed(123) > X <- stabledist::rstable(1000, alpha = 1, beta = 0, 1, 0, pm = 1) > t_grid <- seq(0.1, 1.0, length.out = 50) > phi_vals <- ecf_empirical(X, t_grid) > components <- ecf_components(phi_vals) > y_vals <- components$y_vals > arg_vals <- components$arg_vals > print(y_vals) [1] -2.60186343 -2.37619757 -2.17058999 -1.99652025 -1.86802332 -1.74316387 [7] -1.62302572 -1.55474034 -1.46902282 -1.37020140 -1.28758930 -1.22538101 [13] -1.17547139 -1.12039678 -1.03732094 -1.00782668 -0.95998026 -0.93396745 [19] -0.90078881 -0.86388373 -0.83197813 -0.79970075 -0.77462120 -0.73661452 [25] -0.73069182 -0.68953924 -0.65864355 -0.63038561 -0.57858679 -0.54381410 [31] -0.53143107 -0.48649043 -0.45653985 -0.43754534 -0.40597878 -0.36872446 [37] -0.34610784 -0.32131772 -0.28080330 -0.26485146 -0.22904473 -0.21901130 [43] -0.18163143 -0.15398507 -0.14642906 -0.12129889 -0.10324586 -0.09386809 [49] -0.06941164 -0.03625023 > print(arg_vals) [1] 0.008821858 0.016614452 0.016388042 0.018219473 0.022641592 [6] 0.020947783 0.009882925 0.013498361 0.007193231 0.007672032 [11] -0.001588190 -0.008560958 -0.004757965 -0.008925987 -0.009425465 [16] -0.012042830 -0.019281665 -0.026412909 -0.018409936 -0.023682784 [21] -0.034895020 -0.034354342 -0.037619408 -0.056207287 -0.053601083 [26] -0.057752202 -0.058569216 -0.058761564 -0.062235207 -0.059614042 [31] -0.063810928 -0.054521030 -0.054563449 -0.045558226 -0.049053799 [36] -0.044327881 -0.042916544 -0.049358825 -0.046836499 -0.047199823 [41] -0.048664456 -0.041098818 -0.038150690 -0.037900341 -0.039370411 [46] -0.051246207 -0.049730621 -0.056707077 -0.058219829 -0.080650152 > > proc.time() user system elapsed 2.04 0.35 2.40