Question

import patsy
from patsy import dmatrices, dmatrix, demo_data
dt=pd.DataFrame({'F1':['a','b','c','d','e','a'],'F2':['X','X','Y','Y','Z','Z']})

我知道我可以做到这一点

dmatrix("1+I(F1=='a')",dt)

但是我可以创建一个任意函数patsy吗？我试图在R中模仿公式语言中的相同级别的灵活性，但在python中似乎不能直接实现

def abd(x):
    1 if x in ['a','b','d'] else 0

dmatrix("1+abd(F1)",dt)

Answer 1

IIUC

def abd(x):
    return x.isin(['a','b','d'])
dmatrix("1+abd(F1)",dt)
Out[182]: 
DesignMatrix with shape (6, 2)
  Intercept  abd(F1)[T.True]
          1                1
          1                1
          1                0
          1                1
          1                0
          1                1
  Terms:
    'Intercept' (column 0)
    'abd(F1)' (column 1)

Answer 2

我测试更接近模仿已在R公式系统中建立的内容。以下是接受答案的简单表示。 Python的设计提供了这样的灵活性。 R肯定可以做同样的事情（自定义功能），但更容易被忽略。

import pandas as pd
from patsy import dmatrices, dmatrix, demo_data

dt=pd.DataFrame({'F1':['a','b','c','d','e','a'],'F2':['X','X','Y','Y','Z','Z']})
def xx(x,y):return(np.isin(x,list(y))*1)
dmatrix("1+xx(F1,['a','b'])",dt)

DesignMatrix with shape (6, 2)
  Intercept  xx(F1, ['a', 'b'])
      1                   1
      1                   1
      1                   0
      1                   0
      1                   0
      1                   1
  Terms:
    'Intercept' (column 0)
    "xx(F1, ['a', 'b'])" (column 1)

在Patsy中创建自定义函数

2 个答案: