similar to: Timeseries Data Plotted as Monthly Boxplots

I create a zooreg object that runs from Jan-1-2002 0:00 to Jun-1-2005 0:00... regts.start = ISOdatetime(2002, 1, 1, hour=0, min=0, sec=0, tz="") regts.end = ISOdatetime(2005, 6, 1, hour=0, min=0, sec=0, tz="") regts.zoo <- zooreg( NA, regts.start, regts.end, deltat=3600 ) Upon inspection: > regts.zoo[1:3] 2002-01-01 00:00:00 2002-01-01 01:00:00 2002-01-01 02:00:00

Hi, I have a dataset called 'fish'. fish$Species returns extract 1. When I use fish$Species != c("CRA","PHC"), i.e. I want all species except "CRA" and "PHC", I get extract 2 which is blatantly wrong. Can anyone see what I'm doing wrong? Regards, James EXTRACT 1 > fish$Species [1] ALB ALB ALB ALB ALB ALB ALB ALB ALB ALB ALB ALB ALB

I am new to the world of R/programming so this may be a really easy question. I thank you for your patience and help in advance I would like the characters km^2 to be displayed on the plot subtitle as km squared - two as a superscript. I would also like to have the numbers from the data set for longitude and latitude to be rounded to four decimal places. Thank you. plot (

I KNOW this should be easy, but I'm stuck. My data frame consists of multiple observations from each of a number of stations, and what I would like to do is create another data frame that contains all the variables of the first, but only rows where a certain variable is at its maximum for the station. So, for example: > my.df stn obs v 1 1 1 0.26400396 2 2 1

Hello, I using the fuction auto.arima() from package forecast to predict the values of p,d,q and P,D,Q. My problem is the execution time of this function, for example, a time series with 2323 values with seasonality to the week take over 8 hours to execute all the possibilities. I using a computer with Windows XP, a processor Intel Core2 Duo T7300 and 2Gb of RAM.

Hi all, I am trying to download some tab separated data from the internet. The data is not available directly at the URL that could be known apriori. There is an intermediate form where start and end dates have to be given to get to the required page. For example, I want to download data for a station 03015795. The form for this station is at:

Dear Helpers, I wrote a simple function to standardise variables if they contain more than one value. If the elements of the variable are all identical, then I want the function to return zero. When I submit variables whose elements are all identical to the function, it returns not zero, but NaNs. zt=function(x){if (length(table(x)>1)) y=(x-mean(x))/sd(x) else if (length(table(x)==1)) y=0;

Dear Helpers, I wrote a simple function to standardise variables if they contain more than one value. If the elements of the variable are all identical, then I want the function to return zero. When I submit variables whose elements are all identical to the function, it returns not zero, but NaNs. zt=function(x){if (length(table(x)>1)) y=(x-mean(x))/sd(x) else if (length(table(x)==1)) y=0;

Dear R users, When running the program below I receive the following error message: fit <- optim(parm, objective, yt = tyield, hessian = TRUE) Error in as.vector(data) : no method for coercing this S4 class to a vector I can't figure out what the problem is exactly. I imagine that it has something to do with "tyield" being a matrix. Any help on explaining what's going on

Hi, I would like to know what should I garantee about P and GGt in order to have F = Z %*% P %*% t(Z) + GGt always as a positive definite matrix. Being more precise: I am trying to find minimum likelihood parameters by using the function 'optim' to find the lowest value generated by $LogLik from the function 'fkf' (http://127.0.0.1:27262/library/FKF/html/fkf.html). The

Hello All, I trying to use the function auto.arima(....) from package forecast but I have a problem. My steps after I used the function auto.arima(...) I create the time series like this: >bbrass = scan("C:/Program Files/R/data PTIN/my_file.dat") >regts.start = ISOdatetime(2006, 7, 1, hour=0, min=0, sec=0, tz="GMT") #2006 07 01 00 >regts.end = ISOdatetime(2006, 7,

Hi, I'm trying to estimate the parameters of a state space model of the following form measurement eq: z_t = a + b*y_t + eps_t transition eq y_t+h = (I -exp(-hL))theta + exp(-hL)y_t+ eta_{t+h}. The problem is that the distribution of the innovations of the transition equation depend on the previous value of the state variable. To be exact: y_t|y_{t-1} ~N(mu, Q_t) where Q is a diagonal

Hello, I would like to sync my Linux firewall to my Mac OS X. The Linux is a Debian 4.0 with ext3 filesystem (UTF-8) $ rsync --version rsync version 3.0.2 protocol version 30 The Mac OS X is a 10.5 Leopard with HFS+ (UTF-8) $ /sw/bin/rsync --version rsync version 3.0.3 protocol version 30 (fink version) I have many accents in my filenames: $ ls -l | tee filelist total 12 -rw-r----- 1 tamas

Thanks for everyone's help with zoo -- I think I've got my data set ready. (The data consists of surface weather temperatures, from 2002 to 2005, one observation per hour. Some values are missing... i.e. NA) I have three goals: GOAL #1:Get the data in proper time series form, preserving frequency information: > w4.ts <- as.ts( w3.zoo, frequency=(1/3600) ) I hope that 1/3600

I've been working with an adtran ta 850 hooked to a t100p pretty much all day today, and I haven't gotten past configuring zaptel.conf and zapata.conf. For some reason, when I pick up analog phone hooked up to the first module of a quad fxs card in the second slot of the ta 850, asterisk thinks that all four of the fxs modules in that card are going off hook. If I pick up a phone hooked

I am changing a system from zap to DAHDI. I removed everything zap. when doing the command: sh -x /etc/init.d/dahdi start, I see initlog -q -c 'modprobe wct4xxp' sh: /sbin/ztcfg: No such file or directory FATAL: Error running install command for wct4xxp doing updatedb then, locate zap returns me /var/lib/asterisk/sounds/spy-zap.gsm /usr/share/groff/1.18.1.1/font/devps/zapfdr.pfa

Hi, I'm trying to setup tunnel between a Linux VM (running on a KVM) and my Mac 10.8.4. On the Linux box i am able to start tincd and get interface up for my VPN name, but on the Mac, when i run i get this error ?1 tinc/zt ? sudo /opt/local/sbin/tincd -c /opt/local/etc/tinc/zt -D --debug=4 tincd 1.0.21 (Sep 21 2013 00:13:28) starting, debug level 4 Could not open /dev/tun0: No such file or

Dear r-helpers, Spencer Graves and Manual Morales proposed the following methods to simulate p-values in lme4: ************preliminary************ require(lme4) require(MASS) summary(glm(y ~ lbase*trt + lage + V4, family = poisson, data = epil), cor = FALSE) epil2 <- epil[epil$period == 1, ] epil2["period"] <- rep(0, 59); epil2["y"] <- epil2["base"]

Hello list, I've been working on a package of my own. It works fine on the 32bit linux machines that I've tested it on. Someone using 64bit ubuntu wanted to try it, and surprising segfaults happened. My test code results in no segfault, errors, or leaks from my calls when run under valgrind (I recompiled R with the level 2 valgrind instruments). R and packages are compiled from source,

We (the lme4 authors) are having a problem with doing a proper deep copy of a reference class object in recent versions of R-devel with the LAZY_DUPLICATE_OK flag in src/main/bind.c enabled. Apologies in advance for any improper terminology. TL;DR Is there an elegant way to force non-lazy/deep copying in our case? Is anyone else using reference classes with a field that is an external pointer?