Re: [R] Fwd: Timebased predictions in postgresql.

Hi, hints for your r-related day ahead predictions you find in the second part of Louis Torgo's very usefuel: DataMining with R http://www.liacc.up.pt/~ltorgo/ DataMiningWithR/PDF/DataMiningWithR.pdf And for your PL/SQL intention: http://www.omegahat.org/RSPostgres/ christian Am

Re: [R] How to use the whole dataset (including between events) in Cox model (time-varying covariates) ?

This is the consequence of the use of partial likelihood in the Cox model. You should read the literature on this point (for example, have you read Cox, 1972 and all its discussion, or Anderson, Borgan, Gill Keiding?). It is not an R question. You need to make more assumptions, such as a

[R] Robust R^2

Dear R-users, Do you know if there is any function in R available which calculates an R^2 for a robust MM-regression? With best regards, Carsten Colombier Dr. Carsten Colombier Economist Group of Economic Advisers Swiss Federal Finance Administration Bundesgasse 3 CH-3003 Bern phone +41 31

Re: [R] rgl snapshot command

Joern Kamradt wrote: Hi, I am using the rgl package for 3D display. Unfortunately, I am not able to get the snapshot command running. I tried the following: example(rgl.surface) rgl.sr data(volcano) rgl.sr y - 2 * volcano rgl.sr x - 10 * (1:nrow(y)) rgl.sr z - 10 * (1:ncol(y)) rgl.sr ylim -

Re: [R] truly object oriented programming in R

Dear Jason, Of course you can do almost eveything in R, but I have rarely seen someone prototyping something in Java to then implement it in R. KD-trees are for performance, and interpreted R is not really fast. So interfacing some foreign code might really be the better choice here. Opportunity

[R] how to plot an array with labels

How can i plot an array and instead of having on the x labels the indexes of the array I want to display an other String array of the same length __ [EMAIL PROTECTED] mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting

Re: [R] How to use the whole dataset (including between events) in Cox model (time-varying covariates) ?

