Here's my situation:
In theory, $y = f(x1, x2,...,x5)$, and the exact equation is context dependent.
In practice, I don't have data for $x4$ and $x5$ as yet.
Based on the somewhat promising results of analysis with $y$, $x1$, $x2$, and $x3$, I'd like to make the case for further investment to obtain data for $x4$ and $x5$ and repeat the analysis. If that provides a better fit, it could lead to actionable insights.
Regression results so far on a dataset with 100+ samples look like this:
y ~ x1 has $R^2$ of 42% with p-value less than 5%
y ~ x1+x2 has $adj.R^2$ of 48% with p-value of x1 and overall p-value less than 5% but p-value of x2 is 20%
y ~ x1+x3 has $adj.R^2$ of 45% with p-value of x1 and overall p-value less than 5% but p-value of x3 is 13%
y ~ x1+x2+x3 has $adj.R^2$ of 51% with p-value of x1 and overall p-value less than 5% but p-value of x2 is 20% and x3 is 13%
Residual plots for all of the above show significant bias which I attribute to missing $x4$ and $x5$.
Given that the overall p-values are always less than 5%, is it okay to conclude that the analysis so far explains 51% of the variation in $y$ or should I be mindful of the not-so-great p-values of $x2$ and $x3$ and go with 42%?