Displaying 1 result from an estimated 1 matches for "logtps".
2008 Apr 30
2
wireframe - add data points
...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.1, by=0.25),logwinterp=...