Re: [R] Lattice levelplots and (partial) failure on some devices----Long

2003-06-13 Thread Prof Brian Ripley
This is a Windows run-time quirk (bug): Linux gives nan. It's a problem with what grid/lattice sends the driver. Rather than put protection code in each driver, it needs to be in the central graphics code. We've been here several times before On Thu, 12 Jun 2003, Joel Kincaid wrote: > >

Re: [R] Layout of windows devices

2003-06-13 Thread Prof Brian Ripley
There is already support for this in the R-devel version: it like the current approach assumes one screen, though. From the CHANGES file The initial size and position of the MDI frame can be set in Rconsole: see the comments in .../etc/Rconsole. The initial position of the console window

RE: [R] Layout of windows devices

2003-06-13 Thread Prof Brian Ripley
Oh, I see now you want to do this differently for each instance of a device. That would easy to add via arguments to windows(), as the internal code is already there. On Fri, 13 Jun 2003, David Khabie-Zeitoune wrote: > Duncan -- thanks for your reply. As you point out -- this is more a > "nice-to

Re: [R] Using jpeg() function over cgi

2003-06-13 Thread Prof Brian Ripley
Yes, it has been solved and discussed in the R-help archives many times. The help(jpeg) page is pretty explicit too. One alternative is bitmap(). To use jpeg() under a Unix-alike you need to set up an X server that your R process can use. It's a bit hard to help you do that when you don't even m

Re: [R] problem installing packages from source on win2k

2003-06-14 Thread Prof Brian Ripley
The short answer is that you are not following readme.packages, which says If you want to make compiled html (.chm) files you will need the Microsoft HTML Help Workshop, currently available for download at http://msdn.microsoft.com/library/en-us/htmlhelp/html/hwmicrosofthtmlhelpdownloads.as

Re: [R] stop criterion for stepAIC

2003-06-16 Thread Prof Brian Ripley
You completely misunderstand: have you even *looked* at the help page? Description: Performs stepwise model selection by exact AIC. ! If you don't know what AIC is (and it seems so), please research it (e.g. in MASS the book whose support software you are using without any credit). On Mo

