Hi,

I have three original curves as follows,
n<-seq(20,200,by=10)

t<-c(0.1138, 0.1639, 0.2051, 0.2473, 0.2890, 0.3304, 0.3827, 0.4075, 0.4618, 0.4944,
0.5209, 0.5562, 0.5935, 0.6197, 0.6523, 0.6771, 0.6984, 0.7209, 0.7453)

es<-c(0.3682, 0.4268, 0.5585, 0.6095, 0.7023, 0.7534, 0.8225, 0.8471, 0.8964, 0.9098, 0.9371, 0.9514, 0.9685, 0.9747, 0.9812, 0.9859, 0.9905, 0.9923, 0.9940)

coves<-c(0.3793, 0.4935, 0.5472, 0.6366, 0.6882, 0.7578, 0.8058, 0.8457, 0.8783, 0.9019, 0.9212, 0.9417, 0.9558, 0.9662, 0.9732, 0.9802, 0.9828, 0.9881, 0.9896)

matplot(n,cbind(coves,es,t),type ='l', col=c("red","blue","black"),xlab='n',ylab='power')

I want power curves to be smooth (avoid sawtooth curve), I try the follows,
##Spline
library(modreg)
st<- smooth.spline(t) 
ses<- smooth.spline(es) 
scoves<- smooth.spline(coves) 

plot(scoves,col='red', type='l')
lines(ses, col='blue',type='l')
lines(st,col='black',type='l')

The shapes are what I want, but I want x-axis consistent to the correspond to n. Does anyone know the solution?

Thanks,

Kate 
	[[alternative HTML version deleted]]

