similar to: looking for Variable selections models, techniques, methods

Thank you for your prompt answer. The breathing signal observations are the amplitude values as a function of time and phase. According to our model the hidden states are the different breathing types. Subjects, whose respiratiion process is regular, are likely to breathe, keeping the same cycle pattern/type, for many consecutive cycles. therefore dwelling in the same hidden state. The more

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

How can I make sure the residual signal, after subtracting the trend extracted through some technique, is actually trend-free ? I would greatly appreciate any suggestion about some Stationarity tests. I'd like to make sure I have got the difference between ACF and PACF right. In the following I am citing some definitions. I would appreciate your thoughts. ACF(k) estimates the correlation

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 upgraded R and am currently running the following version: R version 2.10.1 Patched (2010-02-20 r51163) Copyright (C) 2010 The R Foundation for Statistical Computing ISBN 3-900051-07-0 The characteristics of my system are the following: OS: Linux 2.6.27.29-0.1-default x86_64 Current user: mauede@linux-326k System: openSUSE 11.1 (x86_64) KDE: 4.1.3 (KDE 4.1.3) "release

I have attached the signal that causes the error message in this email subject. Only columns 1 and 3 have to be considered. It is the work trajectory of a molecule migrating between two equilibrium conformations. The curve has 2 peaks, as shown in its plot. But I keep missing the 2nd one. Here is my short script: library(wmtsa) setwd("C:/Documents and Settings/Monville/Alanine

I deal with mono-channel breathing signals sampled at 30[Hz] which are non-linear and non-stationary. My goal is to classify the signals according to common breathing patterns Trend remotion is necessary for cluster analysis but quite challenging. In fact, quasi-periodic patterns that span a number of consecutive breathing cycles should not be naively removed as they carry a lot of useful

I would greatly appreciate some example of correct usage of function unz. I have to download and uncompress the following web compressef file: ftp://ftp.sanger.ac.uk/pub/mirbase/targets/v5/arch.v5.txt.homo_sapiens.zip I tried the following command that does not work: Targets.rec <- readLines(zz <- unz("ftp://ftp.sanger.ac.uk/pub/mirbase/targets/v5/arch.v5.txt.homo_sapiens.zip"))

I wonder whether it is possible in R to generate color-coded 3D plots, like the attached example. Basically a function f(x,y) (the 3rd dimension) is rendered through colors intensities. The side color-bar is a guide to the interpretation of the plot. Thank you very much, Maura e tutti i telefonini TIM! Vai su e tutti i telefonini TIM! Vai su

I have recently replaced R-2.9.0 with R-2.10.1 Patched. Apparently the installation completed successfully but right now I realized that the on-line help does not work any more. When I type "?<R-command>" a message pops up warning that "Help will not be available. Path not found" ... regardless of the R-command. How can I get back R on-line man pages ? Thank you very

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

http://gigamail.rossoalice.alice.it/messages/readMessageFrameset.aspx?DeliveryID=ba40cf18-29db-4404-a3ce-af26f760ecf9 Please, paste the website address above shown in your web browser address field. Make sure the whole string is pasted with no space or any other character. Telecom couldn't generate more clumsy website addresses .... Sorry for that. Thank you in advance, Maura tutti i

I am trying to generate a plot whose x-axis values are the following: data_out[,1] [1] 1979 1958 1937 1916 1895 1874 1853 1832 1811 1790 1769 1748 1727 1706 1685 1664 1643 1622 1601 1580 1559 [22] 1538 1517 1496 1475 1454 1433 1412 1391 1370 1349 1328 1307 1286 1265 1244 1223 1202 1181 1160 1139 1118 [43] 1097 1076 1055 1034 1013 992 971 950 929 908 887 866 845 824 803 782 761

I wonder whether R provides an interface to access miRecords data. Particularly, I am looking for extracting humans miRNA and target genes sequences. All such information is stored in there in a set of structured web site pages (http://mirecords.umn.edu/miRecords) I would greatly appreciate any suggestion even about other data bases from where it is possible to get the same sort of data. I had a

Is there any R package which implements non-linear dimensionality reduction (LLE, ISOMAP, GTM, and so on) and/or intrinsic dimensionality estimation ? Thank you, Maura tutti i telefonini TIM! [[alternative HTML version deleted]]

Can biomaRt connect to data base "http://mirecords.umn.edu" or a branch of it ... for instance the validated miRNAs list ..? Thank you very much. Maura tutti i telefonini TIM! [[alternative HTML version deleted]]

Is there any R package that implements the same capability of MatLab toolbox called SimBiology ? We are expecially interested in protein-protein interactions and network analysis. As far as I know SimBiology implements a system of ODEs reflecting the kinetic chemical reactions. We would be more interested in stochastic simulations. Thank you in advance. Maura tutti i telefonini TIM!

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 am assigning subset of a matrix A [n,3] where n>1 to a temporary matrix TMP I do not know how many rows of A will be assigned to TMP because this is established by a run-time test. I expect TMP to be a matrix [m,3], m >=1 But when 1 row only is transferred from A to TMP then TMP becomes [3,1] rather than [1,3] How can I avoid this unwanted transpose operation ? THank you in advance,

I installed package plotrix because reading its vignette it looks like it can help me solve a "legend" problem. The package instaleed correctly on my Mac OS/X 10.5.8 But I cannot reproduce the examples centered on function "lgendg". > library(plotrix) > plot(0.5,0.5,xlim=c(0,1),ylim=c(0,1),type="n", + main="Test of grouped legend function") >