[R] Smooth Spline

From: kate <yhsu6_at_uiuc.edu>
Date: Thu, 05 Jun 2008 00:46:43 -0500


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]]



R-help_at_r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. Received on Thu 05 Jun 2008 - 06:04:26 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Thu 05 Jun 2008 - 07:30:37 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.

list of date sections of archive