我正在尝试使用以下代码,但遇到了指定的错误代码。 谁能告诉我为什么会这样?
代码:
import numpy as np
import xarray as xr
from scipy.stats import gamma
import matplotlib.pyplot as plt
data1 = xr.open_dataset('data1.nc')
avg1 = data1.sel(lat=slice(5,38),lon=slice(60,100))
data1_fit = gamma.fit(avg1['variable1'], floc=0)
错误代码:
FitDataError:
data
中的值无效。使用“gamma”的最大似然估计要求 0 < (x - loc)/scale < inf for each x in data
.
我该如何解决这个问题?
我试图将 netcdf 数据的变量拟合到伽玛分布中。但这就是我陷入困境的地方。