[R] locfit questions/problems
Liaw, Andy
andy_liaw at merck.com
Mon Jan 25 18:45:56 CET 2010
Just replacing preplot() with predict() should be fine.
BTW, it's always a good idea to specify the version of the package
you're using as well.
Best,
Andy
From: mhall at berkeley.edu
>
> Hi,
>
> I'm trying to work through the examples and code in Loader's
> LOCAL REGRESSION AND LIKELIHOOD, and have run into a problem
> with the code for one sided smoothing and change point analysis
> (p. 110-112).
>
> The code, after loading locfit:
>
> midp<-(1945:1988)+0.5
> fitl<-locfit(thickness~left(year), data=penny, alpha=c(0,10),
> deg=1, ev=midp)
> fitr<-locfit(thickness~right(year),data=penny, alpha=c(0,10),
> deg=1, ev=midp)
> plot((preplot(fitr)-preplot(fitl))^2, type="b")
>
> the error message crops up after the plot command and is:
>
> Error in preplot(fitr) - preplot(fitl) :
> non-numeric argument to binary operator
>
> looking at preplot(fitr0 and preplot(fitl) shows a list with
> numeric values.
>
> Any suggestions on how to resolve this problem? I'm using R
> 2.7.1 on a
> Windows XP machine with Service pack 2.
>
> best and thanks, Mark Hall
>
> PS--a similar error case occurs using R 2.3.1 and 2.5.1
>
> ______________________________________________
> 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.
>
Notice: This e-mail message, together with any attachme...{{dropped:10}}
More information about the R-help
mailing list