search for: totalfrost

...oints on which the predicted probabilities are based to the wireframe. I’ve followed some of the previous post on this but get stuck : ## first part yields the wireframe## setwd("C:/Temp") rnp <- read.table("interact.csv",header=T,sep=";") status <-rnp[,1:1] totalfrost <- rnp[,2:2] logtps <- rnp[,3:3] logpd <- rnp[,4:4] logwinterp <- rnp[,5:5] model <- glm(status~totalfrost+logtps+logpd+logwinterp+totalfrost*logtps+totalfrost*logwinterp,binomial) abc <- expand.grid(totalfrost=seq(-1.9,3.6,by=0.25),logtps=seq(-3.3,1.1,by=0.25),logpd=seq(0.95,4...