similar to: Problem with NLSstClosestX; and suggested fix

Dear list, I have been studying R and I would like the aid of more experienced to solve the problems of the analysis below: r = gl(3, 8, label = c('r1', 'r2', 'r3')) e = rep(gl(2, 4, label = c('e1', 'e2')), 3) y = c(26.2, 26.0, 25.0, 25.4, 24.8, 24.6, 26.7, 25.2, 25.7, 26.3, 25.1, 26.4, 19.6, 21.1, 19.0, 18.6, 22.8, 19.4, 18.8, 19.2, 19.8, 21.4,

Hi, I'm wondering why the spectrum() phase of quadrature couple isn't purely +/-pi. But mostly, I'm looking for a recommended way to take a 2-D spectrum and 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(

Dear friends i am trying to fit an Ornstein-Uhlenbeck process by MAXIMUM LIKELYHOOD method. i found these formulas on http://www.sitmo.com/article/calibrating-the-ornstein-uhlenbeck-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)

--- libcelt/Makefile.am | 6 +- libcelt/bands.c | 26 +++++++++- libcelt/bands.h | 2 +- libcelt/celt.c | 23 +++----- libcelt/pgain_table.h | 133 ------------------------------------------------- libcelt/quant_pitch.c | 117 ------------------------------------------- libcelt/quant_pitch.h | 44 ---------------- 7 files changed, 37 insertions(+), 314

Dear R users, My immediate problem is that I do not seem to understand how to use the function scale_x_datetime function in ggplot. My deeper problem is that I don?t understand the use of the POSIXct class. ? My data looks like like the following : T????????? val 2009-11-04 23:59:57?? 972.357117???? 2009-11-05 00:02:21?? 969.389221???? 2009-11-05 00:04:45?? 868.817749???? 2009-11-05 00:07:09??

When plotting over multiple pages in lattice, I'd like to be able to have "same" scales within a page, but "free" scales between pages. In other words, something like: z<-data.frame(x=1:100, y=runif(100), d=rep(1:2,50), p=rep(1:2,each=50)) plot<-xyplot(y~x|d*p, data=z,

Full_Name: Quanah Gibson-Mount Version: 2.3.0, 2.3.1 OS: Linux Submission from: (NULL) (171.66.155.86) libc used to export the symbol: __libc_stack_end; however, newer versions of libc no longer export it. R has some serious problems around this, because the code has: # ifdef linux extern void * __libc_stack_end; # endif and #if defined(linux) R_CStackStart = (uintptr_t)

Hello, nls library provides 6 self-starting models, among them: SSfp, a four parameters logistic function. Its self-starting procedure involves several steps. One of these steps is: pars <- as.vector(coef(nls(y ~ cbind(1, 1/(1 + exp((xmid - x)/exp(lscal)))), data = xydata, start = list(lscal = 0), algorithm = "plinear"))) which assumes an initial value of lscal equal to 0. If lscal

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

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

I've rolled up R-2.12.0.tar.gz a short while ago. This is a development release which contains a number of new features. Also, a number of mostly minor bugs have been fixed. See the full list of changes below. You can get it from http://cran.r-project.org/src/base/R-2/R-2.12.0.tar.gz or wait for it to be mirrored at a CRAN site nearer to you. Binaries for various platforms will appear in

I've rolled up R-2.12.0.tar.gz a short while ago. This is a development release which contains a number of new features. Also, a number of mostly minor bugs have been fixed. See the full list of changes below. You can get it from http://cran.r-project.org/src/base/R-2/R-2.12.0.tar.gz or wait for it to be mirrored at a CRAN site nearer to you. Binaries for various platforms will appear in

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

Full_Name: Ivan the Terrible Version: 2.9.2 OS: Windows XP SP3 Submission from: (NULL) (89.110.13.151) When using the method blackboost_fit of the package mboost appear following error : Error in party:::get_variables(obj at responses) : trying to get slot "responses" from an object (class "boost_data") that is not an S4 object Simple test case that produce bug:

----=_NextPart_ST_21_29_48_Monday_February_23_2004_11097 Content-Type: text/plain; Your message could not be sent. A transcript of the attempts to send the message follows. The number of attempts made: 1 Addressed To: john@spaar.com Mon, 23 Feb 2004 21:29:47 -0600 (2) Failed to send to identified host, [10.0.0.2], 550 5.1.1 User unknown [SMTP server error] --- Message non-deliverable.

Full_Name: Luzi P. Schucan Version: R-base-0.90.1-2.i386.rpm OS: Linux Submission from: (NULL) (141.84.136.129) I just wanted to install the rpm package with rpm --v -i [file], and here is the log: (I DID run it as root!!) (the important thing is in the end: there must be a bug in the post install script. The problem is, that it does install all the files, but it doesn't correctly hang the

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