Re: [R] How to plot two variables using a secondary Y axis

2007-07-11 Thread hadley wickham
On 7/10/07, Felipe Carrillo [EMAIL PROTECTED] wrote: Date Fo Co6/27/2007 57.1 13.96/28/2007 57.7 14.3 6/29/2007 57.8 14.36/30/2007 57 13.97/1/2007 57.1 13.9 7/2/2007 57.2 14.07/3/2007 57.3 14.17/4/2007 57.6 14.2 7/5/2007 58 14.4

[R] Window Version for permtest

2007-07-11 Thread amna khan
Hi Sir There is no Window Binary versin of package permtest.Please provide information in this regard. Thank you -- AMINA SHAHZADI Department of Statistics GC University Lahore, Pakistan. Email: [EMAIL PROTECTED] [EMAIL PROTECTED] [EMAIL PROTECTED] [[alternative HTML version deleted]]

Re: [R] Window Version for permtest

2007-07-11 Thread Prof Brian Ripley
On Wed, 11 Jul 2007, amna khan wrote: Hi Sir There is no Window Binary versin of package permtest.Please provide information in this regard. See http://cran.r-project.org/bin/windows/contrib/2.5/check/permtest-check.log [Sending this once was quite sufficient.] -- Brian D. Ripley,

[R] Plot SpatialLinesDataFrame with xlim ylim

2007-07-11 Thread Edzer J. Pebesma
Michael, The plot method for SpatialLinesDataFrame objects resides in package sp, and questions regarding it are easier noticed on the r-sig-geo mailing list. The reason why they are plotted with aspect ratio 1 is that they are assumed to be spatial (geographical) data, and assume that 1 m

[R] Improved Windows Vista compatibility

2007-07-11 Thread Prof Brian Ripley
Duncan Murdoch and I have been working on Vista compatibility. The first build of R-2.5.1-win32.exe was built with a version of the installer that had some compatibility issues, and it has been rebuilt under the latest version. What now happens is that if you want to install R in an

[R] Installation on Leopard

2007-07-11 Thread Helin Gai
Hi all, I'm wondering whether anyone has installed R on Mac OS X leopard with success. I downloaded the latest version, but couldn't have it installed. Any idea as to how I get around the problem? Thanks! Sincerely, Helin -- Helin (Colin) Gai Class of 2009, Duke University Box 96332 Durham, NC

[R] Making Gehan-Breslow test for Survival data

2007-07-11 Thread Cody Hamilton
Jose, The Gehan-Breslow test provides a generalization of the Kruskal-Wallis test for censored data. As an alternative, try using survdiff with rho=1. This method uses weights w(ti) = S(ti) (where S is the Kaplan-Meier estimate of survival) which yields Fleming and Harrington's version of

[R] Previously saved workspace restored

2007-07-11 Thread Kristi Glover
hi there, i an beginner of R. some one have sent me a file (extension is .Rdata). i have installed R in my computer and i just double clicked the data. then it automatically opened R programme and displayed that [previously saved workspace restored]. the following message was displayed.

Re: [R] ECDF, distribution of Pareto, distribution of Normal

2007-07-11 Thread livia
Thank you very much for your reply. I am afraid I have no idea what is wrong with the pgpg function. The parameters are generated from pre-fitted GPD distribution. 1.544 is the location parameter, 0.4373 is the scale parameter and -0.2398 is the shape parameter. Cound you please give me some

Re: [R] [R-sig-Geo] Plot SpatialLinesDataFrame with xlim ylim

2007-07-11 Thread Roger Bivand
On Wed, 11 Jul 2007, Edzer J. Pebesma wrote: Michael, The plot method for SpatialLinesDataFrame objects resides in package sp, and questions regarding it are easier noticed on the r-sig-geo mailing list. The reason why they are plotted with aspect ratio 1 is that they are assumed to be

Re: [R] type III ANOVA for a nested linear model

2007-07-11 Thread Mark Difford
Indeed! And, apropos of the expression, to be Ripleyed (and so be condemned to eating cookies for a long, long time), what about being Billasted? BestR, Mark. Simon Blomberg-4 wrote: I second the nomination! Simon. On Tue, 2007-07-10 at 10:02 -0600, Greg Snow wrote: I nominate the

