Berechnung des Bestimmungskoeffizienten (R2) und des quadratischen Mittelwertfehlers (RMSE) für die nichtlineare Kurvenanpassung in Python. Der folgende Code gilt bis zur Kurvenanpassung. Wie berechnet man dann R2 und RMSE?

import numpy as np
import matplotlib.pyplot as plt
from scipy.optimize import curve_fit

def func(x, a, b, c):
    return a * np.exp(-b * x) + c

x = np.linspace(0,4,50)
y = func(x, 2.5, 1.3, 0.5)
yn = y + 0.2*np.random.normal(size=len(x))

popt, pcov = curve_fit(func, x, yn)

plt.figure()
plt.plot(x, yn, 'ko', label="Original Noised Data")
plt.plot(x, func(x, *popt), 'r-', label="Fitted Curve")
plt.legend()
plt.show()