Re: [R] binding vectors or matrix using their names

2004-03-25 Thread Patrick Burns
Not at all -- I can obfuscate much better. [EMAIL PROTECTED] wrote: Goodness Patrick, this must surely qualify for the obfuscated R competition finals. I love it! [snip] A more serious, philosophical word on Patrick's solution. It is rarely necessary (in my limited experience, sure) to

[R] mlocal/mtrace inside a loop

2004-03-25 Thread Fred J.
Hello I need some help in figuring Bravington’s debugger out. Ok I have 2 functions, fun1 and fun2 saved in a ASCII file say filename is funs. Fun1 has a loop which calls fun2, fun2 has a loop which fails and I need to find out the value of the variables of the fun2 and fun1 loops at the specific

[R] book recommendation

2004-03-25 Thread Fred J.
Hello right, I need a good book to help me with R. background, programming in matlab, perl, visualBasics. os windows and Linux. some good books to match the need would be appreciated. Thanks __ [EMAIL PROTECTED] mailing list

Re: [R] binding vectors or matrix using their names

2004-03-25 Thread Peter Dalgaard
[EMAIL PROTECTED] writes: A more serious, philosophical word on Patrick's solution. It is rarely necessary (in my limited experience, sure) to have to use parse() like this. Where it provides a quick (kludgy?) solution I often find it a useful exercise to consider alternatives. They often

RE: [R] book recommendation

2004-03-25 Thread Ko-Kang Kevin Wang
Hi, -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Fred J. Sent: Thursday, March 25, 2004 9:25 PM To: r help Subject: [R] book recommendation Hello right, I need a good book to help me with R. background, programming in matlab, perl,

Re: [R] book recommendation

2004-03-25 Thread Prof Brian Ripley
I think you need to learn to use an FAQ: R has one covering books about R. The rest of your questions should be asked in a more appropriate place: the answers really depend on your level of knowledge. On Thu, 25 Mar 2004, Fred J. wrote: right, I need a good book to help me with R. background,

[R] R/ESS/Bioc* mailing lists: 2 hour time out -- 6 hours from now

2004-03-25 Thread Martin Maechler
This is (hopefully) a once in ten years event, happening today, from 17--19 CET (= 16--18 UTC = 11--13 US Eastern == 8--10 US Pacific Time) The main `vault' (security, air condition,..) where our mail server hypatia is located will undergo a major `network re-connection' event that

[R] Ordered Logistic Mixed Model

2004-03-25 Thread Leonard Assis
What is the best solution to fit an Ordered Logistic (Or Probit) Mixed Model in R? []s Leonard Assis Estatístico - CONFE 7439 __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide!

Re: [R] Calling data within a function (was First Variable in lm)

