[R] scatter.smooth() fitted by loess

[David Winsemius]

On Oct 19, 2010, at 5:43 PM, phoebe kong wrote:

> Hi there,
>
> I would like to draw a scatter plot and fit a smooth line by loess.
> Below is the data.
> However, the curve line started from 0, which my "resid" list doesn't
> consist of 0 value.
> It returned some warnings which I don't know if this is the reason
> affecting such problem. Here I also attached the warning messages.
> Please let me know if there is a solution to fix this. Thank you very
> much!

Yes, Increase the span parameter value. (And replace the missing comma  
in line 3 of the offered example.)

-- 
David

>
> YMRS_Sum<-c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0,
> 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 2, 0, 0, 2, 2, 0, 0, 0, 1, 1,
> 4, 0, 0, 2, 0, 1, 0, 2, 0, 0, 0, 0, 0, 0, 0, 3, 1, 4, 2, 1, 2, 0, 1,
> 0, 0, 0, 0, 2, 0, 0, 3, 0, 0, 0, 0, 2, 0, 4, 0, 2, 0, 0)
>
> resid<-c(80.1150, 84.0279, 88.3736, 90.1557, 55.1979, 78.4293,
> 87.4367, 74.0271, 80.8871, 91.5685, 82.4154, 73.3080, 66.7786,
> 70.2486, 82.4971, 77.3792, 70.7731, 66.9593, 85.5515 81.4071, 68.6646,
> 89.8271, 91.6041, 85.1980, 80.6071, 86.4362, 86.2915, 86.4493,
> 87.8664, 84.5150, 64.4975, 79.1246, 84.9350, 89.1608, 92.7546,
> 70.0253, 81.4146, 73.2755, 82.5200, 79.7164, 92.0786, 82.5633,
> 84.4336, 84.0193, 64.8029, 87.4864, 86.3338, 75.6758, 86.8567,
> 85.1077, 88.9533, 81.7240, 84.1713, 80.0400, 77.6050, 81.4436,
> 83.8379, 72.5050, 80.5423, 83.2564, 84.1436, 90.0662, 84.5293,
> 81.6771, 90.6425, 90.3285, 76.2371, 87.3625, 70.7917, 77.0993,
> 88.3608, 89.7200, 79.1031, 79.7421, 84.2469, 83.9371, 73.8800,
> 89.3921, 89.3900, 86.8921, 85.7036, 85.2664, 83.8700, 90.5493)
>
>> scatter.smooth(YMRS_Sum,resid)
>
> Warning messages:
> 1: at  -0.02
> 2: radius  0.0004
> 3: all data on boundary of neighborhood. make span bigger
> 4: pseudoinverse used at -0.02
> 5: neighborhood radius 0.02
> 6: reciprocal condition number  nan
> 7: zero-width neighborhood. make span bigger
> 8: There are other near singularities as well. 1
> 9: at  -0.02
> 10: radius  0.0004
> 11: all data on boundary of neighborhood. make span bigger
> 12: pseudoinverse used at -0.02
> 13: neighborhood radius 0.02
> 14: reciprocal condition number  nan
> 15: zero-width neighborhood. make span bigger
> 16: There are other near singularities as well. 1
> 17: at  -0.02
> 18: radius  0.0004
> 19: all data on boundary of neighborhood. make span bigger
> 20: pseudoinverse used at -0.02
> 21: neighborhood radius 0.02
> 22: reciprocal condition number  nan
> 23: zero-width neighborhood. make span bigger
> 24: There are other near singularities as well. 1
> 25: at  -0.02
> 26: radius  0.0004
> 27: all data on boundary of neighborhood. make span bigger
> 28: pseudoinverse used at -0.02
> 29: neighborhood radius 0.02
> 30: reciprocal condition number  nan
> 31: zero-width neighborhood. make span bigger
> 32: There are other near singularities as well. 1
> 33: In simpleLoess(y, x, w, span, degree, FALSE, FALSE, normalize =
> FALSE,  ... :
>  at  -0.02
> 34: In simpleLoess(y, x, w, span, degree, FALSE, FALSE, normalize =
> FALSE,  ... :
>  radius  0.0004
> 35: In simpleLoess(y, x, w, span, degree, FALSE, FALSE, normalize =
> FALSE,  ... :
>  all data on boundary of neighborhood. make span bigger
> 36: In simpleLoess(y, x, w, span, degree, FALSE, FALSE, normalize =
> FALSE,  ... :
>  pseudoinverse used at -0.02
> 37: In simpleLoess(y, x, w, span, degree, FALSE, FALSE, normalize =
> FALSE,  ... :
>  neighborhood radius 0.02
> 38: In simpleLoess(y, x, w, span, degree, FALSE, FALSE, normalize =
> FALSE,  ... :
>  reciprocal condition number  nan
> 39: In simpleLoess(y, x, w, span, degree, FALSE, FALSE, normalize =
> FALSE,  ... :
>  zero-width neighborhood. make span bigger
> 40: In simpleLoess(y, x, w, span, degree, FALSE, FALSE, normalize =
> FALSE,  ... :
>  There are other near singularities as well. 1
>
>
> Thanks,
> Phoebe
>
> ______________________________________________
> 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.

David Winsemius, MD
West Hartford, CT