Curvas dosis respuesta (daño)

library(readxl)

## Warning: package 'readxl' was built under R version 4.3.3

library(drc)

## Loading required package: MASS

## 
## 'drc' has been loaded.

## Please cite R and 'drc' if used for a publication,

## for references type 'citation()' and 'citation('drc')'.

## 
## Attaching package: 'drc'

## The following objects are masked from 'package:stats':
## 
##     gaussian, getInitial

library(magrittr)
library(dplyr)

## 
## Attaching package: 'dplyr'

## The following object is masked from 'package:MASS':
## 
##     select

## The following objects are masked from 'package:stats':
## 
##     filter, lag

## The following objects are masked from 'package:base':
## 
##     intersect, setdiff, setequal, union

library(gridExtra)

## 
## Attaching package: 'gridExtra'

## The following object is masked from 'package:dplyr':
## 
##     combine

library(broom)

## Warning: package 'broom' was built under R version 4.3.3

library(ggplot2)

 QUINCLORAC

QNC21DR4<- read_excel("D:/Paula/Ensayos invernadero Maestría/DR4/Muestreos DR4 21ddt.xlsx", 
    sheet = "QNC", range = "A1:X304")
QNC21DR4

  Prueba daño QNC 21 ddt

modelQNC21DR4 <- drm(data= subset(QNC21DR4, ID %in% c("FO", "UM", "J2", "VE", "SA", "CU")), 
           formula = Daño21 ~ Dosis, 
           curveid = ID,
           fct = LL.3())
tidy(modelQNC21DR4)

mselect(modelQNC21DR4, fctList = list(LL.4(), LL.3(), W1.4(), W2.4()))

## Warning in sqrt(diag(varMat)): Se han producido NaNs

##         logLik       IC  Lack of fit  Res var
## LL.3 -715.2933 1468.587 4.000085e-04 76.52114
## LL.3 -715.2933 1468.587 4.000085e-04 76.52114
## W2.4 -714.6783 1479.357 4.800676e-05 78.62027
## LL.4 -715.2182 1480.436 3.597917e-05 79.04173
## W1.4 -716.0773 1482.155 2.264842e-05 79.71684

plot(modelQNC21DR4,
     xlab= expression(Dosis (g~i.a~ha^{-1})),
     ylab= "Control (%)",
     lty=1,
     pch = c(1,15,0,5,18,6,7, 10, 3 , 2), cex= 1,
       legendPos = c(90, 105), cex.legend = 0.8,
     col = c("blue", "red", "darkgreen", "purple", "orange", "black"))

ED (modelQNC21DR4,respLev = c(90) , interval = "delta")

## 
## Estimated effective doses
## 
##            Estimate  Std. Error       Lower       Upper
## e:CU:90  2.8387e+04  6.1379e+04 -9.2710e+04  1.4948e+05
## e:FO:90  6.2908e+03  5.0928e+03 -3.7570e+03  1.6339e+04
## e:J2:90  5.5892e+06  1.9859e+07 -3.3591e+07  4.4769e+07
## e:SA:90  3.6261e+05  2.4251e+06 -4.4220e+06  5.1472e+06
## e:UM:90  1.9049e+03  1.5294e+03 -1.1126e+03  4.9224e+03
## e:VE:90  1.4663e+07  6.0132e+07 -1.0397e+08  1.3330e+08

compParm(modelQNC21DR4,"e", "-")

## 
## Comparison of parameter 'e' 
## 
##         Estimate Std. Error t-value   p-value    
## CU-SA   -706.883   2678.169 -0.2639 0.7921196    
## CU-VE  -9422.809  22562.901 -0.4176 0.6767090    
## CU-FO    -92.129    174.820 -0.5270 0.5988318    
## CU-J2 -17359.106  39254.294 -0.4422 0.6588480    
## CU-UM    203.894    157.912  1.2912 0.1982576    
## SA-VE  -8715.926  22720.212 -0.3836 0.7017031    
## SA-FO    614.753   2674.721  0.2298 0.8184727    
## SA-J2 -16652.224  39344.928 -0.4232 0.6726167    
## SA-UM    910.777   2673.669  0.3406 0.7337581    
## VE-FO   9330.680  22562.492  0.4135 0.6796866    
## VE-J2  -7936.298  45276.183 -0.1753 0.8610474    
## VE-UM   9626.703  22562.367  0.4267 0.6701179    
## FO-J2 -17266.977  39254.059 -0.4399 0.6605421    
## FO-UM    296.024     80.504  3.6771 0.0003095 ***
## J2-UM  17563.001  39253.988  0.4474 0.6550984    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

modelFit(modelQNC21DR4)