Python PyMankendall与Pandas Groupby解包结果
我有一个系列,它是使用lambda应用pymannkendall函数的结果。然而,问题是,我需要以某种方式解压该系列的元素,而我不能以其当前形式这样做。数据帧“ncd”的示例如下所示:Python PyMankendall与Pandas Groupby解包结果,pandas,lambda,pandas-groupby,Pandas,Lambda,Pandas Groupby,我有一个系列,它是使用lambda应用pymannkendall函数的结果。然而,问题是,我需要以某种方式解压该系列的元素,而我不能以其当前形式这样做。数据帧“ncd”的示例如下所示: plant_name year season day count mean std min 10% 50% 90% max 0 ARIZONA I 2005 1 1 72.0 9.37 1.36 7.43 7.72 9.26 10.8
plant_name year season day count mean std min 10% 50% 90% max
0 ARIZONA I 2005 1 1 72.0 9.37 1.36 7.43 7.72 9.26 10.86 13.18
1 ARIZONA I 2005 1 2 72.0 9.19 1.00 7.37 7.93 9.22 10.24 12.24
2 ARIZONA I 2005 1 3 72.0 9.28 1.60 6.31 7.10 9.41 11.32 12.63
3 ARIZONA I 2005 1 4 72.0 8.76 1.99 5.59 5.74 9.36 10.80 11.06
4 ARIZONA I 2005 1 5 72.0 8.40 2.42 4.18 5.09 9.33 10.91 13.89
ncd = ncData.groupby(['plant_name','year','season','day'])['wind_speed_ms'].describe(percentiles = [0.9, 0.1]).round(2).reset_index()
nc10= ncd.groupby(['plant_name']).apply(lambda x: mk.original_test(x['10%']))
In [68]: nc10
Out[68]:
plant_name
ARIZONA I (no trend, False, 0.416671903456554, 0.8122086...
CAETITE I (increasing, True, 0.00022557963311364837, 3.6...
CAETITE II (increasing, True, 0.001237301922455858, 3.230...
CAETITE III (increasing, True, 0.00029610116053091495, 3.6...
CALANGO I (no trend, False, 0.06614633155767802, -1.8374...
CALANGO II (decreasing, True, 0.004494175262409472, -2.84...
CALANGO III (no trend, False, 0.08652731348848075, -1.7140...
CALANGO IV (decreasing, True, 0.0203025663568539, -2.3207...
CALANGO V (no trend, False, 0.05860321719835615, -1.8911...
CALANGO VI (decreasing, True, 0.006518506037077154, -2.72...
CANOAS (decreasing, True, 0.04470174299349505, -2.007...
LAGOA I (no trend, False, 0.1335175311613055, -1.50037...
LAGOA II (no trend, False, 0.1395114855149504, -1.47761...
MEL II (no trend, False, 0.09929968197907635, -1.6482...
RIO DO FOGO (no trend, False, 0.5553621888263085, 0.589744...
SANTANA I (decreasing, True, 0.006538153896656684, -2.71...
SANTANA II (decreasing, True, 0.002009326564589742, -3.08...
dtype: object
我创建系列(17,)的代码如下所示:
plant_name year season day count mean std min 10% 50% 90% max
0 ARIZONA I 2005 1 1 72.0 9.37 1.36 7.43 7.72 9.26 10.86 13.18
1 ARIZONA I 2005 1 2 72.0 9.19 1.00 7.37 7.93 9.22 10.24 12.24
2 ARIZONA I 2005 1 3 72.0 9.28 1.60 6.31 7.10 9.41 11.32 12.63
3 ARIZONA I 2005 1 4 72.0 8.76 1.99 5.59 5.74 9.36 10.80 11.06
4 ARIZONA I 2005 1 5 72.0 8.40 2.42 4.18 5.09 9.33 10.91 13.89
ncd = ncData.groupby(['plant_name','year','season','day'])['wind_speed_ms'].describe(percentiles = [0.9, 0.1]).round(2).reset_index()
nc10= ncd.groupby(['plant_name']).apply(lambda x: mk.original_test(x['10%']))
In [68]: nc10
Out[68]:
plant_name
ARIZONA I (no trend, False, 0.416671903456554, 0.8122086...
CAETITE I (increasing, True, 0.00022557963311364837, 3.6...
CAETITE II (increasing, True, 0.001237301922455858, 3.230...
CAETITE III (increasing, True, 0.00029610116053091495, 3.6...
CALANGO I (no trend, False, 0.06614633155767802, -1.8374...
CALANGO II (decreasing, True, 0.004494175262409472, -2.84...
CALANGO III (no trend, False, 0.08652731348848075, -1.7140...
CALANGO IV (decreasing, True, 0.0203025663568539, -2.3207...
CALANGO V (no trend, False, 0.05860321719835615, -1.8911...
CALANGO VI (decreasing, True, 0.006518506037077154, -2.72...
CANOAS (decreasing, True, 0.04470174299349505, -2.007...
LAGOA I (no trend, False, 0.1335175311613055, -1.50037...
LAGOA II (no trend, False, 0.1395114855149504, -1.47761...
MEL II (no trend, False, 0.09929968197907635, -1.6482...
RIO DO FOGO (no trend, False, 0.5553621888263085, 0.589744...
SANTANA I (decreasing, True, 0.006538153896656684, -2.71...
SANTANA II (decreasing, True, 0.002009326564589742, -3.08...
dtype: object
“nc10”系列的外观如下所示:
plant_name year season day count mean std min 10% 50% 90% max
0 ARIZONA I 2005 1 1 72.0 9.37 1.36 7.43 7.72 9.26 10.86 13.18
1 ARIZONA I 2005 1 2 72.0 9.19 1.00 7.37 7.93 9.22 10.24 12.24
2 ARIZONA I 2005 1 3 72.0 9.28 1.60 6.31 7.10 9.41 11.32 12.63
3 ARIZONA I 2005 1 4 72.0 8.76 1.99 5.59 5.74 9.36 10.80 11.06
4 ARIZONA I 2005 1 5 72.0 8.40 2.42 4.18 5.09 9.33 10.91 13.89
ncd = ncData.groupby(['plant_name','year','season','day'])['wind_speed_ms'].describe(percentiles = [0.9, 0.1]).round(2).reset_index()
nc10= ncd.groupby(['plant_name']).apply(lambda x: mk.original_test(x['10%']))
In [68]: nc10
Out[68]:
plant_name
ARIZONA I (no trend, False, 0.416671903456554, 0.8122086...
CAETITE I (increasing, True, 0.00022557963311364837, 3.6...
CAETITE II (increasing, True, 0.001237301922455858, 3.230...
CAETITE III (increasing, True, 0.00029610116053091495, 3.6...
CALANGO I (no trend, False, 0.06614633155767802, -1.8374...
CALANGO II (decreasing, True, 0.004494175262409472, -2.84...
CALANGO III (no trend, False, 0.08652731348848075, -1.7140...
CALANGO IV (decreasing, True, 0.0203025663568539, -2.3207...
CALANGO V (no trend, False, 0.05860321719835615, -1.8911...
CALANGO VI (decreasing, True, 0.006518506037077154, -2.72...
CANOAS (decreasing, True, 0.04470174299349505, -2.007...
LAGOA I (no trend, False, 0.1335175311613055, -1.50037...
LAGOA II (no trend, False, 0.1395114855149504, -1.47761...
MEL II (no trend, False, 0.09929968197907635, -1.6482...
RIO DO FOGO (no trend, False, 0.5553621888263085, 0.589744...
SANTANA I (decreasing, True, 0.006538153896656684, -2.71...
SANTANA II (decreasing, True, 0.002009326564589742, -3.08...
dtype: object
我需要以某种方式解包mk.original_测试(x['%10'])的结果。解包的元素是这种形式的,我不能用我所知道的(17,1)系列“nc10”来解包。也许,我可以用另一种方式将mk.original_测试(数据)函数与groupby('plant_name')一起使用除了lambda函数之外…通过我需要的plant_名称获取解包结果
trend, h, p, z, Tau, s, var_s, slope, intercept = mk.original_test(data)