NOvA is a long-baseline neutrino oscillation experiment that measures oscillations in charged-current <a:math xmlns:a="http://www.w3.org/1998/Math/MathML" display="inline"><a:mrow><a:msub><a:mrow><a:mi>ν</a:mi></a:mrow><a:mrow><a:mi>μ</a:mi></a:mrow></a:msub><a:mo stretchy="false">→</a:mo><a:msub><a:mrow><a:mi>ν</a:mi></a:mrow><a:mrow><a:mi>μ</a:mi></a:mrow></a:msub></a:mrow></a:math> (disappearance) and <d:math xmlns:d="http://www.w3.org/1998/Math/MathML" display="inline"><d:msub><d:mi>ν</d:mi><d:mi>μ</d:mi></d:msub><d:mo stretchy="false">→</d:mo><d:msub><d:mi>ν</d:mi><d:mi>e</d:mi></d:msub></d:math> (appearance) channels, and their antineutrino counterparts, using neutrinos of energies around 2 GeV over a distance of 810 km. In this work we reanalyze the dataset first examined in our previous paper [] using an alternative statistical approach based on Bayesian Markov chain Monte Carlo. We measure oscillation parameters consistent with the previous results. We also extend our inferences to include the first NOvA measurements of the reactor mixing angle <g:math xmlns:g="http://www.w3.org/1998/Math/MathML" display="inline"><g:mrow><g:msub><g:mrow><g:mi>θ</g:mi></g:mrow><g:mrow><g:mn>13</g:mn></g:mrow></g:msub></g:mrow></g:math>, where we find <i:math xmlns:i="http://www.w3.org/1998/Math/MathML" display="inline"><i:mrow><i:mn>0.071</i:mn><i:mo>≤</i:mo><i:msup><i:mrow><i:mi>sin</i:mi></i:mrow><i:mrow><i:mn>2</i:mn></i:mrow></i:msup><i:mn>2</i:mn><i:msub><i:mrow><i:mi>θ</i:mi></i:mrow><i:mrow><i:mn>13</i:mn></i:mrow></i:msub><i:mo>≤</i:mo><i:mn>0.107</i:mn></i:mrow></i:math>, and the Jarlskog invariant, where we observe no significant preference for the <k:math xmlns:k="http://www.w3.org/1998/Math/MathML" display="inline"><k:mi>C</k:mi><k:mi>P</k:mi></k:math>-conserving value <m:math xmlns:m="http://www.w3.org/1998/Math/MathML" display="inline"><m:mi>J</m:mi><m:mo>=</m:mo><m:mn>0</m:mn></m:math> over values favoring <o:math xmlns:o="http://www.w3.org/1998/Math/MathML" display="inline"><o:mi>C</o:mi><o:mi>P</o:mi></o:math> violation. We use these results to examine the effects of constraints from short-baseline measurements of <q:math xmlns:q="http://www.w3.org/1998/Math/MathML" display="inline"><q:msub><q:mi>θ</q:mi><q:mn>13</q:mn></q:msub></q:math> using antineutrinos from nuclear reactors when making NOvA measurements of <s:math xmlns:s="http://www.w3.org/1998/Math/MathML" display="inline"><s:msub><s:mi>θ</s:mi><s:mn>23</s:mn></s:msub></s:math>. Our long-baseline measurement of <u:math xmlns:u="http://www.w3.org/1998/Math/MathML" display="inline"><u:msub><u:mi>θ</u:mi><u:mn>13</u:mn></u:msub></u:math> is shown to be consistent with the reactor measurements, supporting the general applicability and robustness of the Pontecorvo-Maki-Nakagawa-Sakata framework for neutrino oscillations. Published by the American Physical Society 2024