similar to: Relating Spectral Density to Chi-Square distribution

Dear list, What on earth is spec.pgram() normalized too? If you would like to skip my proof as to why it's not normed too the mean squared or sum squared amplitude of the discrete function a[], feel free too skip the rest of the message. If it is, but you know why it's not exact in spec.pgram() when it should be, skip the rest of this message. The issue I refer herein refers only too a

Dear Tom, Short answer, if your using spec.pgram(), use the smoothing kernel to get a better estimate at the frequency centered in the bandwidth. If your frequency bin of interest is wider than the bandwidth of the kernel, average across frequencies (I think). The estimate appears to be normalized already. If you are calculating your PSD independently, then oversample (e.g. 2, perhaps 4 or more

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

Hi everyone Can anyone tell me how I normalise a power spectral density (PSD) plot of a periodical time-series. At present I get the graphical output of spectrum VS frequency. What I want to acheive is period VS spectrum? Are these the same things but the x-axis scale needs transformed ? Any help would be greatly appreciated Tom

Dear Signal Processing Expert, I would like to generate a random stationary signal of gaussian probability density function to simulate narrow band noise at the output of an IF amplifier. I know the receiver's system temperature (Ts) and IF bandwidth (B) therefore I assume that my narrow band noise mean power equals KTsB watts and therefore the power spectral density No=KTs per Hz. Do you

Hello to R uers, I am wondering if there is an easy way to perform a cross-power spectral density estimation of ?two timeseries (x and y) using the Welch's method. Both packages "bspec" and "oce" provide a function to calculate the PSD with the Welch's method, but only for a timeserie. Thank you in advance. Regards, Pascal

I have done a fair bit of spectral analysis, and hadn't finished collecting my thoughts for a reply, so hadn't replied yet. What exactly do you mean by normalize? I have not used the functons periodogram or spectrum, however from the description for periodogram it appears that it returns the spectral density, which is already normalized by frequency, so you don't have to worry about

Can someone explain me this spec.pgram effect? Code: period.6<-c(0,0,0,0,0,10,0,0,0,0,0,10,0,0,0,0,0,10,0,0,0,0,0,10,0,0,0,0,0,10 ,0,0,0,0,0,10,0,0,0,0,0,10,0,0,0,0,0,10,0,0,0,0,0,10,0,0,0,0,0,10) period.5<-c(0,0,0,0,0,10,0,0,0,0,10,0,0,0,0,0,0,10,0,0,0,0,10,0,0,0,0,0,0,10 ,0,0,0,0,10,0,0,0,0,0,0,10,0,0,0,0,10,0,0,0,0,0,0,10,0,0,0,0,10,0) par(mfrow=c(2,1))

Hi to all I found an formular of an ** ***p-Value Calculator for the Chi-Square test* *http://www.danielsoper.com/statcalc/calc11.aspx* *with the formula* *http://www.danielsoper.com/statkb/topic11.aspx* *what's the gamma function of this formula in r?* *df=5* *ch2=25.50878* *the following code does not give the result <0.001 for the values above * *p=

Full_Name: Madeleine Yeh Version: 1.9.1 OS: AIX 5.2 Submission from: (NULL) (151.121.225.1) After compiling R-1.9.1 on AIX 5.2 using the IBM cc compiler, I ran the checks. One of them failed. Here is the output from running the check solo. root@svweb:/fsapps/test/build/R/1.9.1/R-1.9.1/tests/Examples: ># ../../bin/R --vanilla < stats-Ex.R R : Copyright 2004, The R

Hello world, I am actually trying to transfer a lecture from Statistica to R and I ran into problems with spectral analysis, I think I just don't get it 8-( (The posting from "FFT, frequs, magnitudes, phases" from 2005 did not enlighten me) As a starter for the students I have a 10year data set of air temperature with daily values and I try to get a periodogram where the annual

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)

Hi, How do we get the value of a chi square as we usually look up on the table on our text book? i.e. Chi-square(0.01, df=8), the text book table gives 20.090 > dchisq(0.01, df=8) [1] 1.036471e-08 > pchisq(0.01, df=8) [1] 2.593772e-11 > qchisq(0.01, df=8) [1] 1.646497 > nono of them give me 20.090 Thanks, cruz

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, pchisq -> distribution function dchisq -> density function pval is the area under the curve, to calculte it you use distribution function which is the integral of density function. See: http://www.itl.nist.gov/div898/handbook/eda/section3/eda362.htm http://mathworld.wolfram.com/DistributionFunction.html f(x) density function F(x) distribution function =Pr(X<x)= integral(f(x))

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

Dear list, For educational purposes I have been working with the script below. I have a observation: line 31 #CScal[i] = (amo^2) # IT IS WRONG, I KNOW, BUT IT MAKE R TO CRASHES! I'm thinking this is a possible bug in the R! And I have a couple of doubts: 1) line 51 #curve(dchisq(x, n-1), add = T, col = 'red' I think that it is correct, but the function no

Hi There, There are 300 boy students and 100 girl students in a class. One interesting question is whether boy is smarter than girl or not. first given the exam with a difficulty level 1, the number of the student who got A is below 31 for boy, 10 for girl. Then we increase the difficulty level of the exam to level 2, the number of the student who got A is below 32 for boy, 10 for girl. We

Dear R users I have a really simple question (hoping for a really simple answer :-): Having estimated the spectral density of a time series "x" (heart rate data) with: x.pgram <- spectrum(x,method="pgram") I would like to compute the power in a specific energy band. Assuming that frequency(x)=4 (Hz), and that I am interested in the band between f1 and f2, is the

Hi, I am trying to test for pairwise associations between genotypes ( Rows=individuals, Columns =genes, data are up to 4 genotypes per gene, some with 2,3 or 4) where each chisquare comparison is different depending on the genes tested. The test is the observed multilocus (across columns for each individual) genotypes vs the expectation, which is the product of the individual frequency for each