Julia实现OLS时出现UndefVarError: Cnsmr未定义问题求助

错误原因 #

你遇到的UndefVarError: Cnsmr not defined问题，根源在这行代码：

model = lm(eval(Meta.parse(formula_str)), df)

当你把字符串"Cnsmr ~ Mkt"解析并求值时，Julia会尝试寻找名为Cnsmr和Mkt的全局变量，而不是从DataFramedf中读取对应列。即使DataFrame包含这些列，当前作用域中没有同名变量，就会抛出未定义错误。

修复方案 #

方案1：用@formula宏动态构建公式 #

直接通过符号转换和宏插入，避免eval带来的作用域问题：

begin
    industries = ["Cnsmr", "Manuf", "HiTec", "Hlth", "Other"]
    for industry in industries
        if industry ∈ names(df)
            # 将字符串列名转为符号，动态注入公式
            formula = @formula($(Symbol(industry)) ~ Mkt)
            fit_model = lm(formula, df)

            intercept = coef(fit_model)[1]
            p_value = coef_pvalues(fit_model)[1]

            println("Industry: $industry")
            println("Alpha (Intercept): $intercept")
            println("Is statistically significant? ", p_value < 0.05)
            println()
        else
            println("Column $industry does not exist in the DataFrame.")
            println()
        end
    end
end

方案2：用term函数构建公式项（更灵活） #

通过StatsModels包的term函数直接创建公式组件，组合成合法公式：

using StatsModels

begin
    industries = ["Cnsmr", "Manuf", "HiTec", "Hlth", "Other"]
    mkt_term = term(:Mkt)
    for industry in industries
        if industry ∈ names(df)
            industry_term = term(Symbol(industry))
            formula = industry_term ~ mkt_term
            fit_model = lm(formula, df)

            intercept = coef(fit_model)[1]
            p_value = coef_pvalues(fit_model)[1]

            println("Industry: $industry")
            println("Alpha (Intercept): $intercept")
            println("Is statistically significant? ", p_value < 0.05)
            println()
        else
            println("Column $industry does not exist in the DataFrame.")
            println()
        end
    end
end

注意事项 #

• 尽量避免在Julia中随意使用eval，尤其是动态代码场景，不仅容易引发作用域问题，还存在安全风险。
• 确保已提前加载DataFrames、GLM（以及方案2所需的StatsModels）包。

内容的提问来源于stack exchange，提问作者Arusi Berd