similar to: Applying user function over a large matrix

Hi there, I would like to draw a scatter plot and fit a smooth line by loess. Below is the data. However, the curve line started from 0, which my "resid" list doesn't consist of 0 value. It returned some warnings which I don't know if this is the reason affecting such problem. Here I also attached the warning messages. Please let me know if there is a solution to fix this. Thank

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

Hello, I found a little BUG in loess <stats>. It does not receive the iterations parameter. It can be debugged in the following way: THIS IS AN EXCERPT FROM THE CODE: .... fit <- simpleLoess(y, x, w, span, degree, parametric, drop.square, normalize, control$statistics, control$surface, control$cell, iterations, control$trace.hat) Replace argument iterations with

Hi mailinglist members, I’m actually working on a time series prediction and my current approach is to decompose the series first into a trend, a seasonal component and a remainder. Therefore I’m using the stl() function. But I’m wondering how to get the single components in order to predict the particular fitted series’. This code snippet illustrates my problem: series <-

Hi all, I am doing a projection pursuit regression using the ppr() function from modreg. I would also like to use predict.ppr(). However, I cannot find any information about it in the help files. There is a link to predict.ppr in the index for modreg, but that link is to the help for ppr(). Has predict.ppr() not been implemented? If not, does anyone have a suggestion as to how to implement

Hi all, I am just learning R and I am trying to work through the book "Applied Longitudinal Data Analysis" by Singer & Willett. I have some code for this book that supposedly works in S-Plus (I don't have S-Plus so I can't verify that) and I am trying to run the examples in R. Most of the examples run, but I have one plot that gives me an error message. I have

Dear R experts, I have a problem that is a related to the question raised in this earlier post https://stat.ethz.ch/pipermail/r-help/2007-January/124064.html My situation is different in that I have only 2 predictors (coordinates x,y) for local regression but a number of global ("parametric") offsets that I need to consider. Essentially, I have a spatial distortion overlaid over a

hi! i want to change the "robustness weights" used by loess. these are described on page 316 of chambers and hastie's "statistical models in S" book as r_i = B(e_i,6m) where B is tukey's biweight function, e_i are the residulas, and m is the median average distance from 0 of the residuals. i want to change 6m to, say, 3m. is there a way to do this? i cant

Dear List, Having had my appetite sufficiently whetted by Prof. Ripley's email about the new graphics capabilities in Unixes, I wanted to try them out. I updated to svn r44608, configured with the following options: R is now configured for x86_64-unknown-linux-gnu Source directory: .. Installation directory: /usr/local C compiler: gcc -O3 -g -std=gnu99

Skipped content of type multipart/alternative-------------- next part -------------- A non-text attachment was scrubbed... Name: Blank Bkgrd.gif Type: image/gif Size: 145 bytes Desc: not available Url : http://lists.samba.org/archive/samba/attachments/20020911/a89e1b74/BlankBkgrd.gif

I noticed a problem using R 2.7.1 on Windows XP SP2 with the precompiled Atlas Rblas.dll. Running the code below causes R to crash. I started R using Rgui --vanilla and am using the precompiled Atlas Rblas.dll from cran.fhcrc.org dated 17-Jul-2007 05:04 for Core2 Duo. The code that causes the crash: x <- rnorm(100) y <- rnorm(100) z <- rnorm(100) loess(z ~ x * y) loess(z ~ x) does

1. Minimum Mean Square Estimation (for Multi-Dimensional Linear Regression) 2. Two-Sample Kolmogorov-Smirnov Tests 3. Empirical ANOVA Methods All this to be done on empirical data. I am new to the R world, any help would be appreciated. Thank You, Bill Gilliland Lead Sensor Systems Engineer E303, Tactical Communications, Transport, and Sensors Group: Tactical

Dear R-Helpers, I suspect I'm about to ask a FAQ, but I haven't been able to find an answer in the FAQ, AItR or an R Site Search. When I look at the methods of summary (below) it says, "Non-visible functions are asterisked". I looked at the help file for summary.princomp, which did not comment on it being non-visible. I ran its help file example, which printed visible output. I

Hi, I'm not sure what the recommended thing to do here would be: In the lattice library, panel.loess needs to use loess.smooth(), which is in the modreg library. Now should I (1) make lattice depend on modreg, and call require(modreg) in zzz.R OR (2) use autoload("loess", "modreg") I like the second option better, as a particular session might not need panel.loess.

Hello, I have been struggling to do a plot in ggplot(2) that's of lattice equivalent. The following code shows the lattice plot.

I have created an xyplot of a time series with the following code... win.graph(width = 10, height = 7) panel1 = function(x, y) { panel.loess(x, y, lwd=2.5, span=0.5, col="gray") panel.xyplot(x, y, pch=19, col="blue", cex=1.25) } xyplot(oneplusdensity ~ year, data=figdata, aspect="fill", cex=1.5, xlab=NULL, ylab=expression("Crabs per

Is there one can give me some information on how to setup printer in Samba so windows can print to it. My file sharing works fine. In winows can't see the printer. Thanks Matthew Wu School Loans Corp. 10780 Santa Monica Blvd, Ste 225 Los Angeles, ca 90025 Tel. 310-474-7456 ext.235 Email mwu@higheredloans.com

Dear all; I wish to create a graphic object combing an xyplot() and an mtext(). My code looks like following, gmv <- { trellis.device("windows", bg="white", width = 7, height = 7) xyplot(Mvol ~ Age | Nl * Th , data = Hft1, prepanel = function(x, y) prepanel.loess(x, y, span = 1), xlab =list(label = "Age (Years)", font = 2),

Is it possible to allow two Linux computers to browse each other in a manner similar to the way a winbox and browse a Linux box with Samba? Ya know, like the network neighborhood thing, not the smbclient get put thing? Thanks

Hello, I was wondering wether there's a function in R that takes two vectors (of same length) as input and computes mean values for bins (intervals) or even a sliding window over these vectros. I've several x/y data set (input/response) that I'd like plot together. Say the x-data for one data set goes from -5 to 14 with 12,000 values, then I'd like to bin the x-vector in steps of