similar to: R-beta: (0+0i)^2

Why is 19 March, 1947 a little longer than one day? x <- as.POSIXct("1947-04-16") julian(x, origin = as.POSIXct("1947-03-20")) Time difference of 27 days julian(x, origin = as.POSIXct("1947-03-19")) Time difference of 28.04167 days > julian(x, origin = as.POSIXct("1947-03-18")) Time difference of 29.04167 days I am running R-1.7.1 compiled on

Some of you may have seen the pythag() part in the R API definition in "Writing R Extensions" (source = doc/manual/R-exts.texinfo). or followed the report and Prof. Brian Ripley's answer about pythag()'s availability from R's binary. As we say in above manual >> `pythag(A, B)' computes `sqrt(A^2 + B^2)' without overflow or >> destructive

Hi,I am working with a custom built (LFS-based) Linux distro and am attempting to compile Samba4. Samba3 has always compiled without issue. The system has its Python installed in /opt/python2. There are symlinks to put "python" in /usr/local/bin and the "python2.7" lib folder in /usr/include. The configure command used was simply ./configure --enable-fhs --with-quotas The

Hello R-helpers, I was wondering if anyone is porting S lib packages for design of experiments: conf.design, glim, graff(all from statlib)? I was able to change the format of S packages into R (including docs). In conf.design, by W. Venables, sort.list function is missing from R. I just put sort.list <- sort, but I got more error messages. Before doing more work, would like to know if this

Hello R-helpers, I was wondering if anyone is porting S lib packages for design of experiments: conf.design, glim, graff(all from statlib)? I was able to change the format of S packages into R (including docs). In conf.design, by W. Venables, sort.list function is missing from R. I just put sort.list <- sort, but I got more error messages. Before doing more work, would like to know if this

Hi all, I post it on both r-help and r-finance since I don't know where is most appropriate for this topic. Sorry if it bothers you. I did garch fitting on S&P500 monthly returns with garchFit from fSeries. I got same coefficients from all cond.dist except normal. I thought that is probabaly usual for the data. But when I play with it, I got another question. I plot skew normal with

Consider the following scrap of code: > x<- ts(1:50,start=c(1,11),freq=12) > y <- aggregate(x,nfreq=4) > c(y) [1] 6 15 24 33 42 51 60 69 78 87 96 105 114 123 132 141 > y Error in rep.int("", start.pad) : invalid number of copies in rep.int() > tsp(y) [1] 1.833333 5.583333 4.000000 So we can aggregate into quarters, but we cannot print it using

Consider the following scrap of code: > x<- ts(1:50,start=c(1,11),freq=12) > y <- aggregate(x,nfreq=4) > c(y) [1] 6 15 24 33 42 51 60 69 78 87 96 105 114 123 132 141 > y Error in rep.int("", start.pad) : invalid number of copies in rep.int() > tsp(y) [1] 1.833333 5.583333 4.000000 So we can aggregate into quarters, but we cannot print it using

Hi ?atan2 says that atan2(y,x)=atan(y/x) for x and y numeric or complex vectors. Well, I would expect atan2(1,1i) to be equal to atan(-1i), but > atan2(1,1i) Error in atan2(y, x) : Non-numeric argument to mathematical function > R.version _ platform powerpc-apple-darwin8.5.0 arch powerpc os darwin8.5.0 system powerpc, darwin8.5.0

Hi ?atan2 says that atan2(y,x)=atan(y/x) for x and y numeric or complex vectors. Well, I would expect atan2(1,1i) to be equal to atan(-1i), but > atan2(1,1i) Error in atan2(y, x) : Non-numeric argument to mathematical function > R.version _ platform powerpc-apple-darwin8.5.0 arch powerpc os darwin8.5.0 system powerpc, darwin8.5.0

Dear friends, After running the lm() model, we can get summary resluts like the following: Coefficients: Estimate Std. Error t value Pr(>|t|) x1 0.11562 0.10994 1.052 0.2957 x2 -0.13879 0.09674 -1.435 0.1548 x3 0.01051 0.09862 0.107 0.9153 x4 0.14183 0.08471 1.674 0.0975 . x5 0.18995 0.10482 1.812 0.0732 . x6 0.24832 0.10059 2.469 0.0154 * x7

Hi. Background - I am working with a dataset involving around 750K observations, where many of the variables (8/11) are unordered factors. The typical model used to model this relationship in the literature has been a simple linear additive model, but this is rejected out of hand by the data. I was asked to model this via kernel methods, but first wanted to play with the parametric

Today I installed the newest R develepment branch R version 2.14.0 Under development (unstable) (2011-03-18 r54866) on FreeBSD 9.0-CURRENT (amd64). All seems fine so far. After that I updated my R packages with option 'checkBuilt=TRUE'. There are four packages (spatstat, pscl, adehabitatLT, adehabitatHR) which gives an error like this: Error in dyn.load(file, DLLpath = DLLpath, ...) :

Dear All I hope my question is relevant on this forum, else very sorry for the disturbance I want to simply plot a bode diagram of a siso model using the 'bodeplot' command in the control package fail due to error in 'issiso' evaluation... Error in if (issiso(sys)) {: then I try by myself some retro-engineering :) the problem is when I use H <- freqresp(syst, w) mag =

Hello, I am new to non-linear growth modelling in R and I am trying to reproduce an analysis that was done (successfully) in S-Plus. I have a simple non-linear growth model, with no nesting. I have attempted to simplify the call as much as possible (by creating another grouped object, instead of using subset= and compacting the fixed and random expressions.) This is a what the grouped

Deemed colleagues I would appreciate your help with a sentence to calculate the linear distance between 2 geographical points (coordinates in UTM). In advance thnks for your attention, -- Ricardo Bandin Llanos rbandin@udec.cl Estudiante - Magíster Cs. m. Pesquerías Universidad de Concepción, Región del Bio-Bio, Chile Celular: (0056-41) 97949957 [[alternative HTML version deleted]]

Hello, I am trying to perform a logistic regression using counts. For example: cedegren <- read.table("http://www.cloudstat.org/index.php?do=/attachment/download/id_95 /", header=T) attach(cedegren) ced.del <- cbind(sDel, sNoDel) ced.logr <- glm(ced.del ~ cat + follows + factor(class), family=binomial("logit")) This works. However, if I change the family to

Dear All, I would like to be able to plot histograms/densities on a semi-log or log-log scale. I found several suggestions online http://tolstoy.newcastle.edu.au/R/help/05/09/12044.html https://stat.ethz.ch/pipermail/r-help/2002-June/022295.html http://www.harding.edu/fmccown/R/#histograms Now, consider the code snippet taken from http://www.harding.edu/fmccown/R/#histograms # Get a random

https://bugzilla.mindrot.org/show_bug.cgi?id=2085 Bug ID: 2085 Summary: fd_mask and howmany are not defined in Android Classification: Unclassified Product: Portable OpenSSH Version: 6.2p1 Hardware: All OS: Linux Status: NEW Severity: minor Priority: P5 Component: Build system

To win*/NT users: I just recieved a mail from r-help with a virus attachment. I am on an unix system but our mail server deleted the message anyways. here it is: Received: by stat.math.ethz.ch (8.9.1/8.9.1) id RAA27614 for r-help-gang-use; Wed, 10 Jan 2001 17:01:56 +0100 (MET) Received: (from daemon at localhost) by stat.math.ethz.ch (8.9.1/8.9.1) id RAA27412 for