Re: [R] generalized negative binomial

2005-03-13 Thread Prof Brian Ripley
On Sat, 12 Mar 2005, Erin M Simpson wrote: I am looking for code that allows for a more flexible negative binomial model (similar to Stata's gnbreg). Your subject line is not clear to me: Stata appears to fit a negative binomial model, the point being that it is not a glm as fitted by glm.nb.

[R] how to draw the data set processed by hclust?

2005-03-13 Thread XP Sun
hi, all i hava a dataset formated as follow: x1 y1 z1 x2 y2 z2 x3 y3 z3 ... how to cluster it with hclust? and the draw the data with color of cluster? thank you in advance! best! xpsun

[R] Big databases

2005-03-13 Thread Kenneth Cabrera
Hi R users: I am using R 2.01. over a LINUX(Scientific Linux CERN) platform and I got the following problem: It takes too much time (more than 6 hours, because I stop the process) to read a 270MB database. I am using read.table() function. Is there any workaround to read faster a big data base?

Re: [R] XML to data frame or list

2005-03-13 Thread bogdan romocea
I managed to parse more complex XML files as well. The trick was to manually determine the position of the child nodes of interest, after which they can be parsed in a loop. For example: require(XML) doc - xmlTreeParse(file.xml,getDTD=T,addAttributeNamespaces=T) r - xmlRoot(doc) #find the nodes

Re: [R] Multiple histograms in one chart

2005-03-13 Thread Uwe Ligges
Werner Wernersen wrote: Hi! I've searched quite a while and this is probably a very easy question: How do I prevent a new hist() histogram from clearing the chart window before drawing? I would like to have several histograms simply drawn on each other in one chart... See either ?par

Objet: [R] New user of R on Mac OS X - Please help

2005-03-13 Thread Denis Chabot
De: Depiereux Constant [EMAIL PROTECTED] Date: 12 mars 2005 09:44:08 GMT-05:00 À: r-help@stat.math.ethz.ch Objet: [R] New user of R on Mac OS X - Please help Brand new Mac OS X user, I am transfering my R stuffs from my windows machine. When porting some of my functions, I got messages such as

RE: [R] Big databases

2005-03-13 Thread Liaw, Andy
1. Read the R Data Import/Export manual. 2. Check the R-help archive, as this question has been asked and answered many times. Andy From: Kenneth Cabrera Hi R users: I am using R 2.01. over a LINUX(Scientific Linux CERN) platform and I got the following problem: It takes too much time

Re: [R] generalized negative binomial

2005-03-13 Thread Spencer Graves
I whole-heartedly endorse Prof. Ripley's suggestion to write down the log(likelihood) and use optim; I've done that many times with seemingly good results. For confidence intervals, the best procedure is to use 2*log(likelihood ratio) being approximately chi-square. If you are

Re: [R] RODBC, IODBC and Mac OS X

2005-03-13 Thread Charles and Kimberly Maner
Hi. Another, slightly more updated, link to the same info is: http://bioconductor.org/, under the Developmental Packages under Software, contains Rdbi version 1.1.1 and RdbiPgSQL version 1.1.3. I am running them under MacOSX and they work beautifully. I am using it to pull data on PostgreSQL

RE: [R] any book and tutorial about how to manipulate data with R /S+

2005-03-13 Thread Liaw, Andy
From: Wensui Liu I am sorry that I did not state my question clearly. What I mean by data manipulation includes sort, merge, aggregate, transpose, R has functions for doing those: sort(), merge(), aggregate(), and t(), respectively. data export and import, format, date time handle,

Re: [R] Simplex(boot) returning invalid answer

2005-03-13 Thread Duncan Murdoch
On Fri, 11 Mar 2005 09:09:15 -0800, Andrew Stoneman [EMAIL PROTECTED] wrote : In trying to use simplex() from the boot package, I have run into a situation that doesn't seem like it should be possible. It is claiming that it has solved the LP, but returns a vector of all zeros, which does not

[R] cross-validation

2005-03-13 Thread Trevor Wiens
I've been looking at the base and Design libraries and it is unclear to me the best way to approach doing cross-validation. I'm interested in using temporal (I have five years of data), spatial (I've divided my data set up into 5 blocks that make sense and have a block variable attached to my

[R] iterations of nnet training

2005-03-13 Thread Wensui Liu
Does anyone know how to suppress the iterations of nnet training with nnet library? Thank you so much! -- WenSui Liu, MS MA Senior Decision Support Analyst Division of Health Policy and Clinical Effectiveness Cincinnati Children Hospital Medical Center

RE: [R] iterations of nnet training

2005-03-13 Thread Liaw, Andy
From: Wensui Liu Does anyone know how to suppress the iterations of nnet training with nnet library? 1. That's _package_, not _library_. 2. See the `trace' option in ?nnet. Andy Thank you so much! -- WenSui Liu, MS MA Senior Decision Support Analyst Division of Health Policy and

[R] initialising trellis device {lattice}/ postscript

2005-03-13 Thread McClatchie, Sam (PIRSA-SARDI)
Background: OS: Linux Mandrake 10.1 release: R 2.0.0 editor: GNU Emacs 21.3.2 front-end: ESS 5.2.3 - Colleagues I am using: trellis.device(postscript, file=../figures/name.ps) .code to generate the trellis display here graphics.off() to create a

fixed/ [R] initialising trellis device {lattice}/ postscript

2005-03-13 Thread McClatchie, Sam (PIRSA-SARDI)
Thanks Andy out- xyplot() print(out) did the trick nicely. Sam Sam McClatchie, Biological oceanography South Australian Aquatic Sciences Centre PO Box 120, Henley Beach 5022 Adelaide, South Australia email [EMAIL PROTECTED] Telephone: (61-8) 8207 5448 FAX: (61-8) 8200 2481 Research

Re: [R] any book and tutorial about how to manipulate data with R/S+

2005-03-13 Thread Gabor Grothendieck
Wensui Liu liuwensui at gmail.com writes: : transpose, data export and import, format, date time handle, and so Regarding just the date and time classes part of your question, an article in R News 4/1 discusses that. __ R-help@stat.math.ethz.ch

Re: [R] Multiple histograms in one chart

2005-03-13 Thread Adaikalavan Ramasamy
If you use par(new=TRUE), you are overlaying one graph on top of the other. In which case, make sure your xlim and ylim are correctly set. Another way is to split the plotting window. For example par(mfrow=c(2,3)) for(i in 1:6) hist( rnorm(100), main=paste(Histogram, i)) See help(par) for

Re: [R] cross-validation

2005-03-13 Thread Trevor Wiens
On Sun, 13 Mar 2005 15:28:54 -0700 Trevor Wiens [EMAIL PROTECTED] wrote: I've been looking at the stats and Design libraries and it is unclear to me the best way to approach doing cross-validation. I'm interested in using temporal (I have five years of data), spatial (I've divided my data

Re: [R] initialising trellis device {lattice}/ postscript

2005-03-13 Thread Uwe Ligges
McClatchie, Sam (PIRSA-SARDI) wrote: Background: OS: Linux Mandrake 10.1 release: R 2.0.0 editor: GNU Emacs 21.3.2 front-end: ESS 5.2.3 - Colleagues I am using: trellis.device(postscript, file=../figures/name.ps) .code to generate the trellis display here

[R] polymars

2005-03-13 Thread [EMAIL PROTECTED]
Can you help me interpret the output I get with “polymars”? I'd like to find a regression model to predict 2 dependent variables (then called Y1 and Y2) with 2 independent variables (then called X1 and X2). Here is the output: polymars(responses = data[, 13:14], predictors = data[, 2:3]) Model