data(DataNAT) # dataset
Exp <- DataNAT[,2] # identification of the experiments
X <- DataNAT[,3:6] # independent variables (components)
Y <- DataNAT[,11] # dependent variable (response Bitterness)
# effects o the mixtures in the reference mixture
Effects <- rep(c(-0.1,0,0.1,0.2,0.3,0.4,0.5,0.6,0.7),4)
Conc <- as.matrix(DataNAT[,7]) # covariable (process variable)
Res <- Blend(exp = Exp, X = X, Y = Y, conc = Conc, effects = Effects)
print("Predicted and observed values"); Res$MPred
print("Values predicted by components:"); Res$MCPred
print("Design of the experiments:"); Res$MExp
print("Estimates of the linear model parameters:"); Res$Theta
Run the code above in your browser using DataLab