Dalphin, Mark
2006-Sep-08 19:38 UTC
[R] Confidence intervals on Lowess (eg SiZer from J.S.Marron)
Hi, I have some very noisy, relatively sparse data; a biological response of roughly ~8 subjects at ~8 times points). I've been following the data trend using a lowess line, over-plotted with several values of bandwidth, 'f <- seq(0.3, 0.9, by=0.1)'. At this point, we have no models for these data. I wonder if there is any way under R to assign some sort of confidence interval to the lowess line. For example, I have seen a method from the lab of J.S.Marron, which his group has implemented in the Matlab program, SiZer. http://www.stat.unc.edu/faculty/marron/DataAnalyses/SiZer_Intro.html and, the more interesting: http://www.stat.unc.edu/faculty/marron/DataAnalyses/SiZer/SiZer_Basics.html An implementation under R of SiZer would obviously answer this question; does it exist? Suggestions of alternative approaches would also be welcome. My searchs of the R-maillist archive for 'sizer', 'Marron' and 'scale-space' didn't return anything that I recognized as useful here, but I am new to the idea of confidence in a lowess line so I may not be using the appropriate vocabulary. Thank-you, Mark ---------------------- Mark Dalphin Dept Comp Biol, M/S AW2/D3262 Amgen, Inc. 1201 Amgen Court W Seattle, WA 98119 Phone: +1-206-265-7951
Ben Bolker
2006-Sep-08 19:59 UTC
[R] Confidence intervals on Lowess (eg SiZer from J.S.Marron)
Dalphin, Mark <mdalphin <at> amgen.com> writes:> I wonder if there is any way under R to assign some sort of confidence > interval to the > lowess line.bootstrap (e.g. in the simpleboot package)? Ben Bolker
Seemingly Similar Threads
- Regression test failed when building on an "older" Linux cluster
- How to reformat data from database into data.frame?
- Fresh build from source of R-3.2.5 failing "make check" under 64-bit Ubuntu
- Not able to write to PostgreSQL database using "dbWriteTable"
- Core dump using DBI/ROracle