Hi Stephen if you want and uncomplicated graph you can do the following (data is from https://stat.ethz.ch/pipermail/r-help/2008-December/181887.html which has an example of smoothing with a panel function a (follow the thread)) May be un necessary for what you want
dput(test) structure(list(x = 1:200, y = c(-0.211833603383945, -1.0415911337336, -1.1533075572692, 0.321531497643745, -1.5001298759544, -0.445533260997329, 1.73404543257509, 0.511295615763455, 0.0996450365533331, -0.0578911103962697, -1.74278762900362, -1.32495298319389, -0.547933881534948, -1.45638427585301, 0.0826868204409038, 0.927578947155782, -0.716769334937502, 0.962399675187835, 1.54588458200656, -1.00976361360517, 0.557410853590579, 0.168781513971099, 0.155258858301712, 2.36776473870934, -1.58564441351633, -1.10353350945666, 0.902787057776989, 0.0683183515567076, -1.14946795040125, -0.90127498047004, -1.19845352398361, -0.510163823550385, 0.306246343638814, -0.0678770276660271, 0.372202425855059, 1.14449122171645, -1.43961215364919, -0.745407767319761, -1.08426028148107, -0.911810459312264, 0.254473860506602, 0.0545428435691033, 0.86939682821168, -0.495822138795001, 0.788983956205307, -0.0482778318087694, 0.28444147880148, -2.3449130823549, 0.352139869980497, -0.73670929392907, 0.152865554816491, 0.466635540442577, -1.34862166637197, -0.371330904577685, 2.04819500702545, -1.03035368989461, -0.808913708097202, -0.899512468883577, 1.32164919515955, -1.90258984438615, -1.43516196895465, 0.0465901079414443, -0.050190824605222, -1.28753166678936, 1.15411244955901, -1.37573809203985, -1.14829872319585, 1.41625774885229, 0.105763612116132, -0.384659376665841, 1.78448216497794, -1.1152932761564, -0.0721496590810866, -0.654929692334926, 1.23475033091273, -0.355004059316614, 0.0242929525502168, -1.12014689112359, 0.131516716053627, -1.01789380466415, 1.54412508731677, -0.612498366500995, -1.47544725093144, -1.33884056677665, 0.116000294553631, 0.441007318370771, 0.253328166388004, -0.601324732898278, -0.681872432317887, -0.638406595370473, 1.1255146487638, 1.68030106840915, -0.0458152886664351, -0.359702695597414, -1.27276961899765, 1.65150890035352, 0.970733803202561, -0.299164223313285, 0.0863388612193544, -1.56778219134085, -0.934175135343763, -1.75098468827598, -0.350775124642839, -0.497774672576354, 0.771324636489948, -0.69490617657808, -0.776286103585292, 0.0401278987521509, -0.965161644356313, -0.443562752383677, 0.482231363899495, 0.299362730723139, -1.13975716906454, 0.750123380303915, -0.346046820684503, -0.298660657494398, 1.15743160079942, 1.37831825909263, -0.561195292140498, -0.727441974730694, 0.542558262494051, 1.07423749905296, -0.671511488353386, -0.140489880785358, 1.46724871140033, 0.563360651277434, 1.24069952075524, -0.253423940974, 0.196985403588085, -1.29263886475009, -0.700235249467159, 0.386254280503009, -0.489119591847831, -0.863120570714887, 0.21167415527797, -0.104722187456112, 0.885271725582058, -1.33715921875481, -0.419738882938569, 0.114904893946579, 0.579793678083023, 1.3375478692142, 0.00553313645797594, -1.0480172540963, -0.293082520385054, -0.607995018773811, 0.22269595909434, 0.404577306401601, 0.212081041192621, -1.28370971577478, -0.186741150561225, 1.47504624554164, -0.335669248934062, -0.0752650377217833, 0.938248892501711, -1.19915518471386, -0.823130111133665, 1.77585103279627, -0.726506122852284, -0.854585428953646, 0.18721613856176, 0.664751685525484, -0.378660174874064, -0.406279864316906, -0.377953114167246, -0.578277117855972, -0.195896309465962, -0.752879703028692, -1.20718148517937, 0.103036243241138, 0.707569503031958, 0.488185644365067, -0.674354094380548, 1.45780393551308, -0.426823625045982, 0.372084624745127, -0.489364737649147, 2.00664228219939, 0.0918083578741706, -0.563102471737825, -0.0372326923340063, 1.00352069476112, -1.14244200325742, 1.14469060546702, 1.26070182057624, 0.664616631427959, 0.75021683374034, -1.39029533734557, 0.818844192795068, -0.948663137610259, 2.87434018113975, 0.237475434622456, 0.745884070517981, 1.19494964206119, 0.627396315033843, -0.810728351760724, 1.58246498427057, -0.414333063166778, -1.91414027760821, -1.68038800516161), groups = structure(c(1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L), .Label = c("1", "2", "3", "4"), class = "factor")), .Names = c("x", "y", "groups"), row.names = c(NA, -200L), class = "data.frame") xyplot(y~x|groups, test, type = c("p","smooth"), span = 0.8, col.line="red") span and col.line are not the default values You will have to make your own panel function for locfit if you want to use it I have done it in the past - read the help for library(locfit) ?plot.locfit and the links ?lattice::prepanel Regards Duncan Duncan Mackay Department of Agronomy and Soil Science University of New England Armidale NSW 2351 Email: home: mac...@northnet.com.au -----Original Message----- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of Bond, Stephen Sent: Tuesday, 7 October 2014 23:02 To: r-help@R-project.org Subject: [R] lattice add a fit What is the way to add an arbitrary fit from a model to a lattice conditioning plot ? For example xyplot(v1 ~v2 | v3,data=mydata, panel=function(...){ panel.xyplot(...) panel.loess(...,col.line="red") } ) Will add a loess smoother. Instead, I want to put a fit from lm (but not a simple straight line) and the fit has to be done for each panel separately, not one fit for the full data set, so sth like an lm equivalent of panel.locfit (there is no panel.lmfit) Thank you. Stephen B [[alternative HTML version deleted]] ______________________________________________ R-help@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. ______________________________________________ R-help@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.