Question

所以这不是我需要使用的数据集，但这是我正在为一个癌症研究项目使用的庞大数据集（约180万个数据点）的模板，因此我确定是否可以将其用于使用较小的一个，然后我可以将其改编为较大的一个！因此，举例来说，我有以下数据集：

import numpy as np
import pandas as pd
df = pd.DataFrame({
   'cond': ['A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'B', 'B','B', 'B', 'B', 'B', 'B','B','B'],
   'Array':  ['S', 'S', 'TT', 'TT','S', 'S', 'TT', 'TT','S', 'S', 'TT', 'TT','S', 'S', 'TT', 'TT','SS','TT'],
   'X':  [1, 2, 3, 1, 2 , 3, 4, 7.3, 5.1, 3.2, 1.4, 5.5, 9.9, 3.2, 1.1, 3.3, 1.2, 5.4],
   'Y':  [3.1, 2.2, 2.1, 1.2,  2.4, 1.2, 1.5, 1.33, 1.5, 1.6, 1.4, 1.3, 0.9, 0.78, 1.2, 4.0, 5.0, 6.0],
   'Marker':  [2.0, 1.2, 1.2, 2.01, 2.55, 2.05, 1.66, 3.2, 3.21, 3.04, 8.01, 9.1, 7.06, 8.1, 7.9, 5.12, 5.23, 5.15],
   'Area': [3.0, 2.0, 2.88, 1.33,  2.44, 1.25, 1.53, 1.0, 0.156, 2.0, 2.4, 6.3, 6.9, 9.78, 10.2, 15.0, 16.0, 19.0]
})
print(df)

这将产生如下所示的输出：

   cond Array    X     Y  Marker    Area
0     A     S  1.0  3.10    2.00   3.000
1     A     S  2.0  2.20    1.20   2.000
2     A    TT  3.0  2.10    1.20   2.880
3     A    TT  1.0  1.20    2.01   1.330
4     A     S  2.0  2.40    2.55   2.440
5     A     S  3.0  1.20    2.05   1.250
6     A    TT  4.0  1.50    1.66   1.530
7     A    TT  7.3  1.33    3.20   1.000
8     A     S  5.1  1.50    3.21   0.156
9     B     S  3.2  1.60    3.04   2.000
10    B    TT  1.4  1.40    8.01   2.400
11    B    TT  5.5  1.30    9.10   6.300
12    B     S  9.9  0.90    7.06   6.900
13    B     S  3.2  0.78    8.10   9.780
14    B    TT  1.1  1.20    7.90  10.200
15    B    TT  3.3  4.00    5.12  15.000
16    B    SS  1.2  5.00    5.23  16.000
17    B    TT  5.4  6.00    5.15  19.000

好吧，现在我需要做的是基于两个标签“ cond”和“ Array”将它们拆分。我使用

g=df.groupby(['cond','Array'])['Marker']

这将其分成4个较小的集合，分为A-S，A-TT，B-S，B-TT对。现在，我可以使用自定义功能。这是功能的一部分，我将解释其工作原理：

def num_to_delete(p,alpha,N):
    if p==0.950:
        if 1-alpha==0.90:
            if N<=60:
                m=1
            if 60<N<80:
                m=round(N/20-2)
            if 80<=N:
                m=2
        if 1-alpha==0.95:
            if N<=80:
                m=1
            if 80<N<=100:
                m=round(N/20 -3)
            if 100<N:
                m=2
    return m

好的，所以它的工作方式是向我输入一个我选择的“ p”和“ alpha”（实际函数涵盖了p和alpha的更多情况）。馈入其中的N是我的较小数据集的元素数（在这种情况下，对于A-S为5，对于A-TT为4，依此类推）。因此，我试图发生的事情是，对于每个较小的数据集，吐出一些要删除的点（在此示例中，该函数将始终给我们1，但是我尝试使用该函数为应用到超大型数据集）。由于它给出的数字为1，因此我希望它删除该集合的最大1个数据点，然后告诉我剩下的最高点是什么。

因此，例如，对于A-S耦合，我有5个数据点：2.0、1.2、2.55、2.05和3.21。由于有5个数据点，因此我的函数告诉我删除其中1个，因此忽略3.21，然后告诉我剩下的最高数据点是2.55。我想对每个耦合进行此操作，但是在我的实际数据集中，我将拥有不同数量的元素，因此该函数将告诉我为每个耦合删除一个不同的数目。

我的最终目标是拥有一张最终桌，如下所示：

   cond Array   NumDeleted p95/a05  p95/a10       
0     A     S  1.0      2.55   2.55
1     A    TT  1.0      2.01   2.01
2     B     S  1.0      7.06   7.06
3     B    TT  1.0      8.01   8.01

对于较大的集合，最后两列中的值将有所不同，因为在大型数据集中，要删除的值的数量存在更多差异，因此其余值将有所不同。我最终将需要根据我从p95 / a05和p95 / a10获得的值来更改第二个数据集。

无论如何，对不起这么长的解释，但如果有人可以提供帮助，那就太好了！我希望这是一件相当简单的事情，因为我已经坚持了一个多星期。

Answer 1

编辑：更通用的解决方案

首先，创建一个closure来定义您的配置会有所帮助。假设您将来会拥有更多配置：

def create_num_to_delete(p, alpha):
    """Create a num_to_delete function given p and alpha."""
    def num_to_delete(N):
        if p == 0.950:
            if 1 - alpha == 0.90:
                if N <= 60:
                    m = 1
                if 60 < N < 80:
                    m = round(N/20 - 2)
                if 80 <= N:
                    m = 2
            if 1-alpha == 0.95:
                if N <= 80:
                    m = 1
                if 80 < N <= 100:
                    m = round(N/20 -3)
                if 100 < N:
                    m = 2
        return m

    return num_to_delete

然后可以使用此闭包定义配置字典：

configurations = {
    'p95/a05': create_num_to_delete(0.95, 0.05),
    'p95/a10': create_num_to_delete(0.95, 0.10),
}

然后，定义一个汇总数据的函数。此功能应取决于您的配置，以使其保持动态。

def summarize(x):
    # The syntax on the right-hand side is called list comprehension.
    # As you can probably guess, it's essentially a flattened for-loop that
    # produces a list. The syntax starting with "for" is your basic for loop
    # statement, and the syntax to the left of "for" is an expression that
    # that serves as the value of the resulting list for each iteration
    # of the loop.
    #
    # Here, we are looping through the "num_to_delete" functions we defined in
    # our `configurations` dictionary. And calling it in our group `x`.
    Ns = [num_to_delete(len(x)) for num_to_delete in configurations.values()]

    markers = x['Marker'].sort_values(ascending=False)

    highest_markers = []
    for N in Ns:
        if N == len(x):
            highest_markers.append(None)
        else:
            # Since we know that `markers` is already sorted in descending
            # order, all we need to get the highest remaining value is to get
            # the value in the *complete list* of values offset by the 
            # the number of values that need to be deleted (this is `N`).
            #
            # Since sequences are 0-indexed, simply indexing by `N` is enough.
            # For example, if `N` is 1, indexing by `N` would give us
            # the marker value *indexed by* 1, which is,
            # in a 0-sequenced index, simply the second value.
            highest_markers.append(markers.iloc[N])

    # Returning a list from an applied groupby function translates into
    # a DataFrame which the series index as the columns and the series values
    # as the row values. Index in this case is just the list of configuration
    # names we have in the `configurations` dictionary.
    return pd.Series(highest_markers, index=list(configurations.keys()))

最后，apply将该函数添加到您的数据集并重置索引。这会将cond和Array保留为列：

grouped = df.groupby(['cond', 'Array'])
grouped.apply(summarize).reset_index()

输出为：

    cond    Array   p95/a05 p95/a10
0   A   S   2.55    2.55
1   A   TT  2.01    2.01
2   B   S   7.06    7.06
3   B   SS  NaN NaN
4   B   TT  8.01    8.01

希望这会有所帮助。

使用split-apply-combine通过自定义函数删除一些值并结合剩余的内容

1 个答案: