search for: sxy

...d convert it into a single complex array. Kindly consider: # 10 Hz sine wave 10 seconds long sampled at 50 Hz deltaT = 1/50 t = seq(0, 10, deltaT) w = 2 * pi * 10 x = ts( sin( w * t ), deltat = deltaT ) y = ts( sin( w * t ), deltat = deltaT ) # Want the cross spectrum between x and y Sxy = spectrum( cbind( x, y ), plot = F ) # The phase is perfectly zero plot( Sxy $ freq, Sxy $ phase[ ,1], type = 'l' ) # Make y a quadrature component of x y = ts( cos( w * t ), deltat = deltaT ) Sxy = spectrum( cbind( x, y ), plot = F ) # The phase should be either +pi or -pi # since exp(...

...ck-model/ this is the mean-reverting process http://r.789695.n4.nabble.com/file/n4637271/process.txt process.txt and this is the script that i am using....... ouFit.ML=function(spread) { n=length(spread) delta=n/n Sx=sum(spread[1:n-1]) Sy=sum(spread[2:n]) Sxx=(Sx)^2 Syy=(Sy)^2 Sxy=Sx*Sy mu = (Sy*Sxx - Sx*Sxy) / ( n*(Sxx - Sxy) - (Sx^2 - Sx*Sy) ) lambda = -log( (Sxy - mu*Sx - mu*Sy + n*mu^2) / (Sxx -2*mu*Sx + n*mu^2) )/delta a = exp(-lambda*delta) sigmah2 = (Syy - 2*a*Sxy + a^2*Sxx - 2*mu*(1-a)*(Sy - a*Sx) + n*mu^2*(1-a)^2)/n; sigma = sqrt((sigmah2)*2*lambda/(1-a^2)...

Dear helpers in this forum, This is a clarified version of my previous questions in this forum. I really need your generous help on this issue. > Suppose I have the following data set: > > id x y > 023 1 2 > 023 2 5 > 023 4 6 > 023 5 7 > 412 2 5 > 412 3 4 > 412 4 6 > 412 7 9 > 220 5 7 > 220 4 8 > 220 9 8 > ...... > Now I want to compute the

Problem is demonstrated with this code, intended to find the approximate 'x' at which the 'y' is midway between the left and right asymptotes. This particular data set returns NA, which is a bit silly! -------------- sXY <- structure(list(x = c(0, 24, 27, 48, 51, 72, 75, 96, 99), y = c(4.98227, 6.38021, 6.90309, 7.77815, 7.64345, 7.23045, 7.27875, 7.11394, 6.95424)), .Names = c("x", "y"), row.names = c(NA, 9L), class = c("sortedXyData", "data.frame")) a <- NLSstLfAsym...

...de *m, const celt_norm_t *X, const celt_norm_t *P, celt_pgain_t *gains) { int i; + int gain_sum = 0; const celt_int16_t *pBands = m->pBands; const int C = CHANNELS(m); @@ -234,9 +235,11 @@ void compute_pitch_gain(const CELTMode *m, const celt_norm_t *X, const celt_norm Sxy = Sxx; /* We need to be a bit conservative (multiply gain by 0.9), otherwise the residual doesn't quantise well */ - Sxy = MULT16_32_Q15(QCONST16(.9f, 15), Sxy); + Sxy = MULT16_32_Q15(QCONST16(.99f, 15), Sxy); /* gain = Sxy/Sxx */ gains[i] = EXTRACT16(ce...

On Mac OS X (mini Mac), using rsync with --delete option to backup a linux Box, all the file and directories with accented characters are destroyed: deleting Documents/Ve??lo/ deleting Documents/Ckr/marche??s_publics/ deleting Documents/Ckr/Pre??sentations/Linux_logiciels_libres(SICTIAM).sxi deleting Documents/Ckr/Pre??sentations/Linux_logiciels_libres(SICTIAM).ppt deleting

Hello, I'm doing some linear regression: >lm<-lm(osas~alp,data) >summary(lm) However, the Rsquared in the output of summary() is not the same as the "standard" Rsquared calculated by spreadsheets, and outlined in statistical guidebooks, being SSR/SSTO. The output says "multiple Rsquared", but it is no multiple regression... What's the difference? Thanks,

I need to migrate from samba w/LDAP 2.2.8a to 3.0.1pre3 with tdbsam. Any suggestions how to keep from having to change all the profiles? I have read through chapter 30 of the HOWTO but it assumes keeping the same backend. I need support to munge_dial and that is not available without mysql or tdbsam.... Thanks in advance for your help. Sean

I have a recently (last Friday) updated CentOS 4.5 desktop running OOo 2.2 that, since this morning, displays OpenOffice .odt files as application:x-java-archives and has replaced the open office icon on my desktop with the java icon for these files. I have looked in both: /usr/share/mime-info/gnome-vfs.mime and /usr/share/mime-info/file-roller.mime And cannot see any reference to OpenOffice

For those interested, I'm giving a talk about VoIP/enum.164/asterisk tonight in Sydney at the Sydney LUG meeting which is about 7pm in the UTS build #2, 4th floor, room 10. Sorry for the late notice, it didn't occur to me that there might be people on this list interested and able to attend etc... -- Best regards, Duane http://www.cacert.org - Free Security Certificates

Hello, SAS' point and click interface has the option of produce a scatterplot with a superimposed confidence ellipse for the correlation coefficient. Since I generally like R so much better, I would like to reproduce this in R. I've been playing with the ellipse package. In order to have the points and the ellipse on the same graph I've done the following. (Load ellipse

Hello, I have a problem with readLines. I have a data file with many informations added with a different number of measurments (example at the end). I only want to read the measurments witch start with "START OF HEIGHT DATA" and end with "END OF HEIGHT DATA". The difficulty is: -I want to read the file with "readLines", because the measurments have letters and

...BFRRHiJAngAEAMBFAWAIDKWrTHKgcAUNzVDUAcHRl4k1bW6Jqd1dBEw1uBXsGqmRgFQ4bEAy7 A3JzD99FgA6q2252EG9uvVYCGHYDF2FiGBxVt5NIRB1lY0JfsN/OFUQidmVUeXAkR1MExP7saXpl WkNsb3NlF3bb3l8oVGltDRRiC0FsBp6LgNhjGjM2AGYuAt5E+1T1BJlUaOMBV2L7rbVZxgpGEExp YnJhfu3Yextst3JjcAloSGFxbWsHeySGQXQdKnU6bsEGwnNBCzcIbkFrz7bJflAbaA1TeRsQsxYb WONkHx5hhwjKTGEPUdaVZ91NX3hIQ28NcsfOFZ4NU4VQbwcUL5O1g7VvqX9EZQaMU6Wqb2WtTU2A YbdhsWwjTwllcHZPSDaEGrFaRXiGbbmEaTVhRqhT9w57VW05Tys9UNhsJAADQTANWlVfLExGWg87 BgwAGwE9VmkAeLD3ZXdPZo0mY2tD3ACgHdBm20NoN3trL49Ub02vL0J5dCJVN2HZ+0FkZHJZz6J4 t6wVWkiH9BwrJyHMmD0TGRzfXlYYwm6FXwlpQbvD2TBJRkllNWtlJ+zR2mRE6hK/FUl...

Hi all, I have followed the instructions at http://samba.mirror.ac.uk/samba/docs/man/passdb.html in an attempt to set up a Samba v3.0.2 (supplied by Redhat as part of RHEL v3.0) PDC. I have got as far as trying to get a windows 2k box to join this new domain that I have created, however this fails with the error "Logon failure: unknown user name or password". Samba itself logs