[R] question about gamm models

2007-07-11 Thread Matthias an der Heiden
Dear R-Users, I have a question concerning mixed models in R. The strata of my model are the counties of Germany. The differencies between these counties should be modelled as realizations of a normally distributed random variable X. Moreover, the model contains a 0/1 variable A that enters as

Re: [R] variable and value labels

2007-07-11 Thread Julien Barnier
Hi Steve, I am new here - enjoying the power of R compared to SPSS. Looking for sets of tips and tricks for people with old SPSS habits. In particular, would like to know an easy way to set variable labels across a dataset and to set value labels for sets of variables. Maybe you should take

[R] Changing default library

2007-07-11 Thread Schmitt, Corinna
Hallo R experts, I got a question concerning the .libPaths(). if I do the command .libPaths() than the result is /usr/lib/R/library. This is the default folder. I now want to change this one into /home/csc/usr/lib/R/library. I thought it would work with the command

Re: [R] Changing default library

2007-07-11 Thread Julien Barnier
Hi, When I than do the command .libPaths() the result is: /usr/lib/R/library /home/csc/usr/lib/R/library. But if I start R the next time the result of the command .libPaths() is again just /usr/lib/R/library. From libPaths help page : , | The library search path is initialized at

Re: [R] Changing default library

2007-07-11 Thread Achim Zeileis
Corinna: I got a question concerning the .libPaths(). if I do the command .libPaths() than the result is /usr/lib/R/library. This is the default folder. I now want to change this one into /home/csc/usr/lib/R/library. I thought it would work with the command

[R] inquiry about anova and ancova

2007-07-11 Thread Anderson, Mary-Jane
Dear R users, I have a rather knotty analysis problem and I was hoping that someone on this list would be able to help. I was advised to try this list by a colleague who uses R but it is a statistical inquiry not about how to use R. In brief I have a 3x2 anova, 2 tasks under 3

[R] variable and value labels

2007-07-11 Thread Steve Powell
Dear list I am new here - enjoying the power of R compared to SPSS. Looking for sets of tips and tricks for people with old SPSS habits. In particular, would like to know an easy way to set variable labels across a dataset and to set value labels for sets of variables. Grateful for any help,

Re: [R] inquiry about anova and ancova

2007-07-11 Thread Peter Dalgaard
Anderson, Mary-Jane wrote: Dear R users, I have a rather knotty analysis problem and I was hoping that someone on this list would be able to help. I was advised to try this list by a colleague who uses R but it is a statistical inquiry not about how to use R. In brief I have a 3x2

[R] p-value from survreg, library(survival)

2007-07-11 Thread vladosr
dear r experts: I would appreciate advice on how to get the p-value from the object 'sr' created with the function survreg() as given below. vlad sr-survreg(s~groups, dist=gaussian) Coefficients: (Intercept) groups -0.02138485 0.03868351 Scale= 0.01789372 Loglik(model)= 31.1

Re: [R] inquiry about anova and ancova

2007-07-11 Thread Andrew Robinson
Hi Mary, it sounds like you have a split-plot design, or more gruesomely, a split-subject design. The model that I infer from your description of the design can be fit using the lme() function of the nlme() package, along the lines of a similar analysis documented in section 1.6 of Pinheiro

Re: [R] Warning message: cannot create HTML package index

2007-07-11 Thread Prof Brian Ripley
On Fri, 6 Jul 2007, Peter Dalgaard wrote: Leo wrote: On 06/07/2007, Prof Brian Ripley wrote: On Fri, 6 Jul 2007, Leo wrote: I have set R_LIBS=~/R_lib as I don't have root access. The following message shown up every time after installing a package: .. The downloaded packages are

Re: [R] Installation on Leopard

2007-07-11 Thread Douglas Bates
On 7/11/07, Helin Gai [EMAIL PROTECTED] wrote: Hi all, I'm wondering whether anyone has installed R on Mac OS X leopard with success. I downloaded the latest version, but couldn't have it installed. Any idea as to how I get around the problem? Thanks! Simon Urbanek has comments about R for

[R] p-value from survreg(), library(survival)

