Displaying 1 result from an estimated 1 matches for "newax".
Did you mean:
newa
2011 Jan 12
1
graphics: 3D regression plane
...id),min))
t0<-as.Date(d1,origin=as.Date("1970-01-01"));t0
A<-data.frame(id=c("1","2","3"),t=t0)
B<-data.frame(id=id,tempo=time)
C<-merge(A,B);C
rd<-as.vector(C$tempo-C$t);rd #time centered on sbj specific first
occurrence
mod<-lm(y~rd*k)
newax<- expand.grid(
days = giorni<-seq(min(rd),max(rd), length=100),
expl= esplic<- seq(min(k), max(k), length=100)
)
fit <- predict(mod,data.frame(rd=giorni,k=esplic))
graph <- persp(x=giorni, y=esplic,fit,
expand=0.5, ticktype="detailed", theta=-45) #error : z ar...