[R] [R-pkgs] new package on CRAN

2006-05-31 Thread Robin Hankin
Dear List Please find uploaded to CRAN a new package, partitions that enumerates integer partitions using C. A paper describing the package has recently been accepted by JSS (volume 16, codesnippet 1). Integer partitions arise in many branches of combinatorics, and the partition function P(n)

Re: [R] Automate concatenation?

2006-05-31 Thread Martin Maechler
Paul == Paul Roebuck [EMAIL PROTECTED] on Tue, 30 May 2006 16:04:14 -0500 (CDT) writes: Paul On Tue, 30 May 2006, Robert Lundqvist wrote: I have this typical problem of joining a number of vectors with similar names - a1, a2,..., a10 - which should be concatenated into one.

[R] Message deleted

2006-05-31 Thread Mail Admin
Scanner: MailMonitor for SMTP v1.2.2 Problem description: Email data: MessageID: [EMAIL PROTECTED] From: r-help@stat.math.ethz.ch To: [EMAIL PROTECTED] Cc: Subject: Returned mail: see transcript for details Scanning part [] Scanning part [letter.zip] Attachment validity check: passed. Virus

Re: [R] merging

2006-05-31 Thread Gavin Simpson
On Tue, 2006-05-30 at 15:09 -0500, Sundar Dorai-Raj wrote: Gavin Simpson wrote: Dear List, snip / Will this help: rbind.all - function(...) { x - list(...) cn - unique(unlist(lapply(x, colnames))) for(i in seq(along = x)) { if(any(m - !cn %in% colnames(x[[i]]))) {

Re: [R] How to choose columns in data.frame by parts of columns' names?

2006-05-31 Thread Gabor Grothendieck
Use rowMeans: rowMeans(iris[, regexpr(Sepal, colnames(iris)) 0]) gives the means of each row of iris using only the columns whose names contain Sepal. On 5/31/06, Guo Wei-Wei [EMAIL PROTECTED] wrote: Gabor and Peter, Thank you. Both of you give me excellent ways. I have a further problem.

[R] R search plug-in

2006-05-31 Thread arin basu
Hi All: I have written an R search plug-in based on Prof Jonathan Baron's R search site. This Firefox/Netscape plug-in lets you search R sites and R discussion group archives directly from your Mozilla (Seamonkey) or Firefox or Netscape browser. Here is the link to the installer for the

[R] Group Method of Data Handling

2006-05-31 Thread Mike Bocharov
Hello, folk! I need Group Method of Data Handling realization for R... And I can't find it in cran.r-project.org, Do you know it(Group Method of Data Handling) exists for R ? ps Sorry for my english __ R-help@stat.math.ethz.ch mailing list

Re: [R] R search plug-in

2006-05-31 Thread Sundar Dorai-Raj
Hi, Arin, This was already done: http://finzi.psych.upenn.edu/R/Rhelp02a/archive/53597.html And works well. Try firefox search plugin. --sundar arin basu wrote: Hi All: I have written an R search plug-in based on Prof Jonathan Baron's R search site. This Firefox/Netscape plug-in lets

Re: [R] Bug of pmvnorm?

2006-05-31 Thread Torsten Hothorn
On Sat, 27 May 2006, P Ehlers wrote: Looks like a bug (in mvt()?). Note that pmvnorm(lower = c(-Inf, -Inf), upper = c(Inf, 10)) works as expected, as does replacing any of the 4 Infs with a finite value. Function sadmvn() in pkg:mnormt does give 1 with Inf lower/uppers. yes, looks like

Re: [R] Running R

2006-05-31 Thread Pramod Anugu
I have downloaded R-2.3.0.tar.gz for Linux 64 bit processor 1. Untar tar -zxvf R-2.3.0.tar.gz 2. changed the directory to the newly created directory R-2.3.0 3. Typed ./configure 4. Typed make 5. Copied the Script File R from directory R-2.3.0/bin/R to /usr/local/bin/R. 6. Typed R Fatal error:

Re: [R] Running R

2006-05-31 Thread Andrew Perrin
At step 5, don't copy the file. Instead, do: make install -- Andrew J Perrin - andrew_perrin (at) unc.edu - http://perrin.socsci.unc.edu Assistant Professor of Sociology; Book Review Editor, _Social Forces_ University of North

Re: [R] Running R

2006-05-31 Thread Jonathan Baron
On 05/31/06 09:25, Pramod Anugu wrote: I have downloaded R-2.3.0.tar.gz for Linux 64 bit processor 1. Untar tar -zxvf R-2.3.0.tar.gz 2. changed the directory to the newly created directory R-2.3.0 3. Typed ./configure 4. Typed make 5. Copied the Script File R from directory R-2.3.0/bin/R to

Re: [R] Running R

2006-05-31 Thread Pramod Anugu
Still the same problem after changing the step 5 I have downloaded R-2.3.0.tar.gz for Linux 64 bit processor 1. Untar tar -zxvf R-2.3.0.tar.gz 2. changed the directory to the newly created directory R-2.3.0 3. Typed ./configure 4. Typed make 5. Typed make install 6. Typed R Fatal error:

[R] Help on glmmPQL

2006-05-31 Thread Francisco Redelico
Dear R-Users, I am trying to fit a generalized linear mixed models using glmmPQL, mymodel.glmm = glmmPQL(efect~ time+ cov1+ cov2, random = ~ fact1+ fact2+ fact3+ fact1*time + fact2*time+ fact3*- cov1 - cov2- time| subject, family=poisson) I get this error message: Error

Re: [R] Running R

2006-05-31 Thread Liaw, Andy
It's always a good idea to do make check or even make check-all before make install. If the build does not pass the checks, there's no point in installing. (I've run into successful builds that failed the checks several times before.) Andy From: Pramod Anugu Still the same problem

Re: [R] Running R

2006-05-31 Thread Jonathan Baron
On 05/31/06 09:39, Pramod Anugu wrote: 5. Typed make install 6. Typed R Fatal error: unable to open the base package I get the error message. Please advice Other things to try (aside from make check): 1. Make sure you are root before you say make install. 2. Check to see that the

Re: [R] Running R

2006-05-31 Thread Pramod Anugu
Still the same problem after the changes I have downloaded R-2.3.0.tar.gz for Linux 64 bit processor 1. tar -zxvf R-2.3.0.tar.gz 2. changed the directory to the newly created directory R-2.3.0 3. Typed ./configure 4. Typed make 5. Make check 6. make check-all 7. Typed make install 8. Typed R

[R] interpolating a lot of data at once

2006-05-31 Thread David Hugh-Jones
I have a big dataset containing a lot of values for 1970, 1980 and 1990. I want to interpolate values for the years in between, and also if possible to extrapolate to 1968 and 1969. The method doesn't have to be clever but I am looking for a function that will do all the data at once. (Doing

[R] how to produce this pic in R

2006-05-31 Thread Lisa Wang
Hello Everyone, I would like to produce a 3d par chart like the attached. How do I do it in R? thanks a lot Lisa Wang Princess Margaret hospital Toronto,Ca __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE

Re: [R] Running R [Broadcast]

2006-05-31 Thread Liaw, Andy
In the R-2.3.0 directory, what do you get if you type bin/R? Does R starts up OK, or do you get an error? What were the last few lines of output from running configure? Andy _ From: Pramod Anugu [mailto:[EMAIL PROTECTED] Sent: Wednesday, May 31, 2006 11:23 AM To: Liaw, Andy;

[R] FW: Graphical Limitiations of the (D)COM Server AxSTATCONNECTORCLNTLib.AxStatConnectorGraphicsDevice

2006-05-31 Thread Weilbacher, John
A sample picture of this issue: http://hera.pmcl.com/DMap/BeachFX/Bug.gif From: Weilbacher, John Sent: Wednesday, May 31, 2006 10:46 AM To: 'r-help@stat.math.ethz.ch' Subject: Graphical Limitiations of the (D)COM Server

[R] Nesting in Cox proportional hazards survivorship analysis

2006-05-31 Thread Jessica M Pearce
Hello, My advisor and I have been working on some survivorship analyses in R and we are hoping to get some feedback on a particular issue involving nesting. We are interested in patterns of food discovery by ant species. Our observations consist of time to discovery by an ant for three

Re: [R] Running R

2006-05-31 Thread Barry Rowlingson
Jonathan Baron wrote: 5. Copied the Script File R from directory R-2.3.0/bin/R to /usr/local/bin/R. 6. Typed R The instructions say that you should use make install unless you want to run R from the directory into which you unpacked it. Another idea is to have a symlink (here

[R] Graphical Limitiations of the (D)COM Server AxSTATCONNECTORCLNTLib.AxStatConnectorGraphicsDevice

2006-05-31 Thread Weilbacher, John
I am having issues displaying Trellis Graphs using the AxSTATCONNECTORCLNTLib.AxStatConnectorGraphicsDevice. The graph displays correctly in a default R window, but will not draw labels or data in the GraphicsDevice. Is there a workaround for this, or am I doing something incorrect. I can get

[R] 3d par chart

2006-05-31 Thread Lisa Wang
Hello Everyone, I would like to produce a 3d par chart like the attached. How do I do it in R? thanks a lot Lisa Wang Princess Margaret hospital Toronto,Ca__ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do

[R] 3D graph - automated way

2006-05-31 Thread Rado Bonk
Dear R-users, On behalf of my colleague I would like to ask you for a help: I have to produce a 3D graphs from several two-column files. Each file has around 20-30k of lines. EXAMPLE OF TWO_COLUMNS_FILE - col_names:(EPS, UPS) 8 725 8 825 3 1025 10 9950 9 1 8

Re: [R] interpolating a lot of data at once

2006-05-31 Thread Gabor Grothendieck
On 5/31/06, David Hugh-Jones [EMAIL PROTECTED] wrote: I have a big dataset containing a lot of values for 1970, 1980 and 1990. I want to interpolate values for the years in between, and also if possible to extrapolate to 1968 and 1969. The method doesn't have to be clever but I am looking for

[R] X11 device

2006-05-31 Thread Marius Hofert
Hi, I am working with R on a Mac Powerbook G4. As described a couple of days ago, I would like to create an animation consisting of - say - 1000 plots (lattice: xyplot, wireframe etc.). Everything works perfectly fine (thanks to all of you for the suggestions and ideas), but before I am

[R] (what I think are) two easy questions

2006-05-31 Thread Rob Tirrell
Hello, I have a table of data (attachment) from which I am trying to generate a graph. I import the data into a vector (is that the right term? matrix maybe), which seems to work fine except there is an X prepended to each date. This is my first problem. My second is how would I make a graph

[R] plot problem in version 2.3.0 with dates as x

2006-05-31 Thread jun . zang
Hi, I have some old programs running well for older versions 2.2.0 and 2.1.1, but cannot run under 2.3.0. The problem is in plot with dates as x: library(survival) test - data.frame( x = c(1:20), date = as.date(1:20)) plot(test$date, test$x)

[R] trouble with boot ()

2006-05-31 Thread Rogério Rosa da Silva
Dear members, I am trying to use boot () to compute the distributions of a statistic of a data set. The statistic is defined in the following code: eds-function(x) { r-cor(x) paren-1-abs(r)/2 denom-sum(sum(paren)+0.5)

[R] Is this my mistake or a bug?(K-S test and EVT)

2006-05-31 Thread Mu Tian
Hi, I was doing a Kolmogrov test on x, y shown below. They are both 151 long. According to the help file, exact p-value is not available so I set exact = FALSE, but still got the warning. There are no duplicated values in either X or Y. Thank you for your tips. BTW, is there functions in R

[R] follow-up

2006-05-31 Thread Rob Tirrell
Hi, I've been able to get the data _of every individual_ into a graph (there are 16 lines), but I have a few more questions... How do I neatly get the average weights of each group (PC3/CO, PC3/FO, FADU/CO and FADU/FO) and make a graph of four lines. How do I exclude the first column (which is

Re: [R] trouble with boot ()

2006-05-31 Thread Sundar Dorai-Raj
Rogério Rosa da Silva wrote: Dear members, I am trying to use boot () to compute the distributions of a statistic of a data set. The statistic is defined in the following code: eds-function(x) { r-cor(x) paren-1-abs(r)/2

Re: [R] X11 device

2006-05-31 Thread tom wright
Marius, Not quite what your asking but try ?bitmap. This can create png images without requiring X On Wed, 2006-31-05 at 20:55 +0200, Marius Hofert wrote: Hi, I am working with R on a Mac Powerbook G4. As described a couple of days ago, I would like to create an animation consisting of -

[R] follow-up

2006-05-31 Thread Rob Tirrell
Hi, I've been able to get the data _of every individual_ into a graph (there are 16 lines), but I have a few more questions... How do I neatly get the average weights of each group (PC3/CO, PC3/FO, FADU/CO and FADU/FO) and make a graph of four lines. How do I exclude the first column (which is

[R] Installing RMySQL on Solaris

2006-05-31 Thread Ike Nnabugwu
% sudo R CMD INSTALL RMySQL_0.5-7.tar.gz Password: * Installing *source* package 'RMySQL' ... creating cache ./config.cache checking how to run the C preprocessor... cc -E checking for compress in -lz... no checking for getopt_long in -lc... no checking for mysql_init in -lmysqlclient... no

Re: [R] Running R

2006-05-31 Thread Pramod Anugu
Under the below seection it talks about X86_64 archetecture for Linux. 1. tar -zxvf R-2.3.0.tar.gz 2. changed the directory to the newly created directory R-2.3.0 3. Typed ./configure 'r_arch=name' 4. Typed make 5. Make check 6. make check-all 7. Typed make install 8. Typed R I still get the same

Re: [R] Is this my mistake or a bug?(K-S test and EVT)

2006-05-31 Thread Peter Ehlers
Mu Tian wrote: Hi, I was doing a Kolmogrov test on x, y shown below. They are both 151 long. According to the help file, exact p-value is not available so I set exact = FALSE, but still got the warning. There are no duplicated values in either X or Y. But I'll wager that there are some

Re: [R] Installing RMySQL on Solaris

2006-05-31 Thread Prof Brian Ripley
Install libz from www.zlib.net. Note that it is not finding your MySQL libs either, so you need to specify a path to the those and may be it will find libz in the same place. It is even possible that you do not have a developer installation of MySQL. `Solaris' is vague, and MySQL is not a

[R] timeSeq and TimeDate analog in R ?

2006-05-31 Thread markleeds
Hi All : I am attempting tomove a large amount of code from Splus to R and I was hoping that there was an equivalent in R of the Splus functions timeSeq and timeDate ? I did an RSiteSearch but nothing came up ? If the equivalent functions are part of some package, that's fine. Thanks a lot.

Re: [R] sib TDT transmission/disequilibrium test

2006-05-31 Thread David Duffy
Farrel Buchinsky [EMAIL PROTECTED] wrote: Does anyone know if the sib TDT has been implemented in R 1. Spielman, R.S., and Ewens, W.J. (1998) A sibship test for linkage in the presence of association: the sib transmission/disequilibrium test. Am J Hum Genet 62, 450-458 It can be done as a

[R] a problem 'cor' function

2006-05-31 Thread Tao Shi
Hi list, One of my co-workers found this problem with 'cor' in his code and I confirm it too (see below). He's using R 2.2.1 under Win 2K and I'm using R 2.3.0 under Win XP. === R.Version() $platform [1] i386-pc-mingw32 $arch [1] i386 $os [1] mingw32

Re: [R] X11 device

2006-05-31 Thread Paul Roebuck
On Wed, 31 May 2006, Marius Hofert wrote: I am working with R on a Mac Powerbook G4. As described a couple of days ago, I would like to create an animation consisting of - say - 1000 plots (lattice: xyplot, wireframe etc.). Everything works perfectly fine (thanks to all of you for the

Re: [R] Nesting in Cox proportional hazards survivorship analysis

2006-05-31 Thread Greg Snow
The site/station syntax is mainly useful for situations where you have the same station id's withing different sites (e.g. there is a station #1 in site #1 and also a different station still labeled #1 within site #2). It is unclear whether you really need the /station part or not (it generally

Re: [R] follow-up

2006-05-31 Thread Greg Snow
To get the averages, look at ?tapply or ?aggregate To exclude parts of the data look at ?subset To remove the spaces in the lines around points look at ?matplot and ?plot specifically look for type='o' To add a legend look at ?legend. You may also want to look at the lattice package, using

Re: [R] a problem 'cor' function

2006-05-31 Thread Peter Ehlers
Looks like another case of the most F of all FAQs: FAQ 7.31. See if the following makes sense to you: pl - iris[101:150, 3] all.equal(cor(pl,pl), 1) [1] TRUE cor(pl,pl) == 1 [1] FALSE sprintf(%1.22g, cor(pl, pl)) [1] 0.99989 sprintf(%1.22g, pl) [1] 6

Re: [R] plot problem in version 2.3.0 with dates as x

2006-05-31 Thread Prof Brian Ripley
Note, this is about plot.date (not dates, a different class) in the contributed survival package, not R 2.3.0. You need to use log= as plot.date incorrectly uses as its default log=NULL. Please see the posting guide about how to report a problem in a contributed package. On Wed, 31 May

Re: [R] timeSeq and TimeDate analog in R ?

2006-05-31 Thread Roger D. Peng
I think the 'fCalendar' package might be of help. -roger [EMAIL PROTECTED] wrote: Hi All : I am attempting tomove a large amount of code from Splus to R and I was hoping that there was an equivalent in R of the Splus functions timeSeq and timeDate ? I did an RSiteSearch but nothing came

[R] Vertical x axis labels with Trellis dot plots

2006-05-31 Thread Steve Miller
Using R 2.2.1 under Windows, I have developed a very busy Trellis dot plot with 2 columns of 7 panels and dates on the x axes. Is there a way to print the date labels vertically instead of horizontally so I can include more on the graph? Also, can I take some control over which dates are printed

Re: [R] Vertical x axis labels with Trellis dot plots

2006-05-31 Thread Chuck Cleland
See the scales section in the help page for xyplot. This gives you control over the number, location, and labeling of ticks as well as the rotation of the labels. Steve Miller wrote: Using R 2.2.1 under Windows, I have developed a very busy Trellis dot plot with 2 columns of 7 panels and

[R] R News, volume 6, issue 2 is now available

2006-05-31 Thread Paul Murrell
Hi The May 2006 issue of R News is now available on CRAN under the Documentation/Newsletter link. Paul (on behalf of the R News EditorialBoard) -- Dr Paul Murrell Department of Statistics The University of Auckland Private Bag 92019 Auckland New Zealand 64 9 3737599 x85392 [EMAIL PROTECTED]

Re: [R] timeSeq and TimeDate analog in R ?

2006-05-31 Thread Don MacQueen
If you look at the help for seq() in R (?seq), you will find a reference to seq.POSIXt() -Don At 4:40 PM -0500 5/31/06, [EMAIL PROTECTED] wrote: Hi All : I am attempting tomove a large amount of code from Splus to R and I was hoping that there was an equivalent in R of the Splus functions

Re: [R] a problem 'cor' function

2006-05-31 Thread Tao Shi
Hi Peter, Thank you very much for your quick reply! I'm aware of the accuracy issue of the numbers in R. I guess the thing I was puzzled is how come the same function when used in different way produce different results (it's supposed to use the same algorithm). The only explanation for

[R] Some non-fatal startup issues

2006-05-31 Thread Michael Kubovy
Dear r-helpers, I'm running Version 2.3.0 (2006-04-24) under OS X 10.4.6. (1) When I start R from the terminal, it has no complaints. (2) When I start R.app using R for Mac OS X Aqua GUI version 1.15 (3106) I get (I have increased the quote level of what is different from the output of the

Re: [R] Help on glmmPQL

2006-05-31 Thread Lorenz.Gygax
See help file of glmmPQL (the first thing you should be doing): glmmPQL(fixed, random, family, data, correlation, weights, control, niter = 10, verbose = TRUE, ...) ^^ Arguments ... niter maximum number of iterations. ... Thus, you can increase the value of