2007-07-11 Thread Vlado Sremac
dear r experts: It seems my message got spam filtered, another try: i would appreciate advice on how to get the p-value from the object 'sr' created with the function survreg() as given below. vlad sr-survreg(s~groups, dist=gaussian) Coefficients: (Intercept) groups -0.02138485 0.03868351

[R] Some questions about quadratic programming (QP)

2007-07-11 Thread Amir Safari
Dear R Users , As a beginner in QP, I'm trying to solve a Support Vector Machine problem by a QP. In particulare I am using the quadprog package. My questions are here: 1- In the document for the package (The quadprog Package), the inequality constraint is

[R] CDF for pareto distribution

2007-07-11 Thread livia
Hi, I would like to use the following codes to plot the CDF for pareto distribution. Before doing this, I have plot the emperical one. x - seq(1.6, 3, 0.1) lines(x,pgpd(x, 1.544,0.4477557,), col=red) Could anyone give me some advice whether the above codes are correct? Many thanks. -- View

[R] Antigen found FILE FILTER= *.pif file

2007-07-11 Thread Antigen_UM-TSMTPOUT1
Antigen for Exchange found instruction.zip__instruction.doc .pif matching FILE FILTER= *.pif file

[R] RWeka control parameters classifiers interface

2007-07-11 Thread strinz
Hello, I have some trouble in achieving the desired parametrisation for the weka classifier functions, using the package RWeka. The problem is, that the functions result=classifier(formula, data, subset, na.action, control = Weka_control(mycontrol)) do not seem to be manipulated by

Re: [R] Previously saved workspace restored

2007-07-11 Thread john seers \(IFR\)
Hi If you enter the command ls()you will see a list of names that have come with the .Rdata file you double-clicked. If you enter one of these names at the command prompt you will see the data. So, for example if you have some data called mydata: ls() [1] mydata repos mydata [,1]

Re: [R] type III ANOVA for a nested linear model

2007-07-11 Thread Carsten Jaeger
Hello Peter, thanks for your help. I'm quite sure that I specified the right model. Factor C is indeed nested within factor A. I think you were confused by the numbering of C (1..11), and it is easier to understand when I code it as you suggested (1,2,3 within each level of A, as in mydata1 [see

Re: [R] RWeka control parameters classifiers interface

2007-07-11 Thread Achim Zeileis
On Wed, 11 Jul 2007 [EMAIL PROTECTED] wrote: The problem is, that the functions result=classifier(formula, data, subset, na.action, control = Weka_control(mycontrol)) do not seem to be manipulated by the mycontrol- arguments Yes, they are...not all parameter changes have always an

Re: [R] variable and value labels

2007-07-11 Thread Frank E Harrell Jr
Steve Powell wrote: Dear list I am new here - enjoying the power of R compared to SPSS. Looking for sets of tips and tricks for people with old SPSS habits. In particular, would like to know an easy way to set variable labels across a dataset and to set value labels for sets of variables.

[R] 3D plot and interactive PDFs

2007-07-11 Thread Bruno C\.
With version 8 of acrobat reader, it is now possible to have 3D in PDf documents. Does it exist already an R package who manage to produce 3D plots which can be saved as interactive 3D graphs in a PDF file? Best Regards Bruno Cavestro --

[R] elementary statistics with R (rkward?)

2007-07-11 Thread Donatas G.
Hi, I am trying to learn some basic statistics stuff but I cannot find any elementary statistics exercises using R language. Using RKward would be even better... I need that in analysing sociological data, obtained through questionnairres - findind corelations between variables, relations between

Re: [R] elementary statistics with R (rkward?)

2007-07-11 Thread Sarah Goslee
I don't know anything about RKward, but there are many, many tutorials, guides and other documents written for people learning R available online. Try the introduction to R at: http://www.r-project.org/ under manuals, or some of the many fine contributions at:

Re: [R] elementary statistics with R (rkward?)

2007-07-11 Thread Charles Annis, P.E.
Face the music and buy the book: _Introductory Statistics with R_ by Peter Dalgaard. It's perfect for what you need. It's clear and concise and will teach you statistics AND R as painlessly as such a thing can be. It's inexpensive and you can get it on Amazon.com and every other major

Re: [R] CDF for pareto distribution

