11b_rs.R

## ---- include=FALSE-----------------------------------------------------------------------------
library(tidyverse)
library(broom)
library(flextable)
library(knitr)





## ----cement-------------------------------------------------------------------------------------
library(daewr)
cement


## -----------------------------------------------------------------------------------------------
library(rsm)
grout.lin <- rsm(y ~ SO(x1, x2, x3), data = cement, subset = Block == 1)
anova(grout.lin)


## -----------------------------------------------------------------------------------------------
Treb
treb.quad <- rsm(y ~ SO(x1, x2, x3), data = Treb)
summary(treb.quad)


## -----------------------------------------------------------------------------------------------
par(mfrow=c(2,2))
contour(treb.quad, ~ x1+x2+x3)


## -----------------------------------------------------------------------------------------------
par(mfrow=c(2,2))
contour(treb.quad, ~ x1+x2+x3, at = xs(treb.quad))


## -----------------------------------------------------------------------------------------------
par(mfrow=c(2,2))
persp(treb.quad, ~x1+x2+x3, zlab = "Distancia")


## -----------------------------------------------------------------------------------------------
ridge <- steepest(treb.quad, dist = seq(0, 1.412, by=0.1), descent = FALSE)
ridge


## -----------------------------------------------------------------------------------------------
library(daewr)
Design<-DefScreen(m=8,c=0)
Design


## -----------------------------------------------------------------------------------------------
library(daewr)
design<-DefScreen(m=5,c=0,center=2,randomize=FALSE)
design


## -----------------------------------------------------------------------------------------------
library(daewr)
Smeso<-c(241,295,260,338,320,265,275,248,92.5,383,313,305,304)
FitDefSc(Smeso,design,alpha=.05)


## -----------------------------------------------------------------------------------------------
library(daewr)
FitDefSc(Smeso,design,alpha=.1)