Re: [R] thanks

2005-01-12 Thread Lefebure Tristan
see man R example from a shell: echo -e pdf(file=\test.pdf\)\nplot(1:10,11:20)\ndev.off(dev.cur())\ncmd.R R -s cmd.R (write a file of command for R, and than feed R with it) On Tuesday 11 January 2005 15:59, Cserháti Mátyás wrote: Dear all, Thanks to those 3 people who sent me answers to my

[R] (no subject)

2005-01-12 Thread Bigby
Hello, My name is Graham, I am an engineering student and my lecturer wishes us to get a numerical summary of some data. He said use the command numerical.summary(Data), which didnt work, he suggested we try library(s20x) first, which came up with an error on my console. I have version 2.0.1 of R

RE: [R] global objects not overwritten within function

2005-01-12 Thread McGehee, Robert
I would suggest reading the posting guide, (http://www.r-project.org/posting-guide.html) and give a reproducible example, with the error message that you received. As is, I have no idea what you are doing here, and certainly cannot run this code. You use ... as an argument to your functions (why I

Re: [R] Graphical table in R

2005-01-12 Thread Marc Schwartz
On Tue, 2005-01-11 at 14:59 +, Dan Bolser wrote: On 10 Jan 2005, Peter Dalgaard wrote: Dan Bolser [EMAIL PROTECTED] writes: Cheers. This is really me just being lazy (as usual). The latex function in Hmisc allows me to make a .ps file then grab a screen shot of that ps and make a

[R] RODBC package -- sqlQuery(channel,.....,nullstring=0) still gives NA's

2005-01-12 Thread Luis Rideau Cruz
R-help, I'm using the RODBC package to retrieve data froma ODBC database which contain NA's. By using the argument nullstring = 0 in sqlQuery() I expect to coerce them to numeric but still get NA's in my select. I'm running on Windows XP version _ platform

[R] Meeker's SPLIDA Reliability in R

2005-01-12 Thread Cunningham, Colin A
Hi, Would anyone be aware of an R package implementing the functionality found in Meeker's SPLIDA software written for S-Plus? I don't know if anyone has tried to port the s/w to R directly, or if equivalent functions are available within another package. Thanks in Advance. - Colin Colin

Re: [R] thanks

2005-01-12 Thread Jan T. Kim
On Tue, Jan 11, 2005 at 03:59:58PM +0100, Cserh?ti M?ty?s wrote: I would like to run an R script from the Linux prompt. Is there any way possible to do this? The reason is, the calculation that I'm doing takes a few hours, and I would like to automatize it. Or does it mean that I have to

Re: [R] global objects not overwritten within function

2005-01-12 Thread Prof Brian Ripley
On Tue, 11 Jan 2005, bogdan romocea wrote: Dear useRs, I have a function that creates several global objects with assign(obj,obj,.GlobalEnv), and which I need to run iteratively in another function. The code is similar to f - function(...) { assign(obj,obj,.GlobalEnv) } fct - function(...) { for

[R] help on integrate function

2005-01-12 Thread Francisca xuan
here is a function I wrote cdfest=function(t,lambda,delta,x,y){ a1=mean(x t) a2=mean(x t-delta) a3=mean(y1 t) s=((1-lambda)*a1+lambda*a2-a3)^2 s } when I try to integrate over t, I got this message: integrate(cdfest,0,4,lambda=0.3,delta=1,x=x,y=y1) Error in integrate(cdfest, 0, 4,

[R] useR 2005 ?

2005-01-12 Thread Rau, Roland
Dear R-Help-List, are there any plans to organize a useR conference in 2005? Best, Roland + This mail has been sent through the MPI for Demographic Rese...{{dropped}} __ R-help@stat.math.ethz.ch mailing list

Re: [R] thanks

2005-01-12 Thread Uwe Ligges
Cserhti Mtys wrote: Dear all, Thanks to those 3 people who sent me answers to my question. Got the problem solved. Great! Now, another question of mine is: I would like to run an R script from the Linux prompt. Is there any way possible to do this? The reason is, the calculation that I'm doing

[R] getting variable names from formula

2005-01-12 Thread Daniel Almirall
R-list, 1. Given a formula (f) w variables referencing some data set (dat), is there any easier/faster way than this to get the names (in character form) of the variables on the RHS of '~' ? dat - data.frame(x1 = x1 - rnorm(100,0,1), x2 = x2 - rnorm(100,0,1), y = x1 + x2 + rnorm(100,0,1)) f

Re: [R] Graphical table in R

2005-01-12 Thread Dan Bolser
On Tue, 11 Jan 2005, Marc Schwartz wrote: On Tue, 2005-01-11 at 14:59 +, Dan Bolser wrote: On 10 Jan 2005, Peter Dalgaard wrote: Dan Bolser [EMAIL PROTECTED] writes: Cheers. This is really me just being lazy (as usual). The latex function in Hmisc allows me to make a .ps file then

[R] Standard error for the area under a smoothed ROC curve?

2005-01-12 Thread Dan Bolser
Hello, I am making some use of ROC curve analysis. I find much help on the mailing list, and I have used the Area Under the Curve (AUC) functions from the ROC function in the bioconductor project... http://www.bioconductor.org/repository/release1.5/package/Source/ ROC_1.0.13.tar.gz

Re: [R] help on integrate function

2005-01-12 Thread Prof Brian Ripley
On Tue, 11 Jan 2005, Francisca xuan wrote: here is a function I wrote cdfest=function(t,lambda,delta,x,y){ a1=mean(x t) a2=mean(x t-delta) a3=mean(y1 t) s=((1-lambda)*a1+lambda*a2-a3)^2 s } when I try to integrate over t, I got this message:

[R] Breslow Day Test

2005-01-12 Thread Palos, Judit
Breslow-Day test A statistical test for the homogeneity of odds ratios. Homogeneity In javascript:void(0); systematic reviews homogeneity refers to the degree to which the results of studies included in a review are similar. Clinical homogeneity means that, in studies included in a review,

Re: [R] getting variable names from formula

2005-01-12 Thread Douglas Bates
Daniel Almirall wrote: R-list, 1. Given a formula (f) w variables referencing some data set (dat), is there any easier/faster way than this to get the names (in character form) of the variables on the RHS of '~' ? dat - data.frame(x1 = x1 - rnorm(100,0,1), x2 = x2 - rnorm(100,0,1), y = x1 + x2

[R] transcan() from Hmisc package for imputing data

2005-01-12 Thread avneet singh
Hello: I have been trying to impute missing values of a data frame which has both numerical and categorical values using the function transcan() with little luck. Would you be able to give me a simple example where a data frame is fed to transcan and it spits out a new data frame with the NA

Re: [R] getting variable names from formula

2005-01-12 Thread Dimitris Rizopoulos
maybe something like: f - y ~ x1 + x2 attr(terms(f), term.labels) but this wan't work if you have a more complex formula (e.g., f - y ~ x1*x2 + I(x1^2)) and you want only c(x1, x2). I hope it helps. Best, Dimitris Dimitris Rizopoulos Ph.D. Student Biostatistical Centre School of Public

Re: [R] (no subject)

2005-01-12 Thread Kevin Wang
Hi, On Wed, 12 Jan 2005, Bigby wrote: Hello, numerical.summary(Data), which didnt work, he suggested we try library(s20x) first, which came up with an error on my console. I have version 2.0.1 of R library(s20x) is a package written by the Department of Statistics at the University of

Re: [R] Destructor for S4 objects?

2005-01-12 Thread Adam Lyon
Hi Robert, It looks like there is no way to explicitly make an S4 object call a function when it is garbage collected unless you resort to tricks with reg.finalizer. It turns out that Prof. Ripley's reply (thanks!!) had enough hints in it that I was able to get the effect I wanted by using R's

Re: [R] RODBC package -- sqlQuery(channel,.....,nullstring=0) stillgives NA's

2005-01-12 Thread Prof Brian Ripley
PLEASE do read the help page, which says nullstring: character string to be used when reading 'SQL_NULL_DATA' character items from the database. ^^^ so this does not apply to numeric items. You can of course easily change numeric NAs to 0s, if you want to. On Wed,

Re: [R] RODBC package -- sqlQuery(channel,.....,nullstring=0)stillgives NA's

2005-01-12 Thread Luis Rideau Cruz
(1) I do read the posting guide (the fact that I missread o missunderstood something does not imply not reading) (2) I could change NAs to 0 (I know) but I have previously (older versions of R and SQL*Plus) used the same select with the right output (namely with 0s). (3) AFAIK strange is not a

Re: [R] Breslow Day Test

2005-01-12 Thread Tobias Verbeke
On Tue, 11 Jan 2005 10:45:48 -0500 Palos, Judit [EMAIL PROTECTED] wrote: Breslow-Day test A statistical test for the homogeneity of odds ratios. [..some definitions..] Your message was not particularly clear, but if you were looking for R code to do a Breslow-Day test, Google found this

Re: [R] Standard error for the area under a smoothed ROC curve?

2005-01-12 Thread Frank E Harrell Jr
Dan Bolser wrote: Hello, I am making some use of ROC curve analysis. I find much help on the mailing list, and I have used the Area Under the Curve (AUC) functions from the ROC function in the bioconductor project... http://www.bioconductor.org/repository/release1.5/package/Source/

Re: [R] transcan() from Hmisc package for imputing data

2005-01-12 Thread Frank E Harrell Jr
avneet singh wrote: Hello: I have been trying to impute missing values of a data frame which has both numerical and categorical values using the function transcan() with little luck. Would you be able to give me a simple example where a data frame is fed to transcan and it spits out a new data

[R] [survey] R for Reporting - the R Output MAnager (ROMA) project

2005-01-12 Thread Eric Lecoutre
Hi R UseRs, I am interested in providing Reporting abilities to R and have initiated a project called R Output MAnager (ROMA). My starting point was my R2HTML package which provides (rough) HTML exportations. I began with trying to mimic it for LaTeX but fastly did realize that it was a bad

Re: [R] thanks

2005-01-12 Thread Jan T. Kim
On Tue, Jan 11, 2005 at 04:24:11PM +0100, Lefebure Tristan wrote: example from a shell: echo -e pdf(file=\test.pdf\)\nplot(1:10,11:20)\ndev.off(dev.cur())\ncmd.R R -s cmd.R (write a file of command for R, and than feed R with it) This may be on the verge of becoming offtopic, but let me

Re: [R] CUSUM SQUARED structural breaks approach?

2005-01-12 Thread Achim Zeileis
On Tue, 11 Jan 2005 19:33:41 + Rick Ram wrote: Groundwork for the choice of break method in my specific application has already been done - otherwise I would need to rework the wheel (make a horribly detailed comparison of performance of break approaches in context of modelling post

Re: [R] [survey] R for Reporting - the R Output MAnager (ROMA) project

2005-01-12 Thread A.J. Rossini
Your example is sequential, ignoring the tree-like structure of most documents. Why not via a DOM or similar XML-ish structure? While I'd never advocate general purpose XML as a user format, as you note, that is what XSLT is for, and using XML as an electronic internal document representation

Re: [R] useR 2005 ?

2005-01-12 Thread Marc Schwartz
On Tue, 2005-01-11 at 17:39 +0100, Rau, Roland wrote: Dear R-Help-List, are there any plans to organize a useR conference in 2005? Best, Roland As I understand it, no. The next one will be in 2006, so it will be every other year, interleaved with the DSC meeting the odd years.

[R] (no subject)

2005-01-12 Thread nicolas . deig
hi, I am trying to grow a classification tree on some data, but I have a little problem. In order to do so I have to use a function like tree in R and on the internet help(tree) I get the following: The left-hand-side (response) should be either a numerical vector when a regression tree will

Re: [R] What is a factor? [was :(no subject)]

2005-01-12 Thread Marc Schwartz
On Wed, 2005-01-12 at 15:17 +0100, [EMAIL PROTECTED] wrote: hi, I am trying to grow a classification tree on some data, but I have a little problem. In order to do so I have to use a function like tree in R and on the internet help(tree) I get the following: The left-hand-side

[R] (no subject)

2005-01-12 Thread saverio vicario
Dear help desk and R community, I have a problem on how R2.0 handle the RAM, maybe a bug In fact I used R1.9 since january 2004 with large data set using a macosx G5 with 1G of ram without problem . Then I passed to 2.0 and I found myself short in ram using virtual memory. I tried to use

Re: [R] Breslow Day Test

2005-01-12 Thread Thomas Lumley
On Wed, 12 Jan 2005, Tobias Verbeke wrote: On Tue, 11 Jan 2005 10:45:48 -0500 Palos, Judit [EMAIL PROTECTED] wrote: Breslow-Day test A statistical test for the homogeneity of odds ratios. [..some definitions..] Your message was not particularly clear, but if you were looking for R code to do a

Re: [R] (no subject)

2005-01-12 Thread Anne
I think you will find all the doc in the help files ?factor() gets The function 'factor' is used to encode a vector as a factor (the terms 'category' and 'enumerated type' are also used for factors). If 'ordered' is 'TRUE', the factor levels are assumed to be ordered. For

[R] changing langage

2005-01-12 Thread Kurt Sys
Hi all, I've got a small, practical question, which untill now I couldn't solve (otherwhise I wouldn't mail it, right?) First of all, I'm talking about R 2.0.1 on a winxp system (using the default graphical interface being 'Rgui'). When I make plots, using dates on the x-axis, it puts the

[R] defining lower part of distribution with covariate

2005-01-12 Thread Troels Ring
I try again - perhaps it is analysis of covariance with treatment (thio,ultiva) as two categories and antime as covariate. On the basis of such a model, is then the probability of GCS = 12 larger with thio treatment ? Dear friends, forgive me a simple question, possibly related to quantreg but

Re: [R] Kolmogorov-Smirnof test for lognormal distribution with estimated parameters

2005-01-12 Thread Christoph Buser
Hi Kwabena I did once a simulation, generating normal distributed values (500 values) and calculating a KS test with estimated parameters. For 1 times repeating this test I got about 1 significant tests (on a level alpha=0.05 I'm expecting about 500 significant tests by chance) So I think if

[R] model.response error

2005-01-12 Thread Bang
When I installed R 2.0.1 (replacing 1.9.0) for Windows, a code using model.response began acting up. Here are the first several lines of a code I had been tweaking for a spatial model (the code is mostly that of Roger Bivand--I am adapting it to a slightly different data structure and the problem

Re: [R] Kolmogorov-Smirnof test for lognormal distribution with estimated parameters

2005-01-12 Thread Frank E Harrell Jr
Christoph Buser wrote: Hi Kwabena I did once a simulation, generating normal distributed values (500 values) and calculating a KS test with estimated parameters. For 1 times repeating this test I got about 1 significant tests (on a level alpha=0.05 I'm expecting about 500 significant tests by

Re: [R] changing langage

2005-01-12 Thread Peter Dalgaard
Kurt Sys [EMAIL PROTECTED] writes: Hi all, I've got a small, practical question, which untill now I couldn't solve (otherwhise I wouldn't mail it, right?) First of all, I'm talking about R 2.0.1 on a winxp system (using the default graphical interface being 'Rgui'). When I make plots,

Re: [R] Kolmogorov-Smirnof test for lognormal distribution with estimated parameters

2005-01-12 Thread Christian Hennig
For the KS-test of normality with estimated parameters see ?lillie.test in package nortest. Best, Christian On Wed, 12 Jan 2005, Christoph Buser wrote: Hi Kwabena I did once a simulation, generating normal distributed values (500 values) and calculating a KS test with estimated

Re: [R] changing langage

2005-01-12 Thread Gabor Grothendieck
Kurt Sys kurt.sys at pandora.be writes: : : Hi all, : : I've got a small, practical question, which untill now I couldn't solve : (otherwhise I wouldn't mail it, right?) First of all, I'm talking about : R 2.0.1 on a winxp system (using the default graphical interface being : 'Rgui'). : When

Re: [R] changing langage

2005-01-12 Thread Prof Brian Ripley
It uses the language set by LC_TIME: see ?Sys.setlocale and ?format.Date which references it. On Wed, 12 Jan 2005, Kurt Sys wrote: Hi all, I've got a small, practical question, which untill now I couldn't solve (otherwhise I wouldn't mail it, right?) First of all, I'm talking about R 2.0.1 on a

Re: [R] model.response error

2005-01-12 Thread Roger Bivand
On Wed, 12 Jan 2005, Bang wrote: When I installed R 2.0.1 (replacing 1.9.0) for Windows, a code using model.response began acting up. Here are the first several lines of a code I had been tweaking for a spatial model (the code is mostly that of Roger Bivand--I am adapting it to a slightly

RE: [R] defining lower part of distribution with covariate

2005-01-12 Thread Berton Gunter
Troels: It would be best if you discussed this with a local statistician to make sure that the data and analysis are properly addressing the scientific issues. Perhaps that is why no one replied to your previous post. Also, this is primarily a **statistical** issue, not really an ** R-issue **.

Re: [R] transfer function models

2005-01-12 Thread Paul Gilbert
I don't know what SAS does, but transfer functions are essentially MA/AR from an ARMA model, so you should be able to get what you want from the various ARMA estimation tools in R. Paul Gilbert Samuel Kemp (Comp) wrote: Hi, Does anyone know of a function in R that can estimate the parameters of

Re: [R] Standard error for the area under a smoothed ROC curve?

2005-01-12 Thread Frank E Harrell Jr
Dan Bolser wrote: On Wed, 12 Jan 2005, Frank E Harrell Jr wrote: Dan Bolser wrote: Hello, I am making some use of ROC curve analysis. I find much help on the mailing list, and I have used the Area Under the Curve (AUC) functions from the ROC function in the bioconductor project...

Re: [R] global objects not overwritten within function

2005-01-12 Thread bogdan romocea
Apparently the message below wasn't posted on R-help, so I'm sending it again. Sorry if you received it twice. --- bogdan romocea [EMAIL PROTECTED] wrote: Date: Tue, 11 Jan 2005 17:31:42 -0800 (PST) From: bogdan romocea [EMAIL PROTECTED] Subject: Re: [R] global objects not overwritten within

Re: [R] changing langage [SOLVED]

2005-01-12 Thread Kurt Sys
To all that replied, thanks... I have a clue where I can change the settings. tnx, Kurt Sys - Oorspronkelijk bericht - Van : Gabor Grothendieck [mailto:[EMAIL PROTECTED] Verzonden : woensdag , januari 12, 2005 05:26 PM Aan : r-help@stat.math.ethz.ch Onderwerp : Re: [R] changing langage

[R] gbm

2005-01-12 Thread Weiwei Shi
Hi, there: I am wondering if I can find some detailed explanation on gbm or explanation on examples of gbm. thanks, Ed __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide!

[R] Finding seasonal peaks in a time series....

2005-01-12 Thread Dr Carbon
I have a seasonal time series. I want to calculate the annual mean value of the time series at its peak (say the mean of the three values before the peak, the peak, and the three values after the peak). The peak of the time series might change cycle slightly from year to year. # E.g., nPts -

Re: [R] gbm

2005-01-12 Thread Spencer Graves
I just got 25 hits from www.r-project.org - search - R site search. Might one or more of these help you? If they don't solve your problem, I suggest you try the posting guide! http://www.R-project.org/posting-guide.html;. If that still doesn't solve your problem, it should help you

Re: [R] gbm

2005-01-12 Thread Peter Dalgaard
Weiwei Shi [EMAIL PROTECTED] writes: Hi, there: I am wondering if I can find some detailed explanation on gbm or explanation on examples of gbm. What is gbm? Green Belt Movement? Georgie Boy Manufacturing? I'm serious! Well, only sort of, but try Google on gbm and you'll find those two

Re: [R] gbm

2005-01-12 Thread Daniel Almirall
You can also check out: http://www.i-pensieri.com/gregr/gbm.shtml There are reference papers on there, too. HTH, Danny On Wed, 12 Jan 2005, Peter Dalgaard wrote: Weiwei Shi [EMAIL PROTECTED] writes: Hi, there: I am wondering if I can find some detailed explanation on gbm or

[R] Off Topic: Statistical philosophy rant

2005-01-12 Thread Berton Gunter
R-Listers. The following is a rant originally sent privately to Frank Harrell in response to remarks he made on this list. The ideas are not new or original, but he suggested I share it with the list, as he felt that it might be of wider interest, nonetheless. I have real doubts about this, and I

[R] gbm

2005-01-12 Thread Weiwei Shi
Hi, there: Thanks a lot for all people' prompt replies. In detail, I am facing a huge amount of data: over 10,000 and 400 vars. This project is very challenging and interesting to me. I tried rpart which gives me some promising results but not good enough. So I am trying randomForest and gbm now.

Re: [R] Off Topic: Statistical philosophy rant

2005-01-12 Thread Dan Bolser
On Wed, 12 Jan 2005, Berton Gunter wrote: R-Listers. The following is a rant originally sent privately to Frank Harrell in response to remarks he made on this list. The ideas are not new or original, but he suggested I share it with the list, as he felt that it might be of wider interest,

RE: [R] gbm

2005-01-12 Thread Liaw, Andy
From: Weiwei Shi Hi, there: Thanks a lot for all people' prompt replies. In detail, I am facing a huge amount of data: over 10,000 and 400 vars. This project is very challenging and interesting to me. I tried rpart which gives me some promising results but not good enough. So I am

RE: [R] Changing the ranges for the axis in image()

2005-01-12 Thread Mulholland, Tom
There's something that you're not telling me. If you want something other than your data, why use your data. If you have another set of data that you wish to overlay on the image then you are going to have to scale one of the data sources to match the other. I'm not sure where your problem is

[R] multivariate diagnostics

2005-01-12 Thread Yulei He
Hi, there. I have two questions about the diagnostics in multivarite statistics. 1. Is there any diagnostics tool to check if a multivariate sample is from multivariate normal distribution? If there is one, is there any function doing it in R? 2. Is there any function of testing if two

RE: [R] Off Topic: Statistical philosophy rant

2005-01-12 Thread Mulholland, Tom
I have often noted that statistics can't prove a damn thing, but they can be really useful in disproving something. Having spent most of 80s and half of the 90s with the Australian Bureau of Statistics to find out how you collect these numbers, I am disconcerted at the apparent disregard for

RE: [R] Finding seasonal peaks in a time series....

2005-01-12 Thread Mulholland, Tom
You might find breakpoints in strucchange helpful Tom -Original Message- From: Dr Carbon [mailto:[EMAIL PROTECTED] Sent: Thursday, 13 January 2005 6:19 AM To: r-help@stat.math.ethz.ch Subject: [R] Finding seasonal peaks in a time series I have a seasonal time series. I

RE: [R] Please unsubscribe me from you list

2005-01-12 Thread Mulholland, Tom
You stand more chance if you do it yourself https://stat.ethz.ch/mailman/listinfo/r-help -Original Message- From: Kevin Ita [mailto:[EMAIL PROTECTED] Sent: Thursday, 13 January 2005 2:25 AM To: R-help@stat.math.ethz.ch Subject: [R] Please unsubscribe me from you list Please

RE: [R] Finding seasonal peaks in a time series....

2005-01-12 Thread Mulholland, Tom
Sorry I didn't read the question properly. Please disregard, my mind was elsewhere. Tom -Original Message- From: Mulholland, Tom Sent: Thursday, 13 January 2005 10:52 AM To: Dr Carbon; r-help@stat.math.ethz.ch Subject: RE: [R] Finding seasonal peaks in a time series You

[R] [R-pkgs] New package: MatchIt

2005-01-12 Thread Elizabeth Stuart
We would like to announce the release of our software MatchIt, now available on CRAN. MatchIt implements a variety of matching methods for causal inference. Abstract: MatchIt implements the suggestions of Ho, Imai, King, and Stuart (2004) for improving parametric statistical models by