Re: [R] Cannot update some packages after upgrade to 2.1.1

2005-07-16 Thread Prof Brian Ripley
-lf77blas is part of ATLAS, so I do suspect the RPM builder had ATLAS installed. lme4 needs a compatible Matrix installed. I do think installing from the sources would solve this, but probably you need to discuss this with the RPM maintainer as a dependency appears to be missing. On Fri, 15

Re: [R] Coxph with factors

2005-07-16 Thread Kylie-Anne Richards
Thank you for your help. In any case, to specify f.pom You need it to be a factor with the same set of levels. You don't say what the lowest level of pom is, but if it is, say, -3. f.pom=factor(-3, levels=seq(-3,2.5, by=0.5))

[R] cbind a list of matrices

2005-07-16 Thread Mauro Gasparini
Dear users, I have a list of several matrices with the same number of columns, how do I rbind them all with a vectorized command? A related simpler question is, how do I vectorize the instruction that rbinds together several copies of the same matrix? Thanks. -- Mauro Gasparini

Re: [R] R: to the power

2005-07-16 Thread Bernardo Rangel Tura
At 10:11 12/7/2005, [EMAIL PROTECTED] wrote: hi all why does R do this: (-8)^(1/3)=NaN the answer should be : -2 Allan In my computer: (-8)^(1/3) [1] NaN -8^(1/3) [1] -2 -(8^(1/3)) [1] -2 The problem is -8 or the problem is (-8) ? []s Tura -- No virus found in this outgoing

Re: [R] cbind a list of matrices

2005-07-16 Thread Ted Harding
On 16-Jul-05 Mauro Gasparini wrote: Dear users, I have a list of several matrices with the same number of columns, how do I rbind them all with a vectorized command? A related simpler question is, how do I vectorize the instruction that rbinds together several copies of the same matrix?

Re: [R] cbind a list of matrices

2005-07-16 Thread Peter Dalgaard
(Ted Harding) [EMAIL PROTECTED] writes: On 16-Jul-05 Mauro Gasparini wrote: Dear users, I have a list of several matrices with the same number of columns, how do I rbind them all with a vectorized command? A related simpler question is, how do I vectorize the instruction that

Re: [R] Can't get sample function from An Introduction to R to work

2005-07-16 Thread Spencer Graves
When you copied it into MS Word, did you Copy - paste special - unformatted text? This trick sometimes might expose (or eliminate) nonprinting characters or characters with special attributes that might have created the problem. spencer graves David Groos wrote: Thank

Re: [R] rmpi in windows

2005-07-16 Thread Uwe Ligges
Uzuner, Tolga wrote: Hi Folks, Has anyone been able to get rmpi to work under windows ? Rmpi uses the LAM implementation of MPI. See http://www.lam-mpi.org/faq/category12.php3 and read FAQ 2 which implicitly tells us that there is no native port, hence you cannot run it under Windows. The

Re: [R] Proportion test in three-chices experiment

