使用NumPy的 -
a = df.values
out = a[np.abs(a[:,1:] - dfB.values.ravel()).argmin(0),0]
基本上,我們減去dfB
從dfA
的每一行開始,由於我們正在使用NumPy數組(因爲我們已經提取了那些wi .values
),這些在broadcasted manner
中被有效地扣除。然後,我們找到絕對值並沿.argmin(axis=0)
(簡而言之.argmin(0)
)沿着每列查找最小值arg。
如果您也在使用NaNs
,請使用np.nanargmin
來忽略這些。
步步樣品運行,以使事情更容易理解 -
# Extract array from dfA
In [9]: a = dfA.values
# Slice a from col-1 onwards and perform broadcasted differencing with dfB values
In [10]: a[:,1:] - dfB.values.ravel()
Out[10]:
array([[-0.92 , -0.802, -0.642, -0.523],
[-0.02 , -0.042, -0.022, -0.013],
[-0.02 , -0.042, -0.012, -0.003],
[ 0. , -0.022, 0.008, 0.017],
[ 0.02 , -0.002, 0.028, 0.037],
[ 0.04 , 0.018, 0.048, 0.057]])
# Get absolute values
In [11]: np.abs(a[:,1:] - dfB.values.ravel())
Out[11]:
array([[ 0.92 , 0.802, 0.642, 0.523],
[ 0.02 , 0.042, 0.022, 0.013],
[ 0.02 , 0.042, 0.012, 0.003],
[ 0. , 0.022, 0.008, 0.017],
[ 0.02 , 0.002, 0.028, 0.037],
[ 0.04 , 0.018, 0.048, 0.057]])
# Look for argmin along each col
In [14]: idx = np.abs(a[:,1:] - dfB.values.ravel()).argmin(axis=0)
In [17]: idx
Out[17]: array([3, 4, 3, 2])
# First col from a
In [15]: a[:,0]
Out[15]: array([ 0., 1., 2., 3., 4., 5.])
# Index into first col with those indices to select the desired output values
In [16]: a[idx,0]
Out[16]: array([ 3., 4., 3., 2.])
嗯,也許主要是時間上的差異,因爲'NaN'處理......你怎麼看? – jezrael
@jezrael我想它只是它的陣列更快的方式,我猜? Nan處理可以用'np.nanargmin()'來解決。 – Divakar
@jezrael如果你碰巧知道這些減法,大熊貓是否使用'broadcast'? – Divakar