you can always use parametric models and a full likelihood (but you may have to program them yourself). Prof Brian Ripley I started trying this but I could not make the counting process notation work on this. (Andersen, P.K. and Gill, R.D. (1982). Cox's regression model for counting

Re: [R] Robust R^2

On Fri, 13 Aug 2004 [EMAIL PROTECTED] wrote: Do you know if there is any function in R available which calculates an R^2 for a robust MM-regression? It's very easy to do, *but* why would you want a non-robust measure for a robust fit? [Just 1 - sum(residuals(fit)^2)/sum((y-mean(y))^2), for

[R] how to plot an array with labels

How can i plot an array and instead of having on the x labels the indexes of the array I want to display an other String array of the same length Do this: plot(myarray,xaxt=n,xlab=) axis(1,at=1:length(myarray),lab=my.vector.of.names) __ [EMAIL

[R] lapply problem

R-help, I wish to replace NULL elements(or missing) in the following list : z2 $cod mean sd 62.56190 12.65452 $haddock mean sd 36.61490 11.50365 $ling mean sd 86.17949 20.43587 $saithe meansd 50.275847 5.453606 $whiting NULL $norway pout

[R] simtest for Dunnett's test

Hi! I use simtest fonction of multcomp package to compile a Dunnett's test. I have 10 treatments and one control group, so i create a matrix with: m-matrix(0,10,11) m[1,1]--1 m[1,2]-1 m[2,1]--1 m[2,3]-1 m[3,1]--1 m[3,4]-1 m[4,1]--1 m[4,5]-1 m[5,1]--1 m[5,6]-1 m[6,1]--1 m[6,7]-1 m[7,1]--1

RE: [R] simtest for Dunnett's test

Before you do any of that, you should realize the fact that simtest does _not_ do Dunnett's test: It can use `Dunnett' _contrasts_ to do comparisons with a control, but the actual procedure is different. You are unlikely to get the same result as other packages that perform standard Dunnett's

RE : [R] simtest for Dunnett's test

However, on several tests, my results are correct: The t-value seems always to correspond to the q'value calculated by SigmaStat... So,What do I have to use to find the same results as SigmaStat (value and comparaison)? Laurent Houdusse Analyste Programmeur -Message d'origine- De

RE: [R] simtest for Dunnett's test

On Fri, 13 Aug 2004, Liaw, Andy wrote: Before you do any of that, you should realize the fact that simtest does _not_ do Dunnett's test: It can use `Dunnett' _contrasts_ to do comparisons with a control, but the actual procedure is different. You are unlikely to get the same result as

[R] simtest for Dunnett's test

Yes, but my control group is not always in the first colonne. So, I can have to a matrix like this: 1 0 0 -1 0 0 1 0 -1 0 0 0 1 -1 0 0 0 0 -1 0 Laurent Houdusse Analyste Programmeur -Message d'origine- De : Torsten Hothorn [mailto:[EMAIL PROTECTED] Envoyé : vendredi 13

RE: RE : [R] simtest for Dunnett's test

From: Laurent Houdusse However, on several tests, my results are correct: The t-value seems always to correspond to the q'value calculated by SigmaStat... So,What do I have to use to find the same results as SigmaStat (value and comparaison)? Seems like you don't have much of a

RE : RE : [R] simtest for Dunnett's test

Ok, I use simtest because i retrieve the same result But i don't know how to do the comparaison of P Laurent Houdusse Analyste Programmeur -Message d'origine- De : Liaw, Andy [mailto:[EMAIL PROTECTED] Envoyé : vendredi 13 août 2004 14:24 À : Laurent Houdusse; '[EMAIL PROTECTED]'

Re: [R] simtest for Dunnett's test

Yes, but my control group is not always in the first colonne. So, I can have to a matrix like this: 1 0 0 -1 0 0 1 0 -1 0 0 0 1 -1 0 0 0 0 -1 0 see the `base' argument to `contrMat'. Torsten Laurent Houdusse Analyste Programmeur -Message d'origine- De :

[R] Install R for Windows

Hi I'm a Brazilian student and I'd like to work with R. But, I had some problems to install R for windows with shortcut of the statisticals packages in tool bar. I'll be glad if you could help me to install it without miss the additionals packages I had install (more than 50!!) Thanks!

Re: [R] Install R for Windows

Clécio da Silva Ferreira wrote: Hi I'm a Brazilian student and I'd like to work with R. But, I had some problems to install R for windows with shortcut of the statisticals packages in tool bar. What do you mean with shortcut of the statisticals packages in tool bar. What are the statistical

Re: [R] How to use the whole dataset (including between events) in Cox model (time-varying covariates) ?

On Fri, 13 Aug 2004, Mayeul KAUFFMANN wrote: you can always use parametric models and a full likelihood (but you may have to program them yourself). Prof Brian Ripley I started trying this but I could not make the counting process notation work on this. (Andersen, P.K. and Gill, R.D.

[R] barplot and names.arg

R-help Is there any option to get closer the x-axis and names.arg from barplot? Thank you Luis Ridao Cruz Fiskirannsóknarstovan Nóatún 1 P.O. Box 3051 FR-110 Tórshavn Faroe Islands Phone: +298 353900 Phone(direct): +298 353912 Mobile: +298 580800 Fax:

[R] removing observations

Hello R list. I am sure this is an easy question, please forgive my ignorance of R. I have a panel data set (25 years) and I would like to break it up by year, and run a few tests. How is this done in R? In STATA it is quite simple (keep if year =1970). If someone could give me an

RE: [R] coxph question

Thanks Mayeul, I actually would like to test each variable individually and use those have low p-value to build a classifier (not in cox model). Therefore, I need to write a function to subset those low p-value variables, instead of putting them as covariates. Any ideas? Chris -Original

Re: [R] removing observations

Justin G. Gardner [EMAIL PROTECTED] writes: Hello R list. I am sure this is an easy question, please forgive my ignorance of R. I have a panel data set (25 years) and I would like to break it up by year, and run a few tests. How is this done in R? In STATA it is quite simple

Re: [R] truly object oriented programming in R Thank you all

My deep appreciation to everyone who responded to my question. I am digesting your proposals (a little distracted by my governor's dramatic resignation in New Jersey). It seems that I am able to implement the algorithm in R (with the framework kindly provided by Thomas). I will post the code when

Re: [R] barplot and names.arg

On Fri, 2004-08-13 at 09:22, Luis Rideau Cruz wrote: R-help Is there any option to get closer the x-axis and names.arg from barplot? Thank you Using mtext() you can do something like the following: data(VADeaths) # Now place labels closer to the x axis # set 'axisnames' to FALSE so the

Re: [R] How to use the whole dataset (including between events) in Cox model (time-varying covariates) ?

First, I do not know how to specify such a link function in R. Second, if I can specify such alink, I could use (in place of d*j), the smooth baseline estimated after doing a Cox regression. But I don't know how to fit (for instance) a piecewise constant baseline hazard with a Poison glm I found

Re: [R] removing observations

You probably want something like subset(dataset, year == 1970) and loop through the years. Or, perhaps even classier, would be datalist - split(dataset, as.factor(year)) results - lapply(datalist, myAnalysisFunction) -roger Justin G. Gardner wrote: Hello R list. I am sure this is an easy question,

R: [R] removing observations

by is perhaps what you looking for try to write ?by for help I have a panel data set (25 years) and I would like to break it up by year, and run a few tests. How is this done in R? In STATA it is quite simple (keep if year =1970). If someone could give me an example of how to do it for a

Re: [R] coxph question

Chris, I understood the following: you want to try every single covariate in a cox model with only one covariate, then take the best ones according to p-value. Assume your columns look like: stop status event x1 x2 x3 etc You want to add column 3 (x1), then 4, etc. I suggest a for() loop:

[R] Question from Newbie on PostScript and Multiple Plots

Hi, As I'm pretty new to R I hope this question isn't too basic. I am currently looping through my dataset and for each iteration am producing three separate plots. When I output these plots to the screen they are nicely grouped as three plots per page, however, when I try to send it to a

Re: [R] Question from Newbie on PostScript and Multiple Plots

On Fri, 2004-08-13 at 11:28, Johnson, Heather wrote: Hi, As I'm pretty new to R I hope this question isn't too basic. I am currently looping through my dataset and for each iteration am producing three separate plots. When I output these plots to the screen they are nicely grouped as

Re: [R] Question from Newbie on PostScript and Multiple Plots

Have you looked at ?par(mfrow) ? Johnson, Heather [EMAIL PROTECTED] Sent by: [EMAIL PROTECTED] 08/13/2004 12:28 PM To: '[EMAIL PROTECTED]' [EMAIL PROTECTED] cc: Subject:[R] Question from Newbie on PostScript and Multiple Plots Hi, As I'm pretty new

RE: [R] Question from Newbie on PostScript and Multiple Plots

Marc and Partha, Thank you both for your suggestions of using par(mfrow), that's exactly what I needed to do. In fact I already had it in the code (which I didn't write originally) and simply had to put my postscript statement in front of it. :) It works great now! Thanks for the quick

[R] cell counting in R

I have a bunch of jpg files that I need to count cells on. Some cells are green and some cells are red. Is there a function like locator that I could use to manually count the cells? I would like to be able to click on a cell, have it marked so I don't count it again and have a counter keep

Re: [R] cell counting in R

You might find it easier to do this type of thing with an open-source image analysis program such as ImageJ: http://rsb.info.nih.gov/ij/ Rob - Original Message - From: Doug Bourne [EMAIL PROTECTED] To: [EMAIL PROTECTED] Sent: Friday, August 13, 2004 1:10 PM Subject: [R] cell counting

[R] Tutorials on R

Hi !! Is there a good tutorial for the R language ? I really find it hard to find methods in R. Thanks and Cheers ../Murli __ [EMAIL PROTECTED] mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide!

Re: [R] Tutorials on R

An Introduction to R is pretty good (http://cran.r-project.org/doc/manuals/R-intro.pdf). And there are many others available at http://cran.r-project.org/other-docs.html -roger T. Murlidharan Nair wrote: Hi !! Is there a good tutorial for the R language ? I really find it hard to find methods

Re: [R] how to plot an array with labels

Hello to all, There are several functions to perform scatter diagrams with labels in the library 'ade4'. # library(ade4) bob - as.data.frame((array(1:5, c(4,2 row.names(bob) row.names(bob) - paste(row,1:4,sep=.) row.names(bob) par(mfrow=c(1,2)) s.label(bob,

Re: [R] Tutorials on R

There are quite a few helpful free contributions--under contributed go figure-- on the R site. Using R for Data Analysis and Graphics by John Maindonald (PDF [702kB], data sets and scripts are available at JM's homepage). I found this a big help( Thanks to John Maindonald) as well as trying to

[R] numerical accuracy, dumb question

I store an id as a big number, could this be a problem? Should I convert to at string when I use read.table(... example id's 1001001001001 1001001001002 ... 1002001002005 Bigest is probably 1011001001001 Ta, Dan. __ [EMAIL PROTECTED] mailing

RE: [R] Private methods

.class1 is an 'internal' method in the namespace of the package 'methods'. It can be accessed with the ::: operator. A good start for documentation may be Vol 3/1 of R-News. methods:::.class1 function (x) class(x)[[1]] environment: namespace:methods --Matt -Original Message- From:

RE: [R] numerical accuracy, dumb question

If I'm not mistaken, numerics are read in as doubles, so that shouldn't be a problem. However, I'd try using factor or character. Andy From: Dan Bolser I store an id as a big number, could this be a problem? Should I convert to at string when I use read.table(... example id's

Re: [R] Private methods

On Sat, 14 Aug 2004, Matthew Walker wrote: I have read (in a very recent email to r-help) that R's methods package does not support private methods. However, I also looked at the source for the function is: is function (object, class2) { cl - .class1(object) if

[R] SPSS, social science majors, and R

Hi Roland and folks: Roland's and Susanna's comments were very interesting. I think we need to keep in mind that unlike SPSS, R is more than a statistical package. It's a tool. SPSS does indeed make life easy, particularly for data preprocessing (or data cleaning), something that is perhaps

Re: [R] Giving a first good impression of R to Social Scientists

I think I agree with Prof. Lumley: For me, R would likely be the platform of choice for most new statistical algorithm development. It has a steep learning curve, and I would therefore not recommend it for anyone who would not likely be involved in this kind of activity (unless my

RE: [R] numerical accuracy, dumb question

Part of that decision may depend upon how big the dataset is and what is intended to be done with the ID's: object.size(1011001001001) [1] 36 object.size(1011001001001) [1] 52 object.size(factor(1011001001001)) [1] 244 They will by default, as Andy indicates, be read and stored as doubles.

[R] RE: Giving a first good impression of R to Social Scientists

Dear Roland, Have you looked at Zelig (http://gking.harvard.edu/zelig)? Several professors in my department are going to use it to teach R to political science undergraduates and graduate students this fall. We just presented it at the Political Methodology meeting, and will present it again at