similar to: errors appears in my time Series regression fomula

Dear useRs, I'm trying to simply fill in the area under a curve using RGL. Here' the set up: x <- c(0.75,75.75,150.75,225.75,300.75,375.75,450.75,525.75,600.75,675.75, 0.5,50.5,100.5,150.5,200.5,250.5,300.5,350.5,400.5,450.5, 0.25,25.25,50.25,75.25,100.25,125.25,150.25,175.25,200.25,225.25) y <- c(0.05,4.91,9.78,14.64,19.51,24.38,29.24,34.11,38.97,43.84,

Hello, I'm new to this and am trying to teach myself some R by plotting biological data. The growth curve in question is supposed to be fitted to the Verhulst equation, which may be transcribed as follows: f(x)=a/(1+((a-0.008)/0.008)*exp(-(b*x))) - for a known population density (0.008) at t(0). I am trying to rework the example from "An Introduction to R" (p. 72) for my case and

Hi, I am using aggregate() to add up groups of data according to year and month. It seems that the function aggregate() automatically sorts the levels of factors of the grouping elements, even if the order of the levels of factors is supplied. I am wondering if this is a bug, or if I missed something important. Below is an example that shows what I mean. Does anyone know if this is just the way

Hi: I need to plot two time series in the same plot and they cover the same time range and have the same frequency. With RSiteSearch("multiple series plot") i found this post by Gabor Grothendieck: http://finzi.psych.upenn.edu/R/Rhelp02a/archive/42281.html Exactly what i need except for one detail. I want one series to be made of points and the other by a line. When I simply try: xts

Hi: I need to plot two data sets and then I need to change labels of xaxis. I can manage to achieve something like it with this: x<-data.frame(a=c(1,2,3,4,5),b=c(2,4,6,8,10)) y<-data.frame(a=c(3,4,5),b=c(1.5,2,2.5)) xts <- ts(x$b,start=x$a[1]) yts <- ts(y$b,start=y$a[1]) ts.plot(xts,yts,col=c("red","blue")) But ts.plot does not allow to change x labels. Besides,

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 $

What''s a rule of thumb for guesstimating how many Mongrels to use in a cluster for an app? I have an app that gets about 5000 unique visitors per day. I figured I''d give it plenty of Mongrels -- twenty to be specific. After running out of memory and hitting the swap periodically, I scaled it back to five and it still seems to serve up visitors fine. So, is there some super-secret

for unweighted fits using `nls' I compute confidence intervals for the fitted model function by using: #------------------- se.fit <- sqrt(apply(rr$m$gradient(), 1, function(x) sum(vcov(rr)*outer(x,x)))) luconf <- yfit + outer(se.fit, qnorm(c(probex, 1 - probex))) #------------------- where `rr' contains an `nls' object, `x' is the independent variable vector, `yfit'

Hi, don' t understand why the function fomula have this error, i enclose the parameter "a" with the function I() Thank Ruben x<-1:5 y<-c( 2 ,4 , 6 , 8 ,11) formu<-y~I(a*x) form<-formula(formu) dummy<-data.frame(x=x,y=y) fm<-lm(form,data=dummy) Error in unique(c("AsIs", oldClass(x))) : Object "a" not found

Hi, I am fairly new to R so this may be easy for some. I have one column that is the Type="A, B' T" and I have another column that lists Price='894, 895, 896' I wanted to plot three different series based upon the Type...If B than xts..If T yts... How exactly do I do this? Thanks in advance -- View this message in context:

dear list members, I apologize in advance for posting a second time, but probably after one week chances are, the first try went down the sink.. my question concerns computation of confidence intervals in nonlinear fits with `nls' when weigthing the fit. the seemingly correct procedure does not work as expected, as is detailed in my original post below. any remarks appreciated. greetings

This is one of those simple-been-doing-this-forever things that, for some reason, has me stumped today. When I try to NFS (v4) mount a directory, the user/group ownership shows up as user "nobody" even though /etc/passwd has values for the correct user names. How do I get it to mount with the correct user IDs? Hume is the server, running CentOS 6, all updates applied, maybe a week

dear list members, my question concerns computation of confidence intervals in nonlinear fits with `nls' when weigthing the fit. the seemingly correct procedure does not work as (I) expected. I'm posting this here since: (A) the problem might suggest a modification to the `m' component in the return argument of `nls' (making this post formally OK for this list) and (B) I got no

I think this is another one of the same kind of bugs in ts: Version 0.64.1 (May 8, 1999) ... > z <- ts(matrix(1:20,10,2), start=c(1969,1), frequency=12) > (z > 5) | (z < 2) Error: invalid time series parameters specified > Paul -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-devel mailing list -- Read

say we have: > x<-data.frame(a=c(1,2,3,4,5),b=c(1,1,1.5,2,2)) > y<-data.frame(a=c(1,2,3,4,5),b=c(1,2,2,3,3)) How would I plot this so that, with the shared $a as the x-axis values, I have both $b columns plotted together? (a comparison of the two?) thanks. DISCLAIMER: This e-mail message and any attachments are inte...{{dropped}}

Hi, we came across the following unexpected (for us) behavior in terms.formula: When determining whether a term is duplicated, only the order of the arguments in function calls seems to be checked but not their names. Thus the terms f(x, a = z) and f(x, b = z) are deemed to be duplicated and one of the terms is thus dropped. R> attr(terms(y ~ f(x, a = z) + f(x, b = z)),

Dear R-help, This is an example in the {Hmisc} manual under rcorr.cens function: > set.seed(1) > x <- round(rnorm(200)) > y <- rnorm(200) > round(rcorr.cens(x, y, outx=F),4) C Index Dxy S.D. n missing uncensored Relevant Pairs Concordant Uncertain 0.4831 -0.0338 0.0462 200.0000

Hello Kindred R Spirits... I'm trying to get a file (csv) from Dropbox using their shareable link concept. They issue a short URL that goes to a web page where you see a button that says "Download File". They don't really give you the URL of the file itself, just this page. Is there a way to coax R into getting such a file? I don't even really want the file per

Hello samba, I would like an advice on following: is there any way how to tell samba what transfer mode to use? I think I read that samba does only binary and to use text you have to use something else. We have windows machines using text files on linux samba server and i need to transfer those files in text mode. Can anyone help me please? -- Best regards, Martin

Dear Achim, >>>>> Achim Zeileis <Achim.Zeileis at r-project.org> >>>>> on Fri, 10 Mar 2017 15:02:38 +0100 writes: > Hi, we came across the following unexpected (for us) > behavior in terms.formula: When determining whether a term > is duplicated, only the order of the arguments in function > calls seems to be checked but not