2005-07-16 Thread Spencer Graves
Have you considered BTm in library(BradleyTerry)? Consider the following example: cond1 - data.frame(winner=rep(LETTERS[1:3], e=2), + loser=c(B,C,A,C,A,B), + Freq=1:6) cond2 - data.frame(winner=rep(LETTERS[1:3], e=2), + loser=c(B,C,A,C,A,B), +

Re: [R] cbind a list of matrices

2005-07-16 Thread S.O. Nyangoma
Didn't you simply try: A-matrix(c(1.1,1.2,1.3,1.4,1.5,1.6),ncol=3) B-matrix(c(2.1,2.2,2.3,2.4,2.5,2.6),ncol=3) C-matrix(c(3.1,3.2,3.3,3.4,3.5,3.6),ncol=3) A [,1] [,2] [,3] [1,] 1.1 1.3 1.5 [2,] 1.2 1.4 1.6 B [,1] [,2] [,3] [1,] 2.1 2.3 2.5 [2,] 2.2 2.4 2.6

Re: [R] Padding in lattice plots

2005-07-16 Thread Federico Gherardini
On Friday 15 July 2005 17:00, Deepayan Sarkar wrote: On 7/15/05, Federico Gherardini [EMAIL PROTECTED] wrote: On Friday 15 July 2005 14:42, you wrote: Hi all, I've used the split argument to print four lattice plots on a single page. The problem now is that I need to reduce the amount

[R] summary: cbind a list of matrices

2005-07-16 Thread Mauro Gasparini
Thanks to all kind people who answered. Attached is a useful reply, among many others. -- Mauro Gasparini Professore Straordinario di Statistica Dipartimento di Matematica, Politecnico di Torino Corso Duca degli Abruzzi 24 I-10129 Torino, Italy tel: +39 011 564 7546 fax: +39 011 564 7599

Re: [R] Can't get sample function from An Introduction to R to work

2005-07-16 Thread Thomas Lumley
On Fri, 15 Jul 2005, David Groos wrote: Thank you-all very much for your help, your responses and help has been very encouraging. The following doesn't close the case but it tables it... First I copied Ken's code into my R Console and...it worked great! That was baffling as it looked

Re: [R] Coxph with factors

2005-07-16 Thread Thomas Lumley
On Sat, 16 Jul 2005, Kylie-Anne Richards wrote: Thank you for your help. In any case, to specify f.pom You need it to be a factor with the same set of levels. You don't say what the lowest level of pom is, but if it is, say, -3.

Re: [R] R: to the power

2005-07-16 Thread Thomas Lumley
On Sat, 16 Jul 2005, Bernardo Rangel Tura wrote: At 10:11 12/7/2005, [EMAIL PROTECTED] wrote: hi all why does R do this: (-8)^(1/3)=NaN the answer should be : -2 Yes and no. The problem is that the reciprocal of 3 is not exactly representable as a floating point number (it has an

Re: [R] Coxph with factors

2005-07-16 Thread Thomas Lumley
On Sat, 16 Jul 2005, Thomas Lumley wrote: Yes, but you don't need to go via the baseline. The survival curves for any two covariate vectors z1 and z2 are related by S(t; z1)= S(t; z2)^(z1-z2) Actually S(t; z1)=S(t;z2) ^(beta'(z1-z2)) of course. -thomas

Re: [R] Proportion test in three-chices experiment

2005-07-16 Thread Jonathan Baron
I suspect that there are more direct ways to do this test, but it is unclear to me just what the issue is. For example, if there are many subjects and very few stimuli for each, you might want to get some sort of measure of ability for each subject (many possibilities here, then test the measure

[R] topical guide to R packages

2005-07-16 Thread Vivek Rao
I would like to see R packages arranged by topic. CRAN Task Views are at http://lib.stat.cmu.edu/R/CRAN/src/contrib/Views/ , but I'd like something more detailed. For example, the IMSL Fortran library, version 4 is easy to navigate and has procedures arranged according to following topics:

[R] xlabel and ylabel in ROC curve ??

2005-07-16 Thread luk
When I draw a ROC curve, the xlabel is '1-spec:', and ylabel is 'sens:' . My question is: can I specify them to '1-specificity', and 'sensitivity' ? thanks for your help. Lu. __ [[alternative HTML version deleted]]

[R] Confidence Intervals for Arbitrary Functions

2005-07-16 Thread Jeff Newmiller
I have a rather basic background in statistics, and am looking for assistance in solving what I expect is a common type of problem. I have measurements of physical processes, and mathematical models of those processes that I want to feed the measurements into. A simple case is using measurements

Re: [R] topical guide to R packages

2005-07-16 Thread Duncan Murdoch
On 7/16/2005 1:17 PM, Vivek Rao wrote: I would like to see R packages arranged by topic. CRAN Task Views are at http://lib.stat.cmu.edu/R/CRAN/src/contrib/Views/ , but I'd like something more detailed. For example, the IMSL Fortran library, version 4 is easy to navigate and has procedures

[R] PBSmapping and shapefiles

2005-07-16 Thread Denis Chabot
Hi, Is there a way, preferably with R, to read shapefiles and transform them in a format that I could then use with package PBSmapping? I have been able to read such files into R with maptools' read.shape and plot it with plot.Map, but I'd like to bring the data to PBSmapping and plot from

Re: [R] Confidence Intervals for Arbitrary Functions

2005-07-16 Thread Gabor Grothendieck
On 7/16/05, Jeff Newmiller [EMAIL PROTECTED] wrote: I have a rather basic background in statistics, and am looking for assistance in solving what I expect is a common type of problem. I have measurements of physical processes, and mathematical models of those processes that I want to feed

Re: [R] R: to the power

2005-07-16 Thread François Pinard
[Thomas Lumley] It would be nice if R could realize that you meant the cube root of -8, but that requires either magical powers or complicated and unreliable heuristics. The real solution might be a function like root(x,a,b) to compute x^(a/b), where a and b could then be exactly

[R] xfig device - depth

2005-07-16 Thread Thomas Zanon
Hi, I hope this is the right list for my posting, since I've never posted to any R list before. I'm quite extensively using the xfig graphics device and as far as I figured out this device writes all the objects into xfig layer 100 (based on what I saw in the devPS.c file -if this is the file

[R] Time Series Count Models

2005-07-16 Thread Brett Gordon
Hello, I'm trying to model the entry of certain firms into a larger number of distinct markets over time. I have a short time series, but a large cross section (small T, big N). I have both time varying and non-time varying variables. Additionally, since I'm modeling entry of firms, it seems

[R] Time Series Count Models

2005-07-16 Thread Brett Gordon
Hello, I'm trying to model the entry of certain firms into a larger number of distinct markets over time. I have a short time series, but a large cross section (small T, big N). I have both time varying and non-time varying variables. Additionally, since I'm modeling entry of firms, it seems

[R] Is it possible to coerce R to continue proceeding the next command in a loop after an error message ?

2005-07-16 Thread Yimeng Lu
Hello R-users, In a loop, if a function, such as nls, gives an error, is it possible to coerce R to continue proceeding the next command with the same loop? Thanks so much for your advice! Hanna Lu [[alternative HTML version deleted]] __

Re: [R] Is it possible to coerce R to continue proceeding the next command in a loop after an error message ?

2005-07-16 Thread Duncan Murdoch
Yimeng Lu wrote: Hello R-users, In a loop, if a function, such as nls, gives an error, is it possible to coerce R to continue proceeding the next command with the same loop? Yes, see the try() function. The basic usage is something like value - try( some calculation ) if

Re: [R] Is it possible to coerce R to continue proceeding the next command in a loop after an error message ?

2005-07-16 Thread Spencer Graves
?try spencer graves Yimeng Lu wrote: Hello R-users, In a loop, if a function, such as nls, gives an error, is it possible to coerce R to continue proceeding the next command with the same loop? Thanks so much for your advice! Hanna Lu [[alternative HTML

Re: [R] nlme and spatially correlated errors

2005-07-16 Thread Spencer Graves
Have you tried anova(fit1, fit2), where fit1 - lme(one model...) fit2 - lme(a submodel ... ) This anova does about the best that anyone knows how to do -- or at lest did 7 years ago when it was written. If the submodel changes the fixed effects, you should use

[R] printing the name of the arguments passed to a function

2005-07-16 Thread John Sorkin
R2.1.1 Win 2k I have a function, B, within a function, A. I would like to have B print the name of the argument passed to it (not the value of the arguments). i.e., A-function() { B-function(x,y) { fit1-lm(y~x,data=jo) print(summary(fit1) I want B to print

Re: [R] Can't get sample function from An Introduction to R to work

2005-07-16 Thread David Groos
On Jul 16, 2005, at 9:58 AM, Thomas Lumley wrote: On Fri, 15 Jul 2005, David Groos wrote: Thank you-all very much for your help, your responses and help has been very encouraging. The following doesn't close the case but it tables it... First I copied Ken's code into my R Console

Re: [R] 2D contour predictions

2005-07-16 Thread Spencer Graves
It's not so simple, but consider the following: x=rep(1:6, each=2, length=24) y=rep(1:6, each=4) z=rep(0:1, length=24) set.seed(1) tstDF - data.frame(x=x, y=y, z=z, w=(x-3.5)^2+(y-3.5)+z+0.1*rnorm(24)) fit - lm(w~x+y+z+I(x^2)+I(y^2), tstDF) x0 - seq(1, 5, .5) y0 - seq(1, 6,

Re: [R] printing the name of the arguments passed to a function

2005-07-16 Thread Spencer Graves
How about the following: tstFn - function(x){ + xName - deparse(substitute(x)) + cat(xName) + done + } tstFn(Varname) Varname[1] done tstF1 - function(x){ + tstF2 - function(y){ + cat(deparse(substitute(y)), + doesn't work\n) + } + tstF2(x) + xName -

[R] Strange problems with lattice plots

2005-07-16 Thread Wladimir Eremeev
Dear r-help, I need to draw some lattice plots and have stuck in strange problems. I have a data frame of 405 rows: str(dframe) `data.frame': 405 obs. of 4 variables: $ year : num 1990 1990 1990 1990 1990 1990 1990 1990 1990 1990 ... $ month: num 1 2 3 4 5 6 7 8 9 10 ... $ V_A2 :

Re: [R] Strange problems with lattice plots

2005-07-16 Thread Wladimir Eremeev
Yes!!! I have overcome it! The solution is xyplot((V_A3/25)~year|factor(month), strip=strip.custom(factor.levels=c(month.name,Annual mean, Annual integral,Integral (Jan-Mar),Integral (Apr-Jun),Integral (Jul-Sep),Integral (Oct-Dec), and so on...