我有一组数据,我喜欢适合1D-Sersic配置文件(已经在astropy中定义的函数)如果我只考虑没有任何错误条的数据点,我可以很容易地拟合数据:
import numpy as np
from astropy.modeling import models, fitting
from astropy.modeling.models import Sersic1D
##
data_LARS02 =np.loadtxt('/Users/a/data.asc')
##
r_l02 = data[:,0]
flux_l02 = data[:,1]
###########################################Fitting
sers_l02_init = Sersic1D(amplitude=5e40, r_eff=1, n=2)
fit_sers_l02 = fitting.LevMarLSQFitter()
Sers_l02 = fit_sers_l02(sers_l02_init,r_l02,flux_l02)
#####
代码有效,一切都很好。但是,我注意到由于忽略误差条而导致拟合问题。因此,我决定包括错误栏。但我不知道如何实现它。你有什么建议吗?