Código Rstudio

Vista previa del material en texto

A) Código de Rstudio 
 
trsm <- ccd(basis = 2, randomize = F,n0 = c(3,3), coding = list(x1~(RA-163)/15, x2~(CE-250)/35)) 
trsm$Y <- c(199, 370, 292, 549, 356,352, 357, 195, 486, 258,450, 364, 369, 366) 
 
first <- rsm (Y~SO(x1, x2), data = trsm) 
summary (first) 
 
 
second <- rsm(Y~FO(x1,x2)+TWI(x1,x2)+PQ(x1),trsm) 
summary (second) 
 
 
par(mfrow = c(1,2)) 
persp(second, x2 ~ x1, zlab = "Respuesta", contours = list(z = "bottom", col = "colors"), at = 
c(summary(second$canonical$xs)), theta = 15, phi = 20) 
 
 
contour (trsm,~x1+x2, image = T) 
 
 
point <- data.frame(RA=c(159,160,161,162,163),CE=c(rep(260,15))) 
point 
Yp <- predict(second, coded.data(point,x1~(RA-163)/15,x2~(CE-250)/35)) 
time_400 <- data.frame(point,Yp) 
time_400 
 
 
 
 
 
 
 
 
 
 
 
 
 
point <- data.frame(RA=c(165,166,167,168,169),CE=c(rep(260,15))) 
point 
Yp <- predict(second, coded.data(point,x1~(RA-163)/15,x2~(CE-250)/35)) 
time_400 <- data.frame(point,Yp) 
time_400 
 
 
B.- POMEDIO Y DESVIACIÓN ESTANDAR 
 
De acuerdo con las coordenadas otorgadas por el programa se realizaron de manera manual las 
corridas de los datos 
 
distancias <- c(322,339,333,323,331,323.320,314,330,338,325,324,327,325,340) 
promedio <- mean(distancias) 
promedio 
desvstd <- sd(distancias) 
desvstd 
 
 
 
Una vez realizado el procedimiento en el programa de Rstudio, se han obtenido los datos anteriores, 
y de acuerdo con lo revisado en clase, es más complejo predecir los resultados a distancias mayores, 
por lo que es necesario destacar esta información, debido a que no se han obtenido las 400 unidades 
deseadas, aunque en los resultados obtenidos en el analisis de varianza si indica que es un proceso 
predictivo. 
Sin embargo, la gráfica de contornos si predice de mejor manera los valores a utilizar para obtener 
las unidades deseadas. 
Nota: He realizado varias pruebas con la catapulta para verificar, y a distancias menores si, se 
obtienen los valores de predicción. 
Por otro lado, me disculpo por la tardanza, el programa me falló en repetidas ocasiones, adjunto 
una imagen de la evidencia