吸烟行为学术事件图（学术图表）

def solve():
    import pandas as pd
    import matplotlib.pyplot as plt
    import numpy as np
    
    url = 'https://liangdaima.com/static/data/wooldridge/smoke.csv'
    df = pd.read_csv(url)
    
    plt.rcParams['font.size'] = 12
    df['income_q'] = pd.cut(df['income'].rank(method='first'), bins=5, labels=['Q1', 'Q2', 'Q3', 'Q4', 'Q5'])
    grouped = df.groupby('income_q', observed=True)['cigs'].agg(['mean', 'sem']).reset_index()
    plt.figure(figsize=(10,6))
    plt.bar(grouped['income_q'], grouped['mean'], yerr=grouped['sem']*1.96, color='white', edgecolor='black', capsize=5)
    plt.title('收入与吸烟量', fontsize=14)
    plt.xlabel('收入分位数', fontsize=12)
    plt.ylabel('每日吸烟量', fontsize=12)
    plt.grid(axis='y', alpha=0.3)
    plt.tight_layout()
    plt.show()
    return "可视化代码执行完成"