Displaying 1 result from an estimated 1 matches for "ogmod1".
Did you mean:
mod1
2010 Jan 06
1
Lattice Plot formatting problem/s
...for a lattice plot I am
trying to produce.
Here is a pic of the plot as I get it now
http://i235.photobucket.com/albums/ee37/scotrivers/lattice_plot01.jpg
and here is the code I am using:
RN<-read.csv("N:/data.dat",header=T)
DATA<-RN[is.element(RN$age,0:3),] #select data on age
ogmod1<-glm(mig~age+length,family=binomial(link="logit"),data=DATA)
design1<-expand.grid(age=c(0:3),length=seq(30,300,by=10))
design1$pbs<-predict(ogmod1,new=design1,type="response")
design1$migr<-round(design1$pbs)
library(lattice)
wireframe(pbs~age*length,data=design1,asp...