2004-03-25 Thread Gabor Grothendieck
Christian Hoffmann christian.hoffmann at wsl.ch writes: ggrothendieck at myway.com --- data(longley) lm( longley[,7] ~. , data = longley[,-7] ) You cannot call data() inside a function: To call data() within a function: f - function() { data(longley, envir =

Re: [R] binding vectors or matrix using their names

2004-03-25 Thread Rolf Turner
Well, I ***liked*** Patrick's approach. Why? Because it's basically straightforward. You put together a character string by pasting together the names of the objects you want to do things to, and the things (operations) you want to do to them. Then you get R to ``execute'' this string just as

[R] Error in 'legend' help?

2004-03-25 Thread Rau, Roland
Dear all, maybe I have misunderstood something but to me it seems like a minor error in the help for ?legend for the argument 'bg'. There it says: bg: the background color for the legend box. (Note that this is only used if 'bty = n'.) I think, however, that it should be

[R] nlme Crossed Random Effects

2004-03-25 Thread rdlandes
Dear R users: Is it possible to obtain crossed random effects in the nlme function? And if so, how? What follows is an example that pertains. The data set is from a (hypothetical) calibration experiment. There were two devices that were simultaneously calibrated using a reference instrument.

Re: [R] Error in 'legend' help?

2004-03-25 Thread Marc Schwartz
On Thu, 2004-03-25 at 07:17, Rau, Roland wrote: Dear all, maybe I have misunderstood something but to me it seems like a minor error in the help for ?legend for the argument 'bg'. There it says: bg: the background color for the legend box. (Note that this is only used

[R] subsetting based on vector

2004-03-25 Thread Vumani Dlamini
Dear R users; I am trying to write a small program which reads in a data set, and selects observations from certain years before the analysis. I have a problem including the selection criteria in the header of the program. Here is the problem; dataFIT-function(MODEL, MARGINS, yearConsidered){

Re: [R] subsetting based on vector

2004-03-25 Thread Winfried Theis
Hi! Well, you could try a subsetting rule like year %in% yearsconsidered . Hope this helps, Winfried On Thu, 2004-03-25 at 14:33, Vumani Dlamini wrote: Dear R users; I am trying to write a small program which reads in a data set, and selects observations from certain years before the

Re: [R] subsetting based on vector

2004-03-25 Thread Vumani Dlamini
Thanks man. It did the trick. Vumani From: Winfried Theis [EMAIL PROTECTED] To: Vumani Dlamini [EMAIL PROTECTED] CC: [EMAIL PROTECTED] Subject: Re: [R] subsetting based on vector Date: Thu, 25 Mar 2004 16:16:58 +0100 Hi! Well, you could try a subsetting rule like year %in% yearsconsidered .

[R] (no subject)

2004-03-25 Thread boudenne, jb
Hi, Does anybody has a simple example how to use a R function In Cpp programme? If it's possible, a *.cpp-file and what I need (wrapper or what ever). I try several ways, was reading the writing R extension and the windows-FAQ, but failed. I try to display R results in CGRIDControl

[R] downloadable mailing list archive in mbox format

2004-03-25 Thread Oleg Bartunov
Hi there, I'm beginner and would like to have mailing list archive to read it offline and don't disturb mailing list. I'd prefer mbox format. Regards, Oleg _ Oleg Bartunov, sci.researcher, hostmaster of AstroNet,

Re: [R] subsetting based on vector

2004-03-25 Thread Don MacQueen
Since you are interested in subsetting, a little search of the online help such as help.search('subset') will reveal the subset() function. Which you could use, for example, like this: dataFIT-function(MODEL, MARGINS, yearConsidered){ library(foreign

[R] Error : sink stack is full

2004-03-25 Thread ZABALZA-MEZGHANI Isabelle
Hello, I have implemented a method which uses sink to follow the progression status of an iterative process (Below is part of the code) I have already used such kind of code with no problem. Today, I get a sink stack is full error. I wonder if it could be linked with the fact that my .RData has

Re: [R] downloadable mailing list archive in mbox format

2004-03-25 Thread Oleg Bartunov
On Thu, 25 Mar 2004, Roger D. Peng wrote: You can get archives of the list at the URL listed at the bottom of this email. Thanks ! -roger Oleg Bartunov wrote: Hi there, I'm beginner and would like to have mailing list archive to read it offline and don't disturb mailing list.

[R] XEMACS, ESS, on Debian

2004-03-25 Thread Jean Eid
Dear all, This might not be the best place to post this (maybe the ESS help is better) but I will try anyways. I am running ESS, Xemacs, R 8.1 on a debian testing dist (i386 arch). when I envoke the edit(data) command, a data editor appears and all is fine. However executing C-c C-c will only

[R] How to add a top X-axis with a different logarithmic scale?

2004-03-25 Thread Itay Furman
Hi, I am trying to put on one plot two different logarithmic scales, using the bottom and top X-axes. Below there is an example of what I am trying to achieve, using axTicks() -- and fails. I already spent few hours on that, and cannot figure out from ?par and ?axTicks what I am doing wrong.

[R] g-inverse question

2004-03-25 Thread Richard Valliant
I am using the ginv function from MASS and have run across this problem that I do not understand. If I define the matrix A as below, its g-inverse does not satisfy the Moore-Penrose condition A %*% ginv(A) %*% A = A. The matrix A is X'WX in a quadratic regression using some very large dollar

RE: [R] colors, lines, characters .... documentation

2004-03-25 Thread Henrik Bengtsson
Oops, I forgot to check that plotSymbols(TRUE) worked with R --vanilla and it didn't as some of you already noticed; intToHex() and intToOct() were not defined (they're in my R.classes bundle at http://www.braju.com/R/). Below is a self-contained version that should work. Cheers Henrik Bengtsson

Re: [R] downloadable mailing list archive in mbox format

2004-03-25 Thread Roger D. Peng
You can get archives of the list at the URL listed at the bottom of this email. -roger Oleg Bartunov wrote: Hi there, I'm beginner and would like to have mailing list archive to read it offline and don't disturb mailing list. I'd prefer mbox format. Regards, Oleg

Re: [R] Error in 'legend' help?

2004-03-25 Thread Marc Schwartz
On Thu, 2004-03-25 at 09:59, Prof Brian Ripley wrote: On Thu, 25 Mar 2004, Marc Schwartz wrote: So the help _should_ read: bg: the background color for the legend box. (Note that this is only used if 'bty != n'.) and does so read in 1.9.0 beta: I fixed this on Jan 2.

Re: [R] Error in 'legend' help?

2004-03-25 Thread Prof Brian Ripley
On Thu, 25 Mar 2004, Marc Schwartz wrote: So the help _should_ read: bg: the background color for the legend box. (Note that this is only used if 'bty != n'.) and does so read in 1.9.0 beta: I fixed this on Jan 2. Can we remind people that R is beta-testing 1.9.0, and please

Re: [R] How to add a top X-axis with a different logarithmic scale?

2004-03-25 Thread Marc Schwartz
On Thu, 2004-03-25 at 12:24, Itay Furman wrote: Hi, I am trying to put on one plot two different logarithmic scales, using the bottom and top X-axes. Below there is an example of what I am trying to achieve, using axTicks() -- and fails. I already spent few hours on that, and cannot

[R] Aggregating frequency of irregular time series

2004-03-25 Thread Ivan Alves
Dear all, S-Plus has the function AggregateSeries() whose name is self explanatory. For instance one can derive monthly series from daily ones by specifying end-of-period, averages, sums, etc. I looked for a similar function in the packages its and tseries, but found nothing. I also

Re: [R] g-inverse question

2004-03-25 Thread Prof Brian Ripley
Do read the help file! ginv has a tolerance, and your matrix is below it: it is computationally singular to the tolerance of 1e-8. Try changing the tolerance On Thu, 25 Mar 2004, Richard Valliant wrote: I am using the ginv function from MASS and have run across this problem that I do

[R] S+Finmetrics cointegration functions

2004-03-25 Thread Ivan Alves
Dear all, S+Finmetrics has a number of very specilised functions. I am particularly interested in the estimation of cointegrated VARs (chapter 12 of Zivot and Wang). In this context the functions coint() and VECM() stand out. I looked at package dse1, but found no comparable functionality.

Re: [R] Aggregating frequency of irregular time series

2004-03-25 Thread Prof Brian Ripley
R itself has no support for irregular time series, but it does have an aggregate method for regular ones. You need to look into whichever package is handling irregular time series. However, it seems to me that this is not a time series problem at all: you have a set of observations whose

[R] rsync.r-project.org available again

2004-03-25 Thread Douglas Bates
The rsync site is back up. It may take a while for the new location to percolate through the DNS servers. The site should be [EMAIL PROTECTED]:~$ host rsync.r-project.org rsync.r-project.org is an alias for bates4.stat.wisc.edu. bates4.stat.wisc.edu has address 128.105.174.134 and not the

[R] factor based on pattern match ?

2004-03-25 Thread Oleg Bartunov
Hello, is't possible to specify pattern in levels ? y=c(ff,f,m,mm,fm,mf,ffm,mmf,mmm,fff); factor(y) [1] ff f m mm fm mf ffm mmf mmm fff Levels: f ff fff ffm fm m mf mm mmf mmm I want to specify levels using regexp (f.*,m.*) or use some another method. So, I could have 2 levels, say,

RE: [R] factor based on pattern match ?

2004-03-25 Thread Liaw, Andy
Is the following sort of what you want? y = factor(c(ff,f,m,mm,fm,mf,ffm,mmf,mmm,fff)) levels(y) - substring(levels(y), 1, 1) y [1] f f m m f m f m m f Levels: f m Andy From: Oleg Bartunov Hello, is't possible to specify pattern in levels ? y=c(ff,f,m,mm,fm,mf,ffm,mmf,mmm,fff);

Re: [R] Aggregating frequency of irregular time series

2004-03-25 Thread Ivan Alves
Thank you very much Brian. Indeed, by looking at ?tapply() this would do the job for regular time series (whose data-time INDEX renders itself naturally for bundling data in regular groups). However, with irregular time series the story is different, as some careful bundling is necessary

RE: [R] Aggregating frequency of irregular time series

2004-03-25 Thread Liaw, Andy
From: Ivan Alves Thank you very much Brian. Indeed, by looking at ?tapply() this would do the job for regular time series (whose data-time INDEX renders itself naturally for bundling data in regular groups). However, with irregular time series the story is different, as some

RE: [R] mlocal/mtrace inside a loop

2004-03-25 Thread Toby.Patterson
Look at ?mtrace for a start... And probably you need to look at the function go() -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Fred J. Sent: Thursday, March 25, 2004 8:16 PM To: r help Subject: [R] mlocal/mtrace inside a loop Hello I need some

RE: [R] factor based on pattern match ?

2004-03-25 Thread Oleg Bartunov
On Thu, 25 Mar 2004, Liaw, Andy wrote: Is the following sort of what you want? thanks for example, but I need something more :) In real situation, I'd like to use perl like patterns (like in grep), or just implicitly specify levels, for example: level F is (f.*, mmm) y =

Re: [R] How to add a top X-axis with a different logarithmic scale?

2004-03-25 Thread Itay Furman
Thanks, Marc, for your reply. I didn't explain myself well, but thanks to you I think I could phrase my question better. I want to be able to compute new tick-positions for the top X-axis. Why new ticks positions? See in the example below (that uses your contributed minimal code :-) what

[R] xyplot inside a for loop problem

2004-03-25 Thread Vaclav Petricek
Hello everyone Could someone please explain me why are xyplot() calls inside a for loop unsuccessful? Calling plot() is OK but xyplot() just opens the graphics window and that is all. No error, no warning :-( The same xyplot() outside for loop works fine.

[R] Using R's LAPACK Related files in Visual C++

2004-03-25 Thread Greg Tarpinian
I am a relative newcomer to both the R and C/C++ software worlds -- I'm taking a C Programming class currently. I noticed the other day that the C:\Program Files\R1_8_1\src\include\R_ext directory on my WinXP box has the header files BLAS.h Lapack.h Linpack.h RLapack.h I am interested in

[R] Using R's LAPACK Related files in Visual C++

2004-03-25 Thread Greg Tarpinian
I apologize for double-posting; I believe my first post was HTML. I am a relative newcomer to both the R and C/C++ software worlds -- I'm taking a C Programming class currently. I noticed the other day that the C:\Program Files\R1_8_1\src\include\R_ext directory on my WinXP box has the header

Re: [R] Aggregating frequency of irregular time series

2004-03-25 Thread Gabor Grothendieck
Assume r.its is an its time series and assume that chron has been loaded. chron represents dates as days since an origin so the week number can be obtained by dividing by 7. Then tapply the mean or other function. Finally, recreate an its. week - 7 * (as.numeric( chron( format( [EMAIL

[R] model fitting

2004-03-25 Thread J Fresen
I am trying to fit a normal linear model with response y and predictor x and two factors sex and group. I would like each combination of sex and group to have individual slopes and then subsequently have parallel slopes. I tried the model y ~ x*sex*group and it seemed to work for the first