Re: [R] help comparing two median with R

2007-04-18 Thread Prof Brian D Ripley
On Wed, 18 Apr 2007, Greg Snow wrote: For testing, the permutation test may be prefered to the bootstrap (though the bootstrap could be used for a confidence interval). I remember in grad school doing a project on comparing the efficiency of a permutation test on medians compared to the

Re: [R] Error message when building a package

2007-01-29 Thread Prof Brian D Ripley
Please update your Xcode tools. According to http://cran.r-project.org/bin/macosx/RMacOSX-FAQ.html 2.2 or later is needed, and 2.4.1 is current. On Tue, 30 Jan 2007, C. Lillian Yau wrote: I'm trying to build a package. The machine is PowerPC G4 with Mac OS 10.4.8, and I'm using R2.4.1. I

Re: [R] building R on freebsd 6.2 (amd64)

2007-01-18 Thread Prof Brian D Ripley
On Thu, 18 Jan 2007, Hiroyuki Kawakatsu wrote: Hi, I updated my os to freebsd 6.2 and built R-patched with no changes in configure. I pass make check and copy the R executable to /usr/local/bin/R. Then when I start R, I get /libexec/ld-elf.so.1: Shared object libRblas.so not found,

Re: [R] Seeking advice on lattice package in R2.4.0 concerning stripplot

