#==== nonparametric estimation using orthogonal series expansion  ==== 

source("ose_fun.r")

x <- c(rnorm(50),1+rnorm(50))


est <- ose_density(x,type=1) # Kronmal-Tarter t=1
plot(est$tm,est$pdf)

q <- ose_quantile(x,probs=c(0.05,0.25,0.5),type=4)  
print("Kronmal-Tarter t=4")
print(q)

q2 <- ose_quantile(x,probs=c(0.05,0.25,0.5),type=-2) 
print("Diggle-Hall")
print(q2)