DataPlaguicidas

df<-read.csv("Cyproconazole.csv")
df

##    Conc   SolvArea    MatArea
## 1   0.1   8954.362   6653.599
## 2   0.5  14246.035  10396.710
## 3   2.0  31971.220  25477.647
## 4   5.0  55827.935  43897.521
## 5  20.0 161227.412 126715.373
## 6  50.0 295344.068 248193.188
## 7 100.0 408306.020 362952.037

Curvas del ciproconazol: Efecto de la Matriz

ysolv<-df$SolvArea
ymat<-df$MatArea
x<-df$Conc


modelo<-lm(ymat~x+I(x*x))
s<-summary(modelo)
LOD<-s$coefficients[1,2]/s$coefficients[1]*3.3*8
b0mat=s$coefficients[1]
b1mat=s$coefficients[2]
b2mat=s$coefficients[3]


xs <- seq(0, 100, by=0.1)
ycalc1=b0mat+b1mat*xs+b2mat*xs*xs

yexpMat=modelo$fitted.values

modelo<-lm(ysolv~x+I(x*x))
s<-summary(modelo)

b0solv=s$coefficients[1]
b1solv=s$coefficients[2]
b2solv=s$coefficients[3]


xs <- seq(0, 100, by=0.1)
ycalc2=b0solv+b1solv*xs+b2solv*xs*xs

yexpMat=modelo$fitted.values

print(s)

## 
## Call:
## lm(formula = ysolv ~ x + I(x * x))
## 
## Residuals:
##     1     2     3     4     5     6     7 
## -6386 -4120  2362  4227  9562 -7000  1356 
## 
## Coefficients:
##              Estimate Std. Error t value Pr(>|t|)    
## (Intercept) 14582.288   3949.238   3.692 0.020974 *  
## x            7586.805    313.256  24.219 1.72e-05 ***
## I(x * x)      -36.631      3.165 -11.575 0.000318 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 7475 on 4 degrees of freedom
## Multiple R-squared:  0.9985, Adjusted R-squared:  0.9977 
## F-statistic:  1325 on 2 and 4 DF,  p-value: 2.271e-06

n=length(x)



plot(x,ymat,main="Ciproconazol",xlab="[Ciproconazol](ug/Kg)",ylab="Area", ylim=c(0,450000))
points(x,ysolv,col="blue")
lines(xs,ycalc1,lty=1)
lines(xs,ycalc2,lty=1,col="blue")

yr=ycalc1/ycalc2*100

plot(xs,yr,type="l",col="orange",main="Ciproconazol: %Efecto de Matriz",xlab="[Ciproconazol](ug/Kg)",ylab="%EM", ylim=c(0,100))