similar to: Interpolating spline problems and akima

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

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

Hi, I am hired to translate into FORTRAN R functions. It works great and my customer is happy. I am going to translate a function that call aspline in the library AKIMA. This function call 2 FORTRAN routines that must be include (I guess) in the package AKIMA. And I have some questions. Is the fortran code really included in the R package AKIMA ? In other word, do I have a mean to call these 2

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

I have the following function defined as below match.trace <- function(dfobj, distance, day1, day2) { day1 <- substitute(dfobj$day1); day1 day2 <- substitute(dfobj$day2) distance <- substitute(dfobj$distance) xx <- NULL for (i in 0:10) xx[i+1] <- with(dfobj, cor(Lag((day1-day1[1]),i), (day2-day2[1]), use='pair')) i <- match(max(xx), xx) with(dfobj, {

Hello, I am getting output from anova() and summary(aov()) that depends on the order of the factors in the fitted model object, and this has me baffled. I see this dependency with the data.frame below but not with an example (table 6.4) from Montgomery's DOE book. This is with R 1.3.0 on Debian GNU-Linux. Where have I gone wrong? > centerpts run sample CH50mg 1 day1 dev126 0.56 2

I would like to include this in a package. The S3 methods on R CMD check says * checking S3 generic/method consistency ... WARNING window: function(x, ...) window.chron: function(data, day1, hour1, day2, hour2, ...) See section 'Generic functions and methods' of the 'Writing R Extensions' manual. I have looked and can not figure it out. This function is for convience. What

If init doesn't run then disp->outp might not be initialized, resulting in an oops. Signed-off-by: Maarten Lankhorst <maarten.lankhorst at canonical.com> --- diff --git a/drivers/gpu/drm/nouveau/core/engine/disp/base.c b/drivers/gpu/drm/nouveau/core/engine/disp/base.c index c41f656abe64..9c38c5e40500 100644 --- a/drivers/gpu/drm/nouveau/core/engine/disp/base.c +++

Hi, I have a data.frame SAMPLES with columns: Site Site# Season Day1 Day2 Day3 Day1, Day2, Day3 are class "Date", the other columns are numeric or factor. I have a date "mydate" that may or may not be listed in my data.frame and I need to find that out. If "mydate" is there, I want to get the number of the data.frame row where it occurs.

I have this problem with records that have to be alligned on end- and startdate. I came up with the following. but am not convinced that this is the best way to tackle this problem. Anyone care for a braindump? def head_to_tail # Remove days that have startdate >= self.startdate AND enddate <= self.enddate # # before after # # |=====| |=====|

HI all, raw_urine = read.table("Z:\\bruce.9.3.09.sample.stability.analysis\\urine\\mz.spot.sam.dat.new", header = TRUE ) pvalue = read.table("Z:\\bruce.9.3.09.sample.stability.analysis\\urine\\all.urine.features.t.test.result", header = TRUE ) library(compositions) p = function(a,b){ y = pvalue[,a] if(y<0.01){ index = which(y, arr.ind=TRUE) day1 = raw_urine[index,3:7] day2 =

Dear group, Here is my environment : > ls() [1] "l" "PLglobal" "Pos100415" "Pos100416" "Pos100419" "Pos100420" "position" "select" "Trad100415" "Trad100416" "Trad100419" "Trad100420" "trade" "y" With objects : > l [1]

Readers, I have looked at various documents hosted on the web site; I couldn't find anything on interpolation. So I started r and accessed the help (help.start()). (by the way is it possible to configure r to open help in opera instead of firefox?) Initially I read the help for the akima package but couldn't understand it. Next I tried the asplines package help. I tried to copy the

Dear all, I have the following data in excel: day1 y 1 2 3 2 3 x1 0.2 0.3 0.4 0.3 0.2 x2 7 3.4 2 8 6 day2 y 2 4 3 2 2 x1 0.4 0.5 0.3 0.3 0.2 x2 7 8 9.1 6 5 I have the following problems: first of all, when I ask R to read the file (with the package xlsReadWrite and the command read.xls) it has a problem with the fact that the left most corner is labelled the same way, so in order for it

Hello, I'm fitting general linear models using the function lme() from the package nlme. My variables include a number of covariates and a single random factor, because the experiment was laid out in blocks. I'd like to have a statistic to measure the "proportion of explained variance" from my models. With ordinary multiple regression I'd use R-squared, but I can't

== Summary == * Problem: BibTeX entries extracted from R packages via citation() require too much manual editing to be of general use. * Proposal: Date: fields should be made mandatory in package DESCRIPTION files, perhaps beginning with warnings from R CMD check * Proposal: Package authors should be encouraged to use a (new) Contributors: field in the DESCRIPTION file rather than packing all

Hello, I just started learning R and have a very basic question: When I try ar model and extract residuals it returns Null. Could someone explain what's going on here? Does that mean the model is null? But it seems residual has a length=# observation of the original series according to the model summary. I am learning it by myself and have no one to ask here so please allow me to ask this

Hi there, I have searched what seems like everywhere and was unable to find a .NET wrapper so that I can use speex in VB. I did find a speex.NET.dll but this didn't seem to work on what I needed. I therefore decided I could just use DllImport within VB .NET and create my own unmanaged links to libspeex.dll. With the help of this partial c# code

Hi, If I have a data frame A with the following format: Day1 Day2 Day3 Day4 1 1979-11-02 1979-11-03 1979-11-04 <NA> 2 1979-12-06 <NA> <NA> <NA> 3 1979-12-13 1979-12-14 1979-12-15 1979-12-16 4 1979-12-20 <NA> <NA> <NA> And a date "1979-12-14", for

While I had thought I'd tested this before, it looks like this one issue slipped by my original CRC patches. Basically, there seem to be a few rules we need to follow when sending CRC commands to the display controller: * CRCs cannot be both disabled and enabled for a single head in the same flush * If a head with CRC reporting enabled switches from one OR to another, there must be a