On 26/11/2018 7:13 AM, Thanh Tran wrote:
Dear all,
I'm trying to plot a surface over the x-y plane. In my data, the response
is KIC, and four factors are AC, AV, T, and Temp. A typical second-degree
response modeling is as follows
data<-read.csv("2.csv", header =T)
mod <-
lm(KIC~AC+I(AC^2)+AV+I(AV^2)+T+I(T^2)+Temp+I(Temp^2)+AC:AV+AC:T+AC:Temp+AV:T+AV:Temp+T:Temp,
+ data = data)
For two factors, you could use this code:
pred <- function(AC, AV, Temp, T) predict(mod, newdata = data.frame(AC,
AV, Temp, T))
library(rgl)
persp3d(pred, xlim = c(-1, 1), # The range of values for AC
ylim = c(-1, 1), # The range for AV
xlab = "AC", ylab = "AV", zlab = "KIC",
colour = rainbow, # or a fixed colour, or another fn
otherargs = list(Temp = 0, T = 0))
The otherargs list should contain the values of the two factors to your
model that you are holding fixed while plotting the two that are not fixed.
This Stackoverflow answer
https://stackoverflow.com/questions/53349811/how-to-draw-a-response-surface-plot-for-three-factorial-design/53350259#53350259
describes a way to plot the response to 3 factors at once.
Duncan Murdoch
I want to have a response surface of KIC with two factors, i.e., AC and AV
as shown in the attached figure.
When I run the below code, I have a problem which indicates “object 'AC'
not found” even though I added “data = data”
persp(AC,AV,KIC~ AC + I(AC^2) + AV + I(AV^2) + AC:AV,image = TRUE,theta=30,
+ data = data)
Error in persp(AC, AV, KIC ~ AC + I(AC^2) + AV + I(AV^2) + AC:AV,
image = TRUE, :
object 'AC' not found
If anyone has any experience about what would be the reason for error or
how I can solve it? Is there other simple function to plot the response
surface?
I really appreciate your support and help.
Best regards,
Nhat Tran
Ps: I also added a CSV file for practicing R.
______________________________________________
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
______________________________________________
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
