我无法在 R 包 DRC 中绘制多个模型 drm。
我的数据包含三个变量;剂量;除草剂 (Trt) 和时机
当我绘制 Dose 和 Trt 时,剂量响应曲线工作正常,但是当我添加 Timing 时,我收到此错误: pmodelsList [[i]] 中的错误:下标超出范围这是我的代码:
library(pacman)
library(dplyr)
p_load(MASS,
drc,
ggplot2,
reshape2)
t1= read.csv('TimingBB.csv',sep = ",",header=TRUE,strip.white=TRUE)
t2 <- melt(t1,id.vars=c('Dose', 'Rep', 'Timing'), variable.name = "Trt", value.name = 'Germination')
ggplot(t2, aes(x=Dose, y=Germination, color=Timing))+
geom_point()+
facet_wrap(~Trt,nrow = 3)
modl_Full <-drm(Germination~Dose,Trt,Timing, data = t2, fct = LL.3()) # Full model