similar to: Tinc on Mac OSX 10.8.4.

Hello, I'm trying to develop a box plot of time series data to look at the range in the data values over the entire period of record. My data initially starts out as a list of hourly data, and then I've been using this code to make this data into the final ts array. # Read in the station list stn.list <- read.csv("/home/kbennett/fews/stnlist3", as.is=T, header=F) # Read in

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

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

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

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?

>From: "Ivaha C (AT)" <civaha at glam.ac.uk> >Date: 2008/04/10 Thu AM 08:51:14 CDT >To: R Help <r-help at r-project.org> >Subject: [R] Structural Modelling in R-project if you have a univariate time series and you want to break it into its various components, then you get the scalars based on a decomposition. if you have a kalman filter/ basic strucutural

Hello, I ordered the Devel lite kit, and installed it. I am just trying to get the FXO port to work, and am having trouble. To load the card I do the following. modprobe wcfxs modprobe wcfxo ztcfg -vv asterisk -vc My /var/log/asterisk/messages show Nov 3 11:03:39 WARNING[3317]: zt hook failed: Device or resource busy Here is my /etc/zaptel.conf fxoks=1 fxsks=4 loadzone=us defaultzone=us

Hi List I have a series of zip archives each containing several files. One of these files is called goCats.txt and I would like to read it into R from the archive. It's a simple tab delimited text file. pathToZip <- '/home/iain/Documents/Work/Results/bovineMacRNAData/deAnalysis/afInfection/commonNorm/twoHrs/af2hrs.zip' z <- unz(pathToZip, 'goCats.txt', 'r') zT

I'm resending this message because I did not include a subject line in my first posting. Apologies for the inconvenience! Tanja > Hello, > > I'm trying to fit a generalized linear mixed model to estimate diabetes prevalence at US county level. To do this I'm using the glmer() function in package lme4. I can fit relatively simple models (i.e. few covariates) but when

Hi I have a panel dataset with large number of groups and differing number of observations for each group. I want to randomly select say, 20% of the groups or 200 groups, but along with all observations from the selcted groups (with the corresponding data). I guess it is possible to generate a random sample from the groups ids and then match that with the entire dataset to have the intended

I have a question to everybody. After log10 transfprmation, I have done ordinary kriging in gstat in R? I need to back trnasform the prediction values to orgiginal scale. How do I do this in gstat in R? Thanks Zia -- Zia Uddin Ahmed 915 Brad Field Hall Department of Crop and Soil Cornell University Ithaca NY 14850 USA

Hello, we are setting up another system that will run either 1.2.4, the latest version of 1.2 or 1.4. We have not yet decided on the version. Anyhow, this is a higher volume system (dual processor) which will handle 30-50 simultaneous calls with 60 to 100 simultaneous channels lit up. Most calls are g711 with very little g729 and a little gsm mixed in. We have a similar system doing exactly

Hi, Apologies if my problem sound silly! I've been modifying the status.xsl , and now stuck with the problem. I want to disply on the top header (h2) total number of sources mounted on the server . I dont know how/where the webserver is working, thus i dont know how/where to get the total count (if its already stored somewhere at all). Otherwise, what is the easiest way to get the count?