/* llrtest.gpr: program to test llr.g

                John Rust, University of Maryland, October, 2006 */

    new;

    nobs=5000;
    kerneltyp="gaussian";

    x=rndu(nobs,1);

    y=3+sqrt(x)+rndn(nobs,1);


    xplot=seqa(0,.01,101);

    t=hsec;
    yhat=llr(xplot,y,x);

    t=hsec-t;

    rf=3+sqrt(xplot);

    "finished in "$+ftos(t/100,"*.*lf",9,5)$+" seconds ";

    xy(xplot,rf~yhat);