data("Bahrami1")
gender<-Bahrami1$ GENDER
duration <-Bahrami1$ DURATION
y<-Bahrami1$ BMI
X=cbind(gender,duration )
Xstar=cbind(1,X)
p=2;
n=15;
ini=c(0,rep(1,p));
LICurvature(ini,X=X,Xstar=Xstar,y=y,n=n,p=p)
Run the code above in your browser using DataLab