library(readxl)
library(extrafont)
## Registering fonts with R
#font_import()
AO50 <- t(read_xlsx("50AO.xlsx"))
## New names:
## * `` -> ...1
AO25 <- t(read_xlsx("25AO.xlsx"))
## New names:
## * `` -> ...1
AODW <- t(read_xlsx("DWAO.xlsx"))
## New names:
## * `` -> ...1
x = 1:97
plot(x,AO50[2,1:97]-297.95,xlab = "Tiempo (seg)", ylab = expression(Delta~"T(°C)"), family="Times New Roman",type = 'l')
lines(spline(x,AO50[3,1:97]-297.95),n=201,method = 'natural',lty=1,lwd=1, col=2)
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "n" is not a graphical
## parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "method" is not a
## graphical parameter
lines(spline(x,AO50[4,1:97]-297.95),n=201,method = 'natural',lty=1,lwd=1, col=3)
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "n" is not a graphical
## parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "method" is not a
## graphical parameter
lines(spline(x,AO50[5,1:97]-297.95),n=201,method = 'natural',lty=1,lwd=1, col=4)
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "n" is not a graphical
## parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "method" is not a
## graphical parameter
lines(spline(x,AO25[2,1:97]-297.95),n=201,method = 'natural',lty=2,lwd=1, col=1)
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "n" is not a graphical
## parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "method" is not a
## graphical parameter
lines(spline(x,AO25[3,1:97]-297.95),n=201,method = 'natural',lty=2,lwd=1, col=2)
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "n" is not a graphical
## parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "method" is not a
## graphical parameter
lines(spline(x,AO25[4,1:97]-297.95),n=201,method = 'natural',lty=2,lwd=1, col=3)
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "n" is not a graphical
## parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "method" is not a
## graphical parameter
lines(spline(x,AO25[5,1:97]-297.95),n=201,method = 'natural',lty=2,lwd=1, col=4)
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "n" is not a graphical
## parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "method" is not a
## graphical parameter
lines(spline(x,AODW[2,1:97]-297.95),n=201,method = 'natural',lty=3,lwd=1, col=1)
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "n" is not a graphical
## parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "method" is not a
## graphical parameter
lines(spline(x,AODW[3,1:97]-297.95),n=201,method = 'natural',lty=3,lwd=1, col=2)
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "n" is not a graphical
## parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "method" is not a
## graphical parameter
lines(spline(x,AODW[4,1:97]-297.95),n=201,method = 'natural',lty=3,lwd=1, col=3)
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "n" is not a graphical
## parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "method" is not a
## graphical parameter
lines(spline(x,AODW[5,1:97]-297.95),n=201,method = 'natural',lty=3,lwd=1, col=4)
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "n" is not a graphical
## parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "method" is not a
## graphical parameter
labels_image2 = c("Medicion experimental para nanofluido al 0.050%", 'Modelo completo estimado por LM para nanofluido al 0.050%','Modelo RLM estimado por LM para nanofluido al 0.050%', 'Modelo RLM estimado por GN para nanofluido al 0.050%',
"Medicion experimental para nanofluido al 0.025%", 'Modelo completo estimado por LM para nanofluido al 0.025%','Modelo RLM estimado por LM para nanofluido al 0.025%','Modelo RLM estimado por GN para nanofluido al 0.025%',
'Medicion experimental para el agua destilada', 'Modelo completo estimado por LM para el agua destilada','Modelo RLM estimado por LM para el agua destilada','Modelo RLM estimado por GN para el agua destilada')
legend(x = "topleft",c(0, 0.5),legend = labels_image2, lty = c(1,1,1,1,2,2,2,2,3,3,3,3),col=c(1,2,3,4,1,2,3,4,1,2,3,4),bty='n',xpd = TRUE,cex=0.5)

library(readxl)
Y_real<- read_xlsx("y_real_image4.xlsx")
y_pred_200 <- read_xlsx("y_predT_200_image4.xlsx")
error_absoluto = Y_real - y_pred_200
error_absoluto=t(error_absoluto)
plot(error_absoluto[1:97,1])
