search for: cmis

...er" > julian(m) < R crashes> --please do not edit the information below-- Version: platform = i386-pc-mingw32 arch = i386 os = mingw32 system = i386, mingw32 status = major = 2 minor = 2.0 year = 2005 month = 10 day = 06 svn rev = 35749 language = R Bill Venables, CMIS, CSIRO Laboratories, PO Box 120, Cleveland, Qld. 4163 AUSTRALIA Office Phone (email preferred): +61 7 3826 7251 Fax (if absolutely necessary): +61 7 3826 7304 Mobile (rarely used): +61 4 1963 4642 Home Phone: +61 7 3286 7700 mailto:Bill.Venables at...

...lythely assumed that the difference will always be a number of days and others assume it is always seconds. At the very least I think the help information should carry a big red warning about this rather unusual feature. (It may, I suppose, but I couldn't find it.) Comments? Bill Venables, CMIS, CSIRO Laboratories, PO Box 120, Cleveland, Qld. 4163 AUSTRALIA Phone: +61 7 3826 7251 Fax: +61 7 3826 7304 Mobile: +61 4 1963 4642 Home: +61 7 3286 7700 mailto:Bill.Venables@csiro.au http://www.cmis.csiro.au/bill.venables/

...y <- function(x) sin(x + 1) > eval(m) > y function(x) FUN(x+1) However the story doesn't end there. The substitution appears to have been made, even though the printed version, this time, suggests otherwise. > y(pi) [1] -0.841471 > sin(pi+1) [1] -0.841471 > Bill Venables CMIS, CSIRO Laboratories, PO Box 120, Cleveland, Qld. 4163 AUSTRALIA Office Phone (email preferred): +61 7 3826 7251 Fax (if absolutely necessary): +61 7 3826 7304 Mobile (rarely used): +61 4 1963 4642 Home Phone: +61 7 3286 7700 mailto:Bill.Venables at...

...0, Cleveland, 4163 > AUSTRALIA > Office Phone (email preferred): +61 7 3826 7251 > Fax (if absolutely necessary): +61 7 3826 7304 > Mobile: +61 4 8819 4402 > Home Phone: +61 7 3286 7700 > mailto:Bill.Venables at csiro.au > http://www.cmis.csiro.au/bill.venables/ > > -----Original Message----- > From: ripley at stats.ox.ac.uk [mailto:ripley at stats.ox.ac.uk] > Sent: Monday, 2 July 2007 7:55 PM > To: Venables, Bill (CMIS, Cleveland) > Cc: r-devel at stat.math.ethz.ch > Subject: Re: [Rd] termplot with uniform y-li...

I have a big contingency table, approximately of size 60*2*500*500, and I need to count the number of cells containing a count of 1 for each of the factors values defining the first dimension. Here is my attempt: tab1<-with(pisa1,table(CNT,GENDER,ISCOF,ISCOM)) tab2<-apply(tab1,1:4,function(x)ifelse(sum(x)==1,1,0)) tab3<-apply(tab2,1,sum) Computing tab2 is very slow. Is there a faster

...ctorily on a machine running r-2.11.1 under Suse10.3/sles. Thanks Mark Palmer Senior Statistician CSIRO Mathematics, Informatics and Statistics Phone: +61 8 9333 6293 | Fax: +61 8 9333 6121 | Mobile: 0427502353 Mark.Palmer@csiro.au<mailto:Mark.Palmer@csiro.au> | www.csiro.au | www.csiro.au/cmis<http://www.csiro.au/cmis> Address: Private bag 5, Wembley, WA 6913, Australia PLEASE NOTE The information contained in this email may be confidential or privileged. Any unauthorised use or disclosure is prohibited. If you have received this email in error, please delete it immediately and no...

...age:splines, package:methods, package:stats, package:graphics, package:grDevices, package:datasets, package:RBigData, package:mvbutils, mvb.session.info, package:tools, package:utils, package:RBigData, package:RUtilities, package:RBigLibrary, package:g.data, Autoloads, package:base Bill Venables, CMIS, CSIRO Laboratories, PO Box 120, Cleveland, Qld. 4163 AUSTRALIA Office Phone (email preferred): +61 7 3826 7251 Fax (if absolutely necessary): +61 7 3826 7304 Mobile (rarely used): +61 4 1963 4642 Home Phone: +61 7 3286 7700 mailto:Bill.Venables at...

...he call to model.frame, causing pandemonium. As a confirmation of the second change above, unless you do have (expand.dots = FALSE) the last line quoted above has no effect. Regards, Bill Venables. -- ----------------------------------------------------------------- Bill Venables, Statistician, CMIS Environmetrics Project. Physical address: Postal address: CSIRO Marine Laboratories, PO Box 120, 233 Middle St, Cleveland, Queensland Cleveland, Qld, 4163 AUSTRALIA AUSTRALIA Telephone: +61 7 3826 7251...

Hi R, I wanted to know how do we access the elements of a list. In particular, v=list(c(1,2,3,4,5),c(1,2,33,4,5),c(1,2,333,4,5),c(1,2,3333,4,5)) I want to access all the thirds items of the elements of the list. i.e., I want to access the elements, 3,33,333,3333. This can be done through sapply as: sapply(v,function(x) x[3]) But I need to access this without using

...ny authors and teachers adopting R as a standard computational vehicle of great elegance and power, free and here to stay. Well done, Ross, Robert, Kurt, Peter, Doug, Martin, Brian and whoever else I may have missed. In fact very well done. Bill Venables. -- Bill Venables, Statistician, CMIS Environmetrics Project CSIRO Marine Labs, PO Box 120, Cleveland, Qld, AUSTRALIA. 4163 Tel: +61 7 3826 7251 Email: Bill.Venables at cmis.csiro.au Fax: +61 7 3826 7304 http://www.cmis.csiro.au/bill.venables/ -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-....

Full_Name: David Duffy Version: 1.0.1 OS: Linux Submission from: (NULL) (152.98.96.21) Documentation for diag() states that "If x is a vector (or a 1-d array) then diag(x) returns a diagonal matrix whose diagonal is x." > m<-matrix(1:4,ncol=1) > diag(m) [1] 1 diag(as.vector(m)) gives the advertised performance. Fix: Change documentation.

...ubstitute(). > >> > >> ? > >> > >> Bill Venables. > >> > >> -----Original Message----- > >> From: Peter Dalgaard [mailto:p.dalgaard at biostat.ku.dk] > >> Sent: Friday, 22 December 2006 9:47 PM > >> To: Venables, Bill (CMIS, Cleveland) > >> Cc: r-devel at stat.math.ethz.ch; R-bugs at biostat.ku.dk > >> Subject: Re: [Rd] substitute creates an object which prints > >> incorrectly > >> (PR#9427) > >> > >> Bill.Venables at csiro.au wrote: > >> > >&gt...

Dear all, I'm currently use R v0.99 on Windows 98 Second Edition. I have a question on some simple calculations. I wonder that I've done something wrong with the calculation. Here is the imput commands: > a <- 25.01 > b <- 56.08 > a.trunc <- trunc(a) > b.trunc <- trunc(b) > a.tail <- a - a.trunc > b.tail <- b - b.trunc > a.trunc [1] 25 >

Full_Name: Glenn Stone Version: 1.5.1 and 1.6.0 OS: win2000 Submission from: (NULL) (168.140.227.9) The following code produces incorrect fitted values in version 1.5.1 and an error in 1.6.0 Error in "contrasts<-"(*tmp*, value = "contr.treatment") : contrasts apply only to factors In addition: Warning message: variable ihalf is not a factor in:

When doing polynomial regression I believe it is a good idea to use the poly function to generate orthogonal polynomials. When doing this in Splus there is a handy function (transform.poly I think) to convert the coefficients produced by regression with the poly function back to the original scale. Has somebody written something similar for R ? Robert

Hello, I have a date in the format Year-Month Name (e.g. 1990-January) and R classes it as a character. I want to convert this character into a date format, but when I try as.Date(1990-January, "%Y-%B"), I get back NA. The function strptime also gives me NA back. Thanks. [[alternative HTML version deleted]]

For anybody who uses ESS with R, how do you invoke the vsize and nsize options when you call R. I can't find any appropriate variables from an apropos. Thanks, Jord -- Jordan Howarth CSIRO Mathematical and Information Sciences mailto:jordan.howarth at cmis.csiro.au -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "body", not the subject !) To: r-help-reques...

...R-help and I will write more tetchy little notes like this one... :-) Bill Venables. -- Bill Venables, Statistician Tel. +61 7 3826 7251 CSIRO Marine Laboratories, Fax. +61 7 3826 7304 Cleveland, Qld, 4163 Email: Bill.Venables at cmis.csiro.au AUSTRALIA http://www.cmis.csiro.au/bill.venables/ -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subsc...

> On Wed, 2 Feb 2000, Adriane Leal wrote: > > > I'd like to perform a box-cox transformation to a data set and also plot > > lambda versus L(lambda) using R. Does anybody knows how can I do such a > > thing? gnlr3 in my gnlm library does both linear and nonlinear models with Box-Cox transformation. However, it is somewhat nonstandard as it renormalizes to obtain a

...t...but still can't figure out how to do it. > ..Thanks in advance. No, it is not that difficult at all, but unless you are prepared to look at the documentation and spend a little time studying the examples you are going to find it very tough going. Bill Venables. -- Bill Venables, CSIRO/CMIS Environmetrics Project Email: Bill.Venables at cmis.csiro.au Phone: +61 7 3826 7251 Fax: +61 7 3826 7304 Postal: PO Box 120, Cleveland, Qld 4163, AUSTRALIA -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~horni...