I'm reading the paper on Adaptive LASSO estimator (Zou, 2006). In one of the presented numerical simulation examples (Model 0 (Inconsistent lasso path), page 6 (1423)) they claim the following:
To show this numerically [read. inconsistent variable selection], we simulated 100 datasets from the foregoing model for three different combinations of sample size (n) and error variance (σ2). On each dataset, we computed the entire solution path of the lasso, then estimated the probability of the lasso solution path containing the true model.
However, I don't understand what do they mean with "then estimated the probability of the lasso solution path containing the true model". Or, in other words, how would you estimate such probability?
I can guess at least a few different ways to measure how far was the estimated values from the true values, for example:
Repeat one of the following steps for $N$ times:
Given the full LASSO solution path, mark
true
if for any hyperparameter $\lambda$ the estimated solution was $\hat\beta = (\hat\beta_1, \hat\beta_2, \hat\beta_3 , 0)$. That is, for some $\lambda$, the truly insignificant variable got estimated as $\hat\beta_4 = 0$, but the other three as some shrinked value (not necessarily $\hat\beta_j = \beta_j = 5.6$ holds). But this I understand as testing for variable selection, not true model selectionGiven the full LASSO solution path, mark
true
if for any hyperparameter $\lambda$, the estimated solution was $\hat\beta = (\beta_1, \beta_2, \beta_3, 0)$ with true model parameters. This one sounds the most logical for me, but raises too many questions: firstly, no mention of threshold for the estimated parameters ($\hat\beta_j = 5.61 = \text{ or } \neq 5.6 = \beta_j)$, when are they said to be equal (the true model)? Secondly, unless $\lambda$ is close to 0, the estimates should be significantly shrinked, when compared with original true model. Lastly, LASSO surely can't achieve the true model in 50% of the cases, as they show in the paper? At least, not with my simulations..The percentage of the solution path, where only those variables that are present in the true data generating process are selected. But this somewhat ties with (1.) question, and I doubt this to be the case since it would highly depend on $\lambda$ values.
My question: what do the authors actually mean with their presented results? I'm interested in replicating their simulation results.
EDIT: To short this down, I think I am confused by the wording "path containing the true model". My confusion comes from imagining the LASSO path, say like in this random example; their wording seems to suggest that a certain point in the graph would be the true model, which would seem to involve the estimated values, not only which are "nonzero".
On the other hand, it seems they deal with parameter estimation in other examples, so it would make sense this particular example was about variable selection.