I wish to begin a detailed analysis of this problem and welcome comments. I believe the published integral expression for the beta function continuation
$$\left(1-e^{2\pi i\alpha}\right)\left(1-e^{2\pi i\beta}\right)\beta(\alpha,\beta)=\int_P z^{\alpha-1}(1-z)^{\beta-1}dz$$
is misleading in that it is only valid if the correct determination of the integrand is used at the start of integration. However, as I am new to this concept, I could be wrong and if so, I hope someone reading this thread can help me understand it better. This analysis includes what I believe it one of the most difficult concepts in Complex Analysis for new students to master: computing the analytically-continuous argument change of a function around a contour or ##\Delta_C \text{arg}(f)##. So I hope those interested in the topic will find the following helpful.
The following three forms of the pochhammer contour are useful. The four horizontal paths along the real axis over the ##P_2## and ##W## contours are over different determinations of the function:
Consider the computation of ##\beta(\alpha,\beta)## with ##\alpha=3/2## and ##\beta=4/3## using the pochhammer integral
$$\int_P z^{1/2} (1-z)^{1/3} dz=\int_P fdz.$$
Now consider the complex function ##\sqrt[6]{z^3(1-z)^2}##. This function is a six-valued function and therefore the pochhammer contour can traverse the function in six different ways depending upon which determination of the function the integration begins on. In order to unambigously identify these six determinations, consider the algebraic function, ##w(z)## defined as
$$z^3(1-z)^2-w^6=0$$
Arrange the non-zero singular points, ##s_i##, of the function in increasing order of absolute value. Take one-half of the smallest member, ##s_1## and identify the reference point ##z_0=1/2 |s_1|##. Compute ##f(z_0,w)=0## and arrange the values in order of increasing real part then increasing imaginary part to obtain the list ##(w_1,w_2,\cdots, w_n)##. The determinations of this function can now be unambiguously identified by analytic continuation from these reference points. Now consider
$$I_n=\int_{W_n} f dz$$
where the ##n## reflects the particular determination of ##f## used at the start of integration, ##z=1/10## along the blue contour. The following table summarizes numerical results for all determinations.
$$
\begin{array}{|c|c|}
\hline
n & I_n\\
\hline
1 & -1.37669-0.794831 i \\
2 & 1.37669\, +0.794831 i \\
3 & -2.18\times 10^{-7}-1.58966 i \\
4 & -1.37669+0.794831 i \\
5 & 1.37669\, -0.794831 i \\
6 & 2.18 \times 10^{-7}+1.58966 i \\
\hline
\end{array}
$$
We now go back to the ##P_2## contour and begin the analysis of computing ##\Delta \text{arg}## of the function as we proceed around the contour beginning at the point ##z=1/10## and for that analysis, we write the function as:
$$f=|z|^{1/2}e^{i/2(\theta_p+2n\pi)}|1-z|^{1/3}e^{i/3(\theta_q+2k\pi)}, \quad n=0,1,\quad k=0,1,2$$
with ##\theta_p=\text{Arg}(z),\quad \theta_q=\text{Arg}(1-z)## and ##-\pi<\text{Arg}(u)\leq \pi##. In this particular case, we wish to begin the integration of the principle value of the integrand over the blue contour which in the standard order above is the determination identified by ##w_2##. As we proceed along the blue contour, the argument does not change. At the point ##x=9/10##, we begin around the circular arc with ##z=1+\rho e^{it}## as ##t## varies from ##\pi## to ##-\pi##. The argument of the first term does not change, however, the mapping of ##(1-(1+\rho e^{it})## as t goes around this circle is a circle contour around the origin going in a clockwise direction beginning at ##\rho##. When ##t=0## we have for the second term ##e^{i/3(-\pi)}## and for the path to be analytically continuous across the negative real axis, we must join the current branch with the branch ##e^{i/3(\theta_q+4\pi)}=e^{5\pi i/3}=e^{-\pi i/3}## and therefore we have ##e^{4\pi i/3}f## for the value of the integrand along the purple contour. Again, the argument does not change along the horizontal leg of the contour. At the point ##z=1/10##, we now consider the argument change of ##z^{1/2}##. We then have ##e^{i/2(\theta_p+2n\pi)}##. As we go around the circular leg in a counterclockwise direction, we reach the point ##z=-1/10## and have ##e^{i/2(\pi)}## and similar to the above, in order to maintain analytic continuity around the contour, we join the current branch with ##e^{i/2(\theta_p+2\pi)}=e^{i/2(-\pi+2\pi)}=e^{i/2(\theta_p+2\pi)}## and proceed to the orange contour at which point we now have ##e^{i/2(2\pi)}e^{i/3(4\pi)}f##. We do similar analyses along the entire contour and then consider the limit as the radius, ##\rho## around the two singular points approach zero. We are then left with the horizontal legs and have$$\begin{align*}
\lim_{\rho\to 0}\int_{W_2} fdz&=
\int_0^1 f-e^{4\pi i/3}\int_0^1 f+e^{2\pi i/2}e^{4\pi i/3}\int_0^1 f-e^{2\pi i/2}e^{2\pi i/3}e^{4\pi i/3}\int_0^1 f\\
&=\left(1-e^{4\pi i/3}\right)\int_0^1 f+e^{2\pi i/2}e^{4\pi i/3}\int_0^1 f-e^{2\pi i/2}\int_0^1 f\\
&=\left(1-e^{4\pi i/3}\right)\int_0^1 f+e^{2\pi i/2}\left(e^{4\pi i/3}-1\right) \int_0^1 f\\
&=\left(1-e^{2\pi i/2}\right)\left(1-e^{4\pi i/3}\right)\int_0^1 fdz
\end{align*}
$$
but ##e^{4\pi i(1/3)}=e^{4\pi i(4/3-1)}=e^{4\pi i \beta}## and ##e^{2\pi i/2}=e^{2\pi i(3/2-1)}=e^{2\pi i\alpha}## so that
$$\left(1-e^{2\pi i\alpha}\right)\left(1-e^{4\pi i\beta}\right)\int_0^1 z^{1/2}(1-z)^{1/3}dx=\lim_{\rho\to 0}\int_{W_2} z^{1/2}(1-z)^{1/3}dz$$
or:
$$\left(1-e^{2\pi i\alpha}\right)\left(1-e^{4\pi i\beta}\right)\beta(\alpha,\beta)=\int_{W_2} z^{1/2}(1-z)^{1/3}dz$$
which is consistent with numerical results since:
$$\frac{I_2}{\left(1-e^{2\pi i\alpha}\right)\left(1-e^{4\pi i\beta}\right)}\approx 0.458896 - 7.16\times10^{-9}i\approx \beta(3/2,4/3)$$
Note this is not in the form of the standard integral expression of the beta function continuation.