Displaying 1 result from an estimated 1 matches for "sp1m".
Did you mean:
sp1
2007 May 29
1
rgl.postscript
...x3<-seq(10,18)
y1<-seq(8,0)
y2<--1*y1
y3<-rep(0,9)
m1<-cbind(x1,y1,z1)
m2<-cbind(x2,y2,z2)
m3<-cbind(x3,y3,z3)
m3<-rbind(m2[9,],m3)
up1<-m1[,-2]
up2<-m2[,-2]
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....