similar to: how to estimate multidimensional spectral measure of coherence

I deal with a huge amount of Biology data stored in different databases. The databases belongig to Bioconductor organization can be accessed through Bioconductor packages. Unluckily some useful data is stored in databases like, for instance, miRDB, miRecords, etc ... which offer just an interactive HTML interface. See for instance http://mirdb.org/cgi-bin/search.cgi,

Unluckily I dela with miRNA files whose name may contain the character "*". Because of the special meaning of "*" I have to remove it. I found out how to make list.files() extract only those file names which contain a "*" Namely: # list.files(pattern="\\*") Now I have to process all files whose name does NOT contain the character "*". I cannot

Through converting a miRNAs file from FASTA to character format I get a vector which looks like the following: > nml [1] "hsa-let-7a MIMAT0000062 Homo sapiens let-7a" [2] "hsa-let-7b MIMAT0000063 Homo sapiens let-7b" [3] "hsa-let-7c MIMAT0000064 Homo sapiens let-7c" [4] "hsa-let-7d MIMAT0000065 Homo sapiens

I am looking for some explanations about the usage of the poorely documented R paclkage Rwave. Has anyone ever tried out its functions for Wavelet Analysis ? Thank you so much. Maura Alice Messenger ;-) chatti anche con gli amici di Windows Live Messenger e tutti i telefonini TIM! Vai su http://maileservizi.alice.it/alice_messenger/index.html?pmk=footer [[alternative HTML version deleted]]

Greetings, I would like to estimate a spectral coherence between two timeseries. The stats : spectrum() returns a coh matrix which estimates coherence (squared). A basic test which from which i expect near-zero coherence: x = rnorm(500) y = rnorm(500) xts = ts(x, frequency = 10) yts = ts(y, frequency = 10) gxy = spectrum( cbind( xts, yts ) ) plot( gxy $ freq, gxy $

I cannot find sowas package by Douglas Mauran in CRAN packages list- On which platforms does sowas run ? Has anybody used such a package at all ? hHank you very much, Maura tutti i telefonini TIM! [[alternative HTML version deleted]]

I tried to use R version of package I noticed the original MatLab Pckage is much better documented. I had a look at the R demo code "gtm_demo" and found that variable Y is used in advanced of being created: I wrote my own few lines as follows: inDir <- "C:/Documents and Settings/Monville/Alanine Dipeptide/DBP1/DHA" setwd(inDir) T <-

I have attached a text file representing the centralized amplitude of a signal, sampled at 30Hz, whose length N = 6922 My goal is to remove the trend. I am using package "simsalabim". I ran command decompSSA with L = length(Amps)/5 The reason is that I have SSA/MTM toolkit running in Mac/OS. SSA/MTM documentation, relative to SSA, recommends that N/10 <= L <= N/5.

I anticipate lacking of prior experience with dimensionality reduction problems. Some scientists concerned with drug discovery performed several steered Molecular Dynamics simulations of the alanine-dipeptide molecule dragged by a radial force from an equilibrium conformation to another different equilibrium conformation. They sampled at regular intervals 7 dihedral angles, 5 bending angles, and

-----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 The wiki page http://wiki.centos.org/RT_3.4.x_On_CentOS_4.x describes step by step how to install RT. Unluckily the install fails on step 4 because auf a missing dependency: Error: Missing Dependency: mod_perl >= 2.0.0 is needed by package rt Error: Missing Dependency: perl(Apache2::Log) is needed by package perl-Log-Dispatch Error: Missing

Dear R Community, I am currently student at the Vienna University of Technology writing my Diploma thesis on causality in time series and doing some analyses of time series in R. I have the following questions: (1) Is there a function in R to estimate the PARTIAL spectral coherence of a multivariate time series? If yes, how does this work? Is there an test in R if the partial spectral

Hello! I would like to ask if there is in R a function that estimates the spectral density function of a stochastic series at frequency zero by the "plug-in method", advocated by Andrews in his paper "Heteroscedasticity and Autocorrelation Consistent Covariance Matrix Estimation", Econometrica, 59,817-858. I saw R has functions that employ Andrews' plug-in method using an

Hi, I'm using smbclient to replicate a samba share from one Linux machine to another one. Unluckily, the amount of data is quite large, thus I'd like to transfer it with a compression, but I didn't find such as asny option in smbclient. Is it possible? Thanks, Luca P.S. I cannot use, in this case, other tools like rsync, just only smb!

Dear Andrew. You are right, I should have taken a deeper look into the standard output during compilation. I did just assume source4 was the one for Samba4. Anyway, besides the source confusion (really, even if I had found the right one, following the code would have been out of my reach), I don't seem to find how that is related with the documentation issue, or the mput/mask/recurse

I would like to compute the CSD of two finite time series. Is this implemented in R. Spectrum can handle multivariate input, but computes the coherence and the phase, but not the cross spectral density. Are there any packages out there that can do what I am asking? thanks Stephen Sefick -- Let's not spend our time and resources thinking about things that are so little or so large that all

Thanks Jan. I've been trying to obtain an English copy of the Dovecot book for months, prior to starting this project. So far, I just can't find a copy. It's too bad that the author/publisher won't do a second printing or, if they're not interested in making any more money, then release it to the public domain as a PDF. Very frustrating. Michael > -----Original

Hi, I have a bit of a problem with as.numeric or as.double. I read in an excel-file (either xlsx::read.xlsx2 or gdata::read.xls). Select a subset and then try to make it numeric: # read in the excel-file alldata<-read.xlsx2("input.xls",1) # select the subset s<-subset(alldata, select=c("cI","cII","cIII","cIV","cV")) # unluckily

Is there a way to force the number of ticks along an axis ? I read the on-line documentation and tried many combinations of all available parameters from functions par(), axTicks(), axis(), plot() ... but no luck ! Thank you very much, Maura tutti i telefonini TIM! [[alternative HTML version deleted]]

My question is motivated by dealing with pairs of values, like (3,12.5), (16,2.98), and so on that are mapped to a cartesian plain (Time, Frequence) I miss C multidimensional arrays. I am trying to simulate the 3rd dimension by declaring a matrix of lists. R seems to accept the definition of such a matrix. But I'm experiencing some difficulty at storing pairs of data in it. In the following

Dear Jeremy. Really appreciate you took your time to answer. I had already reviewed source4/client/client.c looking for mput (cmd_mput) command, but unluckily saying that my C programming skills are poor woud be really overrating them. I'm sorry to disagree for two reasons 1-. It works for mget, true that the source code is completely different. 2-. If taking a look into smbclient man