similar to: Any package for deconvolution?

Dear useRs, I have a curve which is a mixture of Gaussian curves (for example UV emission or absorption spectrum). Do you have any suggestions how to implement searching for optimal set of Gaussian peaks to fit the curve? I know that it is very complex problem, but maybe it is a possibility to do it? First supposement is to use a nls() with very large functions, and compare AIC value, but it is

Hi Everyone- After searching through posts and my favorite R-help websites I'm still confused about a problem. I have data which is bimodal in nature, but there is no clearly obvious separation between the two peaks. In programs such as Origin, I can deconvolute the two distributions and have it generate a "best guess" as to what the two subpopulations are which make up my

Hi, I'm trying to match peaks between chromatographic runs. I'm able to match peaks when they are chromatographed with the same method, but not when there are different methods are used and spectra comes in to play. While searching I found the ALS package which should be usefull for my application, but I couldn't figure it out. I made some dummy chroms with R, which mimic my actual

The following commands illustrate a problem with graphing relatively large data sets using the "line" option to plot: temp <- runif(200000) # or whatever other numbers you like; # same thing happens with a sine wave plot(1:200000, temp) # everything fine plot(1:200000, temp, type="l") # data gets cut off The number of points isn't crucial. The cutoff point

I would really appreciate it if someone can give suggestions on how to do spectra fitting in R using ordinary least square fitting and non-negativity constraints. The lm() function works well for ordinary least square fitting, but how to specify non-negativity constraints? It wouldn't make sense if the fitting coefficients coming out as negative in absorption spectra deconvolution. Thanks.

Hi, all. I'm not sure whether this is a problem with ESS or R (or even a feature I don't understand!). Say I've got a function named f.count.these.patterns (as I do -- I'm using actual names in case they're somehow important). I previously had a version called f.count.these.patterns.3, but discarded it because it was worse. So the function is no longer there. If I create a

Hi, all. Does anyone have an equivalent to Splus "peaks", which finds local maxima (with locality defined by a parameter "span")? I thought I'd check whether anyone has done it already before trying to put something together myself. Thanks, Matt Wiener -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read

Hello, all. I've got what is probably a simple question about comparison of models using anova, specifically about the situations in which it's valid. I understand, I think, what's going on when the models are strictly nested (as most are in the demo(lm) examples). My question involves what happens when the models aren't strictly nested. In my particular case, I'm doing

Dear R and sweave users A further problem, which I couldn't resolve, using the manual: In R I use the split.screen command to put e.g. two timecourses one above the other into one plot: split.screen(c(2,1)) screen(1) plot(stick,type='h', col="red",lwd=2) screen(2) plot(deconvolution.amplitude,type='h',col="blue",lwd=2) Is there a similar way, doing this

I'm translating a matlab routine to R and I need some equivalent to deconv(): Description: deconv() [q,r] = deconv(v,u) deconvolves vector u out of vector v, using long division. The quotient is returned in vector q and the remainder in vector r such that v = conv(u,q)+r . If u and v are vectors of polynomial coefficients, convolving them is equivalent to multiplying the two polynomials, and

Hi, all. The following does not behave as I think it should, and as it seems to me it has in the past (although I can't check this easily). I know it happens in both R-0.64.2 and R-0.65.0 on an old Power Computing running Linux-PPC 1999, and in R-0.64.2 on an SGI running Irix 6.5. Try the following: t1 <- data.frame(matrix(rnorm(16), nc=4)) > t1 X1 X2 X3 X4 1 -0.7206945

Hi, I wonder if someone could help. I needed to transfer (copy) a workspace file that had been generated in linux (R 2.11) to windows running the same version of R 2.11 (but of course windows binary). Usually, there is no problem in doing this and all objects work as expected. I am often doing this to be able to produce wmf or emf graphic files that I need. This time I had some spectra that I

Hi, all. Using the patched version of R-0.99.0, I cannot reproduce the following example from the strsplit documentation: (Example) unlist(strsplit("a.b.c", ".")) ## [1] "" "" "" "" "" ## Note that `split' is a regexp! ## If you really want to split on `.', use unlist(strsplit("a.b.c", "\."))

Two more details: Someone else who works here, and hasn't used 0.65.0, thinks he's seen the error before. He restarted R, and it went away. Second, I think the relevant code is in SEXP do_function, in eval.c in the main directory. line 695: WrongArgCount("lambda") Beyond that I'm afraid I'm getting in over my head. Matt

Hi. I'm having a problem with disappearing comments when I transfer files (entire .RData files) from one machine to another, and I'm not sure whether it's an R problem or an ESS problem. Since ESS seems more likely, I'm trying this list first. (However, it seems to me this problem only showed up when I upgraded to R-0.65.0 or R-0.65.1.) Functions originally shows up fine -- if

Hi, all. The help file for tapply says that if simplify is true, and the result of the calculation is always a scalar, then tapply will return a vector. Nonetheless: > t1 <- tapply(runif(10), rep(1:5, 2), mean) > is.vector(t1) [1] FALSE > is.array(t1) [1] TRUE > I have found this in version 0.65.1 on an SGI running Irix 6.5, and on a Mac running Linux-PPC. I've also

Hi, I?m trying to use the STLperArea function in the ndvits package. I can run the sample data (?SLPSAs_full?) without any problem. However, when I come to run the function on my own data, I get the following message. Waiting to confirm page change... Error in .checkPath(path) : no TISEAN executables found in that directory. Please set a proper TISEAN executables path using

Hi, all. In trying to compile R v.0.90 on an alpha running Red Hat Linux 6.0, I'm getting the following error: make[3]: Entering directory `/home/mcw/alpha-R/R-0.90.0/src/appl' gcc -I. -I../include -I../../src/include -DHAVE_CONFIG_H -mieee -g -O2 -c cpoly.c -o cpoly.o /tmp/ccPjZNwd.s: Assembler messages: /tmp/ccPjZNwd.s:312: Error: unknown opcode `sqrttsu' /tmp/ccPjZNwd.s:2631:

Why does R CMD INSTALL work for some packages (e.g., lme4) but not others (e.g., nlme)? Thanks, Gang

Dear All, I am a novice when it comes to time-series analysis and at the moment I am actually interested in calculating the Hurst exponent of a time series. This question has already been asked quite some time ago http://bit.ly/98dZsi and I trust some progress has been made ever since. I was able to find some functions in the packages http://cran.r-project.org/web/packages/Rwave/index.html