Re: [R] Sub with and without perl=TRUE

2005-02-14 Thread Prof Brian Ripley
This was a bug which is already fixed in R-2.0.1 patched. From the NEWS file o gsub(perl=TRUE) returned a string which printed with trailing garbage if there was a match at the beginning whose replacement was shorter. (PR#7479) The crucial comment is that it only happened

[R] gam(mgcv) starting values

2005-02-14 Thread Bjrn Stollenwerk
Hi all! Ive got some problems with the function gam (library mgcv). For some models I get the error message : Error: no valid set of coefficients has been found:please supply starting values In addition: Warning message: NaNs produced in: log(x) This is a shortened code I used: gam(y ~ M1 + M3

Re: [R] how can i make my program faster

2005-02-14 Thread Adaikalavan Ramasamy
On Sun, 2005-02-13 at 20:53 -0800, Cuichang Zhao wrote: Hello, right now, i have a program to collect data into a table. right now, my table is table1 - data.frame(trial = NA, x = NA, y = NA) One often uses the term 'table' as in tabulate or cross-tabulate discrete values. for each time

Re: [R] How to solve error : cannot allocate vector of size 1208235 Kb

2005-02-14 Thread Kum-Hoe Hwang
I appreciate all who helped all. My problem was not enough memory definietely. First I tried smaple function that solved my problem nicely. Second, I tried a method=sparse option in lagsarlm model. Comparing results led me to conclude the first trial that estimates better results in the lagsarlm

[R] Re: [Rd] corrupt data frame: columns will be truncated or padded with NAs in: format.data.frame(x, digits = digits)

2005-02-14 Thread Gregor GORJANC
Hello! Sending this also to r-help so anyone can read it also there and maybe also help me with my puzzle if this trivial and I don't see it. Prof Brian Ripley wrote: [... removed some ...] You add a column, not replace part of a non-existent column. Isn't that obvious, given what you wrote? #

Re: [R] gam(mgcv) starting values

2005-02-14 Thread Simon Wood
My guess is that your model is predicting a negative mean for some of your data. Since this is not possible for a Gamma r.v. the deviance calculation returns something non finite, which triggers the error message. This is possible because you have used an identity link. Is it not possible to

[R] R on solaris 10 on a 64-bit ultra60 sparc

2005-02-14 Thread Christoph Lehmann
Hi since I have no experience with solaris and sparc-architecture. We installed the latest Solaris 10 on a 64-bit ultra60 sparc machine. Since the solaris 10 is said to run native linux-applications: can I just download any r-binaries for linux? if yes, for which distribution? or are there any

Re: [R] Re: [Rd] corrupt data frame: columns will be truncated or padded with NAs in: format.data.frame(x, digits = digits)

2005-02-14 Thread Prof Brian Ripley
On Mon, 14 Feb 2005, Gregor GORJANC wrote: Sending this also to r-help so anyone can read it also there and maybe also help me with my puzzle if this trivial and I don't see it. Please don't, and especially do not after having removed the context. So I have replied only on R-devel. Prof Brian

Re: [R] R on solaris 10 on a 64-bit ultra60 sparc

2005-02-14 Thread Prof Brian Ripley
On Mon, 14 Feb 2005, Christoph Lehmann wrote: since I have no experience with solaris and sparc-architecture. We installed the latest Solaris 10 on a 64-bit ultra60 sparc machine. Since the solaris 10 is said to run native linux-applications: can I just download any r-binaries for linux? if

Re: [R] R on solaris 10 on a 64-bit ultra60 sparc

2005-02-14 Thread Uwe Ligges
Christoph Lehmann wrote: Hi since I have no experience with solaris and sparc-architecture. We installed the latest Solaris 10 on a 64-bit ultra60 sparc machine. Since the solaris 10 is said to run native linux-applications: can I just Native for what? 32 bit linux on x86 architecture, 64 bit

Re: [R] missing X11 graphics title bar

2005-02-14 Thread Steven K Thompson
Thanks, Peter. I think that helps narrow down the problem. The change of window manager seems to be a pretty decent solution. If the bug gets fixed or a workaround discovered, though, there might be some advantage to going with the gnome default. Thanks for the keystroke information (on my

[R] 'combinations' in gtools and stack overflow

2005-02-14 Thread Ho-Joon Lee
Dear R-users, Let me ask about the 'stack overflow' error which I got when I used the function 'combinations' in gtools. The following is what I did: - library(gtools) options(expressions=1e5) combinations(500, 3, 1:500) # or combinations(400, 2, 1:400) Error: protect(): stack overflow

RE: [R] combinations without repetition

2005-02-14 Thread davidr
These would more properly be called permutations of a multiset. You can find code on the net by searching for generate permutations multiset (though not in R AFAICS.) David L. Reiner -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] Sent: Sunday, February 13, 2005

Re: [R] 'combinations' in gtools and stack overflow

2005-02-14 Thread Prof Brian Ripley
On Mon, 14 Feb 2005, Ho-Joon Lee wrote: Dear R-users, Let me ask about the 'stack overflow' error which I got when I used the function 'combinations' in gtools. The following is what I did: - library(gtools) options(expressions=1e5) combinations(500, 3, 1:500) # or combinations(400, 2,

Re: [R] graphics - current filename

2005-02-14 Thread Thomas Lumley
On Mon, 14 Feb 2005, Uwe Ligges wrote: Paul Sorenson wrote: I would like to query R for the current (or last used) filename for a graphics device. Eg after png(filename=plot%02d.png) I would like something like the output of dev.cur() but with the %02d expanded to the current name. You cannot,

[R] DCUHRE - Multiple integration

2005-02-14 Thread CYRIL.CAILLAULT
Good morning, Do you know if there exists R package for the DCUHRE fortran routine like for the ADAPT package. Best regards Cyril [[alternative HTML version deleted]] __ R-help@stat.math.ethz.ch mailing list

Re: [R] 64 Bit R Background Question

2005-02-14 Thread Thomas Lumley
On Sun, 13 Feb 2005, Thomas Colson wrote: Hi, I've collected quite a bit of elevation data (LIDAR elevation points) and am looking for a suitable platform to do analysis and modeling on it. The data is sitting in an Oracle database, one table, 200 million rows of x,y, and z. I'm trying to figure

Re: [R] algorithms for matching and Hungarian method

2005-02-14 Thread Kurt Hornik
Martin Olivier writes: Hi all, I would like to match two partitions. That is, if I have exactly the same objects grouped together for the two partitions, the labels may be arbitrarly permuted. and so, i would like to know the correspondances of the groups between the two clusterings. In

Re: [R] graphics - current filename

2005-02-14 Thread Jan T. Kim
On Mon, Feb 14, 2005 at 07:23:23AM -0800, Thomas Lumley wrote: On Mon, 14 Feb 2005, Uwe Ligges wrote: Paul Sorenson wrote: I would like to query R for the current (or last used) filename for a graphics device. Eg after png(filename=plot%02d.png) I would like something like the output

RE: [R] 'combinations' in gtools and stack overflow

2005-02-14 Thread Warnes, Gregory R
The documentation for 'combinations' explicitly describes this problem: Details: Caution: The number of combinations and permutations increases rapidly with 'n' and 'r'!. To use values of 'n' above about 45, you will need to increase R's recursion limit. See the

RE: [R] 'combinations' in gtools and stack overflow

2005-02-14 Thread davidr
combn in package combinat does not use recursion, and so may be simpler to use and may be able to do larger sets. David L. Reiner -Original Message- From: Warnes, Gregory R [mailto:[EMAIL PROTECTED] Sent: Monday, February 14, 2005 11:01 AM To: Ho-Joon Lee Cc: r-help@stat.math.ethz.ch

Re: [R] 'combinations' in gtools and stack overflow

2005-02-14 Thread Ho-Joon Lee
Thank you for your replies, David and Brian. The function 'combn' in the package combinat works really fine without using any particular options. Thanks a lot again for referring me to the package, David. On the other hand, I do not know the meaning of stack overflow technically and am not

RE: [R] 'combinations' in gtools and stack overflow

2005-02-14 Thread Prof Brian Ripley
On Mon, 14 Feb 2005, Warnes, Gregory R wrote: The documentation for 'combinations' explicitly describes this problem: Not so: this is misinformation send in reply to correct and helpful comments. Details: Caution: The number of combinations and permutations increases rapidly with 'n'

Off topic -- large data sets. Was RE: [R] 64 Bit R Background Question

2005-02-14 Thread Berton Gunter
read all 200 million rows a pipe dream no matter what platform I'm using? In principle R can handle this with enough memory. However, 200 million rows and three columns is 4.8Gb of storage, and R usually needs a few times the size of the data for working space. You would likely be

Re: [R] 'combinations' in gtools and stack overflow

2005-02-14 Thread Ho-Joon Lee
Yes, 'combn' is much simpler to use than 'combinations'. Although Greg's 'combinations' is useful as well, 'combn' didn't produce the allocation error about vector size in my particular case. 'combn' works very well even in Windows. Thanks again, David. Lee [EMAIL PROTECTED] wrote: combn in

RE: Off topic -- large data sets. Was RE: [R] 64 Bit R Background Question

2005-02-14 Thread Thomas Colson
The purpose of investigating the entire (200 million record) data set is to investigate several interpolation models for creating gridded elevation data. Most models and algorithms do just that...take a manageable number of points and do the math. My reasoning behind using the entire dataset

[R] roll id

2005-02-14 Thread Omar Lakkis
I am collecting one price for market daily. date1 price1 date2 price2 date3 price3 I have a roll schedule where if date1 is between d1 and d2 then market is roll_id 0 date1 is between d2 and d3 then market is roll_id 1 date1 is between d3 and d4 then market is roll_id 2 and so on for

Re: Off topic -- large data sets. Was RE: [R] 64 Bit R Background Question

2005-02-14 Thread Prof Brian Ripley
On Mon, 14 Feb 2005, Berton Gunter wrote: read all 200 million rows a pipe dream no matter what platform I'm using? In principle R can handle this with enough memory. However, 200 million rows and three columns is 4.8Gb of storage, and R usually needs a few times the size of the data for working

Re: [R] roll id

2005-02-14 Thread Adaikalavan Ramasamy
See help(cut) for general or help(cut.POSIXt) for time related cut. On Mon, 2005-02-14 at 13:08 -0500, Omar Lakkis wrote: I am collecting one price for market daily. date1 price1 date2 price2 date3 price3 I have a roll schedule where if date1 is between d1 and d2 then market is

[R] testing equality of variances across groups in lme?

2005-02-14 Thread Bill Shipley
Hello. I am fitting a two-level mixed model which assumes equality of variance in the lowest-level residuals across groups. The call is: fit3-lme(CLnNAR~CLnRGR,data=meta.analysis, + na.action=na.omit,random=~1+CLnRGR|study.code) I want to test the assumption of equality of variances

[R] Re: Bubble Plot with Pie Charts

2005-02-14 Thread KMarkus
Could somebody provide a basic reference for bubble plots? I've not seen that term used before. Thanks in advance, Keith Markus [[alternative HTML version deleted]] __ R-help@stat.math.ethz.ch mailing list

[R] R News: Call for Papers

2005-02-14 Thread Paul Murrell
Dear useRs and developeRs, the next issue of `R News' is scheduled for the beginning of May and we are now accepting submissions for this first issue in 2005. For more information see http://cran.r-project.org/doc/Rnews/ If you are the author of a package on CRAN and you would like to

[R] a syntax question: $V100

2005-02-14 Thread WeiWei Shi
Hi, there: I have a syntax question. I have a dataset x with 100 variables. I did not set the column name so I used x$V1...x$V100. For my case, I need to put the number (e.g. 20) into another variable, like index so that I can refer to x$V20 by using something like x$V(index) but I don't know how

RE: [R] a syntax question: $V100

2005-02-14 Thread Liaw, Andy
x[,20] or x[[20]] will give you the 20th column of the data frame. Things like these are covered in `An Introduction to R'. If the column names are not necessarily in order, you can always use something like x[[paste(V, 20, sep=)]]. Andy From: WeiWei Shi Hi, there: I have a syntax

[Thank-you] [R] Help: sorting data frames

2005-02-14 Thread Barry Rountree
-BEGIN PGP SIGNED MESSAGE- Hash: SHA1 Hello, Back in 2002, you answered a question posted to the r-help mailing list, something like: Hello everyone: I have a data frame with two columns of data. I want to create a new data frame, in which the first column is sorted in accending

Re: [R] sample

2005-02-14 Thread T. Murlidharan Nair
My R version s 2.0.1 and I am running it under windows. I want to use R directly, but in this case since I could not figure out what went wrong I tried to transpose it in perl. Certainly I understand your concern regarding handling numerical data in Perl. I read the file using read.csv, it

[R] lattice multiple plots per page

2005-02-14 Thread Jeff Jorgensen
Dear R-sters, I was wondering if anyone has encountered the following issues. I've figured out how to get multiple levelplots [library(lattice)] on a single plot. However, when I add text (adding axis labels for the entire four panel plot) the text is missing when I insert the *.eps file I've

RE: [R] using poly in a linear regression in the presence of NA f ails (despite subsetting them out)

2005-02-14 Thread Liaw, Andy
This smells like a bug to me. The error is triggered by the line: variables - eval(predvars, data, env) inside model.frame.default(). At that point, na.action has not been applied, so poly() ended being called on data that still contains missing values. The qr() that issued the error is

RE: [R] testing equality of variances across groups in lme?

2005-02-14 Thread Lorenz . Gygax
I am fitting a two-level mixed model which assumes equality of variance in the lowest-level residuals across groups. The call is: fit3-lme(CLnNAR~CLnRGR,data=meta.analysis, + na.action=na.omit,random=~1+CLnRGR|study.code) I assume that CLnRGR is a factor and thus the groups which might

RE: [R] using poly in a linear regression in the presence of NA f ails (despite subsetting them out)

2005-02-14 Thread Prof Brian Ripley
Andy, I don't think it is a bug. The problem is that poly(x, 2) depends on the possible set of x values, and so needs to know all of them, unlike e.g. log(x) which is observation-by-observation. Silently omitting missing values is not a good idea in such cases, especially if the values are

[R] special symobol / character

2005-02-14 Thread FISCHER, Matthew
Hi all, Is it possible to add a permil (or per mille) symbol to an R plot (I couldn't find this symbol under demo(Hershey) or the plotmath information). In some ascii tables it is symbol no. 137. cheers, Matt. [[alternative HTML version deleted]]

[R] R News: Call for Papers

2005-02-14 Thread Torsten Hothorn
Dear useRs and developeRs, the next issue of `R News' is scheduled for the beginning of May and we are now accepting submissions for this first issue in 2005. For more information see http://cran.r-project.org/doc/Rnews/ If you are the author of a package on CRAN and you would like to

Re: [R] Re: [Rd] corrupt data frame: columns will be truncated or padded with NAs in: format.data.frame(x, digits = digits)

2005-02-14 Thread Gregor GORJANC
James, thanks for the response. I understand now my puzzle. tmp - data.frame(y1=1:4, f1=factor(c(A, B, C, D))) tmp$y2 - NA tmp[1:2, y2] - 2 Does the job. I see that I should add a full column. In my case adding full column of NAs and then adding values, solves my problems. Thanks to all. [EMAIL

Re: [R] special symobol / character

2005-02-14 Thread Eric Lecoutre
Hi Matthew, Most systems allow to enter any ASCII (or extended ASCII) character directly using a key combination. Accessing ANSI charcaters under Windows is possible with: ALT+0xxx (press ALT, hold it down, press 0 and the number of the character, release ALT) Thus: ALT+0137 makes: ‰ The future