2007-07-11 Thread Vincent Goulet
Le 07-07-11 à 07:56, livia a écrit : Hi, I would like to use the following codes to plot the CDF for pareto distribution. Before doing this, I have plot the emperical one. x - seq(1.6, 3, 0.1) lines(x,pgpd(x, 1.544,0.4477557,), col=red) Could anyone give me some advice whether the above

Re: [R] Previously saved workspace restored

2007-07-11 Thread Kristi Glover
hi, thanks for the suggestion, however i am struggling to see the data, here is load(C://R//species.Rdata) ls()[1] abc abc1 abc2 example example1.14 generalislast.warningwhen i wrote ls(), then i saw the the files names I created before in R. there was no

Re: [R] CDF for pareto distribution

2007-07-11 Thread livia
Hi, thank you very much for your reply. The function pgpd() is from the package POT, and the 1.544 is the location parameter, 0.4477557 is the scale parameter and -0.50113 is the shape parameter, which can be both negtive or positive. Vincent Goulet wrote: Le 07-07-11 à 07:56, livia a écrit

Re: [R] p-value from survreg(), library(survival)

2007-07-11 Thread hadley wickham
str(survreg(s~groups, dist=gaussian)) is probably a good place to start. Hadley On 7/11/07, Vlado Sremac [EMAIL PROTECTED] wrote: dear r experts: It seems my message got spam filtered, another try: i would appreciate advice on how to get the p-value from the object 'sr' created with the

[R] Power calculation for the time series experiment

2007-07-11 Thread A Ezhil
Hi All, We are planning to run an experiment, where samples will be taken at different time points (say, 0, 4, 8, 16, 24). If I am interested in the effect size of 1.5 for a reasonably large samples (say 500), what will be the power? Is it a good idea to use F-test (one-way ANOVA) as my test

Re: [R] 3D plot and interactive PDFs

2007-07-11 Thread Duncan Murdoch
On 7/11/2007 9:18 AM, Bruno C. wrote: With version 8 of acrobat reader, it is now possible to have 3D in PDf documents. Does it exist already an R package who manage to produce 3D plots which can be saved as interactive 3D graphs in a PDF file? No, not as far as I know. If you want to help

Re: [R] type III ANOVA for a nested linear model

2007-07-11 Thread Andrew Robinson
Billasted sounds too brutal. How about Billeted, as in what one does to one's breshly-caught bish? Andrew On Wed, Jul 11, 2007 at 01:48:10AM -0700, Mark Difford wrote: Indeed! And, apropos of the expression, to be Ripleyed (and so be condemned to eating cookies for a long, long time),

Re: [R] ECDF, distribution of Pareto, distribution of Normal

2007-07-11 Thread J. R. M. Hosking
livia wrote: Hello all, I would like to plot the emperical CDF, normal CDF and pareto CDF in the same graph and I amusing the following codes. z is a vector and I just need the part when z between 1.6 and 3. plot(ecdf(z), do.points=FALSE, verticals=TRUE,

Re: [R] p-value from survreg(), library(survival)

2007-07-11 Thread Henrique Dallazuanna
Try also: pchisq(summary(sr)$chi, degrees_freedom, lower=FALSE) *You need know your degrees of freedom -- Henrique Dallazuanna Curitiba-Paraná-Brasil 25° 25' 40 S 49° 16' 22 O On 11/07/07, Vlado Sremac [EMAIL PROTECTED] wrote: dear r experts: It seems my message got spam filtered, another

[R] survfit for interval censored data

2007-07-11 Thread Sandra Ellis
pHello,/ppI am a new R-user and would like to use survfit for interval censored data.  Whenever I try I get an error message that states I can only use survfit for right censored or counting process data.  I was wondering if anyone knows if there is an additional package available that can 

Re: [R] p-value from survreg(), library(survival)

2007-07-11 Thread Marc Schwartz
Actually, in this case, looking at the code for: survival:::print.survreg would be better, as the p value is calculate there, rather than being part of the survreg object. As with many R functions, the p value is calculated in the print method for the object. In this case, it is a pretty

Re: [R] p-value from survreg(), library(survival)

2007-07-11 Thread hadley wickham
On 7/11/07, Marc Schwartz [EMAIL PROTECTED] wrote: Actually, in this case, looking at the code for: survival:::print.survreg would be better, as the p value is calculate there, rather than being part of the survreg object. As with many R functions, the p value is calculated in the print

Re: [R] CDF for pareto distribution

2007-07-11 Thread Vincent Goulet
Le 07-07-11 à 10:01, livia a écrit : Hi, thank you very much for your reply. The function pgpd() is from the package POT, and the 1.544 is the location parameter, 0.4477557 is the scale parameter and -0.50113 is the shape parameter, which can be both negtive or positive. The

Re: [R] 3D plot and interactive PDFs

2007-07-11 Thread vito muggeo
This does not answer exactly to your question, anyway..: If you are planning to use latex, the package movie15 allows to include media files in your document (to be processed via pdflatex) vito Bruno C. wrote: With version 8 of acrobat reader, it is now possible to have 3D in PDf

[R] R CMD SHLIB problem with -lg2c and make:

2007-07-11 Thread Richard Zur
Hello, I have 2 problems with R CMD SHLIB. The first problem didn't resolve, it just morphed into the new problem. I first have to admit that even though I've used R CMD SHLIB for a few years, I don't understand C compiling. The problem arose after our linux guy upgraded from Mandriva 2005

[R] installing, removing, upgrading, and downgrading packages

2007-07-11 Thread Christopher W. Ryan
I'm very new to R, trying to learn it. I started with R 2.4, but I have since upgraded to 2.5.0, on WindowsXP. I understand that 2.5.1 is now available. Last night in the course of things I loaded libraries coin and survival. I received warning messages that they had been built under version

Re: [R] elementary statistics with R (rkward?)

2007-07-11 Thread Christopher W. Ryan
As a fellow beginner, I also found Handbook of Statistical Analyses Using R, by Brian Everitt, to be a very useful book. There is an accompanying R package, HSAUR. Also Using R for Introductory Statistics, by John Verzani. There is an accompanying R package, UsingR. Christopher W. Ryan, MD

Re: [R] type III ANOVA for a nested linear model

2007-07-11 Thread Peter Dalgaard
Carsten Jaeger wrote: Hello Peter, thanks for your help. I'm quite sure that I specified the right model. Factor C is indeed nested within factor A. I think you were confused by the numbering of C (1..11), and it is easier to understand when I code it as you suggested (1,2,3 within each

Re: [R] elementary statistics with R (rkward?)

2007-07-11 Thread Chuck Cleland
Christopher W. Ryan wrote: As a fellow beginner, I also found Handbook of Statistical Analyses Using R, by Brian Everitt, to be a very useful book. There is an accompanying R package, HSAUR. Also Using R for Introductory Statistics, by John Verzani. There is an accompanying R package,

Re: [R] p-value from survreg(), library(survival)

2007-07-11 Thread Marc Schwartz
On Wed, 2007-07-11 at 16:41 +0200, hadley wickham wrote: On 7/11/07, Marc Schwartz [EMAIL PROTECTED] wrote: Actually, in this case, looking at the code for: survival:::print.survreg would be better, as the p value is calculate there, rather than being part of the survreg object. As

Re: [R] 3D plot and interactive PDFs

2007-07-11 Thread Bruno C\.
Thanks vito I was aware of movie15 but the point is how to get a U3D or VRLM file out of R :/ I don't know those two standards, nor I know the usual format for R 3D plots ... And unfortunately I am a bit in a rush so no way, right now, to do some reverse engineering about plot file format in

[R] tkfocus issue

2007-07-11 Thread Hao Liu
Dear All: I am stuck with this issue: I have a button on a TK window, once click it, it pops up a individual plot device: individual_plot - function() { tkconfigure(overlay.button, state=normal) options(locatorBell = FALSE) plotfuntype() trellis.focus(panel, 1,

Re: [R] exces return by mktcap decile for each year

2007-07-11 Thread jim holtman
here is one way of doing it using 'ave': dat - read.table(textConnection(mc yrret + 32902.233 01/01/1995 0.426 + 15793.691 01/01/1995 0.024 + 2375.868 01/01/1995 0.660 + 54586.558 01/01/1996 0.497 + 10674.900 01/01/1996 0.405 +859.656 01/01/1996 -0.033 +

Re: [R] Power calculation for the time series experiment

2007-07-11 Thread Greg Snow
The built in power functions are for the fairly straight forward situations. Yours does not appear to fit into any of those. You need to think through your problem a bit more before starting to think about power. What do you mean by effect size of 1.5 (is that 1.5 standard deviations? Or raw

Re: [R] survfit for interval censored data

2007-07-11 Thread Cody Hamilton
Sandra, As far as I am aware, you will have to use a parametric model (survreg) if your survival times are interval-censored. Regards, -Cody Cody Hamilton Edwards Lifesciences -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Sandra Ellis Sent:

[R] Drawing rectangles in multiple panels

2007-07-11 Thread Jonathan Williams
Hi folks, I'm having some trouble understanding the intricacies of panel functions. I wish to create three side-by-side graphs, each with different data-- so far, so good: I rbind() the data, add a column of subscripts as a conditioning variable, load up the lattice package, specify either a

Re: [R] Drawing rectangles in multiple panels

2007-07-11 Thread deepayan . sarkar
On 7/11/07, Jonathan Williams [EMAIL PROTECTED] wrote: Hi folks, I'm having some trouble understanding the intricacies of panel functions. I wish to create three side-by-side graphs, each with different data-- so far, so good: I rbind() the data, add a column of subscripts as a conditioning

Re: [R] error using lp function in linux

2007-07-11 Thread Byran Smucker
Perhaps I should clarify that when I use sample code from the lp help function, I can paste it and run it in my windows-based R with no problems. My question is why the same code won't work on a linux system. Thanks, Byran Smucker On Wed, 11 Jul 2007 12:00:09 +0200 r-help@stat.math.ethz.ch

Re: [R] Drawing rectangles in multiple panels

2007-07-11 Thread Jonathan Williams
Deepayan, Thanks for the clarification. The rectangles are completely external to the panel data, and correspond to 90% confidence intervals built from training data, to be overlaid on these graphs of the test data. - Jonathan At 10:04 AM 7/11/2007, you wrote: On 7/11/07, Jonathan Williams

Re: [R] Drawing rectangles in multiple panels

2007-07-11 Thread Bert Gunter
Deepayan et. al.: A question/comment: I have usually found that the subscripts argument is what I need when passing *external* information into the panel function, for example, when I wish to add results from a fit done external to the trellis call. Fits[subscripts] gives me the fits (or

[R] How to load permtest package

2007-07-11 Thread amna khan
Hi sir Sir I am not understanding how to load permtest on R console for windows? Regards -- AMINA SHAHZADI Department of Statistics GC University Lahore, Pakistan. Email: [EMAIL PROTECTED] [EMAIL PROTECTED] [EMAIL PROTECTED] [[alternative HTML version deleted]]

[R] aov() question

2007-07-11 Thread Leigh E Alexander
Hi all, So I think I have seen some similar questions to mine when I searched the archives, but have not seen any concrete answers and was wondering if any one could help. I have been trying to use R's aov() function to analyze my data. I have a 3 x 4 x 2 repeated measures design. All of

Re: [R] Drawing rectangles in multiple panels

2007-07-11 Thread deepayan . sarkar
On 7/11/07, Bert Gunter [EMAIL PROTECTED] wrote: Deepayan et. al.: A question/comment: I have usually found that the subscripts argument is what I need when passing *external* information into the panel function, for example, when I wish to add results from a fit done external to the trellis

Re: [R] Drawing rectangles in multiple panels

2007-07-11 Thread deepayan . sarkar
On 7/11/07, Jonathan Williams [EMAIL PROTECTED] wrote: Deepayan, Thanks for the clarification. The rectangles are completely external to the panel data, and correspond to 90% confidence intervals built from training data, to be overlaid on these graphs of the test data. Right. So if you

[R] system: Linux vs. Windows differences

2007-07-11 Thread Alberto Monteiro
[I tried to send this messages two days ago, but I guess I mistyped the To: address...] Why system is different in Linux and Windows? Both in R 2.4.1, but in Windows there is an option: system(something, wait = FALSE) while on Linux (Fedora Core 4), there is no such option? Alberto

Re: [R] Drawing rectangles in multiple panels

2007-07-11 Thread hadley wickham
A question/comment: I have usually found that the subscripts argument is what I need when passing *external* information into the panel function, for example, when I wish to add results from a fit done external to the trellis call. Fits[subscripts] gives me the fits (or whatever) I want to

Re: [R] system: Linux vs. Windows differences

2007-07-11 Thread Marc Schwartz
On Wed, 2007-07-11 at 16:16 -0200, Alberto Monteiro wrote: [I tried to send this messages two days ago, but I guess I mistyped the To: address...] Why system is different in Linux and Windows? Both in R 2.4.1, but in Windows there is an option: system(something, wait = FALSE) while

[R] Stepwise GLM selection by LRT?

2007-07-11 Thread Lutz Ph. Breitling
Dear List, having searched the help and archives, I have the impression that there is no automatic model selection procedure implemented in R that includes/excludes predictors in logistic regression models based on LRT P-values. Is that true, or is someone aware of an appropriate function

Re: [R] Stepwise GLM selection by LRT?

2007-07-11 Thread Ravi Varadhan
Check out the stepAIC function in MASS package. This is a nice tool, where you can actually implement any penalty even though the function's name has AIC in it because it is the default. Although this doesn't do an LRT test based variable selection, you can sort of approximate it by using a

[R] How to get weekly Co-Variance

2007-07-11 Thread Felipe Carrillo
Hi: I am trying to migrate from Systat to R but I am facing my first challenge. While I easily can get weekly co-variance for my data, I can't seem to acomplish this with R ( I can't figure out how is done) If interested in looking a sample of my data, please check the Excel attachment. In

[R] make error R-5.1 on sun solaris

2007-07-11 Thread Dan Powers
I hope this is enough information to determine the problem. Thanks in advance for any help. Configure goes ok (I think) ./configure --prefix=$HOME --without-iconv R is now configured for sparc-sun-solaris2.9 Source directory: . Installation directory:/home/dpowers C

Re: [R] system: Linux vs. Windows differences

2007-07-11 Thread Ted Harding
On 11-Jul-07 18:28:19, Marc Schwartz wrote: On Wed, 2007-07-11 at 16:16 -0200, Alberto Monteiro wrote: [I tried to send this messages two days ago, but I guess I mistyped the To: address...] Why system is different in Linux and Windows? Both in R 2.4.1, but in Windows there is an

Re: [R] make error R-5.1 on sun solaris

2007-07-11 Thread Peter Dalgaard
Dan Powers wrote: I hope this is enough information to determine the problem. Thanks in advance for any help. Configure goes ok (I think) ./configure --prefix=$HOME --without-iconv R is now configured for sparc-sun-solaris2.9 Source directory: . Installation directory:

Re: [R] system: Linux vs. Windows differences

2007-07-11 Thread Marc Schwartz
On Wed, 2007-07-11 at 22:32 +0100, [EMAIL PROTECTED] wrote: On 11-Jul-07 18:28:19, Marc Schwartz wrote: On Wed, 2007-07-11 at 16:16 -0200, Alberto Monteiro wrote: [I tried to send this messages two days ago, but I guess I mistyped the To: address...] Why system is different in Linux

Re: [R] Drawing rectangles in multiple panels

2007-07-11 Thread Deepayan Sarkar
On 7/11/07, hadley wickham [EMAIL PROTECTED] wrote: A question/comment: I have usually found that the subscripts argument is what I need when passing *external* information into the panel function, for example, when I wish to add results from a fit done external to the trellis call.

Re: [R] type III ANOVA for a nested linear model

2007-07-11 Thread Richard Rowe
Mark Difford wrote: Indeed! And, apropos of the expression, to be Ripleyed (and so be condemned to eating cookies for a long, long time), what about being Billasted? BestR, Mark. Simon Blomberg-4 wrote: I second the nomination! Simon. On Tue, 2007-07-10 at 10:02 -0600, Greg Snow

[R] rgdal memory error for a small map.

2007-07-11 Thread Milton Cezar Ribeiro
Dear all, I am trying to read a - almost for me - small map on R using rgdal package. This image have dimension of 1701x1401 and are stored on native ArcGis GRID format. When I see the files sizes, it is less then one megabyte. But unfortunately when I try read using readGDAL function it

Re: [R] Drawing rectangles in multiple panels

2007-07-11 Thread Stephen Tucker
Not that Trellis/lattice was entirely easy to learn at first. :) I've been playing around with ggplot2 and there is a plot()-like wrapper for building a quick plot [incidentally, called qplot()], but otherwise it's my understanding that you superpose elements (incrementally) to build up to the

Re: [R] Drawing rectangles in multiple panels

2007-07-11 Thread Stephen Tucker
In the Trellis approach, another way (I like) to deal with multiple pieces of external data sources is to 'attach' them to panel functions through lexical closures. For instance... rectInfo - list(matrix(runif(4), 2, 2), matrix(runif(4), 2, 2), matrix(runif(4), 2, 2))

Re: [R] Drawing rectangles in multiple panels

2007-07-11 Thread Deepayan Sarkar
On 7/11/07, Stephen Tucker [EMAIL PROTECTED] wrote: In the Trellis approach, another way (I like) to deal with multiple pieces of external data sources is to 'attach' them to panel functions through lexical closures. For instance... rectInfo - list(matrix(runif(4), 2, 2),

Re: [R] Drawing rectangles in multiple panels

2007-07-11 Thread Gabor Grothendieck
Your approach of using closures is cleaner than that given below but just for comparison in: http://tolstoy.newcastle.edu.au/R/devel/06/03/4476.html there is a createWrapper function which creates a new function based on the function passed as its first argument by using the components of the

Re: [R] Drawing rectangles in multiple panels

2007-07-11 Thread Stephen Tucker
Regarding this, I meant to imply that lattice was similarly flexible in the sense of handing multiple data sets [IMHO], in regards to other aspects of the 'grammar of graphics' I have no qualifications to justify comment. But the idea and intuitiveness of graph construction in ggplot2 is very

Re: [R] Drawing rectangles in multiple panels

2007-07-11 Thread hadley wickham
On 7/12/07, Deepayan Sarkar [EMAIL PROTECTED] wrote: On 7/11/07, hadley wickham [EMAIL PROTECTED] wrote: A question/comment: I have usually found that the subscripts argument is what I need when passing *external* information into the panel function, for example, when I wish to add

Re: [R] Drawing rectangles in multiple panels

2007-07-11 Thread hadley wickham
On 7/12/07, Stephen Tucker [EMAIL PROTECTED] wrote: Not that Trellis/lattice was entirely easy to learn at first. :) I've been playing around with ggplot2 and there is a plot()-like wrapper for building a quick plot [incidentally, called qplot()], but otherwise it's my understanding that you

Re: [R] Drawing rectangles in multiple panels

2007-07-11 Thread hadley wickham
On 7/12/07, Stephen Tucker [EMAIL PROTECTED] wrote: In the Trellis approach, another way (I like) to deal with multiple pieces of external data sources is to 'attach' them to panel functions through lexical closures. For instance... rectInfo - list(matrix(runif(4), 2, 2),

[R] Subsetting problem

2007-07-11 Thread Cressoni, Massimo (NIH/NHLBI) [F]
I need to perform the Exact Wilcoxon Mann-Whitney on a subset of my database. Assuming that IPPO is my data frame and IPPOBIS is the subset my variable still have 3 different levels and the function wilcox_test (package coin) does not accept it. I do not know how to overcome this problem. ippo -

Re: [R] make error R-5.1 on sun solaris

2007-07-11 Thread Prof Brian Ripley
You are asked for a C99 compiler and configure normally finds one: mixing declarations and code is valid C99. Unless something has been done with environment variables (e.g. in config.site) this gcc is very old. configure should come up with 'gcc -std=gnu99'. Re-ordering the code will help

[R] ggplot2 / reshape / Question on manipulating data

2007-07-11 Thread Pete Kazmier
I'm an R newbie but recently discovered the ggplot2 and reshape packages which seem incredibly useful and much easier to use for a beginner. Using the data from the IMDB, I'm trying to see how the average movie rating varies by year. Here is what my data looks like: ratings -

Re: [R] Subsetting problem

2007-07-11 Thread Prof Brian Ripley
You have three levels of factor 'ippo' and data on two. That is not a two-sample problem, as the error message says. Try IPPOBIS$ippo - IPPOBIS$ippo[drop=TRUE] And please use an informative subject line (see the posting guide). On Thu, 12 Jul 2007, Cressoni, Massimo (NIH/NHLBI) [F] wrote: I