2006-12-15 Thread Prof Brian D Ripley
The problem looks to be that you are specifying the types of symbols in two ways: pch=c(0,3,6,5), points = Rows(trellis.par.get(superpose.symbol),1:4), The simplest thing to do is to run your code directly on the pdf() device, but you could experiment with other ways of setting the key. (We

Re: [R] *** caught segfault *** error

2006-12-01 Thread Prof Brian D Ripley
On Fri, 1 Dec 2006, Juanjo Abellan wrote: [...] Sorry, I didn't understand the issue about the locale and the C locale that Prof Ripley mentioned. Where could I learn about that? In the R-admin manual (but Roger will understand). -- Brian D. Ripley, [EMAIL PROTECTED]

Re: [R] Install bioconductor

2006-11-11 Thread Prof Brian D Ripley
You are missing both the library and headers for what would be called either zlib or libz. (On an RPM-based system it would be zlib-devel that is missing.) However, this is the wrong list for questions anout installing BioC packages: please see the R posting guide. (And probably the packages

Re: [R] Adding predicted values as a new variable in a data frame

2006-09-14 Thread Prof Brian D Ripley
?na.exclude should help you: my guess is that you asked (by using the default) na.action = na.omit) for rows with missing values to be excluded from the residuals. But since you have not mentioned missing values, we have to guess what 'for some reason' was: please note the footer of this

Re: [R] Problem setting options(error=recover) in my Rprofile.site

2006-09-14 Thread Prof Brian D Ripley
On Thu, 14 Sep 2006, Marcus, Jeffrey wrote: I'd like to be able to set options(error=recover) in my Rprofile.site file. When I do this I get the message Error in options(error = recover) : object recover not found I assume this is because the utils package (where recover and dump.frames are

Re: [R] Parameterization puzzle

2006-07-21 Thread Prof Brian D Ripley
On Fri, 21 Jul 2006, Berwin A Turlach wrote: G'day all, BDR == Prof Brian Ripley [EMAIL PROTECTED] writes: BDR R does not know that poly(age,2) and poly(age,1) are linearly BDR dependent. Indeed, I also thought that this is the reason of the problem. BDR (And indeed they only

Re: [R] failed installing rgl

2006-07-21 Thread Prof Brian D Ripley
On FC5: roc% rpm -q --file /usr/include/GL/glu.h mesa-libGLU-devel-6.4.2-6.FC5.3 Do check what the R-admin manual says about -devel RPMs. On Thu, 20 Jul 2006, James Foadi wrote: Dear all, I have tried installing rgl with the usual command: R CMD INSTALL rgl_0.67-2.tar.gz Differently from

Re: [R] How to access results of survival analysis

2006-05-08 Thread Prof Brian D Ripley
The solution is to use a separate function to do this (and call it from the print method). On Mon, 8 May 2006, Heinz Tuechler wrote: At 12:55 07.05.2006 +0100, Prof Brian Ripley wrote: On Sun, 7 May 2006, Heinz Tuechler wrote: Hello Xiaochun Li! Thank you for submitting the function. At

Re: [R] error in color-labeling in plot.lm

2006-05-08 Thread Prof Brian D Ripley
On Mon, 8 May 2006, Ulrich Halekoh wrote: I want to label groups of points in the plot of residuals against predicted values generated by plot.lm(). Labelling the points by color x-1:15 y-rnorm(15) f-factor(rep(c(1:3),each=5)) g-glm(y~x) plot(g,which=1,col=c(1,2,3)[f]) I receive the

Re: [R] Link to useR! 2006 from ww.r-project.org not working

2006-05-03 Thread Prof Brian D Ripley
Yes, this has been reported to the webmasters (not really much point in telling a list). cran.r-project.org and www.r-project.org now point at machines at wu-wien.ac.at (the DNS was changed in the last 24 hours). I have also found that rsync is not being allowed through (e.g. for

Re: [R] R CMD build produces tar error under FreeBSD 5.4

2005-09-25 Thread Prof Brian D Ripley
On Sun, 25 Sep 2005, Andrew Robinson wrote: On Sun, Sep 25, 2005 at 12:07:02PM +0100, Prof Brian Ripley wrote: On Sun, 25 Sep 2005, Andrew Robinson wrote: I have constructed a package using package.skeleton(), when I try $ R CMD build foo * checking for file 'foo/DESCRIPTION' ...

Re: [R] Error in contrasts in step wise regression

2005-06-27 Thread Prof Brian D Ripley
On Mon, 27 Jun 2005, Young Cho wrote: Thanks for the reply. I created a new dataframe and ran step on it. But, still it does not work. detach(dat) attach(ds) dat - ds[,sapply(ds,nlevels)=2] dat$Y - Response detach(ds) attach(dat) fmla - as.formula(paste( ~

Re: [R] Spare some CPU cycles for testing lme?

2004-09-14 Thread Prof Brian D Ripley
As others have said, this needs tools not CPU cycles: gctorture or valgrind. Valgrind found (after a few seconds and on the first pass) ==23057== Invalid read of size 4 ==23057==at 0x3CF4E645: ssc_symbolic_permute (Mutils.c:373) ==23057==by 0x3CF5BF75: ssclme_create (ssclme.c:168)

Re: [R] linking .lib and/or .dll files

2004-07-01 Thread Prof Brian D Ripley
On Thu, 1 Jul 2004, Samuel Kemp wrote: Linking C++ dynamic libraries has become the bane of my R life. I have a piece of C++ code that I would like to make into a .dll (to call into R). However, this piece of C++ code needs to be linked with a .dll or .lib - which I have been created in

Re: [R] how to sort the contents of a list-object?

2004-05-12 Thread Prof Brian D Ripley
On Wed, 12 May 2004, Martin Kerick wrote: I have a list-object RGList named RG containing two vectors (RG$R and RG$G) and a data.frame(RG$genes) with 8 variables. Each of the variables of the data.frame have the same length as the two vectors. I think the data in the RGList is structured in

Re: [R] Size of R user base

2004-04-20 Thread Prof Brian D Ripley
On 20 Apr 2004, Peter Dalgaard wrote: [EMAIL PROTECTED] writes: On 20 Apr 2004 at 10:47, Philippe Grosjean wrote: Of course, this will only work with computers connected to the internet,... but at least, it could be one way to evaluate the number of R users. Would that be an

Re: [R] loading into a list

2004-04-20 Thread Prof Brian D Ripley
Try .saveRDS/.readRDS which enable you to save and load unnamed objects (in your case list(A, B, C)). On Tue, 20 Apr 2004, Tamas Papp wrote: I have the following problem. Use the same algorithm (with different parameters) to generate simulation results. I store these in variables, eg A, B,

Re: [R] New unique name

2004-04-19 Thread Prof Brian D Ripley
On Mon, 19 Apr 2004, Barry Rowlingson wrote: (Ted Harding) wrote: A bit sledgehammer for nut, but as least it meets your needs! You could even use uuidgen to create a universally unique ID, and then use make.names to R-ify it: make.names(system(uuidgen,intern=T)) [1]

Re: [R][S] library question

2004-03-24 Thread Prof Brian D Ripley
In some cases it is possible to port an R *package* to S-PLUS: people have done both that and the reverse. If you have some specific examples in mind we may be able to help you further. Since you are on Windows, it may not be easy unless you have the requisite tools (e.g. Visuall C++ and

Re: [R] rstandard does not produce standardized residuals

2004-02-24 Thread Prof Brian D Ripley
This has laready been changed in the R-devel version (1.9.0 to be). Pre-test versions of that for Windows are available on CRAN. From the NEWS file o rstandard() was wrongly scaled for cases where summary(model)$dispersion != 1. and that ought to mention glm somewhere On

Re: [R] Latin 2 encoding + fonts

2004-02-04 Thread Prof Brian D Ripley
On Wed, 4 Feb 2004, Miha STAUT wrote: Hi, In the FAQ I read about options to specify different fonts than the default ones for the console (in the file Rprofile) and for the graphical output (Rdevga). I would however like to replace Latin 1 with Latin 2 enconding for both (console and

Re: [R] Winedit and R

2003-12-16 Thread Prof Brian D Ripley
You are using mismatched versions of R and the packages, as the warnings say. Please update your version of R. On Tue, 16 Dec 2003, Patrick Giraudoux wrote: Hi all, I am trying to install add-on in R (rw1070) to work with WinEdit. Libraries Swinregistry and Rwinedt have been installed via

Re: [R] Calling C function in R

2003-12-16 Thread Prof Brian D Ripley
On Tue, 16 Dec 2003, Manoj - Hachibushu Capital wrote: Hi All, I am trying to write a c function to optimize loop processing. Having read the R extension and trying out a few samples I was pretty comfortable with the basics. However, I am wondering if there is anyway to call

Re: [R] WinMenus - is there a way of knowing if a WinMenu or WinMenuItem already exists?

2003-12-08 Thread Prof Brian D Ripley
On Mon, 8 Dec 2003, Finn Sando wrote: I am developing a menusystem using the functions WinMenuAdd and WinMenuAddItem etc. I want to be able to shift between different interfaces (ie. different sets of menu-trees). Therefore I would like to be able to ask whether a specific menu already

Re: [R]: GLIM PROBLEMS

2003-12-02 Thread Prof Brian D Ripley
In R you use cbind(successes, failures), not cbind(successes, total) as you appear to have done. And GLIM is a program (with I for interactive): these are GLMs, not GLIMs in th ecommonly accepted terminology (but not that of SAS). On Tue, 2 Dec 2003, allan clark wrote: Hi all I have

Re: [R] Rd Files?

2003-12-01 Thread Prof Brian D Ripley
On Mon, 1 Dec 2003, Wolski wrote: Hi! Seems to me that package.skeleton(...) does not like assignment functions. System: Windows 2000 And that's your problem! It's a problem with your OS, not package.skeleton. And it is known.

Re: [R] howto improve sharpeness of fonts in a jpg-image produced by R ?

2003-11-18 Thread Prof Brian D Ripley
On Tue, 18 Nov 2003, Niels Steen Krogh wrote: I used an expression like this: bitmap(file=barplotx.jpg,type=jpeg) So that was nothing to do with R, as reading the help page would have told you. Your subject line is false: the `jpg-image' was produced by gs and you should look into fixing

Re: [R] confint: which method attached?

2003-11-17 Thread Prof Brian D Ripley
On Mon, 17 Nov 2003, Ulrich Halekoh wrote: the function confint uses the profiling method of the function of the package MASS confint.glm even after the package has been detached! Why the exclamation mark? Note profile.glm is not actually in package:MASS (sic). Try looking for it:

Re: [R] Problem with parser and if/else

2003-11-13 Thread Prof Brian D Ripley
On Thu, 13 Nov 2003, Ben Bolker wrote: In the second case, R stops when it has a syntactically complete clause: if (...) { ... } is complete since an else{} clause is not required. R evaluates it, then moves onto else { ... } which is a syntax error (since it doesn't have an if

Re: [R] GDB under windows

2003-11-06 Thread Prof Brian D Ripley
It's in the rw-FAQ! That describes this as `fraught' with good reason. On Thu, 6 Nov 2003, Marsland, John wrote: Does anybody have some simple instructions to get me going using the GNU debugger GDB under windows with R? The R Extensions manual mentions briefly debugging under unix, but I've

Re: [R] map does not display maps, MacOSX

2003-11-05 Thread Prof Brian D Ripley
Notice the // in the path /Users/glaziou/Library/R/maps/mapdata//world.N Some Windows filesystems do not like that, and my guess is that some MacOS X ones may not either. On Wed, 5 Nov 2003, Philippe Glaziou wrote: Hi, I installed the maps and mapdata libraries on my R-1.8.0 on MacOSX

Re: [R] run R under linux

2003-10-22 Thread Prof Brian D Ripley
On Wed, 22 Oct 2003, Zhen Pang wrote: Thank you for your help. But try() seems to only allow for an expression. My simulation have serveral expressions which possibly have problem. Is there any possibility to include them all in the try()? If not, I seem to use several try(). {line1 line2

Re: [R] png() and/or jpeg(): line missing by using box(which=outer)

2003-10-21 Thread Prof Brian D Ripley
It is probably a bug: does it happen when you copy from the screen in png? I would expect not, hence that may be a workaround for you. When I have both time and access to a Windows machine I may be able it take a closer look: meanwhile you do have access and have the source code so please

Re: [R] ks.test()

2003-08-28 Thread Prof Brian D Ripley
On Thu, 28 Aug 2003, franck allaire wrote: With the Shifted Exponential test, H_0 is data is a sample coming from a Shifted Exponential distribution with shift=30 and lambda = 0.001566907 You got that after looking at the data. H_0 has to be specified before looking at the data. You appear

Re: [R] Marginal (type II) SS for powers of continuous variables ina linear model?

2003-08-14 Thread Prof Brian D Ripley
On Tue, 12 Aug 2003, [iso-8859-1] Bjørn-Helge Mevik wrote: Prof Brian Ripley [EMAIL PROTECTED] writes: drop1 is the part of R that does type II sum of squares, and it works in your example. So does Anova in the current car: I'm sorry, I should have included an example to clarify what I

Re: [R] label chemical formulas

2003-08-04 Thread Prof Brian D Ripley
On Mon, 4 Aug 2003, Roberto Canteri wrote: [This is a repeat of a post earlier today.] Dear r-help, the prcomp matrix (x$rotation) I have, has the columns labels as chemical formulas (i.e: C2H3, 37Cl,...,C10H24N3O3). Plotting them as plain text I do the following with success: xlabs -

Re: [R] rterm

2003-07-18 Thread Prof Brian D Ripley
Well, there's a file README.rterm in the top-level (binary) installation directory, and the section `Invoking R under Windows' in `An Introduction to R' applies equally to RGui and Rterm. All the information is in one of those two places. On Fri, 18 Jul 2003, Manuel Lopez Coello wrote: I

Re: [R] sample function

2003-06-27 Thread Prof Brian D Ripley
R is not S-PLUS, and you need Modern Applied Statistics in S (4th ed) for a description including R. sample in R used a PRNG: see ?RNG in R for the details of PRNGs in R. On Fri, 27 Jun 2003, [iso-8859-1] Ramzi Feghali wrote: i have a question about the sample function used in R, does it work

Re: [R] using R as a script; getting the output to display on screen

2003-06-27 Thread Prof Brian D Ripley
On Fri, 27 Jun 2003, Nick Helseth wrote: So, you are using Windows: please say so up front. I'm trying to run R from the command line using a function as an input; it seems to run ok when I use Rterm --vanilla name_of_R_function.R (by ok, I mean that it will say that it's reading in the

Re: [R] Can't save a graph to pdf in R for MacOS

2003-06-26 Thread Prof Brian D Ripley
On Thu, 26 Jun 2003, p.b.pynsent wrote: I do not have an R solution but use eps2pdf (a Perl script) This can be installed by Fink http://sourceforge.net/projects/fink/. eps2pdf runs GhostScript: so does the R device driver bitmap(). Does the latter work on your system? If so it would save

Re: [R] outlier

2003-06-18 Thread Prof Brian D Ripley
Please do read the help page. which explains this is a random algorithm. In your example you can try cov.rob(ap, cor=TRUE, nsamp=exact) On Wed, 18 Jun 2003, kan Liu wrote: I wrote a .R file (see below)to calculate robust measure of correlation using cov.rob. I got different correlation

Re: [R] dyn.load() function problem help!

2003-06-18 Thread Prof Brian D Ripley
On Wed, 18 Jun 2003, Mu Yunming wrote: Hi! I would like to call a Fortran subroutine within R. For doing this, I first built a shared library for loading into R under Unix and it was successful. How did you do this? Did you use R CMD SHLIB? But when I tried to load a shared library using

RE: [R] filtering ts with arima

2003-04-09 Thread Prof Brian D Ripley
I know that there is no direct equivalent to arima.filt in R. I also know that your R calls are *not* equivalent to the S-PLUS ones. I would use S-PLUS to run S-PLUS code: you are a commercial operation and should be able to afford it, Or you could employ an S programmer to convert the code for

RE: [R] command line limit?

2003-03-25 Thread Prof Brian D Ripley
On Tue, 25 Mar 2003, dederderian wrote: Very insightful. Thanks Martin. Maybe, I need to explain my problem a little better. Yes, please do: you seem to be talking about console input lines, not the command line at all. I am creating an R comfile within a perl script. When I run this R

Re: [R] na.action in model.tables and TukeyHSD

2003-02-28 Thread Prof Brian D Ripley
ndata - na.omit(led1t7sts) and work with ndata. Why is that difficult? On Fri, 28 Feb 2003, CG Pettersson wrote: In 27/2, I got the following answer from Prof. Ripley: (The question is at the bottom) This ia already fixed in R-devel. The answer is the same: don't use na.omit implicitly:

Re: [R] Why does 'exists' need a quoted argument?

2003-02-20 Thread Prof Brian D Ripley
foo - bar exists(foo) asks if the object whose name is bar exists, not that whose name is foo exists. There is a distinction between operating on objects (rm and is.na) and names of objects (get, exists, help). The one exception in your list is help. That should really be help(rm) (which

Re: [R] multivariate sampling question again

2003-02-17 Thread Prof Brian D Ripley
R-devel already contains rmultinom. I don't see why rejection sampling on (m, n) should be slow unless m-n is very small or far out in the tails, in which case why do you want this? In any case, inversion will be a perfectly adequate method as you have pbeta and qbeta. On Mon, 17 Feb 2003, Peng

Re: [R] multilm for simseg/acm

2003-02-10 Thread Prof Brian D Ripley
On Mon, 10 Feb 2003, Christian Schulz wrote: hi, for working with the simseg/acm approach i need multilm and it seems that the last version is 0.1-4.tar ? What's wrong in the description file, if i attempt install the package ? c:\DataMining\rw1062\binRcmd Install multilm Malformed DCF

Re: [R] postgres/R access problems

2003-02-09 Thread Prof Brian D Ripley
On 10 Feb 2003, Brendan Murray wrote: I have a problem that has exhausted my ingenuity and would like pointers to a solution, or at least where to debug. - I am using R and postgresql. And which R package are you using, which version, from which repository? -- Brian D. Ripley,

Re: [R] floating point question

2003-01-31 Thread Prof Brian D Ripley
It's a difference in the `libc'. Asking for more precision than the arithmetic has is asking for fairly random results. The differences are as likely to be in the *printing* as in the computations. On Fri, 31 Jan 2003, Bob Gray wrote: Does anyone know precisely what is different about the

Re: [R] printing reals from C with digits

2003-01-29 Thread Prof Brian D Ripley
On Wed, 29 Jan 2003, Ott Toomet wrote: I want to print real numbers in C code with different values for digits. How to do that? Use Rprintf or PrintValue. You'll need to work hard to convince me that Rprintf is not adequate. As I have understood, what I should do is to call

Re: [R] Logistic regression: At times correlation matrix ofcoefficients gets messed up

2003-01-21 Thread Prof Brian D Ripley
It's not messed up, just someone's idea of a compact display. Options are 1) Use vcov(fit) instead 2) Use print(summary(fit), symbolic.cor=FALSE) Does anyone think that the current arrangement (use this scheme for more than 4 coefficients) is sensible? Surely the abbreviations are not (( for

Re: [R] read.table(file=clipboard,...) for R?

2003-01-20 Thread Prof Brian D Ripley
On Mon, 20 Jan 2003, Prof Brian D Ripley wrote: On Mon, 20 Jan 2003, David Parkhurst wrote: I use read.table(file=clipboard,...) a lot in s-plus (under windows 2000), but it does not seem to work in R (and is not in the help screen for read.table). Am I missing something? Would