Re: [R] wildcards and removing variables

2005-10-10 Thread Prof Brian Ripley
On Mon, 10 Oct 2005, Marc Schwartz (via MN) wrote: On Mon, 2005-10-10 at 10:37 -0400, Afshartous, David wrote: All, Is there are a wildcard in R for varible names as in unix? For example, rm(results*) to remove all variable or function names that begin w/ results?

Re: [R] passing char to Fortran routine

2005-10-10 Thread Prof Brian Ripley
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html Programming questions are appropriate for R-devel, as it says there. The issue is OS-specific, but we don't see what your Fortran code is. If you send a reproducible example to the R-devel list, people may be able

Re: [R] greek symbols using pch

2005-10-11 Thread Prof Brian Ripley
This is now well off the topic of the subject line, but I am afraid some misinformation has been propagated (and that is the `bug'). There _are_ bugs in the code shown: the postscript fonts support 32:255, not 1:256, and pch:0:31 are not taken from the font. It seems an uninformed

Re: [R] Writing to a file with fixed precision

2005-10-11 Thread Prof Brian Ripley
On Mon, 10 Oct 2005, Marc Schwartz wrote: On Mon, 2005-10-10 at 19:50 -0400, Richard Hedger wrote: Hi, I'm trying to ouput to a filled with a fixed precision: eg. if I have data x=c(1.0,1.4,2.0), I want to be able to ouput the following to a file: 1.00 1.40

Re: [R] knn

2005-10-11 Thread Prof Brian Ripley
Please do read the posting guide and supply a reproducible example as it asks. Here it matters critically what you mean by `fix l' (to what value?). Is this knn in package class? Have you read the help page? Have you read the references (especially the first)? On Tue, 11 Oct 2005, Haleh

Re: [R] Q: Suggestions for long-term data/program storage policy?

2005-10-11 Thread Prof Brian Ripley
On Tue, 11 Oct 2005, Alexander Ploner wrote: we are a statistical/epidemiological departement that - after a few years of rapid growth - finally is getting around to formulate a general data storage and retention policy - mainly to ensure that we can reproduce results from published

Re: [R] probs in installing packages with R 2.2.0

2005-10-11 Thread Prof Brian Ripley
I expect you need to set a proxy. This is covered in the rw-FAQ that we do ask you to read before posting. See the item The internet download functions fail. On Tue, 11 Oct 2005, Giannitrapani, Marco GSUK-GSSC wrote: Dear R users. I was wondering if you could help me with this

Re: [R] Under-dispersion - a stats question?

2005-10-11 Thread Prof Brian Ripley
On Mon, 10 Oct 2005, Martin Henry H. Stevens wrote: Hello all: I frequently have glm models in which the residual variance is much lower than the residual degrees of freedom (e.g. Res.Dev=30.5, Res.DF = 82). Is it appropriate for me to use a quasipoisson error distribution and test it with

Re: [R] step.gam and number of tested smooth functions

2005-10-12 Thread Prof Brian Ripley
step.gam is a tricky function to use correctly. You will need to consult the original documentation (in Chambers Hastie ca 1992) or ask the package author for help. BTW, it uses loess not lowess. On Wed, 12 Oct 2005 [EMAIL PROTECTED] wrote: Hi, I'm working with step.gam in gam package.

Re: [R] adding 1 month to a date

2005-10-12 Thread Prof Brian Ripley
On Wed, 12 Oct 2005, bogdan romocea wrote: Simple addition and subtraction works as well: as.Date(1995/12/01,format=%Y/%m/%d) + 30 If you have datetime values you can use strptime(1995-12-01 08:00:00,format=%Y-%m-%d %H:%M:%S) + 30*24*3600 where 30*24*3600 = 30 days expressed in seconds.

Re: [R] running JPEG device on R 1.9.1 using xvfb-run on Linux

2005-10-13 Thread Prof Brian Ripley
On Wed, 12 Oct 2005, David Zhao wrote: Does anybody have experience in running jpeg device using xvfb-run on linux? I've been having sporadic problem with: /usr/X11/bin/xvfb-run /usr/bin/R --no-save Rinput.txt, with error saying: error in X11 connection. Especially when I run it from a perl

Re: [R] About Qusi-Monte carlo program

2005-10-13 Thread Prof Brian Ripley
I suspect you mean `Quasi-Monte Carlo', but that is used to do integration not simulation. Using help.search led ne to LowDiscrepancy(fOptions) Low Discrepancy Sequences QUnif(sfsmisc) Quasi Randum Numbers via Halton Sequences both of which generate the

Re: [R] Optim with two constraints

2005-10-13 Thread Prof Brian Ripley
This is actually quadratic programming, so why do you want to use optim()? There are packages specifically for QP, e.g. quadprog. A more general approach is to eliminate one variable, which gives you an inequality constrained problem in n-1 variables to which you could apply contrOptim().

Re: [R] drawing against a date

2005-10-13 Thread Prof Brian Ripley
On Thu, 13 Oct 2005, [EMAIL PROTECTED] wrote: hy all I wish to draw a graph against a date, I have a set of date like this DD/MM/, corresponding to it a set of integer , i wish to draw on x side the dates (the space between the dates have to be constant, not based on the time between

Re: [R] Help with Matrix package

2005-10-13 Thread Prof Brian Ripley
The first thing is to ensure that you are using an optimized BLAS. On Windows, use Goto's BLAS if you have it (is not currently available and redistribution is not allowed) or one of the pre-built ATLAS-based Rblas.dll on CRAN or (best of all) optimize your own build of ATLAS. The Matrix

Re: [R] Do Users of Nonlinear Mixed Effects Models Know Whether Their Software Really Works?

2005-10-13 Thread Prof Brian Ripley
On Thu, 13 Oct 2005, dave fournier wrote: Internal Virus Database is out-of-date. Talk about not being careful! -- Brian D. Ripley, [EMAIL PROTECTED] Professor of Applied Statistics, http://www.stats.ox.ac.uk/~ripley/ University of Oxford, Tel: +44 1865 272861

Re: [R] arima: warning when fixing MA parameters.

2005-10-14 Thread Prof Brian Ripley
It's a bug: the code has 1:arma[1], i.e. 1:0. Replace by seq(length=arma[1]). On Thu, 13 Oct 2005, John Maindonald wrote: I am puzzled by the warning message in the output below. It appears whether or not I fit the seasonal term (but the precise point of doing this was to fit what is

Re: [R] fast and stable way to read data in R

2005-10-14 Thread Prof Brian Ripley
On Fri, 14 Oct 2005, Wensui Liu wrote: Dear useRs, I am wondering what is the most fast and stable way to read data (pretty large) into R. Right now, the methods I could think of are: 1) use read.table to read *.csv or *txt 2) use RODBC to read excel or access But I don't know which is

Re: [R] How to find S4 generics?

2006-07-19 Thread Prof Brian Ripley
. After I find something with methods(class=...), then I have better ideas for what to look for in the documentation. Best Wishes, Spencer Graves Prof Brian Ripley wrote: On Tue, 18 Jul 2006, Spencer Graves wrote: * * methods * * You have

Re: [R] error when compiling stats library in R-2.3.1 on Solaris x86

2006-07-19 Thread Prof Brian Ripley
How did 'cc -xtarget=generic64' get there? AFAIK R does not know about it, so presumably you specified it for CC. You need to do the same thing for *all* the compilers, that is CC, CXX, F77 and FC. The INSTALL file asked you to read the R-admin manual: there you will find a very similar

Re: [R] Automating package building packages and repository uploading

2006-07-20 Thread Prof Brian Ripley
On Thu, 20 Jul 2006, Carlos J. Gil Bellosta wrote: Dear Rusers, I have developed two packages for a client of mine. After new features are added or bugs corrected, I upload them to my own web repository. I create both source and binary versions. binary Linux packages, it seems. The latter

Re: [R] Automating package building packages and repository uploading

2006-07-20 Thread Prof Brian Ripley
://www.datanalytics.com http://www.data-mining-blog.com Quoting Prof Brian Ripley [EMAIL PROTECTED]: On Thu, 20 Jul 2006, Carlos J. Gil Bellosta wrote: Dear Rusers, I have developed two packages for a client of mine. After new features are added or bugs corrected, I upload them

Re: [R] Help on making code faster-would C help ?

2006-07-20 Thread Prof Brian Ripley
Without details it is hard to say, but if you do the WLS via lm.wfit, probably not. If you are using lm() and taking the overhead 1 times, that is the first thing to avoid. As an example of this sort of thing, look at the C code for lqs{MASS} which does thousands of regressions quite

Re: [R] how to print table with more columns per row?

2006-07-20 Thread Prof Brian Ripley
This is controlled by options(width): characters not columns. On Thu, 20 Jul 2006, [EMAIL PROTECTED] wrote: When printing a table it is broken at some point (depending how long are the associated names) see example below. Is there a way to control number of columns being printed for a

Re: [R] Q. regarding optim()

2006-07-21 Thread Prof Brian Ripley
On Fri, 21 Jul 2006, N. Goodacre wrote: Dear R mailing group, The second parameter for the function optim()is a function whose parameters are to be optimized. The description of this function given in the help file is the following: fn: A function to be minimized (or maximized),

Re: [R] Parameterization puzzle

2006-07-21 Thread Prof Brian Ripley
R does not know that poly(age,2) and poly(age,1) are linearly dependent. (And indeed they only are for some functions 'poly'.) I cannot reproduce your example ('l' is missing), but perhaps glm(deaths ~ poly(age,2) + poly(age,1)*Smoke + offset(l), poisson) was your intention? On Fri, 21 Jul

Re: [R] from character to numeric over multiple columns

2006-07-21 Thread Prof Brian Ripley
Are the columns factors or character? I'll try to write code that copes with both: nm - unique(c(as.character(col1), as.character(col2), as.character(col3))) DF[] - lapply(DF, function(x) match(x, nm)) On Fri, 21 Jul 2006, Federico Calboli wrote: Hi All, I have a data frame of three

Re: [R] glm cannot find valid starting values

2006-07-21 Thread Prof Brian Ripley
On Fri, 21 Jul 2006, Dan Bebber wrote: glm(S ~ -1 + Mdif, family=quasipoisson(link=identity), start=strt, sdat) gives error: Error in glm.fit(x = X, y = Y, weights = weights, start = start, etastart = etastart, : cannot find valid starting values: please specify some strt

Re: [R] : Arial font for text in lattice plots under Linux

2006-07-21 Thread Prof Brian Ripley
This is a question about devices, not lattice per se, so applies to all forms of plotting. For jpeg, you can specify the fonts via the 'fonts' argument: see ?X11. For postscript, you can convert the fonts to Type 1 via e.g. ttf2pt1, or you can use ttf2afm to make .afm files and a postscript

Re: [R] compile R with ACML support | RHEL 4

2006-07-23 Thread Prof Brian Ripley
I doubt if /opt/acml3.5.0/gnu/lib is in your library path (it might be in your ldcache paths). So you need to set LD_LIBRARY_PATH or supply -L Look in config.log to find out what actually happened. BTW: this was more of an R-devel question than R-help. On Sat, 22 Jul 2006, Evan Cooch

Re: [R] running jobs in background

2006-07-25 Thread Prof Brian Ripley
On Tue, 25 Jul 2006, Uwe Ligges wrote: Nair, Murlidharan T wrote: Can I run jobs in the background and the check the status of it from time to time in Windows version of R? Depends on your version of Windows. I am running the automatical building of R packages in the background on a

Re: [R] deparse - width.cutoff

2006-07-25 Thread Prof Brian Ripley
On Mon, 24 Jul 2006, johan Faux wrote: I have a question about deparse function in R What is the reason that deparse use an argument like width.cutoff ? Why the maximum cutoff is 500? I was manipulating an R formula and used deparse. Since the length of user's formula was greater then

Re: [R] [Way OT] New hardware

2006-07-26 Thread Prof Brian Ripley
On Wed, 26 Jul 2006, Martin Maechler wrote: Sean == Sean Davis [EMAIL PROTECTED] on Tue, 25 Jul 2006 17:16:02 -0400 writes: Sean Can anyone share experience with opteron versus the Sean xeon (woodcrest) for R under linux? I am looking at Sean using 16-32Gb of ram in a

Re: [R] RODBC on linux

2006-07-26 Thread Prof Brian Ripley
On Wed, 26 Jul 2006, Marc Schwartz wrote: On Wed, 2006-07-26 at 17:52 -0400, Armstrong, Whit wrote: Anyone out there using Linux RODBC and unixODBC to connect to a Microsoft SQL server? If possible can someone post a sample .odbc.ini file? I saw a few discussions on the archives a

Re: [R] how to skip certain rows when reading data

2006-07-28 Thread Prof Brian Ripley
On Thu, 27 Jul 2006, [EMAIL PROTECTED] wrote: Dear all, I am reading the data using read.table. However, there are a few rows I want to skip. How can I do that in an easy way? Suppose I know the row number that I want to skip. Thanks so much! The easy way is to read the whole data frame

Re: [R] HTTP User-Agent header

2006-07-28 Thread Prof Brian Ripley
What OS is this? If Windows, see the rw-FAQ Q2.19. Otherwise, see ?download.file and choose a different download method, or look at the source code (src/modules/internet/nanohttp.c) and submit a patch. On Fri, 28 Jul 2006, James P. Howard, II wrote: I am using R in an environment where the

Re: [R] Running R on a 64 bit processor

2006-07-28 Thread Prof Brian Ripley
Be careful not to use clock speed as a measure of computer performance. Pentium 4s (and the comparable Xeons) were intended to be run very fast, but never managed it. So a 2.4GHz P4 proved to be slower than a 1GHz PIII. Unless you are running Windows 64, the chip having some 64-bit

Re: [R] Non-interpreted strings

2006-07-29 Thread Prof Brian Ripley
On Thu, 27 Jul 2006, [EMAIL PROTECTED] wrote: I am new to R, so please forgive me if there is an obvious answer to this question. I have done fairly extensive searching through R docs, google and a few R users and have not found an answer to my question. Is there a way to create a

Re: [R] Please HELP: Problem with BUILD command

2006-07-31 Thread Prof Brian Ripley
As we have said to you before, we need a reproducible example, and it is very likely that you are speculating. There is no 'BUILD' command in R: do you mean 'R CMD build'? If so, that does not parse any `R program'. On Mon, 31 Jul 2006, Zajd, John wrote: Greetings, I am unable to

Re: [R] glmmNQ

2006-07-31 Thread Prof Brian Ripley
On Mon, 31 Jul 2006, Maria Salomé Esteves Cabral wrote: Hi! Can anyone let me know where is the function glmmNQ? It's said that it is in the MASS library but I can not find it. Where is it said to be in the MASS *package*? Not in MASS the book, for sure. The function of that name used in

Re: [R] Monospaced fonts in legends

2006-07-31 Thread Prof Brian Ripley
On Mon, 31 Jul 2006, Erich Neuwirth wrote: Is there a way of using monospaced fonts in legends in a plot, but still using standard proportionally spaced fonts for all the titles? Yes, just use the family= argument as appropriate, e.g. par(family=mono) before calling legend. -- Brian D.

Re: [R] Sweave error in example code

2006-07-31 Thread Prof Brian Ripley
On Mon, 31 Jul 2006, Sundar Dorai-Raj wrote: LL wrote: Hi.. I am running R version 2.3.1 on a Windows XP machine with the latest Miktex 2.5 installed. I get no errors from R when running the Sweave example, testfile - system.file(Sweave, Sweave-test-1.Rnw, package = utils)

Re: [R] help on fitting negative binomial distribution with MLE

2006-08-01 Thread Prof Brian Ripley
This is a special case of MASS::fitdistr MASS::glm.nb and the first will be easiest for you. On Tue, 1 Aug 2006, zhijie zhang wrote: Dear friends, Anybody knows how to fit the negative binomial distribution with MLE using R or other software? I can't find the solution, any suggestions or

Re: [R] open DLL in R

2006-08-01 Thread Prof Brian Ripley
On Tue, 1 Aug 2006, qian li wrote: I have downloaded a DLL file. I want to look at the contents in the DLL file. How can I do it in R? You need a disassembler such as VC++'s DUMPBIN, but looking at compiled code you did not write is not an easy task. (Or objdump from the MinGW toolset.)

Re: [R] libRmath

2006-08-01 Thread Prof Brian Ripley
On Wed, 2 Aug 2006, TAN Chong Hui wrote: Hi I'm trying to use the standalone library libRmath. According to the documentation, I need to builld it in src/nmath/standalone But this directory does not exist in the R I installed. What's the problem here, anyone? Did you install a

Re: [R] R an AIX

2006-08-01 Thread Prof Brian Ripley
Did you look in config.log for more information about the error? It looks like you cannot even run your compiler. (This is not really an R question, and off-topic for R-help: see the posting guide.) On Tue, 1 Aug 2006, kamil Marcinkowski wrote: Hello All, I am having trouble with

Re: [R] open DLL in R

2006-08-01 Thread Prof Brian Ripley
to know what are the exported functions and constants and imported functions and constants in a easy way. Yes, see the R-admin manual, using non-R tools. thanks -gaurav. Prof Brian Ripley [EMAIL PROTECTED] Sent by: [EMAIL PROTECTED] 01-08-06 11:18 PM To qian li [EMAIL

Re: [R] Support vector in lcrabs example

2006-08-02 Thread Prof Brian Ripley
The examples in the book were run in 2001, using S-PLUS (as it says). The R package e1071 has changed since then, and hence the results it gives have changed. However, the S-PLUS version (which has not been updated) still gives the results in the book. The `problem' is your expectation that R

Re: [R] Summary method needed?

2006-08-02 Thread Prof Brian Ripley
On Wed, 2 Aug 2006, Christian Hennig wrote: Hi list, I'm updating my fpc package at the moment and will add some new functions. I learned that there should be print and summary methods for the key functions. for 'classes', I think. The purpose of the summary methods seems to be to

Re: [R] Summary method needed?

2006-08-02 Thread Prof Brian Ripley
On Wed, 2 Aug 2006, Christian Hennig wrote: Thank you Brian! I'm updating my fpc package at the moment and will add some new functions. I learned that there should be print and summary methods for the key functions. for 'classes', I think. Yes. But in some cases the print

Re: [R] From 2.2.1 to 2.3

2006-08-03 Thread Prof Brian Ripley
On Wed, 2 Aug 2006, [EMAIL PROTECTED] wrote: Hello everyone. Currently I am running R 2.2.1 (windows), and I will like to update to 2.3. I wanted to ask if it is possible to update without having to removed 2.2.1; or do I first need to delete 2.2.1? There is no such version as `2.3', as

Re: [R] question about dll crashing R

2006-08-03 Thread Prof Brian Ripley
On Thu, 3 Aug 2006, Thomas Lumley wrote: On Thu, 3 Aug 2006, Benn Fine wrote: I have ported some R code to C to make it faster. I can perform .Call(foobar,) once and it works fine. Absolutely correct answer. If I put a loop inside foobar and run the main code routine more

Re: [R] User input from keyboard

2006-08-04 Thread Prof Brian Ripley
On Fri, 4 Aug 2006, Petr Pikal wrote: Hi cat(\n,Enter x,\n) # prompt y-scan(n=1) prompts for user imput and scans 1 line from console. (that scans one number: use readLines(n=1) to get a string). But readline() is probably easier. HTH Petr On 4 Aug 2006 at 17:06, chiya sharma

Re: [R] meta characters in file path

2006-08-04 Thread Prof Brian Ripley
On Fri, 4 Aug 2006, Li,Qinghong,ST.LOUIS,Molecular Biology wrote: Thanks. I tried them, it works for most of those characters except * and ?. Those are not valid characters in Windows file paths (/ * : ? \ | are invalid in file or dir names). Does regular expression work in file names

Re: [R] Error when loading odesolve

2006-08-04 Thread Prof Brian Ripley
Looking at odesolve, src/Makevars is PKG_LIBS=$(BLAS_LIBS) Now, the documentation says that if you have $(BLAS_LIBS) you must also have $(FLIBS), so please change this to PKG_LIBS=$(BLAS_LIBS) $(FLIBS) and take this up with the package maintainer (which is what the posting guide asked you to

Re: [R] expression() - Superscript in y-axis, keeping line break in string

2006-08-04 Thread Prof Brian Ripley
On Fri, 4 Aug 2006, Marc Schwartz (via MN) wrote: On Fri, 2006-08-04 at 09:47 -0600, Andrew Kniss wrote: I've tried several different ways to accomplish this, but as yet to no avail. My y-axis for a plot has a rather long label, and thus I have been using /n to break it into two lines.

Re: [R] Postscript fonts

2006-08-05 Thread Prof Brian Ripley
On Sat, 5 Aug 2006, Erich Neuwirth wrote: How can I find out what fonts are available for par(family= for the postscript device? This is dynamic: for the current value names(postscriptFonts()) [1] serifsans mono [4] symbol AvantGarde

Re: [R] AIC for lognormal model

2006-08-05 Thread Prof Brian Ripley
On Sat, 5 Aug 2006, Andreas Beyerlein wrote: Dear all, I want to compare some different models for a dataset by QQ plots and AIC. I get the following AICs: - linear model: 19759.66 - GAMLSS model: 18702.7 - linear model with lognormal response: -7862.182 The QQ plots show that the

Re: [R] Kmeans - how to display results

2006-08-05 Thread Prof Brian Ripley
There are worked examples in MASS (the book, see the FAQ). On Sat, 5 Aug 2006, Ffenics wrote: I'm very new as regards to R. I have managed to work out how to use dist and kmeans but am now wondering how best to display the results from kmeans in a graphical form. If anyone has any general

Re: [R] paired t-test. Need to rearrange data?

2006-08-06 Thread Prof Brian Ripley
On Sun, 6 Aug 2006, Henrik Parn wrote: Dear all, I have received some data arranged like this: # original data id - rep(letters[1:6], each=2) time - as.factor(rep(1:2, 6)) y - as.vector(replicate(6, c(rnorm(n=1, mean=1), rnorm(n=1, mean=2 test.data - data.frame(id, time, y)

Re: [R] Source installation error: gfortran and gcc disagree on int and double ...

2006-08-07 Thread Prof Brian Ripley
First, you do not need -fPIC: it is the same as -fpic which R selects on your platform. Second, please look at config.log to find the exact problem: it well be that your compilers are not properly installed (as was the case in the reference you quote below). Finally, your compilers are pretty

Re: [R] Backquote in R syntax

2006-08-07 Thread Prof Brian Ripley
?Quote tells you all about quotes in R (and so does ?backquote in R-devel, but many people call it backtick and that is all 2.3.1 has). On Mon, 7 Aug 2006, [EMAIL PROTECTED] wrote: Hi Folks, Can someone satisfy my curiosity (well, at least about the following query!) Reading the draft

Re: [R] CPU Usage with R 2.1.0 in Windows (and with R 2.3.1)

2006-08-07 Thread Prof Brian Ripley
On Mon, 7 Aug 2006, Singu wrote: You could compile R using specific libraries if threading/dual core is the issue Look eg here: http://tolstoy.newcastle.edu.au/R/devel/05/12/3355.html Hmm, quoting me. None of those libraries for Windows is multithreaded, nor does it say anthing about

Re: [R] Capturing stderr from system()

2006-08-07 Thread Prof Brian Ripley
On Mon, 7 Aug 2006, McGehee, Robert wrote: Hello, I have a system command that occasionally fails and writes output to standard error, which R will print to the screen when ignore.stderr = FALSE. No, R does not print to the screen, your shell does. You can capture shell and command errors

Re: [R] Source installation error: gfortran and gcc disagree on int and double ...

2006-08-08 Thread Prof Brian Ripley
that is not being found or is broken. One possibility is a missing symlink that is in a -devel RPM. This is not an R issue. On Mon, 7 Aug 2006, Mike wrote: Prof Brian Ripley wrote: First, you do not need -fPIC: it is the same as -fpic which R selects on your platform. Right, I

Re: [R] Netiquette, was Re: ... gfortran and gcc...

2006-08-08 Thread Prof Brian Ripley
On Tue, 8 Aug 2006, Peter Dalgaard wrote: Prof Brian Ripley [EMAIL PROTECTED] writes: First, you replied to the list and not to me, which was discourteous. You mean that he replied to the list *only*, I hope. Yes, and it was written as if to me, and was a reply to an email from me. I

Re: [R] prefixing list names in print

2006-08-08 Thread Prof Brian Ripley
On Tue, 8 Aug 2006, Laurent Deniau wrote: With print(list(A=a,B=b)) it displays $A [1] a $B [1] b I would like to add a common prefix to all the list tags after the $. `prefix' ... `after'? You seem to want to prefix component names: why? What do you want for component $a$b$c?

Re: [R] problem with tseries (was unable to restore saved data)

2006-08-09 Thread Prof Brian Ripley
Hmm. 'time' is in package stats, and you seem to have saved an object that is pulling in the namespace 'tseries'. Loading the workspace is done before loading the standard set of packages. This appears to be a bug in tseries, for if its namespace needs 'time', it should import it from 'stats',

Re: [R] How to draw the decision boundaries for LDA and Rpart object

2006-08-09 Thread Prof Brian Ripley
On Wed, 9 Aug 2006, Am Stat wrote: Hello useR, Could you please tell me how to draw the decision boundaries in a scatterplot of the original data for a LDA or Rpart object. There are examples in MASS (the book). For example: library(rpart) fit.rpart - rpart(as.factor(group.id)~.,

Re: [R] missing documentation entries

2006-08-09 Thread Prof Brian Ripley
This is discussed in `Writing R Extensions', which both points you to the 'internal' keyword, and (in the current version) mentions using name spaces to hide such functions. This really was a question for R-devel: please do study the posting guide. `R-devel is intended for questions and

Re: [R] Error message when using optim

2006-08-09 Thread Prof Brian Ripley
(Subject changed to something less perjorative. This is not `optim error'.) On Wed, 9 Aug 2006, Frank Black wrote: Dear all, There have been one or two questions posted to the list regarding the optim error non-finite finite-difference value [4]. The error apparently means that the

Re: [R] How to draw the decision boundaries for LDA and Rpart object

2006-08-10 Thread Prof Brian Ripley
in MASS (the book). Thanks very much! Leon - Original Message - From: Prof Brian Ripley [EMAIL PROTECTED] To: Am Stat [EMAIL PROTECTED] Cc: r-help@stat.math.ethz.ch Sent: Wednesday, August 09, 2006 3:04 AM Subject: Re: [R] How to draw the decision boundaries for LDA

Re: [R] decimal accuracy in pnorm( )

2006-08-10 Thread Prof Brian Ripley
On Thu, 10 Aug 2006, Chronopoulos Dimitris wrote: Dear R users Is there any way to increase the decimal accuracy for the normal probability distribution? When one needs an accurate p-value for instance this is provided by pnorm(10,lower.tail=F) [1] 7.619853e-24 However, what

Re: [R] extractAIC using surf.ls

2006-08-11 Thread Prof Brian Ripley
Roger, thank you for looking into this. However, the posting guide asked the poster to contact the maintainer. Had (s)he done so, I would have pointed out that spatial 7.28-2 (the current version for R-devel) has this corrected (in a slightly simpler way). On Thu, 10 Aug 2006, Roger Bivand

Re: [R] - factanal scores correlated?

2006-08-11 Thread Prof Brian Ripley
This is a Heywood case, and you don't have a valid fit: myfac Call: factanal(x = m1, factors = 3, scores = regression) Uniquenesses: v1v2v3v4v5v6 0.005 0.101 0.005 0.224 0.084 0.005 notice no less than 3 very small uniquenesses. On Fri, 11 Aug 2006, Christian Montel

Re: [R] bug in interaction order when using drop?

2006-08-11 Thread Prof Brian Ripley
On Thu, 10 Aug 2006, Petr Pikal wrote: Ooops, my first suggestion reorders factor itself but if (drop) factor(ans) else ans instead of whole drop construction shall preserve levels order without changing order of factor Even easier would be to return ans[,drop=drop]. It seems to me

Re: [R] anova.mlm for single model (one-way repeated measured anova)

2006-08-12 Thread Prof Brian Ripley
On Sat, 12 Aug 2006, takahashi kohske wrote: Dear list members: I'd like to one-way repeated measured anova by using mlm. I'm using R-2.3.1 and my code is: dat-matrix( c(9,7,8,8,12,11,8,13, 6,5,6,3,6,7,10,9, 10,13,8,13,12,14,14,16, 9,11,13,14,16,12,15,14),

Re: [R] Changing R network connections

2006-08-12 Thread Prof Brian Ripley
Since you mention 'preferences', it seems you might be using Windows: the information you require is in the rw-FAQ. Since the posting guide asked you to state your platform and read the relevant FAQ, that is the place to start: see the footer below. On Fri, 11 Aug 2006, Alexandre Fabrice

Re: [R] Puzzling warning using 2.3.1...

2006-08-13 Thread Prof Brian Ripley
On Sun, 13 Aug 2006, Brian Lunergan wrote: Greetings folks: Stepped away from a win/lin dual boot system and spent the better part of the past week resetting a windows only arrangement. Before the switch I was running 2.3.1 without any strange behaviours. When I reached the point of

Re: [R] left-justified fixed-width format

2006-08-14 Thread Prof Brian Ripley
On Mon, 14 Aug 2006, roger bos wrote: I need to save data in fixed-width format without headers and reading the help archive leads me to believe that sprintf is pretty much the only way to do this. My question is, is there anyway to change the output so the text in each column is left

Re: [R] Attempt to access unmapped memory

2006-08-14 Thread Prof Brian Ripley
On Mon, 14 Aug 2006, Alexandre Aguiar wrote: Hi, I am usiing R 2.3.1 under Linux kernel 2.6.11 with libreadline/libhistory 5.1. That is known to be buggy. Either patch it or downdate to 5.0. Bothe R and readline were compiled without quircks with gcc 3.3.3 and g77. Every time I try

Re: [R] Fast way to load multiple files

2006-08-14 Thread Prof Brian Ripley
On Mon, 14 Aug 2006, Berton Gunter wrote: A reproducible example here would help (please see posting guide). A guess: is your filelist a list of (quoted) character strings? Correct pathnames to the files with correct separators for your OS? I think the issue is (from the help page) Usage:

Re: [R] help: cannot allocate vector of length 828310236

2006-08-15 Thread Prof Brian Ripley
Does it make any statistical sense to do polr or probit regression (not the same thing) on `really huge data'? There are few regression-like problems in which model inadequacy does not swamp estimation uncertainty for as few as a 1000 cases. If you want to do that sort of thing, by all means

Re: [R] Protection stack overflow

2006-08-15 Thread Prof Brian Ripley
R has a command-line option to set the ppstack size, --max-ppsize=NSet max size of protect stack to N Looks like you need to supply this (and it can be done with embedded R) if the problem persists with current R. You could also try arima0 or even ar to do the fitting. On Tue, 15 Aug

Re: [R] nls

2006-08-15 Thread Prof Brian Ripley
You problem is x^c for x = 0. If you intended only c 1, try a starting value meeting that condition (but it seems that the optimal c is about 0.27 is you increase x slightly). Why have you used ~~ ? (Maybe because despite being asked not to, you sent HTML mail?) On Tue, 15 Aug 2006,

Re: [R] Aliases for arguments in a function

2006-08-15 Thread Prof Brian Ripley
foo - function(arg1, this=that, that) ... works. On Tue, 15 Aug 2006, [EMAIL PROTECTED] wrote: Hi all. I have a function that I would like to use either the argument name as originally defined or another name. Within the function (and other functions) use the argument name as

Re: [R] nls

2006-08-15 Thread Prof Brian Ripley
On Tue, 15 Aug 2006, Peter Dalgaard wrote: Prof Brian Ripley [EMAIL PROTECTED] writes: You problem is x^c for x = 0. If you intended only c 1, try a starting value meeting that condition (but it seems that the optimal c is about 0.27 is you increase x slightly). Surely you mean c

Re: [R] question re: summarry.lm and NA values

2006-08-15 Thread Prof Brian Ripley
On Tue, 15 Aug 2006, Petr Pikal wrote: Hi On 15 Aug 2006 at 7:01, r user wrote: Date sent:Tue, 15 Aug 2006 07:01:13 -0700 (PDT) From: r user [EMAIL PROTECTED] To: rhelp r-help@stat.math.ethz.ch Subject: [R] question re:

Re: [R] read.csv issue

2006-08-16 Thread Prof Brian Ripley
Set allowEscapes = FALSE when reading. See the help page for more details. There is perhaps an argument for changing the default for allowEscapes under read.csv, especially as people have now changed that for comment.char (in R-devel). On Wed, 16 Aug 2006, Doran, Harold wrote: I'm trying to

Re: [R] read.csv issue

2006-08-16 Thread Prof Brian Ripley
On Wed, 16 Aug 2006, Prof Brian Ripley wrote: Set allowEscapes = FALSE when reading. See the help page for more details. There is perhaps an argument for changing the default for allowEscapes under read.csv, especially as people have now changed that for comment.char (in R-devel). Oops

Re: [R] Plots Without Displaying

2006-08-16 Thread Prof Brian Ripley
Yes, see ?jpeg ?bitmap and as you didn't tell us your OS we don't know if these are available to you. jpeg(file=test.jpg) boxplot(sample(100)) dev.off() may well work. 'An Introduction to R' explains about graphics devices, including these. On Wed, 16 Aug 2006, Lothar Botelho-Machado

Re: [R] Plots Without Displaying

2006-08-17 Thread Prof Brian Ripley
:[EMAIL PROTECTED] On Behalf Of Lothar Botelho-Machado Sent: Wednesday, August 16, 2006 4:49 PM To: r-help@stat.math.ethz.ch Subject: Re: [R] Plots Without Displaying Prof Brian Ripley wrote: Yes, see ?jpeg ?bitmap and as you didn't tell us your OS we don't know

Re: [R] Maximum length of R GUI input line?

2006-08-18 Thread Prof Brian Ripley
On Fri, 18 Aug 2006, Philippe Grosjean wrote: You should put your SQL query in a variable, and use this variable in your call. Something like: MyQuery - bla bla bla MyQuery - paste(MyQuery, more bla bla) # doMyRequest(MyQuery) Indeed. For the record, R has an limit of 1024 chars

Re: [R] as.data.frame(cbind()) transforming numeric to factor?

2006-08-18 Thread Prof Brian Ripley
On Fri, 18 Aug 2006, Tom Boonen wrote: Dear List, why does as.data.frame(cbind()) transform numeric variables to factors, once one of the other variablesused is a character vector? # x.1 - rnorm(10) x.2 - c(rep(Test,10)) Foo - as.data.frame(cbind(x.1)) is.factor(Foo$x.1) Foo -

Re: [R] R-update - what about packages and ESS?

2006-08-18 Thread Prof Brian Ripley
Without knowing your OS this is hard to answer (and is the wrong list for the ESS question). For Windows users, the packages part is covered in the rw-FAQ. For Unix-alikes, it all depends how the update was done, but normal package update mechanisms (such as RPM) will not wipe out previously

Re: [R] string-to-number

2006-08-19 Thread Prof Brian Ripley
On Sat, 19 Aug 2006, Marc Schwartz wrote: On Sat, 2006-08-19 at 07:58 -0400, Charles Annis, P.E. wrote: Greetings, Amigos: I have been trying without success to convert a character string, repeated.measures.columns [1] 3,6,10 into c(3,6,10) for subsequent use.

Re: [R] problem with Rcmd check and fortran95, makefile

2006-08-19 Thread Prof Brian Ripley
First, this is off-topic for R-help (see the posting guide, including the non-use of HTML, please). Second, the appropriate list needs to see exactly what the output was when you did this. Remember R CMD SHLIB is a *make* facility, and quite possibly the objects were not remade after changes.

Re: [R] need to find (and distinguish types of) carriage returns in a file that is scanned using scan

2006-08-19 Thread Prof Brian Ripley
On Sat, 19 Aug 2006, Quicke, Donald L J wrote: Hope this is not too trivial I am reading a large file using scan. Why scan? In one part of this file there is a chunk of text within which i need to know the positions of line breaks. But scan seems only only what? An example of the file

Re: [R] question about cbind()

2006-08-20 Thread Prof Brian Ripley
On Sat, 19 Aug 2006, [EMAIL PROTECTED] wrote: Dear all, I have a question about how to get a matrix by combining a large number of columns from a data file. Suppose I read a file which have 1000 columns by: test = read.table(dat.txt, header=F) I know I could use cbind(). It's easy to

Re: [R] issues with Sweave and inclusion of graphics in a document

2006-08-20 Thread Prof Brian Ripley
savePlot is just an internal version of dev.copy, part of the support for the menus on the windows() graphics device. It is described in `An Introduction to R' (the most basic R manual). On Sat, 19 Aug 2006, Thomas Harte wrote: the problem is a little hard to explain; the .Rnw files (below)

<    1   2   3   4   5   6   7   8   9   10   >