similar to: A Question on lowess() function

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

Hi, I have a scatter plot of the variables GNI and Lifeexp (Gross National Income and Life Expectancy, both metric). So I plotted them and I want to add a regression line and a lowess line. I use lowess and not loess because I have missing values. My code: plot(GNI,Lifeexp) abline(lm(Lifeexp~GNI), col="red") y.loess<-loess(Lifeexp~GNI,na. action = na.exclude)

As I've mentioned a number of times. I find it very useful to have lowess() become a generic function so that a lowess.formula() can be defined. Below is a patch that makes both changes, as well as updating the corresponding help documentation. Gregory R. Warnes Manager, Non-Clinical Statistics Pfizer Global Research and Development Tel: 860-715-3536 ? DESCRIPTION ? Makefile ?

Colleagues ---------------------------------- System info: R version rw1022 on NT ESS v. 5.1.18 using emacs ver. 20.4 ---------------------------------- I would like to access the residuals from a lowess fit (function lowess() ) as an aid to assessing whether I am over-smoothing or under-smoothing a dataset (as in Cleveland's book "Visualizing data", section 3.3) Unless I

Dear all, I am a newbie in R. I encounter a problem as follows. I have 2 vectors X and Y that have a equal length of several thousand. I see Y as the function of X. Both of them are random. X is not arrranged in any order. Of course, I do plot(X,Y). Now, I want to use a sliding narrow window to run over each X, then calculate the variances within that window. Anyone knows easy way in R to do

Hi, Can anybody tell me how to open new a new window for plot()? Thanks. Minghua

Hi, Anyone can tell me how to generate "A01", "A02", ..., "A99"? paste("A", 1:99, sep="") generates "A1", "A2",..., "A99". This is not what I want. Thanks for the help. -MY [[alternative HTML version deleted]]

code: x <- rnorm(32) y <- rnorm(32) plot(x,y) lines(lowess(x,y),col='red') Now I need to sample the lowess function into matrix where one series will be X and other will be values of lowess at particular X. -- View this message in context: http://r.789695.n4.nabble.com/how-to-sample-lowess-loess-into-matrix-tp3053458p3053458.html Sent from the R help mailing list archive at

Hi there, I'm plotting some glass RI values just by plotting plot(x) then I put on my lowess smoother lines(lowess(x)) now I want to put on some 95% Confidence Interval bands of the lowess smoother, but don't know how?? Thanks -- Gareth Campbell PhD Candidate The University of Auckland P +649 815 3670 M +6421 256 3511 E gareth.campbell@esr.cri.nz gcam032@gmail.com [[alternative

Dear UseRs: Recently, I have read an article regarding the association between age and lymph node metastases. http://jco.ascopubs.org/content/27/18/2931.long In statistical analysis, the authors stated "Because a nonlinear relationship between age and lymph node involvement was expected based on existing literature, lymph node involvement was also regressed on age using nonparametric

Dear all, please do you have any recommendation about a more advanced function in R for lowess/loess regression ? the basics lowess() or loess() do not perform as well as I would expect. thanks very much, Bogdan [[alternative HTML version deleted]]

Hi Every Body, I don't know why some regression functions have no related prediction function. For example lowess, loess, lokerns, lpridge, lpepa, and ksmooth. What could help? Is there any global or wrapper function so that can help? Regards, Amir Safari --------------------------------- [[alternative HTML version deleted]]

Hi there, I cant seem to find a way for the lowess smoothing function to handle "NA" values. Can anyone help?? Regards, Wayne Dr Wayne R. Jones Statistician / Research Analyst KSS Group plc St James''s Buildings 79 Oxford Street Manchester M1 6SS Tel: +44(0)161 609 4084 Mob: +44(0)7810 523 713 KSS Ltd A division of Knowledge Support Systems Group plc Seventh Floor St

I'm using this code to plot a smoothed line. These two columns of data really represent 4 groups and I'd like to plot a separate line for each group but have them all in the same plot. The R-Docs for lowess do not seem to indicate some type of "GROUPS=var_name" option. What would be the syntax for this? plot(AWGT ~ lipid ) lines(lowess(lipid , AWGT, f=.8)) -- Dean

Hi All, I've been providing a generic version of lowess() in my 'gregmisc' library, along with a 'lowess.formula' and 'lowess.default'. Can/should these be moved into base? -Greg LEGAL NOTICE Unless expressly stated otherwise, this message is confidential and may be privileged. It is intended for the addressee(s) only. Access to this E-mail by anyone else

I am trying to change the attributes of the lowess lines fit to an xyplot command, but have been unsuccessful in my search of the online help. Right now, both the points and lowess line come out in the same color (blue). I am unsure how I can change the properties of the lowess line separately. xyplot ( SnowLineElevation ~ Year | Model, data = data, ylim = c(0,1800), type =

Hi there, Just out of curious, is there any difference between loess() and lowess() in R (and Splus in fact). Which one is more often used? Thanks, Ko-Kang Wang ------------------------------------------------------------------------------ Ko-Kang Kevin Wang Statistical Analysis Division Leader Software Developers' Klub (SDK) University of Auckland New Zealand

Hi folks. If I have the following in my "data" event pH1 pH2 1 4.0 6.0 2 4.3 5.9 3 4.1 6.1 4 4.0 5.9 and on and on..... for about 400 events Is there a way I can get R to plot event vs. pH1 and event vs. pH2 and then do a loess or lowess line for each?? Thanks in advance David [[alternative HTML version deleted]]

Hi, I understand that the core of lowess() function is written in C. Is the C source code available and if yes, how can I get to it? I want to write a multidimensional extension of lowess() if you know of anybody who has implemented this already I would greatly appreciate any references or pointers. thank you very much, pawel

I was trying to plot some data in R. I used the following code to draw a loess fit and got the output as >?lines(lowess(log(abs(t(res))), log(abs(t(synthesised)))), col="red") Error in lowess(log(abs(t(res))), log(abs(t(synthesised)))) :?? NA/NaN/Inf in foreign function call (arg 1) Then I thought to use your Limma package for background correction. Do you think it's a right