similar to: debug

R Users: How can I catch R errors and depend on this error call other R functions in Perl? The foolowing is the error message when I use perl call R. Error in nls(modout ~ exp(-b1 * modin)/(b2 + b3 * modin), trace = F, start = c(b1 = 0.005, : singular gradient Segmentation fault Thanks, Zhongming

Hi, I use lm or loess to make smoothing. After smoothing I need "Residual Standard Error" in my script. Could you please tell me how can I get this information? Thanks,

Hi, I am trying to use piecewise linear regression to approximate a nonlinear function. Actually, I don't know how many linear functions I need, therefore, I want build an array of regression models to automate the approximation job. Could you please give me any clue? Attached is ongoing code: rawData = scan("c:/zyang/mass/data/A01/1.PRN", what=list(numeric(),numeric())); len =

An embedded message was scrubbed... From: "Zhongming Yang" <Zhongming.Yang at cchmc.org> Subject: how to make exprSet Date: Mon, 14 Jul 2003 12:16:33 -0400 Size: 683 Url: https://stat.ethz.ch/pipermail/r-help/attachments/20030714/6e4061d3/attachment.mht

Dear R developers! Considering I got no response or comments in the general r-help forum so far, perhaps my question is actually better suited for this list? I have added some more hopefully relevant technical details to my original post (edited below). Any comments gratefully received! Best regards, David Kreil. ---------- Dear R experts, I have a problem that is a related to the question

Hi, I am segfaulting when using predict.loess() (checked with r62092). I've traced the source with the help of valgrind (output pasted below) and it appears that this is due to int overflow when allocating an int work array in loess_workspace(): liv = 50 + ((int)pow((double)2, (double)D) + 4) * nvmax + 2 * N; where liv is an (global) int. For D=1 (one x variable), this overflows at

Dear All: I need draw some figure through trellis.device and save them as pdf files. How can I close trellis.device (something like dev.off() in nonlattice figure)? Many thanks Zhongming Yang --------------------------------- [[alternative HTML version deleted]]

I have a failure in 'make check' and I would appreciate assistance regarding what to do about it. Thanks -Don Installing R-1.2.3 on build[99]% uname -sr SunOS 5.8 (for those who may recall my previous request for help, I have updated GCC to 2.95.3, and configured with gcc, GCC's c++, and f77 [not g77]) Running 'make check' got this: creating `no-segfault.R'

I would like to experiment with testing the fit of a loess model against the fit from an ordinary linear regression. The 1988 JASA paper by Cleveland and Devlin *appears* to indicate that this can be done, at least ``approximately''. They, as I read it, advocate the use of an ANOVA type test with degree of freedom chosen to make the ``F ratio'' have an approximate F distribution

Hi Folks, I'm looking to do Confidence bands around LOESS smoothing curve. If found the older post about using the Standard error to approximate it https://stat.ethz.ch/pipermail/r-help/2008-August/170011.html Also found this one http://www.r-bloggers.com/sab-r-metrics-basics-of-loess-regression/ But they both seem to be approximations of confidence intervals and I was wonder if there was

Hi all, I have a vector, most of which is missing. The data is always increasing, but may do so in jumps. I would like to interpolate the NAs with 'best guesses', using something like filter(), which doesn't work due to the NAs. Here is an example: > x <- c(2,3,NA,NA,NA,3.2,3.5,NA,NA,6,NA) > x [1] 2.0 3.0 NA NA NA 3.2 3.5 NA NA 6.0 NA I would like a function that

I am 100% certain that there is an easy way to do this, but after experimenting off and on for a couple of days, and searching everywhere I could think of, I haven't been able to find the trick. I have this piece of code: ... attach(d) if (ORDINATE == 'ds') { lo <- loess(percent ~ ncms * ds, d, control=loess.control(trace.hat = 'approximate')) grid <-

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

Dear All: In the attached file, I have 3 group patients, and there are 5 in each group (the groups are decided by the prefix of the idno). I want draw a repeat measurement comparison figure. My goal is to list 5 patients from same group on one horizontal line. But xyplot sounds pick them randomly (or I was confused?). Could you please help me modify the following code to accomplish this?

Hi, R-devel causes me still trouble. It dies with Error: couldn't find function "predict.loess" Execution halted for the package affy, but bash-2.04$ grep predict.loess ~/R-devel/src/library/modreg/R/*.R /misc/homes/laurent/R-devel/src/library/modreg/R/loess.R:predict.loess <- function(object, newdata = NULL, se = FALSE, ...) bash-2.04$ ..the function seems to be there...

I've a situation where I need to perform different regressions on different subsets of my data and then combine the resulting fitted values back into the original dataset (partly this is because loess doesn't do grouping, but I'd also like to approximate the shrinkage estimators of mixed models by giving low weights to observations in the groups other than the one I'm

Dear All: Could you please tell me how I can draw figure formatted like figure 4.18 of MASS (4th) with the attached data set? Thanks Zhongming Yang --------------------------------- -------------- next part -------------- An embedded and charset-unspecified text was scrubbed... Name: sample.txt Url: https://stat.ethz.ch/pipermail/r-help/attachments/20050316/abfdb85e/sample.txt

Hello, I was hoping that someone familiar with the implementation details of the loess algorithm might be able to help me resolve some difficulties I am having. I am attempting to reproduce some of the functionality of the loess() function in C++. My primary motivation is that I would like to understand the algorithm in detail. So far I have managed to create a working port in C++ for the

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

I ran the following script from xyplot Examples using Tin-R on Windows and saw no plot produced. EE <- equal.count(ethanol$E, number=9, overlap=1/4) xyplot(NOx ~ C | EE, data=ethanol, prepanel = function(x,y) prepanel.loess(x, y, span=1), xlab="Compression Ratio", ylab="NOx (micrograms/J)", panel = function(x,y) { panel.grid()(h = -1, v=2)