Second order Response Surface


#1

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?


#3

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.


#5

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