I would like to approximate # in terms of $\pi(x)$.
The nearest I have found is by the answer of @draks:
Interpolating the primorial $p_{n}\#$
It seems to use asymptotic integration:
$\begin{eqnarray} \sum_{k=1}^n \log p_n &=& \int_2^n \log k\; d\pi(k)\\ &=& \log(k)\pi(k)\biggr|_{2}^{n}+\int_{2}^{n}\frac1k \pi(k)dk. \end{eqnarray}$
using $\pi(n)\sim \frac{n}{\log n}$
$\log p_n\# \sim \log(k)\frac{k}{\log k}\biggr|_{2}^{n}+\int_{2}^{n}\frac1k \frac{k}{\log k}dk = (x-1)+\text{Li}(x) \;. \tag{$*$}$
then exponentiate*.
1) No error bound is given and I assume (maybe wrongly) it wound be whatever error there is in the estimation of $\pi(n)$?
2) No reference is given so I am not sure if this method is widely accepted? Would it work? What are the bounds?
3) Is there something better for approximating primorial # in terms of π(x)?.
EDIT: Regarding question 2) There is at least one relevant reference (page 4-5) to using asymptotic integration for finding θ(x) in terms of π(x):
http://math.tufts.edu/faculty/rlemkeoliver/teaching/250/02-RiemannStieltjes.pdf
Hopefully this might help somebody looking at the stackexchange link given above.
Please also note correction to equation above given in @mixedmath answer.