The R(D (*)) anomalies observed in B → D (*) τ ν decays have attracted much attention in recent years. In this paper, we study the B → D (*) τ ν, b → c τ ν, B c → (J/ψ, η c)τ ν, B → X c τ ν, and B c → τ ν decays, all being mediated by the same quark-level b → cτ ν transition, in the Standard Model Effective Field Theory. The most relevant dimension-six operators for these processes are Q (3) lq , Q ledq , Q (1) lequ , and Q (3) lequ in the Warsaw basis. Evolution of the corresponding Wilson coefficients from the new physics scale = 1 TeV down to the characteristic scale μ b m b is performed at three-loop in QCD and one-loop in EW/QED. It is found that, after taking into account the constraint B(B c → τ ν) 10%, a single C (3) lq 3323 () or C (3) lequ 3332 () can still be used to resolve the R(D (*)) anomalies at 1σ , while a single C (1) lequ 3332 () is already ruled out by the measured R(D (*)) at more than 3σ. By minimizing the χ 2 (C i) function constructed based on the current data on R(D), R(D *), P τ (D *), R(J/ψ), and R(X c), we obtain eleven most trustworthy scenarios, each of which can provide a good explanation of the R(D (*)) anomalies at 1σ. To further discriminate these different scenarios, we predict thirty-one observables associated with the processes considered under each NP scenario. It is found that most of the scenarios can be differentiated from each other by using these observables and their correlations.