similar to: nonlinear regression

Hi It is preferable to echo your posts to r-help, you usually get more answers and some definitelly superb to mine. It is also better to start a new mail if your question has nothing to do with original subject "Maura E Monville" <maura.monville at gmail.com> napsal dne 01.10.2007 17:44:43: > Unluckily I do not have the privilege of practising with R all day > long. I

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 have 316 files storing a number of cycles of the respiratory signal from 316 different patients. Each file record is made up of the following data: amplitude,phase,timestamp,validflag,ttlin,mark,ttlout amplitude: is the position in cm relative to an arbitrary reference (signal amplitude) phase: is the phase value for the current sample (the time point at which the breathing signal was

I am reading the comprehensive on-line documentation about msm. The positive side is that it seems it has been designed for biomedical statistics, like Clinical Trials. The bad side is that it does not seem to model observations sequences that are not independent but instead are autocorrelated, as it is my case. I did not find any mention to correlated observations therefore I assume the authors

I am reading some documentation about Cross Spectrum Analysis as a technique to compare spectra. My understanding is that it estimates the correlation strength between quasi-periodic structures embedded in two signals. I believe it may be useful for my signals analysis. I was referred to the R functions that implement this type of analysis. I tried all the examples which generated a series of

*I would like to plot some semi-periodic signals as a function of the phase expressed as a value in [0,2PI]* *The problem is that the phase data is reported as the residual of the division by 2PI.* *For instance if the phase is 10.359 rd then:* *15.359 / 6.283185 = 2*6.283185 + 2.79263* *then only the residual 2.79263 is stored. * *In order to plot a number of consecutive cycles I have to add

I have 316 files. Each file represents a patient's breathing track (respiratory signal recorded for a variable number of cycles). All files have the same are made up of a header followed by a variable number of records. Each record contains 7 comma separated fields. The patient ID is recorder in the header which is stripped off when reading the file into a R data.frame. Since I need to keep

We are trying to build a human respiration model. Preliminary analysis of some breathing signals has shown that humans breathe through switching among a finite number of patterns. Hidden Markov seems to be the right approach. Since most of our code is written in R scripting language, finding an R package implementing an HMM that we can use for our prototype would be very helpful. I have been

I would appreciate some help with the followiing error ... what did I do wrong and what is the woraroud. > library(nlme) > glmod <- gls(A ~ T +cosP +cos2P +cos4P +cos5P +sin3P +sin5P +cosP2 +sinP3 +P2, correlation=corAR1(),method="ML") Error: cannot allocate vector of size 491.3 Mb I have a bunch of humans' breathing signals. The recording frequency was 30 samples per

I tried to install "segmented" from three different repositories but I keep getting the following message on R window: ... Please, advice where to find this package. Thank you in advance. Maura > utils:::menuInstallPkgs() --- Please select a CRAN mirror for use in this session --- Warning: package 'segmented' is in use and will not be installed > utils:::menuInstallPkgs()

What is wrong in the following nested if-else statements: if (Condition_1) { # begin IF_1 statement_1 statement_2 statement_3 if (Condition_2) { # begin IF_2 a<- a +1 } # end IF_2 statement_4 statement_5 statement_6 statement_7 if (Condition_3) {

Here is another occurrence of wmTSA internal error. My time series is a short breathing cycle (2425-Cyle_9.txt). Since wmtsa functions that extract extrema seem to expect longer series than I have, I tried the folowing two tricks: 1) I prolong the 1-cycle series on both ends through duplicating the first value (on the left end( and the last value on the right end as any ties as the 1-cycle

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

Can R handle graphic text ? I mean something like a pop-up message window or a text widget ? Thank you so much. -- Maura E.M

I saw there exists an R version for Mac/OS. I'd like to hear from someone who is running R on a Mac/OS before venturing on getting the following computer system. I am in the process of choosing a powerful laptop 17" MB PRO 2.6GHZ(dual-core) 4GBRAM .... Thank you so much, -- Maura E.M [[alternative HTML version deleted]]

I have downloaded the source code for Linux. the ./configure procedure can find gcc but cannot find f77 or f2c I have looked for the above compiler with Test and installed all what pertains to Fortran. Still not even the man pages for f77 anf fc2 work >. Where is f77 and f2c with linux 10.3 ?? Thank you do much -- Maura E.M [[alternative HTML version deleted]]

Hi all, I'm trying to fit a nonlinear (logistic-like) regression, and I'd like to get some recommendations for which package to use. The expression I want to fit is something like: y ~ A * exp(X * Beta1) / (1 + exp(-(x + X * Beta2 - xmid)/scal)) Basically, it's a logistic function, but I want to be able to modify the saturation amplitude by a few parameters (Beta1) and shift the

I recently installed R 2.6 on Linux/SuSE When I was running the previous version on Windows I was able to save my plots from a script. There was a command "savePlot" that is no more retrieved in the last version. In this scenario, how can I tell R to save the currently displayed plot to a file ? The "save" command is generic and do not know how to tell R that the plot is to be

I'd like to know why I cannot get a plot and the QQnorm in the same sheet. The commands are simple but: library(nlme) glmod1 <- gls(upfmla,correlation=corAR1(),method="ML") summary(glmod1) par(mfrow = c(2,1)) plot(glmod1, main="GLS Residuals vs. GLS Fitted") qqnorm(glmod1) No matter what (I tried different permutations of the plotting commands) the second drawing

I have a loop that processes all the raw data files in agiven directory. It should allocate only the memory to accomodate the filenames list. Each file is loaded into RAM one at a time and saved to another directory after processing. The R script runs fine but it had some memory problems so it was terminated before getting to the end of the list. This is a Windows XP desktop with ~ 1GB RAM In the