Displaying 1 result from an estimated 1 matches for "bavg".
Did you mean:
bang
2007 May 29
1
rgl.postscript
...lp<-m3[,-2]
p1<-rbind(up1, lp[-1,])
p2<-rbind(up2, lp[-1,])
sp1<-spline(p1)
sp2<-spline(p2)
sp1m<-cbind(sp1$x,sp1$y)
sp2m<-cbind(sp2$x,sp2$y)
ge9<-(sp1$x>=9)
ge9recs<-seq(1,length(ge9))[ge9]
b1<-sp1m[ge9recs,]
b2<-sp2m[ge9recs,]
b1b2<-cbind(b1[,2],b2[,2])
bavg<-apply(b1b2,1,mean)
blow<-cbind(sp1m[ge9recs,1],bavg)
path.one<-rbind(sp1m[-ge9recs,],blow)
path.two<-rbind(sp2m[-ge9recs,],blow)
uy1<-9-path.one[-ge9recs,1]
ly1<-rep(0,length(ge9recs))
y1<-c(uy1,ly1)
uy2<--1*(9-path.two[-ge9recs,1])
ly2<-rep(0,length(ge9recs))
y2<-...