search for: woriking

I have two pairs of related vectors x1,y1 and x2,y2 I wish to do a test for differences in means of x1 and y1, ditto x2 and y2. I am getting odd results. I am not sure I am using 'pt' properly... I have not included the raw vectors as they are long. I am interested if I am using R properly... > c(length(x1), length(y1), length(x2), length(y2)) [1] 3436 1619 2677 2378 First

I have trouble with this: as.Date("Sep-1981", format="%b-%Y") Returns "NA" >From documentation for strftime '%b' Abbreviated month name in the current locale. (Also matches full name on input.) '%Y' Year with century. What am I doing wrong? cheers Worik [[alternative HTML version deleted]]

It seems obvious to me that the empty string "" is length 0. cheers Worik [[alternative HTML version deleted]]

Is it possible to get a line number with an error report? I have a long script and an error: Error in `[.xts`(x, xsubset) : subscript out of bounds It would be very helpful, and save a lot of time, if there was some indication in the error message which line the error was. I can find it using binary search but that is a painful process. cheers Worik [[alternative HTML version deleted]]

I am interested in the package 'rgp' But.... > install.packages("rgp") Installing package(s) into ‘/home/worik/R/x86_64-pc-linux-gnu-library/2.12’ (as ‘lib’ is unspecified) Warning message: In getDependencies(pkgs, dependencies, available, lib) : package ‘rgp’ is not available It is in cran: http://cran.r-project.org/web/packages/rgp/index.html How do I install it?

When I say: > Sys.time() [1] "2012-11-12 21:30:14 NZDT" > But that is not what my clock on the wall and my system say. Cannot show you my clock but... worik@lemy:/tmp$ date Tue Nov 13 10:32:20 NZDT 2012 Sys.time() is returning GMT $version.string [1] "R version 2.14.1 (2011-12-22)" > Sys.timezone() [1] "NZDT" > I'm a little lost! Worik

I think I am whistling in the wind, but is there a modern symbolic debugger for R programmes? I am working through some one else's code, thousands of lines, that has the occasional bug in it, and a lot in my understanding of it. I cannot make setBreakpoint or findLineNum work. I get "No source refs found.". I am starting to loose my mind! A debugger where I could set

Friends I am trying to format a number to a string so 2189.745 goes to "2,189.35" and 309283.929 goes to "309,283.93" I have tried to use formatC(X, big.mark=",",drop0trailing=FALSE, format="f") but it does not get the number of decimals correct. Specifying digits does not work as that is significant digits. I could use a switch statement switching on

Friends I have an xts that I wish to access. Browse[2]> DATA.ba[[p]]["2012-03-20 00:59:57","bid"] bid 2012-03-20 00:59:57 1.4993 So far so good. Now putting the index into a variable: Browse[2]> Time [1] "2012-03-20 00:59:57 NZDT" Browse[2]> DATA.ba[[p]][Time, "bid"] bid Where has it gone? Looking closer....

If I create a vector thusly > v1 <- runif(20, min=0, max=1) > v1 [1] 0.9754443 0.6306228 0.3238158 0.3175769 0.6791534 0.6956507 0.3840803 [8] 0.1421328 0.8592398 0.4388306 0.9472040 0.4727435 0.5645302 0.7391616 [15] 0.6116199 0.2727754 0.2657867 0.5261744 0.8764804 0.2032126 And I want to create a logical vector the same length that is true if v1<.3 or v1 > .7 how do I do it

How can I increment the value of a Date class? I want to add a day, month or year to a date. cheers Worik [[alternative HTML version deleted]]

This programme for(T in 1:3){ for(j in 1:(5-1)){ for(k in (j+1):5){ for(l in (j+2):5){ print(paste("1 JKL:", j,k,l,sep=" ")) } } } } Prints out (among other things) [1] "1 JKL: 4 5 6" That is for(l in (j+2):5) sets l to 6 one more than the upper limit. cheers Worik [[alternative HTML version deleted]]

I have a large programme that after running half an hour or so fails with an error Error in x[value] <- NA : only 0's may be mixed with negative subscripts How can I find out where that error occurs? If I have to do a binary search using error messages it will take a long time! Is there some way I can generate a stack trace? cheers Worik [[alternative HTML version deleted]]

Friends If I have a matrix such as... [,1] [,2] [1,] 7 7 [2,] 7 9 [3,] 9 2 [4,] 7 9 And I want to find the row number that has the minimum value of column 2 (row 3 in this case) how can I do it? Is there a simple way? cheers Worik [[alternative HTML version deleted]]

Friends I am looking at Rcpp and I am a bit stuck on a simple matter. (I am calling R from c++, if there is a better way...) Given this simple example using the TTR package and the SMA function which returns a simple moving average.... Rcpp::NumericVector rv; for(int i = 0; i < 100; i++){ rv.push_back(rand()); } Rcpp::Environment TTR("package:TTR");

I have time series data in named vectors. They are all the same length for the same dates. The dates are in a separate vector. I want to create a vector of numeric data for every named series, associated with the dates in a data.frame. So if... Names <- c("a", "b", "c") d <- data.frame(dates=Dates, a=vector(mode="numeric", length=length(Dates),

Friends I am new to R (and statistics) so am struggling a bit. Briefly... I am interested in getting the P-Value from cor(X) where X is a matrix. I have found cor.test. Verbosely... I have 4 vectors and can generate the corellation matrix... > cor(cbind(X1, X2, X3, X4)) X1 X2 X3 X4 X1 1.00000000 -0.06190365 -0.156972795 0.182547517 X2

Friends I clearly donot understand how sapply and vapply work. What I have is a function that returns a matrix with an indeterminate number of rows (some times zero) but a constant number of columns. I cannot reliably use an apply function to assemble the matrices into a matrix. I am not sure it is possible. I can demonstrate the core of my confusion with this simple code. A.f <-

Friends I am starting on a (section of the) project where I need to build a matrix with on the order of 5 million rows and 200 columns I am wondering if I can stay in R. I need to do rollapply type operations on the columns, including some that will be functions of (windows of) two columns. I have been looking at the ff and bigmemory packages but am not sure that they will do. Before I get

Hi! I have just added login feature to my rails application (followed the instruction from (Agile Web Development with Rails) but now I am getting strange errors (looks to me like some kind of routing problem) under the passenger/apache. The error is: 500 Internal Server Error The server encountered an internal error or misconfiguration and was unable to complete your request. From the log: