使用drm拟合Weibull 1.4模型时upperl参数报错的解决咨询

问题描述 #

尝试使用Weibull 1.4模型构建剂量反应模型分析百分比数据，但始终触发上限参数（upperl）错误。调整上限参数值后问题依旧，求错误原因及修复方法。可复现代码如下：

Percent.adj <- c(".002",".005", "003",".5", ".87", "1",".001",".025", "005",".5", ".90", "1")
TempC <- c("0", "2","4","6","8","10","0", "2","4","6","8","10")
Genotype<-c("A","A", "A","A", "A","A", "B","B","B","B","B","B")

df<-data.frame(Percent.adj, TempC, Genotype)


drm(Percent.adj ~ TempC, Genotype, data=df,
    fct=W1.4(),
    upperl=c(NA,2,NA)) -> Leaf.drm

错误原因分析 #

1. 数据类型不匹配：Percent.adj（响应变量）和TempC（剂量变量）均为字符型（带引号的字符串），而drm函数要求这两类变量必须是数值型。字符型数据会导致模型在拟合参数时出现数值计算错误，进而触发参数约束相关的报错。
2. upperl参数维度不匹配：Weibull 1.4（W1.4()）模型包含4个参数，顺序为：下限（lower）、斜率（slope）、ED50、上限（upper）。你传入的upperl=c(NA,2,NA)仅包含3个元素，与模型参数数量不匹配，导致参数约束逻辑失效。

修复方法与修正代码 #

步骤1：转换数据类型为数值型 #

将字符型的响应变量和剂量变量转换为数值型，确保模型能进行正常数值计算。

步骤2：匹配upperl参数维度 #

根据W1.4()模型的4个参数顺序，传入长度为4的upperl向量，NA表示不对该参数设置上限约束。

修正后的代码如下：

# 转换响应变量和剂量变量为数值型
Percent.adj <- as.numeric(c(".002",".005", "003",".5", ".87", "1",".001",".025", "005",".5", ".90", "1"))
TempC <- as.numeric(c("0", "2","4","6","8","10","0", "2","4","6","8","10"))
Genotype <- c("A","A", "A","A", "A","A", "B","B","B","B","B","B")

df <- data.frame(Percent.adj, TempC, Genotype)

# 适配W1.4模型的4个参数，设置upperl约束（示例：仅限制斜率上限为2，其余参数无约束）
drm(Percent.adj ~ TempC, Genotype, data=df,
    fct=W1.4(),
    upperl=c(NA, 2, NA, NA)) -> Leaf.drm

# 查看模型拟合结果
summary(Leaf.drm)

补充说明 #

• 数据类型转换是解决问题的核心，字符型变量无法参与剂量反应模型的非线性拟合计算，这是触发参数错误的根本原因。
• 若需要对其他参数设置上限约束，只需在upperl向量对应位置替换NA为目标数值即可，例如要限制上限参数不超过1（因为你的响应变量是百分比，最大为1），可设置upperl=c(NA, 2, NA, 1)。

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