similar to: Loess Fit

Aloha all, I love using ggplot. It took a while to get used to the grammar of graphics, but it is starting to get easy now that I am thinking in a more structured way. A question. I'm making a wind rose that I'd like to be oriented with due north straight up. I've discovered that the orientation is sensitive to how north is represented. When north is represented as 0,

Given two numeric vectors of possibly unequal length, I'd like to pair each element of the shorter vector with an element of the longer vector such that the sum of squared differences between the pairs is minimized. Can someone point me to an R function or an algorithm for accomplishing this? All the best, Tom Thomas S. Dye, Ph.D. T. S. Dye & Colleagues, Archaeologists, Inc.

Aloha all, I have an adjacency matrix for an acyclic digraph that contains transitive relations, e.g. (u,v), (v,w), (u,w). I want a DAG with only intransitive relations. Can someone point me to an R function that will take my adjacency matrix and give me back one with only intransitive relations? In the example, I'd like to get rid of (u,w) and keep (u,v) and (v,w). All the best, Tom --

Aloha all, I have a data frame with 4 columns. The first three are factors (f1, f2, f3) and the fourth is numeric. I'd like to explore these data using median polish. To do that I plan to use medpolish() on the matrix[f1,f2xf3], then medpolish on the resulting matrix[f2,f3]. This approach is described by Cook on page 141 of Exploring Data Tables, Trends, and Shapes. split()

Hi! How can I fit a loess curve to an array (384 x 2). How can I obtain the equation for thi fi? Thanks in advance. -- View this message in context: http://r.789695.n4.nabble.com/How-to-it-a-loess-curve-and-obtain-the-equation-in-R-tp3927622p3927622.html Sent from the R help mailing list archive at Nabble.com.

This message is in MIME format. The first part should be readable text, while the remaining parts are likely unreadable without MIME-aware tools. Send mail to mime@docserver.cac.washington.edu for more info. --1465796871-645734172-974620647=:18687 Content-Type: TEXT/PLAIN; charset=US-ASCII Hi all, I've discovered and (hopefully) fixed a bug in the stars function I wrote some time ago.

Hi, I wonder why my attempt to extend an existing loess fit to a new data set is producing error. I was trying the following: dat = read.csv(choose.files()) x = dat[,2]; y = dat[,1] x.sort = sort(x) y.loess = loess(y~x, span=0.75) # For testing the above fit with a new dataset: test = read.csv(choose.files()) # test data new_x = test [,1]; new_y = test[,2] new_x.sort = sort(new_x) predicted

May I inquire about how to obtain the coefficients from a loess fit of multiple covariates? I have tried terms and coef, but it does not provide me with the coefficients to identify the terms of the entire regression equation. Many thanks Paul -- Paul R. Hutson, Pharm.D. Associate Professor UW School of Pharmacy 777 Highland Avenue Madison WI 53705-2222

I'm at my wit's end, and have searched all of my sources. I need to generate a relatively large number of individual LOESS fits each month of data (I have about 16 months of data). Fitting the polynomial is not my problem, figuring out what the formula that describes that polynomial is. apologies in advance if this is so simple, but I need a hand here. Thanks. -- View this message in

Hi, I'm trying to normalize data by fitting a line through the highest density of points (in a 2D plot). In other words, if you visualize the data as a density plot, the fit I'm trying to achieve is the line that goes through the "crest" of the mountain. This is similar yet different to what LOESS does. I've been using loess before, but it does not exactly that as it takes

I'm updating the loess routines to allow for, among other things, arbitrary local polynomial degree and number of predictors. For now, I've given the updated package its own namespace. The trouble is, panel.loess still calls the original code in package:stats instead of the new loess package, regardless of whether package:loess or package:lattice comes first in the search list. If I

Dear R community, I am trying to understand how the predict function, specifically, the predict.loess function works. I understand that the loess function calculates regression parameters at each data point in 'data'. lo <- loess ( y~x, data) p <- predict (lo, newdata) I understand that the predict function predicts values for 'newdata' according to the loess regression

For some reason the following code is not plotting as I want it to. I want to plot a "loess" line plotted over a scatter plot. I get a jumble, with lines connecting all the points. I had a similar problem with "lowess". I solved that by dropping "NA" rows from the data columns. Please help. library(stats) attach(gini_pci_wdi_narm) plot(ny_gnp_pcap_pp_kd, si_pov_gini)

Dear R-community, maybe someone can help me with this: I've been using the loess() smoother for quite a while now, and for the matter of documentation I'd like to resolve the acronym LOESS. Unfortunately there's no explanation in the help file, and I didn't get anything convincing from google either. I know that the predecessor LOWESS stands for "Locally Weighted

Thanks, John. However, loess.smooth() is producing a very different curve compared to the one that results from applying predict() on a loess(). I am guessing they are using different defaults. Correct? On Thu, 23 Mar 2023 at 20:20, John Fox <jfox at mcmaster.ca> wrote: > Dear Anupam Tyagi, > > You didn't include your data, so it's not possible to see exactly what >

I have a question about running an optimization function on an existing LOESS function defined in R. I have a very large dataset (1 million observations) and have run a LOESS regression. Now, I want to run a Newton-Raphson optimization to determine the point at which the slope change is the greatest. I am relatively new to R and have tried several permutations of the maxNR and nlm functions with

Dear Anupam Tyagi, You didn't include your data, so it's not possible to see exactly what happened, but I think that you misunderstand the object that loess() returns. It returns a "loess" object with several components, including the original data in x and y. So if pass the object to lines(), you'll simply connect the points, and if x isn't sorted, the points

If a run a LOESS model and then produce a smoothed surface: Is there any way to determine the coordinates of the local maxima on the surface? [[alternative HTML version deleted]]

Hello Masters, I run the loess() function to obtain local weighted regressions, given lowess() can't handle NAs, but I don't improve significantly my situation......, actually loess() performance leave me much puzzled.... I attach my easy experiment below #------SCRIPT---------------------------------------------- #I explore the functionalities of lowess() & loess() #because I have

Gurus, I used predict(loess(Y~X)); where Y and X are of the same length. But there are same NA's in both Y and X. Those NA's are in the same locations in Y and X. The following is the error messageI got: Error in "[<-"(*tmp*, , i, value = predict(loess(Y~X))) : number of items to replace is not a multiple of replacement length If I replaced the NA's with a