RE: [R] Error in .C("open_netcdf"

2003-06-16 Thread Prof Brian Ripley
Then the remedy is simple: do as the rw-FAQ and ReadMe's suggest and compile it from the sources yourself. Something is obviously different between your setup and mine. On Mon, 16 Jun 2003, antonio rodriguez wrote: > Hi, > > Well, I didn't compile from sources R_1.7.0 just downloaded the execu

Re: [R] Re: R 1.7.1 is released

2003-06-16 Thread Prof Brian Ripley
On Mon, 16 Jun 2003, Dan E. Kelley wrote: > Peter Dalgaard BSA wrote: > > >I've rolled up R-1.7.1.tgz a short while ago. > > > > > Thanks, Peter. This is terrific! I find that this R installation is > only accessible to the root. In case it helps other folks, I've > attached a patch on th

Re: [R] R1.7.0 on Solaris 9

2003-06-16 Thread Prof Brian Ripley
On Mon, 16 Jun 2003, A.J. Casadevall wrote: > On Solaris 9 and gcc, R 1.7.0, configure fails. > What to do ? Give us a hint as to what the problem is? Also, it will be easier to help you with the current version, 1.7.1. -- Brian D. Ripley, [EMAIL PROTECTED] Professor of Applied

Re: [R] cut.POSIXct problem

2003-06-17 Thread Prof Brian Ripley
This works in R 1.7.1: please upgrade. [That's assuming that your sequence actually covers the data: you need to go beyond max(cas.param) to be sure, and it looks to me that with right=FALSE and include.lowest-FALSE the breaks never will cover the data.] On Tue, 17 Jun 2003, Petr Pikal wrote:

Re: [R] hist density...

2003-06-17 Thread Prof Brian Ripley
A density integrates to one, so the total area of the bins is one. It's your `Thought' which is incorrect. For hist(freq=FALSE) the area of each rectangle (not its height) represents the proportion of the data falling into the base of the rectangle. On Tue, 17 Jun 2003, wolski wrote: > Hi! > D

Re: [R] help.start

2003-06-17 Thread Prof Brian Ripley
On Tue, 17 Jun 2003, W. C. Thacker wrote: > Dear R experts, > > I would greatly appreciate your help with R's hypertext online > documentation, which I can pass on to our system manager. We are > running version 1.7 on sparc-sun-solaris2.9 using netscape7 browser > for displaying the help. The

Re: [R] vector maximum length

2003-06-17 Thread Prof Brian Ripley
Depends on the R platform, but much larger than that. On most systems the maximum length will be 2^31-1 (over 2 billion), but on 32-bit systems the storage needed for the vector will hit its limit first. You can expect to be able to use a small number of hundred million elements, depending on

Re: [R] kernel smoothing for ordinal data

2003-06-17 Thread Prof Brian Ripley
They do exist: the term has a somewhat different meaning for categorical data. Titerington, D. M. (1980) Technometrics 22, 259-268 might be a good start. On 17 Jun 2003, Peter Dalgaard BSA wrote: > Meike Gebel <[EMAIL PROTECTED]> writes: > > > Hi there, > > during my work I have to use kernel

Re: [R] outlier

2003-06-17 Thread Prof Brian Ripley
On Tue, 17 Jun 2003, kan Liu wrote: > I want to calculate the R-squared between two variables. Can you advice > me how to identify and remove the outliers before performing R-squared > calculation? Easy: you don't. It make no sense to consider R^2 after arbitrary outlier removal: if I remove a

Re: [R] make check and postscript output in R 1.7.1

2003-06-17 Thread Prof Brian Ripley
It's telling us that your machine gives different output from *all* the other machines that were checked in the beta-test period. So what machine is this that no one beta tested? (In future please do participate in the beta test so your unusual system gets tested.) As Rob G has reported thi

Re: [R] outlier

2003-06-17 Thread Prof Brian Ripley
e indicate very clearly that the assumptions are > > incorrect. See > > Anscombe (1973) "Graphs in Statistical Analysis", > > The American > > Statistician, 27: 17-22, reproduced in, e.g., du > > Toit, Steyn and Stumpf > > (1986) Graphical Exploratory

Re: [R] Query: Sd2Rd and nroff macros in S docs

2003-06-18 Thread Prof Brian Ripley
On Wed, 18 Jun 2003 [EMAIL PROTECTED] wrote: > Documentation for S3 functions is apparently written in troff markup > with macro tags like > > .BG .FN .TL .DN .CS ... > > Inspection of S3 documentation source files gives a pretty clear idea > of what these mean, semantically (and Sd2Rd is

Re: [R] downloading packages and AntiVirus program

2003-06-18 Thread Prof Brian Ripley
I've just run that file past two different virus scanners (Sophos and Norton, both fully updated) wiht no hits. I think you can reasonably ask your sysadmins to demonstrate to you what the virus is and in which file in the zip it is. I suspect they do not check ftp downloads, so you could try t

Re: [R] Forward stepwise procedure w/ stepAIC

2003-06-18 Thread Prof Brian Ripley
It starts with the initial model you gave it, which contains all your explanatory variables already. You presumably wanted to use a simpler model to start with. As in: Arguments: object: an object representing a model of an appropriate class. This is used as the initial model in t

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

2003-06-18 Thread Prof Brian Ripley
On Wed, 18 Jun 2003, Thomas Lumley wrote: > On Wed, 18 Jun 2003, Mu Yunming wrote: > > > > > > My source code is writen in fortran 90 and i complied the source code using > > command f90. How do you think I should complie my source file and build the > > shared library? > > > > You probably need

Re: [R] 3-way Interactions w/ stepAIC

2003-06-18 Thread Prof Brian Ripley
On Wed, 18 Jun 2003 [EMAIL PROTECTED] wrote: > I'm attempting to use stepAIC to select a model through a forward procedure. > I want to consider up to all 3-way interactions. > > > > I've attempted to use the following code: > > > > m2.Fwd3way <- stepAIC(m1.Ionly, direction="forward", > >

Re: [R] Ltsreg and nsamp="exact"

2003-06-18 Thread Prof Brian Ripley
Are subsets of your data collinear? On Wed, 18 Jun 2003, Christopher Adolph wrote: > > I'm trying to use least trimmed squares using ltsreg with nsamp="exact". > When I use the following: > > rg <- ltsreg(x,y,nsamp="exact") > > I get: > > Error in lqs.default(x, y, nsamp = "exact", method = "

Re: [R] Multiple graph sheets

2003-06-18 Thread Prof Brian Ripley
On Wed, 18 Jun 2003, Paul, David A wrote: > In Splus the code > > test.lm <- lm(y ~ x, data = test.data) > plot(test.lm) > > generates a graphics window that contains > multiple graph sheets that one may choose > from via the "page" tabs at the bottom of > the window. > > Is there a way to do

Re: [R] What's wrong with ar for my data?

2003-06-18 Thread Prof Brian Ripley
You are apparently fitting a series for which the selected order is zero, and ar.burg is not designed to cope with that (and would in any case tell you nothing useful). The default method does cope, from your output. Why are you fitting an AR model to a series with apparently no correlation? O

Re: [R] Grouping binary data

2003-06-19 Thread Prof Brian Ripley
On Thu, 19 Jun 2003, Henric Nilsson wrote: > Dear all, > > I'm analyzing a binary outcome using glm() with a binomial distribution and > a logit link, and have now reached the point where I'd like to do some > model checking. Since my data are in binary form I'd like to collapse over > the cro

Re: [R] clines loading error

2003-06-19 Thread Prof Brian Ripley
I think your R is too old. Packages can indicate they need a particular version of R, but it is unreasonable to expect them to be tested on versions predating the package (which as clines is not on CRAN I can't check but suspect to be the case). On Thu, 19 Jun 2003, Karim Elsawy wrote: > I'm

Re: [R] saving window size of a graphics device in windows possible?

2003-06-19 Thread Prof Brian Ripley
On Fri, 20 Jun 2003, Masayoshi Hayashi wrote: > In MDI mode of Windows, is it possible to save window size and other > parameter settings of a graphics device when it is activated, much like > Rconsole file? No, but you don't need to. Size is part of the arguments of the windows() device, and yo

Re: [R] namespaces not available

2003-06-20 Thread Prof Brian Ripley
I think you have saved in your workspace an object whose environment is a namespace, probably mva, that is not available when .RData is loaded but is once the session is running. On Fri, 20 Jun 2003, Göran Broström wrote: > This happened to me with R-1.7.0, Linux RH9: > > --

Re: [R] Spedd: R vs S-plus

2003-06-20 Thread Prof Brian Ripley
On Fri, 20 Jun 2003, Uwe Ligges wrote: > (Ted Harding) wrote: > > Hi Folks, > > > > Sorry to raise what has probably been discussed before, > > but I an repeatedly struck by the comparative slowness > > of S-plus for Windows compared with R for Linux when doing > > much the same thing. > > > > I

Re: [R] Silhouette question

2003-06-20 Thread Prof Brian Ripley
?try is your friend here. On Sat, 21 Jun 2003, Jonck van der Kogel wrote: > Hi all, > I am momentarily experimenting with Silhouette from the cluster library > but I am getting some errors. Since Silhouette can be seen as a quality > measure for a clustering what I want to do is run a series of

Re: [R] error functions

2003-06-21 Thread Prof Brian Ripley
Yes, statisticians call the natural versions of the cumulative normal distribution, pnorm. As I recall erf(x) = 2 * pnorm(x * sqrt(2)) - 1 erfc(x) = 2 * pnorm(x * sqrt(2), lower=FALSE) On Sat, 21 Jun 2003, Salvatore Barbaro wrote: > does anybody know if R contains error functions like > "erf

Re: [R] optim with contraints

2003-06-21 Thread Prof Brian Ripley
Adelchi, R is a volunteer project, and we would need a volunteer to look at this. `The developers of optim' (who are R-core) used published code for this method. Since you are getting different answers on different platforms this might be a bug in your compiler or run-time rather than R. Please

Re: [R] core dump on solaris 2.9

2003-06-22 Thread Prof Brian Ripley
[I presume this is sparc-Solaris, but you did not say.] Please do as the INSTALL file asks and read the R-admin manual: this *is* documented there! That version of gcc is obselete, and unusable on sparc-Solaris. You will have to find a working compiler: if all else fails, try turning optimizatio

Re: [R] clicking the stop button

2003-06-22 Thread Prof Brian Ripley
This is PR#3285, fixed in R-patched and R-devel. Search on R-bugs for the full details. On Sun, 22 Jun 2003, syed gillani wrote: > Clicking the stop button results in failure to accept any further > keyboard input.Is it just my machine? > > RGui: 1.7.0 & 1.7.1 > WindowsXp Prof. > Pentium4 2.5

Re: [R] Cross-compiling R packages

2003-06-22 Thread Prof Brian Ripley
On Mon, 23 Jun 2003, Laurent Gautier wrote: > I tried to use the 'make' file discussed in the last R-news. > > The step 'make R' dies with: > make[2]: Entering directory > `/home/laurent/these/R/RCrossBuild/WinR/R-1.7.1/src/gnuwin32' > sed -e s/@RVER@/`cut -d' ' -f1 ../../VERSION | sed -n 1p`/g

Re: [R] right assignment ("->") and functions

2003-06-22 Thread Prof Brian Ripley
On Mon, 23 Jun 2003, Robin Hankin wrote: > Hi everyone > > check this out [R-1.7.0]: > > R> f1 <- function(x){x^2} > R> f1 -> f2 > R> f2(4) > [1] 16 > R> > R> function(x){x^2} -> f3 > function(x){x^2} -> f3 > R> f3(4) > Error: couldn't find function "f3" > > Why does right assignment "->" work

Re: [R] How can I do a spinning plot in R?

2003-06-22 Thread Prof Brian Ripley
There are interfaces to XGobi (the xgobi package on CRAN) and GGobi (RSGGobi from Omegahat, and I have a variation on xgobi()) to allow spinning of point clouds. You can spin surfaces in the R<-->GL packages, e.g. Duncan Murdoch's at http://www.stats.uwo.ca/faculty/murdoch/software/default.htm a

RE: [R] "cloud" function doesn't plot

2003-06-23 Thread Prof Brian Ripley
On Sun, 22 Jun 2003, Tony Ross wrote: > >From: "Simon Blomberg" <[EMAIL PROTECTED]> > >Date: Mon, 23 Jun 2003 16:26:39 +1000 > > >If you type ?cloud, you can see that you need to provide a formula as well > >as a dataset. So if your variables are x,y,z in dataframe dat, then > >cloud(z~x*y, da

Re: [R] Problem with installation under R 1.7.1 and unzip

2003-06-23 Thread Prof Brian Ripley
On Mon, 23 Jun 2003, Ramzi Feghali wrote: > Dear all, > > i was working under R 1.6.2 but i wanted to install bioconductor under WinXP, and > because R 1.7.1 contains the TclTk package and automatic installation of > bioconductor packages, so i install it. > > Everything is fine except this

Re: [R] read.spss

2003-06-23 Thread Prof Brian Ripley
The message is quite correct. See the rw-FAQ, Q2.14, or the FAQ Q7.10. On Mon, 23 Jun 2003, Harold Doran wrote: > I have loaded the foreign package and am still having problems with an > import. I get a message that reads, unable to open file. Whe I try > different files I get the same message. H

Re: [R] erase.screen bug?

2003-06-23 Thread Prof Brian Ripley
Is your background colour set to transparent? There is no way to erase screens in the R/S model: all you can do is overpaint with the background. And the default on-screen background colour is ... transparent. This described in the Warning (and elsewhere) on the help page. With a solid backgrou

Re: [R] ?plot problem

2003-06-23 Thread Prof Brian Ripley
On Mon, 23 Jun 2003, Paul, David A wrote: > R1.7.0, Win2k: > > When I use plot( ) on a groupedData object, > if I have 165 subjects I'm supposed to be able > to use plot(..., layout = c(5,3,11)) to get > all 165 on 11 separate sheets. The graphics > window is only displaying the first 10. If I

Re: [R] Lwd ignored when printing on Windows

2003-06-23 Thread Prof Brian Ripley
What printer driver are you using? I've just tried this and it works exactly as one would expect on my HP 970CXi, as well as cut-and-paste into other applications. It also worked printing to Acrobat Distiller (although all the lines were thinner there than on-screen and on the 970CXi, the rati

Re: [R] Warnings using MASS Library

2003-06-23 Thread Prof Brian Ripley
On Mon, 23 Jun 2003, Williams, Brian wrote: > I am running R version 1.7.1. In the process of checking the results from > some code against the results obtained from version 1.5.1, I encountered a > set of warning messages upon exiting from the session: > > Save workspace image? [y/n/c]: y > War

Re: [R] help on R programming.

2003-06-23 Thread Prof Brian Ripley
On Mon, 23 Jun 2003, Murad Nayal wrote: > - what is the correct way to -remove- a component from a list. this > seems to do the trick: list[[1]] = NULL, however, you'd think this > should simply attach a NULL object at the first component position? This is in the FAQ, section 3.3.3, and is an S

Re: [R] Creating Surface Plots in R

2003-06-23 Thread Prof Brian Ripley
?persp library(lattice) ?wireframe both have several examples. On Tue, 24 Jun 2003, Neil Osborne wrote: > I have three columns of data in afile that I would like to analyse by means > of a surface plot. I would be very grateful if anyone could show me how to > create a surface plot from this

Re: [R] excel files and R

2003-06-24 Thread Prof Brian Ripley
On Tue, 24 Jun 2003, Victor H. Marím wrote: > I am new at R. My questions is rather basic. Looking R manuals looks > like there should be a way to read MS excel files into R. Could > somebody tell me which library should I use for that? Several ways are in the R Data Import/Export Manual (the

Re: [R] ?plot problem

2003-06-25 Thread Prof Brian Ripley
ot is recorded if you move away from it or start a new plot. On Tue, 24 Jun 2003, Paul, David A wrote: > Thanks to Dr. Thomas Hotz, Prof. Brian Ripley, Dr. Dennis > Murphy, and Dr. David Scott for their replies. > > Trying an idea: > > > x <- c(1,2,3,4) > > y<-c(2,4.2

RE: [R] excel files and R

2003-06-25 Thread Prof Brian Ripley
f, but not for batch > processing. > Please someone tell me what I missed - how can I open the connection > within > an R script? > (Windows 98)(not my fault) > > TIA > (sorry for long disclaimer, can't switch it off) > > -----Original Message- > From: Pr

RE: [R] excel files and R

2003-06-25 Thread Prof Brian Ripley
On Wed, 25 Jun 2003, Simon Fear wrote: > Many many thanks. I did look at this but I have absolutely no idea of > the > background so got completely lost. Can you recommend a gentle > introduction/overview to this area, based on the assumption that my > current > knowledge equals zero? Indeed, coul

Re: [R] Execution of R code

2003-06-25 Thread Prof Brian Ripley
I am not sure I fully understand the Qs. There are two phases. 1) The source code is parsed. 2) The parsed code is evaluated. If you run code from source() or a file or the command line, it is parsed and evaluated. However, evaluating a function assignment makes an function object containing t

Re: [R] creating R help page

2003-06-25 Thread Prof Brian Ripley
See the `Writing R Extensions' manual. On Wed, 25 Jun 2003, Vincent Spiesser wrote: > Does anybody know how to create an R help page wich can be opened from R > console (like help(glm)) ? > Particularly, I would like to know : > - what kind of file the help file is ? > - Where does it take place

Re: [R] Help on using read.table with files containing dates

2003-06-25 Thread Prof Brian Ripley
On Wed, 25 Jun 2003, Chriss, Neil wrote: > I am a relatively new user of R and have a question about using the > read.table command for reading in .csv files that contain dates. Unfortuantely them appear to contain perversions of dates, not the dates recognised by the ISO standard. I read you

Re: [R] Mac carbon - foreign - read.spss

2003-06-25 Thread Prof Brian Ripley
On Wed, 25 Jun 2003, Jean-Pierre Muller wrote: > Can someone confirm that "read.spss" doesn't work with > usual (old) mac paths ("disk:dir:dir:file"), but only with > *nix path (Volumes/disk/...)? > > rm171 (carbon) - MacOSX 10.2.6 The exact code used to open files is (variations on) FIL

Re: [R] Help

2003-06-25 Thread Prof Brian Ripley
On Wed, 25 Jun 2003, Tomer Maymon wrote: > my name is Tomer,im using R version 1.4 As the current version is 1.7.1 (and there never was a 1.4), you need to upgrade. > i've encontered a couple of problems with my my R programing. > 1. i cant understand what is the whights in the nnet$wts command

Re: [R] rw1062

2003-06-25 Thread Prof Brian Ripley
The problem is that you have not put the rw1071/bin directory in your path, and so your DCOM-using application is unable to find the dependent DLLs of R. This is a problem with the R (D)COM instructions. On Wed, 25 Jun 2003, Søren Merser wrote: > I need 'rw1062.zip' since i can't get excel, R1

Re: [R] logLik.lm()

2003-06-25 Thread Prof Brian Ripley
Your by-hand calculation is wrong -- you have to use the MLE of sigma^2. sum(dnorm(y, y.hat, sigma * sqrt(16/18), log=TRUE)) Also, this is an inappropriate use of AIC: the models are not nested, and Akaike only proposed it for nested models. Next, the gamma GLM is not a maximum-likelihood fit un

Re: [R] robust regression

2003-06-25 Thread Prof Brian Ripley
You can use the quantreg package. However, neither l1fit nor that do `robust regression', so you need to think more carefully about what you really want. There are almost always better alternatives than L1 fits. On Wed, 25 Jun 2003, Rafael Bertola wrote: > Is there a command in R that make th

Re: [R] logLik.lm()

2003-06-25 Thread Prof Brian Ripley
On Wed, 25 Jun 2003, Spencer Graves wrote: > Dear Prof. Ripley: > > I gather you disagree with the observation in Burnham and Anderson > (2002, ch. 2) that the "complexity penalty" in the Akaike Information > Criterion is a bias correction, and with this correction, they can use > "den

Re: [R] probelem of function inside function

2003-06-25 Thread Prof Brian Ripley
On Wed, 25 Jun 2003, array chip wrote: > Actually, It is the question I encountered in S-Plus. > Sorry that I forgot to mention. As Simon just pointed > out, the function works fine in R (Thanks!). But in > any case, does anyone know how to solve the problem in > S-Plus? Assign to frame 1: see `S

Re:[R] extracting effects estimates from aov model with Error()

2003-06-25 Thread Prof Brian Ripley
I think you are applying effects to the summary, but you would need to apply it to each stratum in turn, e.g. effects(fit[[2]]) Consider example(aov) effects(npk.aov[[3]]) On Wed, 25 Jun 2003, Peter Ho wrote: > Hi all, > > Does anyone know how I might be able to extract effects estimates fr

RE: [R] Change default parameters of panel.smooth

2003-06-26 Thread Prof Brian Ripley
On Thu, 26 Jun 2003, Dowkiw, Arnaud wrote: > Thanks Simon, > > in fact, problems occur when I use panel.smooth inside pairs, not when I > use panel.smooth alone : > > > pairs(Fingaroy.F3.cross3.Streetonout.df[,c(5:9)],main="Kingaroy -- F3 -- > > Cross3",lower.panel=panel.smooth,upper.panel=pane

Re: [R] DB2 ODBC

2003-06-26 Thread Prof Brian Ripley
On 26 Jun 2003, Steffen Neumann wrote: > Babu Prathap R <[EMAIL PROTECTED]> writes: > [...] > > UDB 8.1 using R. I understand that I need to use R ODBC. But I would like > > What OS ? > > On Linux/Solaris you need unixODBC > (or maybe alternatively iODBC) installed. > > Once that works you in

Re: [R] Plots using POSIX

2003-06-26 Thread Prof Brian Ripley
On Thu, 26 Jun 2003, Shawn Way wrote: > Is there a reason that the bottom axis changes color when POSIX data is used > in plot function? It's not the same plot function, that's why. > For example: > > > timedata <- c("2/3/2003","3/4/2003","5/4/2003") > > timedata2 <- strptime(timedata,format="%

Re: [R] Plots using POSIX

2003-06-26 Thread Prof Brian Ripley
On Thu, 26 Jun 2003, Duncan Murdoch wrote: > On Thu, 26 Jun 2003 07:59:00 -0500, Shawn Way <[EMAIL PROTECTED]> wrote : > > > > >Is there a reason that the bottom axis changes color when POSIX data is used > >in plot function? > > It's the old problem of too much of ... being passed onwards. Her

Re: [R] problems with library in 1.7.1

2003-06-26 Thread Prof Brian Ripley
Direct use of .lib.loc was deprecated in 1.6.0. Use .libPaths() instead. .Library is looked for package:base, and your assigning a copy in .GlobalEnv will make no difference. On Thu, 26 Jun 2003, R. Heberto Ghezzo wrote: > Hello. > I am using R 1.7.1 just downloaded on Win98. With the old 1.6.

Re: [R] lm diagnostics and qr (fwd)

2003-06-26 Thread Prof Brian Ripley
On Thu, 26 Jun 2003, Jean Eid wrote: > I have been struggling to find some informaation on what lm exactly does. > I know it uses the QR decomp. However, I was recently faced with a > somewhat badly scaled matrix and summary(lm) said > Coefficients: ( 4 not defined because of singularities)

RE: [R] assignment in lists

2003-06-26 Thread Prof Brian Ripley
Philippe, as.list(NULL) is the same as list(), and that is what I think you should be using in both cases. However, I do think that either both or neither of your examples should work: my preference would be `neither' but as S allows both it should be `either'. Brian On Thu, 26 Jun 2003, Philip

RE: [R] dropping factor levels in subset

2003-06-26 Thread Prof Brian Ripley
A more transparent solution is old.factor[1:3, drop = TRUE] That has worked for a long time, but apparently not been documented in R until 1.7.1 (docs added a couple of hours before release). So you could do (probably, since there are some bugs prior to 1.8.0) crb[] <- lapply(crb, function(x) x[

RE: [R] assignment in lists

2003-06-26 Thread Prof Brian Ripley
On Fri, 27 Jun 2003, Philippe Grosjean wrote: > Prof. Brian Ripley wrote: > >Philippe, > > >as.list(NULL) is the same as list(), and that is what I think you should > >be using in both cases. > > OK, thank you. > > >However, I do think that either

Re: [R] NA points in loess function

2003-06-27 Thread Prof Brian Ripley
On Fri, 27 Jun 2003, Yao, Minghua wrote: > Gurus, > > I used > predict(loess(Y~X)); > > where Y and X are of the same length. But there are same NA's in both Y and > X. Those NA's are in the same locations in Y and X. The following is the > error messageI got: > > Error in "[<-"(*tmp*, ,

Re: [R] feffects in model.tables

2003-06-27 Thread Prof Brian Ripley
On Fri, 27 Jun 2003, Peter Ho wrote: > Is there is there an equivalent of S-plus "feffects" argument in > model.tables() in R for calculating effects of a factorial design? No (as reading the help page would have told you). > Or is > there another way they can be calculated from an aov fit.

Re: [R] connection DB2 and R (R under linux)

2003-06-27 Thread Prof Brian Ripley
On Fri, 27 Jun 2003, laurent buffat wrote: > I would like to access to DB2 on R on a linux machine (R and DB2 are on 2 pc > linux ) > I have read on the R-archive that RODBC work with DB2, but it is not clear > for me, if it's only for R on a windows machine or if it's work also for R > under lin

RE: [R] dropping factor levels in subset

2003-06-27 Thread Prof Brian Ripley
n.) On Fri, 27 Jun 2003, Marc Schwartz wrote: > >-Original Message- > >From: [EMAIL PROTECTED] > >[mailto:[EMAIL PROTECTED] On Behalf Of Prof > >Brian Ripley > >Sent: Friday, June 27, 2003 1:35 AM > >To: Marc Schwartz > >Cc: [EMAIL PROTECTED];

RE: [R] NA points in loess function

2003-06-27 Thread Prof Brian Ripley
On Fri, 27 Jun 2003, Yao, Minghua wrote: > Thank you, Prof. Ripley. By adding na.action=na.exclude to loess function, > the program runs great. > > I did omit some trivial things. But I still think it was NA's that caused > troubles in predict or loess. *NO* the errors message is from your usag

Re: [R] make

2003-06-27 Thread Prof Brian Ripley
As the R-admin manual says, you can use either the GNU make or the Sun native make which is in /usr/ccs/bin. The phrase `does not work' is dreaded by all technical support personnel: it is maximally uniformative. So please try again to explain the actual problem! On Fri, 27 Jun 2003, Roger D.

Re: [R] logLik.lm()

2003-06-27 Thread Prof Brian Ripley
; model from the "true" model, it is not clear that the models have to be > nested. > > Any thoughts or comments on this issue? > > Best, > Ravi. > > > - Original Message - > From: Prof Brian Ripley <[EMAIL PROTECTED]> > Date: Wedn

Re: [R] histogram cells

2003-06-27 Thread Prof Brian Ripley
?hist says breaks: one of: ... * a single number giving the number of cells for the histogram, ... In the last three cases the number is a suggestion only. On Fri, 27 Jun 2003, Tommy E. Cathey wrote: > Why does the following code generate a Histogram

Re: [R] question

2003-06-28 Thread Prof Brian Ripley
On Sat, 28 Jun 2003, Jan de Leeuw wrote: > Is it correct that > > lsfit(matrix(0,10,1),1:10,intercept=FALSE) > > returns zero residuals ? No. I would use lm.fit for such problems, and that seems to behave (probably rather coincidentally). -- Brian D. Ripley, [EMAIL PROTECTE

Re: [R] How could plot a chart with some Chinese's words text?

2003-06-29 Thread Prof Brian Ripley
On Mon, 30 Jun 2003, Ken Lee wrote: > Dear all, > I want to plot a chart title with some Chinese words at unix system, but I do > not how to do. > first, it is no problem when I use write.table like: > > a<-"中文" > write.table(a,file="z.txt",sep="\t") Not, that is not Chinese for me and

Re: [R] ./configure needs /sw/lib?

2003-06-30 Thread Prof Brian Ripley
On Sun, 29 Jun 2003, Tony Marlboro wrote: > Hello, > > I am using R on Mac OS X. I have tried to install the package > "netCDF", but have run into a problem. The install.packages command > fails during a call to the "configure" shell script in the package > build directory, because it can

Re: [R] Fitting inter-arrival time data

2003-06-30 Thread Prof Brian Ripley
On Sun, 29 Jun 2003, M. Edward Borasky wrote: > I have a collection of data which includes inter-arrival times of requests > to a server. What I've done so far with it is use "sm.density" to explore > the distribution, which found two large peaks. However, the peaks are made > up of Gaussians, and

Re: [R] repeatedly applying function with matrix-rows as argument

2003-06-30 Thread Prof Brian Ripley
apply(mymatrix, 1, function(x) myfun(x[1], x[2], x[3])) On Mon, 30 Jun 2003, Maarten Speekenbrink wrote: > Suppose I have a function which takes three arguments. I would like to > repeatedly apply the function, using a matrix N*3 in which each row > supplies the three argements for the function.

Re: [R] ./configure needs /sw/lib?

2003-06-30 Thread Prof Brian Ripley
On Mon, 30 Jun 2003, Tony Marlboro wrote: > > Thank you for your reply, Professor Ripley. > > > Yes, to use R CMD INSTALL directly: use R CMD INSTALL --help to find > > out how. You need to supply --with-netCDF=/sw, I believe (but then > > the headers may not be found). > > R CMD INSTALL --con

Re: [R] R CMD check

2003-06-30 Thread Prof Brian Ripley
You have an error in the R files of your package. Before you even do R CMD check, do try loading the package in R. On Mon, 30 Jun 2003, John Bjørnar Bremnes wrote: > when using R CMD check mypkg I get the error message > ... > * checking R files for library.dynam ... OK > * checking generic/meth

Re: [R] make 1.7.1 crashes in Mac OS X

2003-06-30 Thread Prof Brian Ripley
This is described in the MacOS X section of the R-admin manual. The file INSTALL does ask you to read it: please don't expect R-help readers to do so for you. On Mon, 30 Jun 2003, Michael Herron wrote: > I trying to compile R-.1.7.1 under OS X. > > Configure worked fine; then make crashes with:

Re: [R] date to age

2003-06-30 Thread Prof Brian Ripley
On Mon, 30 Jun 2003, William Briggs wrote: > I have files which have columns of data that look like this: > > DOB > 9/27/1964 > ... > > That is, dates in month/day/year format. When variables like DOB are read > in, they are converted to factors. Why? You can avoid that, of course. > So, I w

Re: [R] masked objects

2003-07-01 Thread Prof Brian Ripley
On Tue, 1 Jul 2003, ATHANASIA KAMARIOTIS wrote: > When opening the software R it appears this message: > > > > - > Attaching package 'methods': > > > The following object(s) are masked _by_ .GlobalEnv : > > new

Re: [R] namespaces and saved objects

2003-07-01 Thread Prof Brian Ripley
On Tue, 1 Jul 2003, Laurent Gautier wrote: > Hi, > > While saving (function 'save') a 'lmList' object (pack nlme), > R issues the warning messages: > 1: namespaces may not be available when loading > 2: names in persistent strings are currently ignored > 3: namespaces may not be available when

Re: [R] X-emacs and R for Windows

2003-07-01 Thread Prof Brian Ripley
On Tue, 1 Jul 2003, A.J. Rossini wrote: > Philippe Glaziou <[EMAIL PROTECTED]> writes: > > > Chriss, Neil <[EMAIL PROTECTED]> wrote: > >> Is there any way to make XEmacs (or any other editor) the default editor for > >> R so that when I type > >> > >> > sample <- function(x,y) { > >> + z<-x+y

Re: [R] Warning message in scatter.smooth (modreg)

2003-07-01 Thread Prof Brian Ripley
On Tue, 1 Jul 2003, Gavin Simpson wrote: > Dear list, > > In using the scatter.smooth() function (modreg) on a small data set (100 > obs) the following error was produced: > > > scatter.smooth(Na, S) > Warning message: > k-d tree limited by memory. ncmax= 200 > > I haven't used scatter.smooth

Re: [R] Memory problem:Failing to increase dynamic memory

2003-07-01 Thread Prof Brian Ripley
On Tue, 1 Jul 2003, Tapan Mehta wrote: > Hello, > > I am trying to use R1.7 on Linux. I having some > problem with memory. The task has to handle 100 files > of 10MB (each file is .CEL file) and is related to > microarrays. I am trying to run this task on a 2GB or > a 4 GB node of a Bewoulf Linux

Re: [R] Matrix to Vector

2003-07-01 Thread Prof Brian Ripley
On Tue, 1 Jul 2003, Arun Horne wrote: > I have a matrix of values that I want to convert to a vector, but remove > certain entries in the matrix first, i.e. in the vector I want every value > in the matrix that is greater than 1000. I would also like to find out how > many values are excluded in t

Re: [R] How long is a day?

2003-07-01 Thread Prof Brian Ripley
On Wed, 2 Jul 2003, Laimonis Kavalieris wrote: > Why is 19 March, 1947 a little longer than one day? A change from Daylight Savings Time in your time zone (unstated), according to your OS. Note, 0.04167 is exactly one hour, not `a little longer'. > x <- as.POSIXct("1947-04-16") > > julian(

Re: [R] Swedish characters in data frames

2003-07-01 Thread Prof Brian Ripley
On Wed, 2 Jul 2003, Neil White wrote: > I have some data that was prepared while I was Sweden. The columns > labelled using ä or å are fine, but for some reason ö is changed to a > period. The reason is in ?data.frame and ?read.table > so Blåbär stay the same, but Dödved becomes D.dved > > I'

Re: Unit of legend() coordinates (was: Re: [R] lines and legend)

2003-07-02 Thread Prof Brian Ripley
On Wed, 2 Jul 2003, Paul Lemmens wrote: > I would like to make a suggestion for an improvement (IMHO) of ?legend. > > Please add a remark that the x,y positioning coordinates are in the units > of the plot() itself, *not* in pixels or anything the like. That's true of all x-y coordinates of all

Re: [R] using [i] to plot & label all vector elements

2003-07-02 Thread Prof Brian Ripley
On Wed, 2 Jul 2003, Robert Schick wrote: > I'm using 1.6.2 on a win 2k box. (I know I'm due for an upgrade.) > > I've used brute force to label a series of vectors, and I'm wondering if > there's a better way to do this. I wanted to create a biplot of the 2nd > & 3rd PCs, and did this: You seem

Re: [R] fractional seconds from format.POSIXct

2003-07-02 Thread Prof Brian Ripley
On Wed, 2 Jul 2003, Dutky, Steve wrote: > Is there a format that yields fractional seconds from format.POSIXct and/or > related methods? No, by definition of the class "POSIXct". See my article on this in an early R Newsletter. > I'm attempting to use irts with millisecond events. -- Brian D

Re: [R] AID and CHAID: news about?

2003-07-03 Thread Prof Brian Ripley
AID and CHAID are (as I understand the terms, e.g. Kass, 1980, Applied Statistics)) more than 20 years old procedures with known pitfalls. They appeared to have been completely eclipsed by other techniques (e.g. rpart). On Thu, 3 Jul 2003, Daniele Medri wrote: > I just search if is there somethi

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