Re: [R] assign values

2005-11-13 Thread Prof Brian Ripley
Please do read the posting guide. The answer to your question can be found by help(-) or help(=). The '-' can be used anywhere, but the '=' is only allowed at the top level (that is, in the complete expression typed by the user) or as one of the subexpressions in a braced

Re: [R] voronoi

2005-11-13 Thread Renaud Lancelot
See package tripack. You can find this information with RSiteSearch(voronoi, restrict = functions) Best, Renaud 2005/11/13, Robert [EMAIL PROTECTED]: Is there any pure r code to do delaunay or voronoi diagrams? Thanks! - [[alternative HTML version

Re: [R] computation on a table

2005-11-13 Thread Claus Atzenbeck
On Sun, 13 Nov 2005, Peter Dalgaard wrote: jim holtman [EMAIL PROTECTED] writes: [...] index - match(colnames(table2), colnames(table1), nomatch=0) t(t(table1[,index]) / table2[index != 0, drop=FALSE]) [...] or even sweep(table1, 2, table2[colnames(table1)], /) Perfect. I was not

Re: [R] voronoi

2005-11-13 Thread Jari Oksanen
On 13 Nov 2005, at 10:17, Renaud Lancelot wrote: See package tripack. You can find this information with RSiteSearch(voronoi, restrict = functions) So the answer is no. The question was about pure r code (with this capitalization), and voronoi.mosaic() in 'tripack' uses .Fortran. cheers,

[R] Legend

2005-11-13 Thread Mark Miller
I use the following to plot two graphs over each other and then insert a legend, but the two items in the legend both come up the same colour x = seq(0,30,0.01) plot(ecdf(complete), do.point=FALSE, main = 'Cummlative Plot of Monday IATs for Data and\n Fitted PDF over Entire 15 Weeks') lines(x,

Re: [R] Legend

2005-11-13 Thread Sundar Dorai-Raj
Mark Miller wrote: I use the following to plot two graphs over each other and then insert a legend, but the two items in the legend both come up the same colour x = seq(0,30,0.01) plot(ecdf(complete), do.point=FALSE, main = 'Cummlative Plot of Monday IATs for Data and\n Fitted PDF over

Re: [R] Legend

2005-11-13 Thread Adaikalavan Ramasamy
And you want to have different colored lines but black texts, try legend(x = 5, y = 0.2, legend = c(Data Set, Fitted PDF), col = c(black, red), lty=1) The advantage of this is that you can use dotted (lty option) or lines with different weights (lwd option). Regards, Adai On Sun,

[R] How to show numerical values on boxplots

2005-11-13 Thread Yue Li
Hi, dear all, I want to show numerical values with decimal points on the boxplots. Here is what I did: For example; x1-rnorm(100,2,2); x2-rexp(100); label-rep(1:2, rep(100,2)) median-round(c(median(x1), median(x2)),3) boxplot(c(x1, x2)~label, medpch=paste(median), medcex=1.2) It only shows

Re: [R] How to show numerical values on boxplots

2005-11-13 Thread Uwe Ligges
Yue Li wrote: Hi, dear all, I want to show numerical values with decimal points on the boxplots. Here is what I did: For example; x1-rnorm(100,2,2); x2-rexp(100); label-rep(1:2, rep(100,2)) median-round(c(median(x1), median(x2)),3) boxplot(c(x1, x2)~label, medpch=paste(median),

Re: [R] About: Error in FUN(X[[1]], ...) : symbol print-name too long

2005-11-13 Thread Uwe Ligges
gsmatos1 wrote: Hi, I´m trying to use the Win2BUGS package from R and I have a similar problem that reurns with the message: Error in FUN(X[[1]], ...) : symbol print-name too long I took a first look. The above looks like a bug in package rbugs which does some really strange things

Re: [R] computation on a table

2005-11-13 Thread Peter Dalgaard
Claus Atzenbeck [EMAIL PROTECTED] writes: On Sun, 13 Nov 2005, Peter Dalgaard wrote: jim holtman [EMAIL PROTECTED] writes: [...] index - match(colnames(table2), colnames(table1), nomatch=0) t(t(table1[,index]) / table2[index != 0, drop=FALSE]) [...] or even sweep(table1, 2,

Re: [R] assign values

2005-11-13 Thread Peter Dalgaard
Prof Brian Ripley [EMAIL PROTECTED] writes: Please do read the posting guide. The answer to your question can be found by help(-) or help(=). The '-' can be used anywhere, but the '=' is only allowed at the top level (that is, in the complete expression typed by the user)

[R] Julien Ruiz est absent.

2005-11-13 Thread julien . ruiz
Je serai absent(e) du 12/11/2005 au 16/11/2005. Je répondrai à votre message dès mon retour. I will be out of the office from 14-MAR-2005 until 18-MAR-2005 I will reply to your message on my return. Julien Ruiz L'acces immediat aux meilleurs tarifs Air France et au billet

[R] simulation of compound Poisson and Cox process

2005-11-13 Thread Barbora Kocúrová
Hello. I have this problem. It is modeling high-frequency financial data. The gamma OU process X (t), BDLP compound Poisson with intensity h 0 and E(a) exponential (a) distribution of jump. Lévy density w of Z (1): w(x) =

[R] problem with grid animation

2005-11-13 Thread Zepu Zhang
I'm trying to do animation with grid. Basically it's a vector field, like what 'quiver' in Matlab creates. I need to update it with grid.edit(). It seems grid erases the whole thing first, then redraws. Therefore the evident 'flash' between frames. Any way to avoid this flash? Thanks.

[R] selection of missing data

2005-11-13 Thread [EMAIL PROTECTED]
Hi i'm a french medical student, i have some data that i import from excel. My colomn of the datafram are the localisations of metastasis. If there is a metatsasis there is the symbol _. i want to exclude the row without metastasis wich represent the NA data. so, i wrote this mela is the data

Re: [R] selection of missing data

2005-11-13 Thread Adaikalavan Ramasamy
I do not quite follow your post but here are some suggestions. 1) You can the na.strings argument to simplify things df - read.delim(file=lala.txt, na.strings=- ) 2) If you can count the number of metastasis per row first, then find the rows with zero sum. met.cols -

[R] R for reliability analysis (with censored samples)?

2005-11-13 Thread Wolfgang Keller
Hello, first: I am a newbie to this list, so if this isn't the right place to ask, thanks for pointing me to the right direction... I'm currently working on adding a degree in RAMS engineering to my general engineering education, and consequently I'm looking for software to get my

Re: [R] assign values

2005-11-13 Thread Robert
Professor, What does if ((x=2)) 2 mean? Thanks, Robert Peter Dalgaard [EMAIL PROTECTED] wrote: Prof Brian Ripley writes: Please do read the posting guide. The answer to your question can be found by help(-) or help(=). The '-' can be used anywhere, but the '=' is only allowed at the top

Re: [R] voronoi

2005-11-13 Thread Robert
It uses FORTRAN code and not in pure R. Renaud Lancelot [EMAIL PROTECTED] wrote:See package tripack. You can find this information with RSiteSearch(voronoi, restrict = functions) Best, Renaud 2005/11/13, Robert : Is there any pure r code to do delaunay or voronoi diagrams? Thanks!

Re: [R] problem with grid animation

2005-11-13 Thread Paul Murrell
Hi Zepu Zhang wrote: I'm trying to do animation with grid. Basically it's a vector field, like what 'quiver' in Matlab creates. I need to update it with grid.edit(). It seems grid erases the whole thing first, then redraws. Therefore the evident 'flash' between frames. Correct. Any

Re: [R] voronoi

2005-11-13 Thread Roger Bivand
On Sun, 13 Nov 2005, Robert wrote: It uses FORTRAN code and not in pure R. The same applies to deldir - it also includes Fortran. So the answer seems to be no, there is no voronoi function only written in R. Renaud Lancelot [EMAIL PROTECTED] wrote:See package tripack. You can find this

Re: [R] assign values

2005-11-13 Thread Peter Dalgaard
Robert [EMAIL PROTECTED] writes: Professor, What does if ((x=2)) 2 mean? Thanks, Robert it assigns 2 to x, then uses the result (i.e. whether it is non-zero) as the condition for the if() construct. So in this case, the condition is always TRUE and the result is always 2. Peter Dalgaard

Re: [R] assign values

2005-11-13 Thread Robert
Peter Dalgaard [EMAIL PROTECTED] wrote: Robert writes: Professor, What does if ((x=2)) 2 mean? Thanks, Robert it assigns 2 to x, then uses the result (i.e. whether it is non-zero) as the condition for the if() construct. So in this case, the condition is always TRUE and the result is

[R] correlating irregular time series

2005-11-13 Thread paul sorenson
I have some time stamped events that are supposed to be unrelated. I have plotted them and that assumption does not appear to be valid. http://metrak.com/tmp/sevents.png is a plot showing three sets of events over time. For the purpose of this exercise, the Y value is irrelevant. The series

[R] open source and R

2005-11-13 Thread Robert
Roger Bivand [EMAIL PROTECTED] wrote: On Sun, 13 Nov 2005, Robert wrote: It uses FORTRAN code and not in pure R. The same applies to deldir - it also includes Fortran. So the answer seems to be no, there is no voronoi function only written in R. I am curious about one thing: since the

Re: [R] assign values

2005-11-13 Thread Roger Bivand
On Sun, 13 Nov 2005, Robert wrote: Peter Dalgaard [EMAIL PROTECTED] wrote: Robert writes: Professor, What does if ((x=2)) 2 mean? Thanks, Robert it assigns 2 to x, then uses the result (i.e. whether it is non-zero) as the condition for the if() construct. So in this case, the

Re: [R] open source and R

2005-11-13 Thread Roger Bivand
On Sun, 13 Nov 2005, Robert wrote: Roger Bivand [EMAIL PROTECTED] wrote: On Sun, 13 Nov 2005, Robert wrote: It uses FORTRAN code and not in pure R. The same applies to deldir - it also includes Fortran. So the answer seems to be no, there is no voronoi function only written in R.

Re: [R] open source and R

2005-11-13 Thread Robert
If I do not know C or FORTRAN, how can I fully understand the package or possibly improve it? Robert. Roger Bivand [EMAIL PROTECTED] wrote: On Sun, 13 Nov 2005, Robert wrote: Roger Bivand wrote: On Sun, 13 Nov 2005, Robert wrote: It uses FORTRAN code and not in pure R. The same

[R] Robust Non-linear Regression

2005-11-13 Thread Vermeiren, Hans [VRCBE]
Hi, I'm trying to use Robust non-linear regression to fit dose response curves. Maybe I didnt look good enough, but I dind't find robust methods for NON linear regression implemented in R. A method that looked good to me but is unfortunately not (yet) implemented in R is described in

Re: [R] open source and R

2005-11-13 Thread Roger Bivand
On Sun, 13 Nov 2005, Robert wrote: If I do not know C or FORTRAN, how can I fully understand the package or possibly improve it? By learning enough to see whether that makes a difference for your purposes. Life is hard, but that's what makes life interesting ... Robert. Roger Bivand

Re: [R] Robust Non-linear Regression

2005-11-13 Thread roger koenker
you might consider nlrq() in the quantreg package, which does median regression for nonlinear response functions url:www.econ.uiuc.edu/~rogerRoger Koenker email [EMAIL PROTECTED] Department of Economics vox:217-333-4558

Re: [R] open source and R

2005-11-13 Thread roger koenker
On Nov 13, 2005, at 3:24 PM, Robert wrote: I am curious about one thing: since the reason for using r is r is a easy-to-learn language and it is good for getting more people involved. Why most of the packages written in r use other languages such as FORTRAN's code? I understand some

Re: [R] open source and R

2005-11-13 Thread Mike Miller
On Sun, 13 Nov 2005, Roger Bivand wrote: On Sun, 13 Nov 2005, Robert wrote: If I do not know C or FORTRAN, how can I fully understand the package or possibly improve it? By learning enough to see whether that makes a difference for your purposes. Life is hard, but that's what makes life

Re: [R] open source and R

2005-11-13 Thread Ted Harding
On 13-Nov-05 Roger Bivand wrote: On Sun, 13 Nov 2005, Robert wrote: If I do not know C or FORTRAN, how can I fully understand the package or possibly improve it? By learning enough to see whether that makes a difference for your purposes. Life is hard, but that's what makes life

[R] poker package -- comments?

2005-11-13 Thread Duncan Murdoch
Over the weekend I wrote a small package to evaluate poker hands and to do some small simulations with them. If anyone is interested in looking at it, I'd appreciate comments and/or contributions. The package is available at http://www.stats.uwo.ca/faculty/murdoch/software. (Look at the

Re: [R] R for reliability analysis (with censored samples)?

2005-11-13 Thread Thomas Lumley
On Sun, 13 Nov 2005, Wolfgang Keller wrote: So here are the questions: How suitable is R for this kind of work (reliability analysis)? Does it handle (right-/left-/interval-)censored samples by default? The survreg() function in the survival package fits accelerated failure models to

Re: [R] open source and R

2005-11-13 Thread Robert
Thanks for all the nice discussions. Though different users have various needs from R, It's always good to stand on the shoulders of giants (as roger said). How far we will see depends our ability to understand what have been done by other languages. The package written in pure R might be

Re: [R] open source and R

2005-11-13 Thread Prof Brian Ripley
On Sun, 13 Nov 2005 [EMAIL PROTECTED] wrote: [...] There is one aspect though where R users are in the cold when it comes to C and FORTAN. If you want to understand the function 'eigen', say, then you can ?eigen to learn about its usage. You can enter eigen to see the R code, and indeed that