search for: taper

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 symbo...

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[, i]'. However, the program actually...

...triangle of pgram is necessary and pgram is not returned by the function, therefore, calculating the lower left portion (ie when j < i) is not required. We change two nested loops to index from i:ncol(x) instead of 1L:ncol(x) : newspec.pgram <- function (x, spans = NULL, kernel = NULL, taper = 0.1, pad = 0, fast = TRUE, demean = FALSE, detrend = TRUE, plot = TRUE, na.action = na.fail, ...) { ## Estimate spectral density from (smoothed) periodogram. series <- deparse(substitute(x)) x <- na.action(as.ts(x)) xfreq <- frequ...

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.

...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 = FALSE , fast = FALSE )$spec, ty = 'l' ) lines(spec.pgram(x, taper = 0 , detrend = FALSE ,plot = FALSE , fast = TRUE )$spec, col = 'red' ) lines(((abs(fft(x))^ 2 )[ 2 : 50 ])/ 99 , col = 'blue' , lty = 2 ) Best r...

...the extremes, values seem to result from kernel application to collections of high and low frequencies. Related to this, shouldn't the frequency range shorten according to the size of kernel used? My doubt came from reading Diggle(1990) p105. e.g. data(lh) x<-spec.pgram(ldeaths, detrend=T, taper=0)$freq y<-spec.pgram(ldeaths, kernel("modified.daniell", c(6,6)), detrend=T, taper=0)$freq x==y Thanks in advance, Nuno Prista

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 works fine, but if I use some taper, e.g. S=SDF(data, method="lag window", tap...

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 the log of the periodogram with per and that with spec.pgram. This vertical shift amounts t...

Great info. I found some descriptions of st5 (md5 fp) and ffp, where I assume that "fp" is just a Taper abbreviation for "fingerprint" - or signature as it's called in the flac header. However, I could not find these utilities or source code. What I found looked like instructions for a gui-based program. I think it would be easier to support Mac if st5 were available as an op...

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.i...

...sub.dat$ymf <- subdat$yf : subdat$mf > sub.datm <- aggregate(subdat$DA,list(subdat$ymf),mean) > > serie.day <- ts(subdat$DA,start(sub.dat$Year[1],sub.dat$Day[1],frequency=365)) > n.d <- length(serie.day) > espec.day <- SDF(serie.day,method="lag window",window=taper(type="parzen",n.sample=n.d,cutoff=(2*sqrt(n.d))),npad=2*n.d) > > serie.month <- ts(sub.datm,start(sub.dat$Year[1],sub.dat$Month[1]),frequency=12) > n.m<-length(serie.month) > n.m [1] 322 > espec.month<-SDF(serie.month,method="lag window",window=taper(type...

...'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 call spec.pgram by spec.pgram (h, log="no", taper=0.5). So, I have some peaks that appear and I would like to interpret them but I do not know how they are computed and so what a peak with a value of 10000 represents in comparison with a peak of value 600 with another histogram. I looked at the source code of the function spec.pgram to better unde...

...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_...

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...

...tes(mfdeaths.spc) $names [1] "freq" "spec" "coh" "phase" "kernel" "df" [7] "bandwidth" "n.used" "orig.n" "series" "snames" "method" [13] "taper" "pad" "detrend" "demean" $class [1] "spec" Thank you so much. Yours Faithfully, -- Maura E.M [[alternative HTML version deleted]]

I realize that it isn't much of an improvement, but AIFF supports 4GB recordings, and flac is compatible with this. Being an avid "taper" myself, I have, on many occasions recorded up to this limit, and I always back up my original recordings using flac. W64 support is more than welcome, but AIFF support gets you twice the length right away. Brian Willoughby Sound Consulting On May 29, 2007, at 23:00, Justin Frankel w...

shntool (http://etree.org/shnutils/, http://shnutils.freeshell.org/) is Windows-only, but the source is available. > Great info. I found some descriptions of st5 (md5 fp) and ffp, where I > assume that "fp" is just a Taper abbreviation for "fingerprint" - or > signature as it's called in the flac header. However, I could not find > these utilities or source code. What I found looked like instructions > for a gui-based program. I think it would be easier to support Mac if > st5 were availa...

...er side of a single data point (thus making an irregularly spaced time series), a spectrum with a very large peak at the highest frequencies is produced. An example of this is shown below: a <- runif(100) x <- 1:length(a) a <- a[-c(4,6)] x <- x[-c(4,6)] spec.ls(x=x, y=a, fast=FALSE, taper=0) Does anyone know why this is? Thank you very much in advance for any help you can give me with this issue. Best regards, Bernardo