similar to: (no subject)

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 all; I have a function written in R that returns as a list of values as output that has associated some user defined attributes to it. How can hide these attributes when printing the output on screen? I'm using R-2.8.1 on WinXP....it's like hiding the attr of the output from the scale function.... Thanks in advance PM

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 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

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?

>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

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

I am not certain how nlme works so I followed an example from the web ( http://www.menne-biomed.de/gastempt/gastempt1.html). I was able to successfully reproduce the example. However, when I modified my the example to use my data and with my formula, I get a set of errors having to do with the log() function. I get 10 of them (all exactly the same) and there are 10 levels in my factor variable.

Hi Duncan, On 17.09.2014 21:10, Duncan Sands wrote: > Hi Oleg, > > On 17/09/14 18:45, Oleg Ranevskyy wrote: >> Hi, >> >> Thank you for all your helpful comments. >> >> To sum up, below is the list of correct folding examples for fadd: >> (1) fadd %x, -0.0 -> %x >> (2) fadd undef, undef -> undef