0

我无法在 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
4

0 回答 0