similar to: Release of new version of caMassClass package and new package caT ools

I am pretty pleased with baselines I fit to chromatograms using the runquantile() function in caTools(v1.6) when its probs parameter is set to 0.2 and its k parameter to ~1/20th of n (e.g., k ~ 225 for n ~ 4500, where n is time series length). This ignores occasional low- side outliers, and, after baseline subtraction, I can re-adjust any negative values to zero. But runquantile's

Hi, Lately I run into a problem that my code R code is spending hours performing simple moving window statistical operations. As a result I did searched archives for alternative (faster) ways of performing: mean, max, median and mad operation over moving window (size 81) on a vector with about 30K points. And performed some timing for several ways that were suggested, and few ways I come up

Hi All, I did a google search and could not find the answer. Thus I post this message. I found runmax only work for positive numbers. x = rep(-1,10) runmax(x,3) x = rep(0, 10) runmax(x,3) for 32-bit R, i got very small numbers: 2.121996e-314 for 64-bit R, i got NaN. Is it a bug in runmax? or should this be included in the documentation that runmax only for positive numbers? however,

Hi, First, let me thank Jaroslaw for making this survey. I find it quite illuminating. Now the questions: * the #1 solution below (based on cumsum) is numerically unstable. Specifically if you do the runmean on a positive vector you can easily get negative numbers due to rounding errors. Does anyone see a modification which is free of this deficiency? * is it possible to optimize the

Hi Everyone, So I've been working a little with mzxml files generated from proteomics data. The current packages available to read these are not sufficient for what I am doing and I would like to build my own. The trouble I am having is how do I read in an mzxml file without the available packages (like reading in a table). I've been searching through Google for a couple of days and am

Hi, Consider the following code: require(XML) xmlFile = paste( "<?xml version=\"1.0\" encoding=\"ISO-8859-1\"?>\n", "<mzXML xmlns=\"a\" xmlns:xsi=\"b\" xsi:schemaLocation=\"c\">\n", "<parentFile a=\"a\" b=\"b\" />\n",

I am getting following error when running this script :------------Error in length(runsum) : object 'runsum' not found ##error function dat$fit <- NULL dat$fit[dat$majorlandmarks] <- dat$Close[dat$majorlandmarks] run <- rle(dat$majorlandmarks) runvalue <- run$values runsum <- cumsum(run$lengths) run <- run$lengths for(i in 1:(length(runsum)-1)){ if

Hi, I am hoping some one can help me. I am learning to use C and would like to learn how to call c code in R. I have look at Writing R Extensions and I tried to copy the example on page 38 void convolve(double *a, int *na, double *b, int *nb, double *ab) { int i, j, nab = *na + *nb - 1; for(i = 0; i < nab; i++) ab[i] = 0.0; for(i = 0; i < *na; i++) for(j = 0; j <

Dear Mr. Tuszynski, I would like to request what I believe would be a beneficial update / patch to the runmean() function in the caTools package. Consider the following R>> x = 1:100 R>> is.integer(x) [1] TRUE R>> library(caTools) R>> head(runmean(x, 5, alg="exact")) [1] 8.487983e-314 1.060998e-313 1.273197e-313 1.697597e-313 2.121996e-313 2.546395e-313

Dear Dr Hornik / R help list, I am using the mda package and in particular the fda routine to classify a set of 162 neurons falling in to 11 neuronal cell types according to 17 morphological variables. I would like to use a cross-validation approach in which I split the data, train with one part amd then test the predictive accuracy of the discriminant functions with the remaining part. However

Hi, Lately my Latex engine used for checking packages stop working, and I can not figure out how to fix it. When creating my package (caMassClass) I get the * checking caMassClass-manual.tex ... ERROR LaTeX errors when creating DVI version. This typically indicates Rd problems. message. I look in caMassClass-manual.log file and find ! LaTeX Error: File `Rd.sty' not found.

Dear all, you can find below my solution for sliding a window. Please find below the code for the two alternatives and the benchmarks.     install.packages('caTools') require(caTools) do_sliding_for_a_window_duty_cycle <- function(DataToAnalyse,  windowSize) {   data<-DataToAnalyse   out <- numeric()   elements<- numeric()   if (length(data[,1]) >= windowSize){       for

R gurus, On Aug 20, 2002, I asked in R-help about calculating a running 5-day median on a large matrix. Thanks to Martin Maechler <maechler@stat.math.ethz.ch> and Ray Brownrigg <Ray.Brownrigg@mcs.vuw.ac.nz> for responding. I ended up writing C code (and an R interface) to do it, which is about 1000x faster than the naive method! (72s became .09s on a 223 x 520 matrix). I added a

Hi I am trying to define a function fu() in the following way but when I try to run I get the error that ma1 is not found. I am not sure where I am going wrong? Does the scope of ma1 not extend to an expr.frame object? expr.frame() is under library tradesys. function (y,ma1,ma2) { x <- y[, c("Open","Close")] d <- expr.frame(x, list(MAf=quote(SMA(Close, ma1)),

Hi, My package "caMassClass" depends on several other packages, one of them "PROcess" residing on "Bioconductor" website. Bioconductor repository is not, listed in default set of repositories in current version of R (Windows R 2.1.1). As a result people installing my package have to change that default, before installing the package. Is there a way of including that

Hi Everyone, I'm currently using the latest build of R and R-Studio server (both are amazing products) I'm still very new to this but I came across this issue: I'm trying to do a select from postgres and put the data into and xts object like so: # Libs library('RPostgreSQL') # http://code.google.com/p/rpostgresql/ library('quantmod') library('TTR')

Dear list, My name is Juan Fernandez-Tajes and I´m trying to use MALDIquant for analysing MS data from ABI 4700 series, I´ve exported a spot set to T2Dextractor in order to convert in mzXML format. I´ve imported this into R with mqReadMzXml function and have followed the suggested flowchart, however when I´m in the step for peak detections I obtained the following error: peaks <-

Hi, I'm very beginner for R but I think it is a time to start as it is very useful. I have a coverage read file (illusmp454merCbed) for whole genome ~ 450 Mbp. This is head of this file. Scaffold sca_position coverage Scaffold1 1 0 Scaffold1 2 0 Scaffold1 3 0 Scaffold1 4 0 Scaffold1 5 0 Scaffold1 6 0 Scaffold1 7 1 Scaffold1 8 3 Scaffold1 9 3 I would like to plot everage coverage for every 1

Hello, I'm new to R and would like to know how to create a vector of "rolling sums". (I have seen the Rmetrics package and the rollMean function and I would like to do the same thing except Sum instead of Mean.) I imagine someone has done this, I just can't find it anywhere. Example: x <- somevector #where x is 'n' entries long #what I would like to do is: x1

Dear all, I have a problem while working with hourly data of fx rates. I've read from a csv file, the following way: csv-file like: Date,Open,High,Low,Close,Volume 2011-08-11 03:00:00,1.41758,1.42205,1.41625,1.42174,8974 ... 2011-08-12 04:00:00,1.42175,1.42413,1.42067,1.42172,7229 ... 2011-12-30 05:00:00,1.42173,1.42341,1.42062,1.42171,6703 ... raw<-