# NOT RUN {
T<-60
change<-35
x<-rnorm(n=T, m=0, sd=1)
e<-scale(rt(n=T,3), scale=FALSE)
y1<-5*x[1:(change-1)]+e[1:(change-1)]
y2<--2*x[change:T]+e[change:T]
y<-c(y1,y2)
model <- changingModel(x=as.data.frame(x),
y=as.data.frame(y),
c(change))
# }
Run the code above in your browser using DataLab