search for: tapers

I am trying to load a .so file and get the following error message: > dyn.load("taper.so",local=F) Error in dyn.load(x, as.logical(local), as.logical(now)) : unable to load shared library "/home/hpc1367/runs/taper/taper.so": ld.so.1: /usr/local/lib/R/bin/R.bin: fatal: relocation error: file /home/hpc1367/runs/taper/taper.so: symbol f90_init: referenced symbol not

Full_Name: Michael Stein Version: Version 2.1.1 OS: linux Submission from: (NULL) (128.135.149.112) The documentation for spec.taper says p: The total proportion to be tapered, either a scalar or a vector of the length of the number of series. Details: The cosine-bell taper is applied to the first and last 'p[i]/2' observations of time series 'x[,

Hello, I propose two small changes to spec.pgram to get modest speedup when dealing with input (x) having multiple columns. With plot = FALSE, I commonly see ~10-20% speedup, for a two column input matrix and the speedup increases for more columns with a maximum close to 45%. In the function as it currently exists, only the upper right triangle of pgram is necessary and pgram is not returned by

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

Hi list, Does anyone know if there is a library in R that does MTM-SVD method for spectral analysis? Thanks ----- Yasir H. Kaheil Columbia University -- View this message in context: http://www.nabble.com/Spectral-analysis-with-mtm-svd-Multi-Taper-Method-Combined-with-Singular-Value-Decomposition-tp21671934p21671934.html Sent from the R help mailing list archive at Nabble.com.

Dear all, This is potentially a bug in spec.pgram, when the number of samples is odd,spec.pgramreturns a different result withfast = TRUE, the example below contains the two varieties with a reference spectrum calculated manually. the number of returned spectra is also larger (50 compared to 49) whenfast = TRUE x <- rnorm( 99 ) plot(spec.pgram(x, taper = 0 , detrend = FALSE , plot =

Hi, I am far from experienced in both R and time series hence the question. The code for spec.pgram() seems to involve a circularity of the kernel (see below) yielding new power estimates to all frequencies computed by FFT. " if (!is.null(kernel)) { for (i in 1:ncol(x)) for (j in 1:ncol(x)) pgram[, i, j] <- kernapply(pgram[, i, j], kernel, circular = TRUE)

...he decibel() function expectantly throws error. f=attr(vysl,"frequency")[-1] s=decibel(as.vector(S)[-1]) Error in decibel(as.vector(vysl)[-1]) : cannot convert negative values to decibels Is it possible to use some other function than decibel, or how else can estimate spec.density with tapers from this package? Thank you very much. [[alternative HTML version deleted]]

Linux RH 7.3, Samba 2.2.7a NT 4.0 SP6 - Backup Exec 8.5 The linux box is running the Backup Exec client for Unix ('http://seer.support.veritas.com/docs/192093.htm') Occasionally, when doing a backup the entire linux system freezes. I have to do a hard reboot without a shutdown, etc. There is 40 plus gigs on the backup, so I dont think a tar and copy would be bery efficient. I am

hello, I have been using the per function in package longmemo to obtain a simple raw periodogram. I am considering to switch to the function spec.pgram since I want to be able to do tapering. To compare both I used spec.pgram with the options as suggested in the documentation of per {longmemo} to make them correspond. Now I have found on a variety of examples that there is a shift between

...shn md5) files. But the original poster and I want a command-line solution to create an st5 for comparison against metaflac's report. I think I have xACT, but I don't want to be grabbing the mouse and clicking dialogs when I want to check a whole directory of file signatures. Any Tapers around? Maybe someone can clarify this for us newcomers to st5. Brian On Feb 8, 2008, at 00:01, rappard@dds.nl wrote: > I have a wav file and would like to see what the flac fingerprint > would be. > > To do this I run flac to encode the wav file and write the flac > file t...

I am interested in tapping into LAPACK functionality directly from R. Using R-2.4.1 for Windows, I was able to do so ala: dyn.load("bin/Rlapack.dll") is.loaded("dstebz") # returns TRUE N <- 100 NW <- 4 n.tapers <- 5 tpW <- (2 * pi * NW)/N otNmo <- 1:N D <- as.double(cos(tpW) * ((N - 1 - 2 * (0:(N - 1)))/2)^2) E <- as.double((otNmo * (N - otNmo))/2) z <- .Fortran("dstebz", "I", "B", as.integer(N), double(1), double(1), as.integer(N - n.tapers + 1), as.in...

Hi everybody, I'm a beginner R user and I'm trying to use the package "sapa" to estimate the spectral density function of several time series using the SDF function. For each time series, I want to calculate the density function at two temporal resolutions (daily and monthly). The monthly values calculated as a mean of the daily values. I first create a ts object for both series

Hi everyone, first of all, I would like to say that I am a newbie in R, so I apologize in advance if my questions seem to be too easy for you. Well, I'm looking for periodicity in histograms. I have histograms of certain phenomenons and I'm asking whether a periodicity exists in these data. So, I make a periodogram with the function spec.pgram. For instance, if I have a histogram h, I

...crude test from one guest over nfs, dd if=/dev/zero of=/nfs/test bs=32768 count=32768 This yields almost always 95-100MB/sec When I run two simultaneously, I cannot seem to get above 25MB/sec from each. It starts off with a large burst like each can do 100MB/sec, but then in a couple of seconds, tapers off to the 15-40MB/sec until the dd finishes. Things I have tried (installed on the host and the guests) net.core.rmem_max = 16777216 net.core.wmem_max = 16777216 net.ipv4.tcp_rmem = 4096 87380 16777216 net.ipv4.tcp_wmem = 4096 65536 16777216 net.ipv4.tcp_no_metrics_save = 1 net.ipv4.tcp_m...

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

...format to get around the 2GB .wav file limitation. W64 was invented by Sonic Foundry, and is an open format as far as I know. The only programs I know about using the .w64 format at the moment are Sound Forge and Steinberg Nuendo, although there may be others out there. With increasing number of tapers recording concerts in 24bit, which takes up 50% more disk space than 16bit files at the same sample rate (about 1GB per hour), it is easy to blow past 2GB for a single file, which is roughly 2 hrs of audio. Some audio editors will stop after the 2GB limit is reached, or start recording a new file,...

I am reading some documentation about Cross Spectrum Analysis as a technique to compare spectra. My understanding is that it estimates the correlation strength between quasi-periodic structures embedded in two signals. I believe it may be useful for my signals analysis. I was referred to the R functions that implement this type of analysis. I tried all the examples which generated a series of

....wav file limitation. W64 was invented by Sonic Foundry, and > is an open format as far as I know. The only programs I know about > using the .w64 format at the moment are Sound Forge and Steinberg > Nuendo, although there may be others out there. With increasing > number of tapers recording concerts in 24bit, which takes up 50% > more disk space than 16bit files at the same sample rate (about 1GB > per hour), it is easy to blow past 2GB for a single file, which is > roughly 2 hrs of audio. Some audio editors will stop after the 2GB > limit is reached, or...

...the > original poster and I want a command-line solution to create an st5 for > comparison against metaflac's report. I think I have xACT, but I don't > want to be grabbing the mouse and clicking dialogs when I want to check > a whole directory of file signatures. > > Any Tapers around? Maybe someone can clarify this for us newcomers to st5. > > Brian > > > On Feb 8, 2008, at 00:01, rappard@dds.nl wrote: >> I have a wav file and would like to see what the flac fingerprint would be. >> >> To do this I run flac to encode the wav file and wri...