similar to: multiple replications of a ts

There is probably an elegant and efficient way to do this -- perhaps already implemented, and some of you knows it. Myself, I'm getting stuck in long loops ... Here is the question: I have a long (80,000 + obs) binary time series. For each of its elements I'd like to compute and save in a vector "the number of steps 'til the first switch of state happens". For example,

Hi all I have a dataset whose response is a categorical variable of the ordinal scale type (6 levels). I'm interested in building classification models, and I'm wondering if there is something implemented in R (or its packages) that I'm not aware of, to treat the ordinal scale measurements straightforwardly. I can think of the alternative of building a conditional (hierarchical)

I'm running R-1.4.1 on a linux debian (version 2.2) box. I've tried to upgrade to R-1.6.0, and the execution of make halts with the following error messages: make[4]: Entering directory `/home/tebaldi/R/R-1.6.0/src/library/methods' dumping R code in package 'methods' Error in testRversion(descfile) : This package has not been installed properly See the Note in ?library

Hi all I'm running lm()/leaps() kind of computations within a loop, and I submitted the entire thing by sourcing a file of R code in interactive mode. I'd like to interrupt the evaluation of this (it's taking way too long) but I am not able to. I'm running R-1.3.0 from an emacs window (ESS-5.0) and I'm on a linux box, (Debian). I tried Control-c Control-c (nothing happens),

Hi I'm trying to print several plots to a multi-page postscript file. Some plots are output of basic commands, some are produced by "trellis" commands (from the package lattice). Right now I'm not able to get a straightforward black and white color scheme for the latter kind: when I open the postscript file the "trellis" plots are invisible, or parts of them are

Hello I'm upgrading from R-1.4.1 to R-1.5.1, locally on my Linux PC running Debian. after running 'configure' this is what I get (just in case it helps figuring out what follows later...) R is now configured for i686-pc-linux-gnu Source directory: . Installation directory: /usr/local C compiler: gcc -D__NO_MATH_INLINES -mieee-fp -g -O2 C++

Dear R-helpers, I'm plotting geophysical data in the form of contours using "filled.contour". The display would be much more effective if the areas with negative values could be color coded by -- say -- "cold colors" in the range of blue to green, and conversely the areas with positive values got plotted with "warm colors", from yellow to red. Right now if I use

I'd like to apply the ideas of Pinheiro et al. (1997): "Robust Estimation in LInear Mixed Effects Models using the Multivariate t-Distribution" and I'm wondering if the EM-based estimation described in that paper has been implemented by somebody in R already (I'm assuming it's not yet available within the package nlme -- whose subtitle is "Fit and compare

I thought I had plotted a map of oregon and washington and overlayed some points using R. Is that possibe since I can't seem to find the source file I used or commands for plotting maps and points? If not, sorry for the bother... Thanks, Jeff. Jeff D. Hamann Hamann, Donald & Associates, Inc. PO Box 1421 Corvallis, Oregon USA 97339-1421 Bus. 541-753-7333 Cell. 541-740-5988 jeff_hamann

Hello Is there a simple way of fitting a 'smooth curve' to a time series with the additional constraint that the first derivative be not negative? I know that if I choose a "large enough" window for a moving average sort of filter on my data I end up with a non-decreasing curve anyway, but I'm wondering if there is a pre-built method that can include the constraint from

Here is the answer to an old question of mine, regarding maps in R Hope this help. claudia ---------- Forwarded message ---------- Date: Mon, 25 Jun 2001 13:48:05 -0400 From: Kieran Healy <kjhealy at Princeton.EDU> To: Claudia Tebaldi <tebaldi at ucar.edu> Subject: Re: [R] maps in R Hi Claudia -- you can get the maps() library from here:

Hi all Just this morning I upgraded to R 2.12.0 (for Mac OS X 10.6.4). All went well until I needed to run a help() or help.search() in my session, which I'm running within Emacs (ESS 5.3.7). Say I need help with the command 'density'. When I type help(density) or ?density the ESS help buffer opens, it is titled *help[R](density)* but it contains only a couple of lines saying,

Greetings, I hope this isn't obvious but I've read through the FAQ and the various other R/S books I have and I cannot, for the life of me, figure out how to sort a data frame by a specific variable. Suppose I have a four variable data.frame x and the variables are var1, var2, var3, var4. I want to sort on say var3 simply for display purposes. How can this be accomplished? Thanks

Dear R-helpers, I have lat/lon coordinates of regularly spaced grid points, about 4Km apart, covering the entire US continental region. I would like to mask this rectangular grid in order to extract all and only the grid points within a specific region. Today I want to extract Montana, say, from this grid, and I am hoping to somehow use the returned value of the function

Hello I'm trying to render in 3D what I usually plot by image(), or image.plot() from the library fields, followed by a map("world",add=TRUE) type of command. More concretely, I have a field of temperature values, for a given geographic area, and I would like to plot a 3D surface whose x and y axes consist of longitude and latitude values and whose height and color-coding correspond

Hi all I'm plotting colored contour maps using filled.contour. My levels are very unevenly spaced, with, say, high resolution in the small numbers but ranges that can be an order of magnitude or two larger in absolute value compared to where the action takes place. Aside from transforming the data, is there a way to control the color spacing in the key to the right of the map? Right now I

Hello I'm having trouble figuring out how to use the output of "segmented()" with a new set of predictor values. Using the example of the help file: ??set.seed(12) xx<-1:100 zz<-runif(100) yy<-2+1.5*pmax(xx-35,0)-1.5*pmax(xx-70,0)+15*pmax(zz-.5,0)+rnorm(100,0,2) dati<-data.frame(x=xx,y=yy,z=zz) out.lm<-lm(y~x,data=dati) o<-## S3

Hi I have integer values representing unix time, stored into a simple numeric vector, no POSIX* class inheritance of any kind. Something like crazylookingtimes_c(993340800.00, 993254400.00, 994809600.00) I'd like to convert them into a format like YYDDMMHHMMSS but I cannot find out how, at least not in the documentation on DateTimeClasses that I've skimmed. Can anybody help, please?

Muchas gracias por la explicación Isidro! El 10 de septiembre de 2013 08:19, Isidro Hidalgo <ihidalgo@jccm.es>escribió: > Cuando utilizas "==" estás comparando el vector "graph$NCar" (longitud 892) > con otro de longitud 2: c("160s","180s"). > Es decir, comparas el primer valor de "graph$NCar" con el primer valor de >

Estimados, tengo el siguiente problema: estoy trabajando con un data.frame que tiene 892 filas y 9 columnas, el cual adjunto para que puedan probar. Cuando consulto por las filas que cumplen una cierta condicion, por ejemplo, ... graph[graph$NCar==c("160s"),] obtengo el siguiente conjunto de datos: sol con dia NCar NIns isom area rep tipo 1 con 0.001 1 160s 0 s