similar to: Problem with sapa package and spectral density function (SDF)

Hi, I want to apply a function (in my case SDF; package ?sapa?) repeatedly over discrete sections of a daily time series object by sliding a time window of constant length (e.g. 10 consecutive years or 1825 days) over the entire ts at increments of 1 time unit (e.g. 1 year or 365 days). So for example, the first SDF would be calculated for the daily values of my variable recorded between years 1

Hello, I'm new to R. In Sapa package there's an example for function SDF like data <- as.numeric(sunspots) methods <- c("direct","wosa","multitaper","lag window") S <- lapply(methods, function(x, data) SDF(data, method=x), data) x <- attr(S[[1]], "frequency")[-1] y <- lapply(S,function(x) decibel(as.vector(x)[-1])) which

Hello R-sig-debian and (hopefully) Dirk: On Debian wheezy, I have the R packaging that CRAN (you) provide. I run into a little trouble while trying to fiddle with alternative BLAS. I know you and I went around on this last year and I think perhaps I've found something wrong in the framework, or I've just done something wrong. I installed the packages openblas-base and openblas-dev, and

Hello, I am new to R and am having difficulty formatting numerical output from a regression analysis. My code iteratively performs linear regression on a dataset while excluding certain data ranges. My code: rm(list = ls(all = TRUE)) sink("outfile") dat <- read.table("testdat", sep="\t", header=TRUE) int = 0.2 for (x in c(0:20)) { subdat <- subset(dat, time

I suspect Prof. Ripley's response suffices. However, there *is* a Parzen kernel for kernel smoothing: Parzen K(z) = { 4/3 - 8z^2 + 8|z|^3 if |z| <= 1/2 8(1 -|z|)^3/3 if 1/2 < |z| <= 1 0 otherwise If I'm not mistaken, this appeared in Parzen's original 1962 paper on kernel density estimation. I also seem to recall

I am getting a difference in results when running some analysis using by and tapply compare to using a for loop. I've tried searching the web but had no luck with the keywords I used. I've attached a simple example below to illustrates my problem. I get a difference in the mean of yvar, diff and the p-value using tapply & by compared to a for loop. I cannot see what I am doing wrong.

Dear R-user and -programmer, has one R-package the ability to compute smoothed periodograms of time series using the Tukey-window and/or the Parzen-window? In the ts- and tseries-packages I have found only Daniell-smoothers. With many thanks in advance for any hint Albrecht Kauffmann -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read

Dear All, I am looking for: A software to select the lag length for a lag window spectral estimator. Also, I have a small query in the reprex given below. Background for the above, from the book by Percival and Walden: 1. We are given X_1,...,X_n which is one realization of a stochastic process. 2. We may compute the periodogram using FFT, for example by the function spectrum in R. 3. The

Dear all, What is the usage of the ".Call()"? What is the meaning of the follows: S=.Call("RS_fractal_spectral_density_function_direct",x,as.vector(taper.),as.logical(center),as.logical(recenter),TRUE,as.integer(npad),COPY=rep(FALSE,6),CLASSES=c(rep("matrix",2),rep("logical",3),"integer"),PACKAGE="ifultools") And is a function

Good afternoon! I need to evaluate the goodness-of-fit (aka calibration) for survival probability estimates from a Cox model. I tried to use 'calibrate' in the Design package but I'm not sure if it should/would produce what I need (ie a chi-sq type statistic with a table of expected vs observed probabilities). Any other functions I should be aware of? Also, has anybody come across

I thought this would be trivial, but I can't find a package or function that does this. I'm hoping someone can guide me to one. Imagine a simple case with two survival curves (e.g. treatment & control). I just want to calculate the difference in KM estimates at a specific time point (e.g. 1 year) plus the estimate's 95% CI. The former is straightforward, but the estimates not

Hi all, First I dont have much experience with R so be gentle. OK, I am dealing with a dataset (~ tens of thousand lines, each line ~ 10 columns of data). I have to create some subset of this data based on some certain conditions (for example, same first column with another dataset etc...). Here is how I did it: # import data dat <- read.table( "test.txt", header=TRUE,

ActiveRecord::Base.logger = Logger.new(STDOUT) -- View this message in context: http://old.nabble.com/SQL-Queries-in-Rails-Console-%3Ca-href%3D%22sdf%22%3E%3C-a%3E%3Cobject%3E%3C-object%3E-tp30901356p30901356.html Sent from the RubyOnRails Users mailing list archive at Nabble.com. -- You received this message because you are subscribed to the Google Groups "Ruby on Rails: Talk"

When I try to use the wineoss.drv sound driver, Wine hangs up my whole system (as described by me in Wine Bug #2579). From what I have learned, this is caused by a kernel bug, but I can't fix it as I only found Andreas Mohr's patch for SB16-compatible cards (http://lkml.org/lkml/2002/6/16/45), while I'm using a Yamaha sound card (ymfpci.c kernel driver). When I disable the sound in the

Een ingesloten tekst met niet-gespecificeerde tekenset is van het bericht gescrubt ... Naam: niet beschikbaar Url: https://stat.ethz.ch/pipermail/r-help/attachments/20061106/756382a8/attachment.pl

Hello, I am trying to reformat some data so that it is organized by group in the columns. The data currently looks like this: group X3.Hydroxybutyrate X3.Hydroxyisovalerate ADP 347 4 4e-04 3e-04 5e-04 353 3 5e-04 3e-04 6e-04 359 4 4e-04 3e-04

I've been working with R for all of about 8 hours, so anyone with more experience than this should be able to help me. General comments about my methods of work are also welcomed. I have a table that I've imported thusly: > w <- read.table("woodford.data", header=T) > w start therms gas KWHs elect temp days 1 10-Jun-98 9 16.84 613 63.80 75 40 2

Dear colleagues, I have the following code. This code is to 'filter' the data set. It works on the data frame 'whole' with four numeric columns: a,b,d, and c. Every row in the data frame is considered as a point in 3-D space. Variables a,b, and d are the point's coordinates, and c is its value. This code looks at every point, builds a cube 'centered' at this

sent this to the list yesterday but didn't see it listed in the daily summary ... apologies if you receive it twice ... ________________________________ From: Afshartous, David Sent: Tuesday, June 27, 2006 10:02 AM To: 'r-help@stat.math.ethz.ch' Subject: sapply question All: I'm trying to use sapply to break up data within another function. (tapply doens't seem to work

Hello there. I use Online Xanto ups'es so I did some googling for their support in nut. Turns out that... - networkupstools.org talks about an 'optiups' driver in the development version (http://www.networkupstools.org/compat/dev.html - by the way, the web page says 'stable' in the title and headings just like the actual stable one) - Andreas Thienemann seems to have a