pythonを用いてのlogistic regression analysis

pythonを用いてのlogistic regression analysis

ここに実現したいことを箇条書きで書いてください。

risk factorのOdds比　信頼区間　ｐ値を出したい

ここに質問の内容を詳しく書いてください。

１２個ほど危険因子を入れてます。すべてカテゴリ変数化
欠損値なし

Dtype : int64と一つだけfloat64

で

logistic = smf.glm(formula = "rec1 ~ gender + age + location + obstruction + type + patho. + p壁深達度 + ly + v + LN_12",
data = df,
family = sm.families.Binomial()).fit()

print(logistic.summary())

と入力すると

エラーメッセージ

Traceback (most recent call last):

File "/usr/local/lib/python3.8/dist-packages/IPython/core/interactiveshell.py", line 3326, in run_code
exec(code_obj, self.user_global_ns, self.user_ns)

File "<ipython-input-53-e20be902c944>", line 1, in <module>
logistic = smf.glm(formula = "rec1 ~ gender + age + location + obstruction + type + patho. + p壁深達度 + ly + v + LN_12",

File "/usr/local/lib/python3.8/dist-packages/statsmodels/base/model.py", line 169, in from_formula
tmp = handle_formula_data(data, None, formula, depth=eval_env,

File "/usr/local/lib/python3.8/dist-packages/statsmodels/formula/formulatools.py", line 63, in handle_formula_data
result = dmatrices(formula, Y, depth, return_type='dataframe',

File "/usr/local/lib/python3.8/dist-packages/patsy/highlevel.py", line 309, in dmatrices
(lhs, rhs) = _do_highlevel_design(formula_like, data, eval_env,

File "/usr/local/lib/python3.8/dist-packages/patsy/highlevel.py", line 164, in _do_highlevel_design
design_infos = _try_incr_builders(formula_like, data_iter_maker, eval_env,

File "/usr/local/lib/python3.8/dist-packages/patsy/highlevel.py", line 66, in _try_incr_builders
return design_matrix_builders([formula_like.lhs_termlist,

File "/usr/local/lib/python3.8/dist-packages/patsy/build.py", line 689, in design_matrix_builders
factor_states = _factors_memorize(all_factors, data_iter_maker, eval_env)

File "/usr/local/lib/python3.8/dist-packages/patsy/build.py", line 354, in _factors_memorize
which_pass = factor.memorize_passes_needed(state, eval_env)

File "/usr/local/lib/python3.8/dist-packages/patsy/eval.py", line 478, in memorize_passes_needed
subset_names = [name for name in ast_names(self.code)

File "/usr/local/lib/python3.8/dist-packages/patsy/eval.py", line 478, in <listcomp>
subset_names = [name for name in ast_names(self.code)

File "/usr/local/lib/python3.8/dist-packages/patsy/eval.py", line 109, in ast_names
for node in ast.walk(ast.parse(code)):

File "/usr/lib/python3.8/ast.py", line 47, in parse
return compile(source, filename, mode, flags,

File "<unknown>", line 1
patho.
^
SyntaxError: invalid syntax

該当のソースコード

python ### 試したこと 色々しました ### 補足情報（FW/ツールのバージョンなど） ここにより詳細な情報を記載してください。