Thanks a lot Bert! Are z, y and x are of the same length?
And after the loess I just "reshape" the "stacked" outputs back into a 2d matrix? Thanks again! On Wed, Jan 11, 2012 at 5:33 PM, Bert Gunter <gunter.ber...@gene.com> wrote: > You have to set the data up as 3 columns: > > z = Values > y = row > x = column > > (This is trivial to do) > > and then use: > > fit.lo <- loess(z ~ x + y,...) > > For the other stuff, check references and documentation. If this is > too difficult for you, get local help or look elsewhere (e.g. > Bioconductor) for something that better fits your background and > amount of effort your willing to make. > > -- Bert > > On Wed, Jan 11, 2012 at 3:26 PM, Michael <comtech....@gmail.com> wrote: > > I looked the "loess" up... > > > > But I am confused... what should be my "formula" if I just want to smooth > > the numbers of my matrix x? > > And those lots of parameters are daunting... what are the essential ones > > that I should definitely set up? > > I remember for 2d Gaussian filter, I just need to set one parameter > which is > > a scaling lambda.... > > > > Could you please help me? > > > > Thanks a lot! > > > > > > loess(formula, data, weights, subset, na.action, model = FALSE, > > span = 0.75, enp.target, degree = 2, > > parametric = FALSE, drop.square = FALSE, normalize = TRUE, > > family = c("gaussian", "symmetric"), > > method = c("loess", "model.frame"), > > control = loess.control(...), ...) > > > > > > > > > > > > > > On Wed, Jan 11, 2012 at 5:18 PM, Bert Gunter <gunter.ber...@gene.com> > wrote: > >> > >> ?loess ## in base R > >> does 2d filtering , though it's not a Gaussian one > >> > >> -- Bert > >> > >> On Wed, Jan 11, 2012 at 3:00 PM, Rolf Turner <rolf.tur...@xtra.co.nz> > >> wrote: > >> > On 12/01/12 09:17, Michael wrote: > >> >> > >> >> Hi all, > >> >> > >> >> I am looking for a command for doing 2D filtering (rectangular or > >> >> Gaussian) in R... > >> >> > >> >> I have looked at ksmooth, filter and convolve but they seem to be > 1D... > >> >> > >> >> Any thoughts? > >> >> > >> >> Thanks a lot! > >> > > >> > > >> > I'm not sure what you really need/want to do, but you might have a > look > >> > at "smooth.ppp" (and possibly "density.ppp") from the "spatstat" > >> > package. > >> > > >> > cheers, > >> > > >> > Rolf Turner > >> > > >> > ______________________________________________ > >> > 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<http://www.r-project.org/posting-guide.html> > >> > and provide commented, minimal, self-contained, reproducible code. > >> > >> > >> > >> -- > >> > >> Bert Gunter > >> Genentech Nonclinical Biostatistics > >> > >> Internal Contact Info: > >> Phone: 467-7374 > >> Website: > >> > >> > http://pharmadevelopment.roche.com/index/pdb/pdb-functional-groups/pdb-biostatistics/pdb-ncb-home.htm > > > > > > > > -- > > Bert Gunter > Genentech Nonclinical Biostatistics > > Internal Contact Info: > Phone: 467-7374 > Website: > > http://pharmadevelopment.roche.com/index/pdb/pdb-functional-groups/pdb-biostatistics/pdb-ncb-home.htm > [[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.