I tried to create a response surface with 56 design variables. But I can't generate a second order response surface. I am successfully able to generate first order surface. How can I generate a response surface with second order accuracy?
Hi there,
Adding some more information to your post, including code with relevant data and some context, may help you garner useful replies. Can you please add a reprex (reproducible example) with a little context to your question? This will ensure we're all looking at the same data and code. A guide for creating a reprex can be found here.
Thank you for the message. How can I upload my excel file for the example? The code I used to implement the data is here.
library(rsm)
library(readxl)
Data_82DV_RSM <- read_excel("D:/Internship and Thesis/RSM Tutorials and Software/Data for Trials/82 Variables/Data_82DV_RSM.xlsx", sheet = "DV_RSM _Input1")
View(Data_82DV_RSM)
CODED_82DV <- coded.data(Data_82DV_RSM, Y1 ~ (x1 - 0.95)/0.2, Y2 ~ (x2 - 11)/2, Y3 ~ (x3 - 0.4)/0.1, Y4 ~ (x4 - 17.25)/3.25, Y5 ~ (x5 - 1.13)/0.2, Y6 ~ (x6 - 9)/2, Y7 ~ (x7 - 0.2)/0.1, Y8 ~ (x8 - 11)/2, Y9 ~ (x9 - 0.95)/0.15, Y10 ~ (x10 - 11)/2, Y11 ~ (x11 - 0.4)/0.1, Y12 ~ (x12 - 19)/5, Y13 ~ (x13 - 0.46)/0.07, Y14 ~ (x14 - 7.25)/1.5, Y15 ~ (x15 - 0.165)/0.035, Y16 ~ (x16 - 8.75)/2, Y17 ~ (x17 - 0.645)/0.145, Y18 ~ (x18 - 8.5)/2.5, Y19 ~ (x19 - 0.25)/0.05, Y20 ~ (x20 - 14)/4, Y21 ~ (x21 - 1.015)/0.185, Y22 ~ (x22 - 11)/2, Y23 ~ (x23 - 0.4)/0.1, Y24 ~ (x24 - 18.25)/4.25, Y25 ~ (x25 - 1.166)/0.2, Y26 ~ (x26 - 12)/2, Y27 ~ (x27 - 0.3)/0.1, Y28 ~ (x28 - 11.25)/1.75, Y29 ~ (x29 - 1.15600000000005)/0.19, Y30 ~ (x30 - 14)/3, Y31 ~ (x31 - 20)/5, Y32 ~ (x32 - 0.4)/0.2, Y33 ~ (x33 - 1.16600000000005)/0.2, Y34 ~ (x34 - 13.25)/1.75, Y35 ~ (x35 - 0.4)/0.1, Y36 ~ (x36 - 18.5)/4.5, Y37 ~ (x37 - 1.65)/0.35, Y38 ~ (x38 - 11)/2, Y39 ~ (x39 - 2.4)/0.4, Y40 ~ (x40 - 9.5)/2, Y41 ~ (x41 - 0.65)/0.15, Y42 ~ (x42 - 9)/2, Y43 ~ (x43 - 0.25)/0.1, Y44 ~ (x44 - 14)/3, Y45 ~ (x45 - 0.45)/0.15, Y46 ~ (x46 - 7)/2, Y47 ~ (x47 - 0.1)/0.05, Y48 ~ (x48 - 9)/2, Y49 ~ (x49 - 0.445)/0.145, Y50 ~ (x50 - 7)/2, Y51 ~ (x51 - 0.1)/0.05, Y52 ~ (x52 - 9)/2, Y53 ~ (x53 - 1.65)/0.25, Y54 ~ (x54 - 2.4)/0.4, Y55 ~ (x55 - 1.595)/0.195, Y56 ~ (x56 - 2.345)/0.345)
CODED_56DV_RSM <- rsm(W_HULL ~ FO (Y1,Y2,Y3,Y4,Y5,Y6,Y7,Y8,Y9,Y10,Y11,Y12,Y13,Y14,Y15,Y16,Y17,Y18,Y19,Y20,Y21,Y22,Y23,Y24,Y25,Y26,Y27,Y28,Y29,Y30,Y31,Y32,Y33,Y34,Y35,Y36,Y37,Y38,Y39,Y40,Y41,Y42,Y43,Y44,Y45,Y46,Y47,Y48,Y49,Y50,Y51,Y52,Y53,Y54,Y55,Y56), data = CODED_82DV)
summary(CODED_56DV_RSM)
I couldn't put a second order polynomial in the rsm code.
Thank you