Re: [R] Help: lda predict

2005-08-26 Thread Prof Brian Ripley
methods. It does, in detail, in sections 2.4 and 2.5. Thanks you, Shengzhe On 8/26/05, Prof Brian Ripley [EMAIL PROTECTED] wrote: On Fri, 26 Aug 2005, Shengzhe Wu wrote: Thanks for your reply. Actually I called function as below. p1 = predict(object, newdata, dimen=1) p2 = predict(object

Re: [R] Shutting down a trellis plot (was R-help Digest, Vol 30, Issue 26)

2005-08-26 Thread Prof Brian Ripley
I suspect you have not print()-ed your graphics, see FAQ Q7.22. It is then possible to include dev.off() within the function. E.g. testit - function(fn = test.eps) { trellis.device(postscript, file=fn, onefile = FALSE, horizontal=FALSE) print(stripplot(voice.part ~ jitter(height), data =

Re: [R] passing arguments from nnet to optim

2005-08-26 Thread Prof Brian Ripley
On Fri, 26 Aug 2005, Tarca, Adi wrote: Hi everyone, According to R reference manual, the nnet function uses the BFGS method of optim to optimize the neural network parameters. What the help page says is ...: arguments passed to or from other methods. That means methods of nnet().

Re: [R] Matrix oriented computing

2005-08-26 Thread Prof Brian Ripley
Try profiling. Doing this many times to get an overview, e.g. for sapply with df=1:1000: % self% total self secondstotalsecondsname 98.26 6.78 98.26 6.78 FUN 0.58 0.04 0.58 0.04 unlist 0.29 0.02

Re: [R] problem with certain data sets when using randomForest

2005-08-26 Thread Prof Brian Ripley
Look at ?[.factor: finaldataset$Species - finaldataset$Species[,drop=TRUE] solves this. On Fri, 26 Aug 2005, Martin Lam wrote: Hi, Since I've had no replies on my previous post about my problem I am posting it again in the hope someone notice it. The problem is that the

Re: [R] parts of data frames: subset vs. [-c()]

2005-08-26 Thread Prof Brian Ripley
Are there NAs in the variable? SYNTAX==Ditrans and SYNTAX!=Ditrans are not mutually exclusive. On Fri, 26 Aug 2005, Stefan Th. Gries wrote: Dear all I have a problem with splitting up a data frame called ReVerb: » str(ReVerb) `data.frame': 92713 obs. of 16 variables: $ CHILD: Factor

Re: [R] class 'named'

2005-08-26 Thread Prof Brian Ripley
That is one of the S4 vs R differences. See the complements. On Fri, 26 Aug 2005, Phineas Campbell wrote: I'm working through the examples in Venables and Ripley in the 'New-style Classes' chapter. On a call to representation, in the lda example, it is unable to find the class named. Is

Re: [R] Help in Compliling user -defined functions in Rpart

2005-08-26 Thread Prof Brian Ripley
On Fri, 26 Aug 2005, Luwis Tapiwa Diya wrote: I have been trying to write my own user defined function in Rpart.I imitated the anova splitting rule which is given as an example.In the work I am doing ,I am calculating the concentration index(ci) ,which is in between -1 and +1.So my deviance

Re: [R] Creating factors from continuous variables

2005-08-26 Thread Prof Brian Ripley
?cut This is in `An Introduction to R', the manual which ships with R and basic reading. On Fri, 26 Aug 2005, David James wrote: What is the quickest way to create many categorical variables (factors) from continuous variables? This is the approach that I have used: # create sample data

Re: [R] zoo, zooreg, ISOdatetime

2005-08-27 Thread Prof Brian Ripley
On Fri, 26 Aug 2005, David James wrote: I create a zooreg object that runs from Jan-1-2002 0:00 to Jun-1-2005 0:00... regts.start = ISOdatetime(2002, 1, 1, hour=0, min=0, sec=0, tz=) regts.end = ISOdatetime(2005, 6, 1, hour=0, min=0, sec=0, tz=) regts.zoo - zooreg( NA, regts.start,

Re: [R] survival parametric question

2005-08-27 Thread Prof Brian Ripley
The survival package is a recommended package in R and contains survreg() which uses the AFT definitions for Weibull survival. This is well documented, and MASS (the book) has comparisons of PH and AFT parametrization for a Weibull example. I think you mean Frank Harrell's `Design' package.

Re: [R] Logistic regression with constrained parameter ??

2005-08-28 Thread Prof Brian Ripley
On Sat, 27 Aug 2005, Benn Fine wrote: I want to fit a logistic regression model under a specified null hypothesis, say Ho:Beta_k=1 Is there a way to constrain a parameter in glm.fit ?? You should be calling glm(). Then you can use offset() in your formula. (You can also use it as an

Re: [R] ttda on R 2.1.1: error

2005-08-29 Thread Prof Brian Ripley
Looks like the package is not written for current R. native.enc was used prior to R 2.1.0. Please take the advice of the posting guide and ask the maintainer. BTW, there is no package ttda on CRAN, and you have not told us where you found it or what version you used. If you mean

Re: [R] prediction intervals

2005-08-29 Thread Prof Brian Ripley
What is an `mle object'? If it is something produced by mle() in stats4, in what sense is `prediction' relevant? On Mon, 29 Aug 2005 [EMAIL PROTECTED] wrote: Hello. I am looking for a reference to compute prediction intervals from a mle object. Thank you. Tony. [[alternative

Re: [R] memory

2005-08-30 Thread Prof Brian Ripley
On Mon, 29 Aug 2005, Ferran Carrascosa wrote: Hi, I have a matrix with 700.000 x 10.000 cells with floating point data. I would like to work with the entire table but I have a lot of memory problems. I have read the ?memory I work with Win 2000 with R2.1.0 The only solution that I have

Re: [R] question about custom contrasts in ANOVA

2005-08-30 Thread Prof Brian Ripley
I was puzzled as to what the question actually was. If you set options(contrasts=c(contr.sum, contr.poly)) the interaction contrasts are precisely those created manually. Nothing fancier is required. But I am not sure what you want to do with them once you have them. On Tue, 30 Aug 2005,

Re: [R] Font Encodings --- some work, some don't

2005-08-31 Thread Prof Brian Ripley
Are you talking about the pdf() device? You never say so, but you do say The font encoding error does not appear in the postscript device driver, and it works fine. This report is far too vague: what did you do and what reported the file was corrupt? (Note that Acrobat reader is famous for

Re: [R] R binaries

2005-08-31 Thread Prof Brian Ripley
On Wed, 31 Aug 2005, Berwin A Turlach wrote: G'day Nam-Ky, NKN == Nam-Ky Nguyen [EMAIL PROTECTED] writes: NKN I intend to burn some R CDs to colleagues in Vietnam. I want NKN to put all binary files for base as well as contributed NKN packages (for both Windows and Linux). I

Re: [R] Installation help

2005-08-31 Thread Prof Brian Ripley
On Tue, 30 Aug 2005, b wrote: Solved the problem with the 'R' installation! One of the processes that loads up on startup was conflicting with it. It was a program called Gameutil.exe which seems to enhance games by altering display modes. For some reason a conflict existed. Anyhow I dont

Re: [R] tcl/tk return problem

2005-08-31 Thread Prof Brian Ripley
On Wed, 31 Aug 2005, Philippe Grosjean wrote: deggle wrote: Hello, I'm very new in working with tcl/tk in R and have a problem which will probably sound silly to most of you. Here is the code I have problems with: readcelfiles - function() { require(tcltk) tt -

Re: [R] tcl/tk return problem

2005-08-31 Thread Prof Brian Ripley
On Wed, 31 Aug 2005, Duncan Murdoch wrote: Prof Brian Ripley wrote: On Wed, 31 Aug 2005, Philippe Grosjean wrote: deggle wrote: Hello, I'm very new in working with tcl/tk in R and have a problem which will probably sound silly to most of you. Here is the code I have problems

Re: [R] Splitting the string at the last sub-string

2005-09-15 Thread Prof Brian Ripley
On Thu, 15 Sep 2005, Barry Rowlingson wrote: Prof Brian Ripley wrote: substring(str, c(1, 26), c(25,length(str))) nchar(str) surely? Yes, or anything larger: I actually tested 1. regexps can be rather slow though. Here's two functions: But that's not the way to do this repeatedly

Re: [R] make check FAILS - Error code 1...comparing d-p-q-r-tests.Rout

2005-09-15 Thread Prof Brian Ripley
The problem is a known one and solved in R 2.1.1-patched. On your OS the result is not as accurate as most, but the tolerance set was too tight so the test failure is not something to worry about. Please install R-patched instead, as it has many bug fixes in place. On Thu, 15 Sep 2005,

Re: [R] Coefficients from LM

2005-09-16 Thread Prof Brian Ripley
On Fri, 16 Sep 2005, Dieter Menne wrote: Charles Annis, P.E. Charles.Annis at StatisticalEngineering.com writes: Here's an example with a glm; lm() works the same way but has fewer internal objects. mod3 - glm(tree ~ altitude, family = binomial) You can use names() to find out what's

Re: [R] About princomp

2005-09-16 Thread Prof Brian Ripley
You send HTML mail (please see the posting guide and do as we ask) and it seems to have not produced a working text version. My guess is that you have a changed version of USArrests around. In any case, there is no problem in R 2.1.1 (sic), so this is local to you. On Fri, 16 Sep 2005,

Re: [R] Integrate functions with loops

2005-09-16 Thread Prof Brian Ripley
Try to give a vector result with one element for each of element of x, e.g. integrandtotest - function(x) colSums(outer(t, x, ^)) works, although in fact this integration can be done analytically (it is a sum of exponentials). On Fri, 16 Sep 2005, A.Brennan wrote: Hi i am having a

Re: [R] How to update R from ver 1.9.1 to 2

2005-09-19 Thread Prof Brian Ripley
See the rw-FAQ ..., Q2.8 in the latest version. On Sat, 17 Sep 2005, Srinivas Iyyer wrote: Dear group, apologies if this is a stupid question. I searched CRAN sites. I am afraid I missed it. Can any one help me if I can update my windows version of 1.9.1 to 2 or higer. -- Brian D.

Re: [R] library (tree)- which samples belong to each terminal branch?

2005-09-19 Thread Prof Brian Ripley
On Mon, 19 Sep 2005, Chris Buddenhagen wrote: What I would like to know is if there is a way to know which specific samples fall within a terminal branch created by tree? The classification summarizes data, but I want to know which specific samples fall within each classification (branch in

Re: [R] Interpretation of csplit from rpart.object

2005-09-21 Thread Prof Brian Ripley
Your message *was* received, and you can check the archives to see it at https://stat.ethz.ch/pipermail/r-help/2005-September/077889.html You need to read the code to answer the question for yourself. There is lots of code interpreting csplit in the rpart package. These lines might be a

Re: [R] user ignored an informative warning message (was ts.intersect bug?)

2005-09-21 Thread Prof Brian Ripley
On Wed, 21 Sep 2005, Antonio, Fabio Di Narzo claimed: This code gives an error: a - ts(1:10, start=0, freq=10) b - ts(1:10, start=1, freq=10) ts.intersect(a,b) In reality, it gives a helpful warning: ts.intersect(a,b) NULL Warning message: non-intersecting series in: .cbind.ts(list(...),

Re: [R] size of subplots with par() / layout()

2005-09-21 Thread Prof Brian Ripley
On Wed, 21 Sep 2005, Robin Hankin wrote: If I do this: par(mfrow=c(2,2)) persp(matrix(1:4,6,6),box=F,phi=33,theta=33) persp(matrix(1:4,6,6),box=F,phi=33,theta=33) persp(matrix(1:4,6,6),box=F,phi=33,theta=33) persp(matrix(1:4,6,6),box=F,phi=33,theta=33) ( or indeed

Re: [R] cdecl and stdcall

2005-09-21 Thread Prof Brian Ripley
This really is the wrong list. Please do read the posting guide as we do ask (before posting). In so far as this is relevant to R at all, it is within the description of the R-devel list. On Wed, 21 Sep 2005, Bogner, Konrad (LfU) wrote: I'm trying to load a dynamic link library and it seems

Re: [R] Fwd: ts.intersect bug?

2005-09-21 Thread Prof Brian Ripley
The problem was that a tolerance in window.default was abs(start)*ts.eps, which is no tolerance at all if start = 0. So the special feature here was that the intersection of the series started at zero. I had already tested and committed a fix. window.default does use a tolerance for en, just

Re: [R] source(file) = file becomes readonly

2005-09-21 Thread Prof Brian Ripley
This was not really accurate. Being `read-only' is a Windows phenomenon on open files (even files open for reading). In 2.0.1 and 2.1.1 the file was open whilst it was parsed and so `read-only' during that time. In 2.0.1 it was then closed, but in 2.1.1 it was (potentially) opened several

Re: [R] controlling usage of digits scientific notation in R plots; postscript margins

2005-09-22 Thread Prof Brian Ripley
Concise answers to verbose questions: 1) Use options(scipen) (and probably change the margin sizes). Or something like options(scipen=10) par(mar=c(5,8,4,2)+0.1) plot(x, y, axes=FALSE) axis(2, las=2) axis(1, labels=FALSE) axis(1, at = c(-2e8, 2e8), labels = expression(-2 %*% 10^8, 2 %*% 10^8))

Re: [R] accessing source code in R packages

2005-09-22 Thread Prof Brian Ripley
The original reply was deliberately (I guess) vague. (I've removed the history, as attributions had already been removed, in violation of copyright law. If you cite someone, you MUST credit the author.) Sometimes a little knowledge is a dangerous thing, and we have had a number of partially

Re: [R] Rf_initEmbeddedR in Windows

2005-09-23 Thread Prof Brian Ripley
How to generate an import library for R.dll for Visual Studio is discussed in README.packages: please do read it. Do note the recommendation not to install R in a path with spaces - your tool is one of many that does not work with spaces. But you don't want a Borland import library (any more

Re: [R] as.character on OS X tiger

2005-09-23 Thread Prof Brian Ripley
This is just a question of rounding error. You are computing 11/999, which is not exactly representable. Most computers use IEC60559 arithmetic, which gives a precision of .Machine$double.eps ~ 2e-16. So you can expect rounding error around 11/999 * .Machine$double.eps ~ 2e-18. However,

Re: [R] Strange behaviour of as.Date function

2005-09-23 Thread Prof Brian Ripley
datai[1,1] appears to be a 'POSIXct' object. Which date that is depends on the locale, and as.Date uses UTC (see ?as.Date). For me: d - as.POSIXct(2005-07-01) d [1] 2005-07-01 BST format(d, tz=GMT) [1] 2005-06-30 23:00:00 as.Date(d) [1] 2005-06-30 Use as.Date(d + 23.99*3600) to avoid

Re: [R] Problem with read.spss() and as.data.frame(), or: alternative to subset()?

2005-09-24 Thread Prof Brian Ripley
On Fri, 23 Sep 2005, Thomas Lumley wrote: On Wed, 21 Sep 2005, Martin Maechler wrote: Dirk == Dirk Enzmann [EMAIL PROTECTED] on Wed, 21 Sep 2005 13:18:32 +0200 writes: Dirk The selection problem can be solved by Dirk dr2000=read.spss('myfile') Dirk

Re: [R] summary nls output

2005-09-25 Thread Prof Brian Ripley
On Sun, 25 Sep 2005, Chun-Ying Lee wrote: Dear R user: I bulid a package, and in the package I use the function nls to solve some questions. If I have two sets of data, and I want to summary these two data's nls output, I write the command in the package source code like: { ..

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

2005-09-25 Thread Prof Brian Ripley
On Sun, 25 Sep 2005, Andrew Robinson wrote: Hi R-helpers, I am trying to build a package under FreeBSD 5.4-RELEASE #0 using R Version 2.1.1. I have constructed a package using package.skeleton(), when I try $ R CMD build foo * checking for file 'foo/DESCRIPTION' ... OK * preparing

Re: [R] nls and na/Nan/Inf error

2005-09-26 Thread Prof Brian Ripley
This works if you omit the deriv() step. Use R's options(error=dump.frames) and debugger(). This gives Browse[1] rhs [1] 0.433 0.4272571 0.3994105 0.3594037 0.3270730 0.3104752 0.3000927 [8] 0.2928445 0.2874249 0.2831787 attr(,gradient) VrVm alphalamda

Re: [R] k-d tree for loess

2005-09-26 Thread Prof Brian Ripley
First a warning: loess in R is only loosely related to loess in S, being derived from a C implementation (by the same authors). In R I don't think you can do this. Those details are never exposed, and are hidden in an undocumented C/Fortran workspace. On Mon, 26 Sep 2005, Carlisle Thacker

Re: [R] anova on binomial LMER objects

2005-09-26 Thread Prof Brian Ripley
On Mon, 26 Sep 2005, Martin Henry H. Stevens wrote: Hello all, 1. Does Matrix 0.98-7 fix any of this? 2. Assuming no, how does one acquire Matrix 0.95-13? It is in the Archive on CRAN, e.g. http://cran.r-project.org/src/contrib/Archive/M/Matrix_0.95-13.tar.gz Cheers, and thank you kindly

Re: [R] reading SAS data files

2005-09-26 Thread Prof Brian Ripley
On Mon, 26 Sep 2005, Dean Sonneborn wrote: I am attempting to read in a SAS 9.1 data file. After starting R I change to the directory containing the sas data file and use the dir command to confirm that it is there. Then I run the following R-code: library(foreign) sashome - /Program

Re: [R] Make check fails on d-p-q-r-tests.R...

2005-09-27 Thread Prof Brian Ripley
You will have to show us the error! It will be shown in d-p-q-r-tests.Rout.fail (unless this was a segfault or similar). It is not OK to skip the test, but note that this test is random and does fail about 1 in 50 times, so you could just try rerunning it. On Mon, 26 Sep 2005, Jeff Ross

Re: [R] Error in make check-all

2005-09-27 Thread Prof Brian Ripley
This is what happens if you don't have a usable X11 display. Did you perhaps use a root account on a console owned by a normal user? On Tue, 27 Sep 2005, Fernando Mayer wrote: Dear R-users, i'm a very newbie in linux, but decided to build R from source. Following the R Installation and

Re: [R] quick points question

2005-09-27 Thread Prof Brian Ripley
On Tue, 27 Sep 2005, Richard Nixon wrote: Hi Just one of those niggles... I've just been trying to plot a filled circle. I thought that this would do it plot(1,1,type=n) points(1,1,pch=1,bg=blue,cex=5) #bg: background (fill) color for open plot symbols But I need to do this instead

Re: [R] Error in make check-all

2005-09-27 Thread Prof Brian Ripley
you need the console running the user that is going to be running R. (I see this when running tests remotely, as then I have no X11 display set to save bandwidth. I just ignore it.) Thanks, Fernando Mayer. Prof Brian Ripley escreveu: This is what happens if you don't have a usable X11

Re: [R] Make check fails on d-p-q-r-tests.R...

2005-09-27 Thread Prof Brian Ripley
Prof Brian Ripley [EMAIL PROTECTED] wrote: You will have to show us the error! It will be shown in d-p-q-r-tests.Rout.fail (unless this was a segfault or similar). It is not OK to skip the test, but note that this test is random and does fail about 1 in 50 times, so you could just try

Re: [R] Question on lm(): When does R-squared come out as NA?

2005-09-28 Thread Prof Brian Ripley
I've not seen a reply to this, nor ever seen it. Please make a reproducible example available (do see the posting guide). On Sun, 25 Sep 2005, Ajay Narottam Shah wrote: I have a situation with a large dataset (3000+ observations), where I'm doing lags as regressors, where I get: Call:

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

2005-09-28 Thread Prof Brian Ripley
On Tue, 27 Sep 2005, Eric van Gyzen wrote: On Sun, Sep 25, 2005 at 12:07:02PM +0100, Prof Brian Ripley wrote: On Sun, 25 Sep 2005, Andrew Robinson wrote: Hi R-helpers, I am trying to build a package under FreeBSD 5.4-RELEASE #0 using R Version 2.1.1. I have constructed a package using

Re: [R] anova on binomial LMER objects

2005-09-28 Thread Prof Brian Ripley
On Wed, 28 Sep 2005, Robert Bagchi wrote: Hi Patrick thanks for your advice. I have now tried glmmPQL, and it worked fine - I'm getting consistent results between plots and models fitted by glmmPQL. Plus it allows predict() and resid() which is another advantage over lmer at present. quick

Re: [R] Errors in odbcConnectExcel()

2005-09-28 Thread Prof Brian Ripley
This is entirely a function of the database format the ODBC Excel driver expects. It is *not* an error in odbcConnectExcel as your subject accuses. To use a spreadsheet as a database you need to have column names. On Thu, 29 Sep 2005, [ISO-2022-JP] ÅÄÃæ Áï wrote: Dear R-help I would like

Re: [R] Change console language ?

2005-09-28 Thread Prof Brian Ripley
On Wed, 28 Sep 2005, Sérgio Nunes wrote: Hi, I'm pretty new to R, I've just installed version 2.1.1 on Windows. I have a simple (it seems) doubt - how do I change the Console default Language ? It's set to Portuguese but I would like to view it in English. Do read the administration manual

Re: [R] how to fix the level-1 variances in lme()?

2005-09-29 Thread Prof Brian Ripley
On Thu, 29 Sep 2005, Mike Cheung wrote: Dear all, Edmond Ng (http://multilevel.ioe.ac.uk/softrev/reviewsplus.pdf) provides an example to fit the mixed effects meta-analysis in Splus 6.2. The syntax is: lme(fixed=d~wks, data=meta, random=~1|study, weights=varFixed(~Vofd),

Re: [R] is it possible to form matrix of matrices...and multiple arrays

2005-09-29 Thread Prof Brian Ripley
On Thu, 29 Sep 2005 [EMAIL PROTECTED] wrote: booop booop a écrit : 1...Kindly tell me is it possible to form a matrix which contains a no of matrices.. for eg.. if a,b,c,d are matrices and e is a matrix which contains a,b,c,d as rows and columns.. I don't think you can use

Re: [R] priceIts

2005-09-29 Thread Prof Brian Ripley
On Thu, 29 Sep 2005, Mulholland, Tom wrote: Well I downloaded the data using the link in your message which suggests that the code is right. I don't have its loaded (I assume it's from the irregular time series package) The package is called its, as the message correctly said. so I can't

Re: [R] R-2.1.1 on AIX 5.*

2005-10-01 Thread Prof Brian Ripley
Note your first error message is connections.c, line 2580.18: 1506-052 (S) Duplicate case label for value 4. Labels must be unique. and the relevant line is #if SIZEOF_LONG == 8 case sizeof(long): so you have somehow got options to give you 4-byete longs although configure

Re: [R] update MASS

2005-10-01 Thread Prof Brian Ripley
On Sat, 1 Oct 2005, Ramón Casero Cañas wrote: I'd like to update MASS from version 7.2-11 to version 7-2.19. I am running R 2.0.1 on ubuntu, that installs MASS with the package r-cran-vr. First, MASS is part of VR, so it is VR you update. Second, please do *READ THE HELP PAGE*. The first

Re: [R] Placing axes label strings closer to the graph?

2005-10-01 Thread Prof Brian Ripley
Use mtext() to put labels in the margin line you want. title() writes xlab and ylab in line 3 AFAIR. So you could use something like par(mai=c(0.6,0.6,0.2,0.2)) plot(D$dates, D$inflation, type=l, lwd=2, col=cadetblue4, cex.axis=0.6) mtext(X label string, 1, line=2, cex = 0.6) mtext(Y

Re: [R] Regression slope confidence interval

2005-10-01 Thread Prof Brian Ripley
On Thu, 29 Sep 2005, Christian Hennig wrote: ?confint Thank you to all of you. As far as I see this is not mentioned on the lm help page (though I presumably don't have the recent version), which I would suggest... and I would suggest that you study a good book on the subject. (confint

Re: [R] ACCESS R and dates

2005-10-01 Thread Prof Brian Ripley
On Fri, 30 Sep 2005, sloan jones wrote: I have used the RODBC package to read in data I have stored in an Access file. When I am using data from files other than ACCESS I have no problem using the survival package to work with dates; however, with the ACCESS data the dates are reading-in

Re: [R] dec2bin?

2005-10-01 Thread Prof Brian Ripley
On Sat, 1 Oct 2005, Haiyong Xu wrote: I just want to ask if there is any function that can convert decimal number to binary vector. What do you mean by a binary vector? Function intToBits might be what you are looking for, e.g. x - as.integer(1234) y - intToBits(x) options(width=50) y

Re: [R] Y-Axis Label Overwritten By TickMark Values

2005-10-01 Thread Prof Brian Ripley
On Sat, 1 Oct 2005, Paul Roebuck wrote: On Sat, 1 Oct 2005, Marc Schwartz wrote: On Sat, 2005-10-01 at 06:07 -0500, Paul Roebuck wrote: Can someone tell me how to fix the left margin of plot region such that the tick values don't overwrite the x-axis label? I haven't been able to set the

Re: [R] User error (was arima.sim bug?)

2005-10-02 Thread Prof Brian Ripley
On Sun, 2 Oct 2005, Kemp S E (Comp) wrote: Hi, I am using the arima.sim function to generate some AR time series. However, the function does not seem to produce exactly the same time series when I specify the innov parameter. For example r - rnorm(300) x - arima.sim(300,

Re: [R] modeling language for optimization problems

2005-10-02 Thread Prof Brian Ripley
On Sun, 2 Oct 2005, Paolo Cavatore wrote: Does anyone know whether R has its own modeling language for optimization problems (like SIMPLE in NuOPT for S-plus)? No. Note that SIMPLE is the language of NUOPT, not of S-PLUS. There is an (extra-cost) interface module S+NUOPT, but it is an

Re: [R] Memory management on Windows (was Size of jpegs/pngs)

2005-10-02 Thread Prof Brian Ripley
I think this an issue about the amount of graphics memory. You are asking for an image of about 17*2*3 = 102Mb, and you need more than that. From the help page: Windows imposes limits on the size of bitmaps: these are not documented in the SDK and may depend on the version of

Re: [R] Local install of a contributed package under Linux (was: R-code for binormla distribution)

2005-10-03 Thread Prof Brian Ripley
On Sun, 2 Oct 2005, Jonathan Baron wrote: Try R CMD INSTALL -l lib pkgs The help file is in the utils package. I'm sure this is documented in the manual too. But you still need to ensure that library() can find the packages. I think the simplest way is that we have set up by default for

Re: [R] Problem building/checking library that requires input from user

2005-10-03 Thread Prof Brian Ripley
Either use \dontrun in a \examples{} section of a help page (see Writing R Extensions), or something like num_reps - if(interactive()) readline(How many reps do you have... ) else 500 On Mon, 3 Oct 2005, Ken Termiso wrote: Hi all, I've got a package i've written that i am trying to check,

Re: [R] modeling language for optimization problems

2005-10-03 Thread Prof Brian Ripley
] On Behalf Of Prof Brian Ripley Sent: Sunday, October 02, 2005 10:38 AM To: Paolo Cavatore Cc: r-help@stat.math.ethz.ch Subject: Re: [R] modeling language for optimization problems On Sun, 2 Oct 2005, Paolo Cavatore wrote: Does anyone know whether R has its own modeling language for optimization

Re: [R] dec2bin?

2005-10-03 Thread Prof Brian Ripley
On Mon, 3 Oct 2005, Tuszynski, Jaroslaw W. wrote: It is unclear what you are trying to do, but check bin2raw in caTools package: print(x - (1:5)*pi) [1] 3.141593 6.283185 9.424778 12.566371 [5] 15.707963 print(y - bin2raw(x)) [1] 18 2d 44 54 fb 21 09 40 18 2d 44 54 fb 21 19

Re: [R] gnomeGUI installation

2005-10-03 Thread Prof Brian Ripley
On Mon, 3 Oct 2005, Daniel Pick wrote: I have successfully downloaded the sources and built R as a shared library on a Red Hat Enterprise Level 3 box. I am now trying to build the GNOME GUI, but configure is barfing on glade. According to the system logs, the RPM for

Re: [R] modeling language for optimization problems

2005-10-03 Thread Prof Brian Ripley
help.search(convex programming) No help files found with alias or concept or title matching 'convex programming' using fuzzy matching. ... Thanks, Roger On 10/3/05, Prof Brian Ripley [EMAIL PROTECTED] wrote: On Mon, 3 Oct 2005, Huntsinger, Reid wrote: Have you looked at the R interface

Re: [R] Getting eps into Word documents.

2005-10-04 Thread Prof Brian Ripley
On Mon, 3 Oct 2005, Marc Schwartz (via MN) wrote: On Mon, 2005-10-03 at 22:00 +0100, Ted Harding wrote: Rolf ( Marc) On 03-Oct-05 Marc Schwartz (via MN) wrote: On Mon, 2005-10-03 at 16:31 -0300, Rolf Turner wrote: A student in one of my courses has asked me about getting R graphics output

Re: [R] missing handling

2005-10-04 Thread Prof Brian Ripley
On Tue, 4 Oct 2005, Weiwei Shi wrote: Hi, Jim: I tried your code and get the following error: trn1-read.table('trn1.svm', header=F, na.string='.', sep='|') Med-apply(trn1, 2, median, na.rm=T) Ind-which(is.na(trn1), arr.ind=T) trn1[Ind]-Med[Ind[,'col']] Error in [-.data.frame(`*tmp*`, Ind,

Re: [R] prediction with nnet function for regression in R

2005-10-05 Thread Prof Brian Ripley
On Wed, 5 Oct 2005 [EMAIL PROTECTED] wrote: Hi, Good Day Everybody. Can anybody help me I am working with R ...and i am using nnet() function for regression, but after trained model, i m getting all predictions as 0.99 or 1., why is so, i dont understand? As a guess, you

Re: [R] Fisher's discriminant functions

2005-10-05 Thread Prof Brian Ripley
On Wed, 5 Oct 2005, Leonardo Lami wrote: Hi, I have the same problem. I found a solution but I think there is something of more simple and correct. I take a group and I put 0 the values of the other group, after I use the function glm to abtain the Fisher's discriminant function for this

Re: [R] eliminate t() and %*% using crossprod() and solve(A,b)

2005-10-05 Thread Prof Brian Ripley
On Wed, 5 Oct 2005, Robin Hankin wrote: I have a square matrix Ainv of size N-by-N where N ~ 1000 I have a rectangular matrix H of size N by n where n ~ 4. I have a vector d of length N. I need X = solve(t(H) %*% Ainv %*% H) %*% t(H) %*% Ainv %*% d and H %*% X. It is possible to

Re: [R] transparent surface in rgl

2005-10-05 Thread Prof Brian Ripley
On Wed, 5 Oct 2005, Duncan Murdoch wrote: On 10/5/2005 9:31 AM, Prof. Paul R. Fisher wrote: Hi all I am a complete newbie to this list (just subscribed) and a newcomer to R (an S user from olden times). I have been using scatter3d to create a 3d scatter plot with surface. The graphic is

Re: [R] --gui none problem

2005-10-06 Thread Prof Brian Ripley
On Thu, 6 Oct 2005, Luc Vereecken wrote: Version 2.x no longer has a -gui=none option. Running BATCH jobs supposedly disables a gui though. There is no --gui=none option in 2.1.x (but there is in 2.0.x). The gui is not disabled in 2.1.x either according to the help page nor in practice

Re: [R] multiple line plots

2005-10-06 Thread Prof Brian Ripley
# Not sure why this did not work #sapply(tw, setlen, len) It probably did, but you discarded the result. Try tw - sapply(tw, setlen, len) On Thu, 6 Oct 2005, sosman wrote: Marc Schwartz wrote: On Wed, 2005-10-05 at 22:19 +1000, sosman wrote: I have some data in a CSV file:

Re: [R] isdst

2005-10-06 Thread Prof Brian Ripley
On Thu, 6 Oct 2005, Omar Lakkis wrote: Can someone, please, explain the difference is results below (notice the isdst value) Timezones. (That is carefully documented on the help page for as.POSIXlt.) unlist(as.POSIXlt('2005-7-1')) sec min hour mday mon year wday yday isdst 0

Re: [R] blank graphics window while R is working

2005-10-06 Thread Prof Brian Ripley
Your R is an obselete version. There are some changes in 2.2.0 which affect the hinting given to X11 WMs, so please try the current version of R. (I know it had only recently been released, but the posting guide refers you to the development versions.) I don't know if this will solve it (I

Re: [R] R-2.2.0 compilation problem [on MacOS X]

2005-10-07 Thread Prof Brian Ripley
On Fri, 7 Oct 2005, Robin Hankin wrote: Hi I tried to compile R-2.2.0 just now. On MacOS X. configure worked fine, but compilation stopped with Please ask MacOS-specific questions on the r-sig-mac list. (Missing RestFP, SaveFP.) I believe this is well-known and indicates an incompatible

Re: [R] builiding R from sources

2005-10-07 Thread Prof Brian Ripley
Please tell us which version of R. If this is 2.2.0 (I think it is) you need to review the tools you are using, as mingw-runtime-3.8 is required. (It says so explicitly in that manual, and that is the current version of that package.) For R-patched you will also need a recent enough version

Re: [R] missing values in step procedure

2005-10-07 Thread Prof Brian Ripley
On Fri, 7 Oct 2005, Andreas Cordes wrote: I have the problem that for the step procedure stops due to missing values. There are no options in Step or stepAIC to handle missing values. Is there any way to run stepwise modelselection in R in an automated way in this case? Try the hint it gives

Re: [R] problems with loess

2005-10-07 Thread Prof Brian Ripley
On Fri, 7 Oct 2005, Dimitri Giunchi wrote: Hi all, I was unable to obtained a smoothed line using the loess function.I used the following code reported in the examples of R documentation: cars.lo - loess(dist ~ speed, cars) Then I tried to plot both the data and the smoothed line

Re: [R] sscanf equivalent

2005-10-07 Thread Prof Brian Ripley
Why not use a text connection? On Fri, 7 Oct 2005, Paul Roebuck wrote: I have a data file from which I need to read portions of data but data location/quantity can change from file to file. I wrote some code and have a working solution but it seems wasteful to have to do it this way. Here's

Re: [R] Applying a function to each element of an array

2005-10-07 Thread Prof Brian Ripley
If A is the matrix the answer is 2*(1 - 2^(-A)), which took about 10secs for an example of your size. From \sum_{i=1}^n x^{1-i} = (1-x^{-n})/(1-x), E OE. On Fri, 7 Oct 2005, Tim Smith wrote: I have a 7000x7000 matrix, and each element is an integer. For each element, I want to apply the

Re: [R] Applying a function to each element of an array

2005-10-07 Thread Prof Brian Ripley
On Fri, 7 Oct 2005, Prof Brian Ripley wrote: If A is the matrix the answer is 2*(1 - 2^(-A)), which took about 10secs for an example of your size. From \sum_{i=1}^n x^{1-i} = (1-x^{-n})/(1-x), E OE. NB: I have assumed n = 1 here, as people nornally do when using 1:count. On Fri, 7 Oct

Re: [R] Applying a function to each element of an array

2005-10-07 Thread Prof Brian Ripley
On Fri, 7 Oct 2005, Patrick Burns wrote: If there weren't an analytic solution to your problem, then you could build a vector of the answers from 1 to the maximum in the matrix. Call that 'wtvec'. Then: ans - array(NA, dim(A), dimnames(A)) ans[] - wtvec[as.vector(A)] should get you what

Re: [R] how to control ticks in plots with yasp or xasp

2005-10-08 Thread Prof Brian Ripley
On Sat, 8 Oct 2005, Marc Schwartz wrote: On Sat, 2005-10-08 at 09:28 -0400, Denis Chabot wrote: Hi, A few times I tried to control the number and position of tick marks in plots with the yasp or xasp parameters. For example, a y axis was drawn by default with tick marks at 0, 20, 40, 80 and

Re: [R] add leading 0s to %d from png() {was Automatic creation of file names}

2005-10-08 Thread Prof Brian Ripley
The example on the png help page, Rplot%03d.png, _is_ what you want. (More details are on ?postscript.) formatC() is an S/R peculiarity: sprintf() is the cross-language way to do this sort of thing. On Sat, 8 Oct 2005, bogdan romocea wrote: Dear useRs, Is there a way to 'properly' format %d

Re: [R] how to control ticks in plots with yasp or xasp

2005-10-08 Thread Prof Brian Ripley
On Sat, 8 Oct 2005, Denis Chabot wrote: Hi, sorry about the bad syntax, though the right syntax would not have worked either, according to your tests (Mark, Brian, Peter). It DOES work according to my tests! (Do give us the credit for testing our advice: we would appreciate your showing

Re: [R] Rpart -- using predict() when missing data is present?

2005-10-08 Thread Prof Brian Ripley
On Sat, 8 Oct 2005, Ajay Narottam Shah wrote: I am doing library(rpart) m - rpart(y ~ x, D[insample,]) D[outsample,] y x 8 0.78391922 0.579025591 9 0.06629211 NA 10 NA 0.001593063 p - predict(m, newdata=D[9,]) Error in model.frame(formula,

Re: [R] lm.ridge

2005-10-09 Thread Prof Brian Ripley
On Sat, 8 Oct 2005, Erin Hodgess wrote: Dear R People: I have a question about the lm.ridge function, please. In the example, there is one set of output values in the select function but another in the comment section. Am I missing something please? The values in the examples were

Re: [R] sqlFetch on MySQL-DB

2005-10-10 Thread Prof Brian Ripley
On Mon, 10 Oct 2005, Bernd Weiss wrote: I successfully set up a local MySQL-database. Connecting via RODBC is not problem, the same in fetching 3 of 4 tables. But trying to connect to table 4 fails. author-sqlFetch(test,author) Error in fromchar(unclass(x)) : character string is not in a

Re: [R] Show Progress in loop

2005-10-10 Thread Prof Brian Ripley
On Mon, 10 Oct 2005, Philippe Grosjean wrote: Hello, You must explicitly use print(), show() on an object -here, use print(i)- in a loop or alternatively, use cat() to display string like: cat(loop, i, \n) With RGui under Windows, there is another subtility: if you have turn on 'Misc -

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