[R] Extract data from edit chart

2005-09-20 Thread Chun-Ying Lee
Dear R users: I wonder if it is possible to use data.frame and edit to show the form like: parameter upper lower A B C when I type the command edit(abc), assume the data.frame named abc. And then I want to extract the number form it, like abc[1,2] which mean the upper level of

Re: [R] Extract data from edit chart

2005-09-20 Thread Gabor Grothendieck
By specifying space like that you are specifying that the columns be factors. If you want them to be character columns use abc - data.frame(Parameter = 1:10, lo = I(), hi = I()) and if you want them to be numeric columns specify something that is numeric like 0, NaN, Inf, -Inf or something

Re: [R] waveform filtering

2005-09-20 Thread Uwe Ligges
Huntsinger, Reid wrote: You should probably have a look at the sound packages for R, tuneR and sound, I believe, on http://cran.r-project.org. Applying a filter can be done with filter(), but you need to come up with filter coefficients. High-pass and low-pass have simple descriptions in

[R] Problem installing packages on Windows

2005-09-20 Thread John Field
I've recently had a series of similar errors trying to install packages to R 2.1.1 running under Windows XP. utils:::menuInstallLocal() package 'randomForest' successfully unpacked and MD5 sums checked updating HTML package descriptions Warning message: unable to move temporary installation

[R] predicting residual expected survival times

2005-09-20 Thread Anne
DeaR-Helpers Is there an implemented method to predict residual expected survival times for parametric/Cox PH models ? (I have modelled my data using the survival library) I would like to predict for a given subject (with a given profile ) having survived up to time Ts the expected residual

[R] scale function within a for-loop?

2005-09-20 Thread Frank Schmid
I have a data matrix containing around 180 variables and more than 470 observations for each and no missing values. Within a for-loop, the first step of calculations is to standardize each column, such that the mean of each column is zero and the sd is one. The for-loop starts with a subset of

Re: [R] Xgird, R, parallel computing

2005-09-20 Thread Christophe Pouzat
Hi, Can start by checking out these two addresses: http://adm.wustl.edu/rcluster.php http://adm.wustl.edu/xgrid.php Christophe. Mahdi Osman wrote: Hi, list, Sorry if I am bothering you. I am interested in using xgrid with R for distributed computing. I am using MacOS X and my R 2.1. 1 is

Re: [R] scale function within a for-loop?

2005-09-20 Thread Petr Pikal
Hi I maybe mistaken but scale(your.matrix) gives you matrix scaled in the way you want. apply(scale(as.matrix(kalcin[,3:7])), na.rm=T,2,sd) vodofe stroz l a b 1 1 1 1 1 apply(scale(as.matrix(kalcin[,3:7])), na.rm=T,2,mean) vodofe

Re: [R] Are least-squares means useful or appropriate?

2005-09-20 Thread Felipe
-BEGIN PGP SIGNED MESSAGE- Hash: SHA1 Hi. My question was just theoric. I was wondering if someone who were using SAS and R could give me their opinion on the topic. I was trying to use least-squares means for comparison in R, but then I found some indications against them, and I wanted

Re: [R] column-binary data

2005-09-20 Thread David Barron
Thanks for the replies. That's not quite what I meant. These data are multipunched to allow more than one variable to be coded in the same column. For example, the first 7 columns of the first card of the data I'm trying to read contain the following: Column Rows Description

Re: [R] Extended Hypergeometric Distribution

2005-09-20 Thread Martin Maechler
Narcyz == Narcyz Ghinea [EMAIL PROTECTED] on Mon, 19 Sep 2005 12:38:27 +1000 writes: Narcyz Dear R Users, Narcyz There exists a non-central hypergeometric Narcyz distribution function in the (MCMCpack) package, and Narcyz a hypergeometric distribution function in the

Re: [R] Teaching R - In front of the computer?

2005-09-20 Thread Christian Hennig
Dear Roland, I wrote an eight-pages manuscript with basic commands and instructions to use the help system and gave the students series of many small exercises. For further concepts such as matrix computations, index manipulations, coercion, I gave short presentations (15 minutes or so), again

Re: [R] Possible bug in lmer nested analysis with factors

2005-09-20 Thread Yan Wong
On 18 Sep 2005, at 16:04, Douglas Bates wrote: You are correct that good documentation of the capabilities of lmer does not currently exist. lmer is still under active development and documentation is spread in several places. The vignette in the mlmRev package explores some of the

Re: [R] Possible bug in lmer nested analysis with factors

2005-09-20 Thread Yan Wong
On 18 Sep 2005, at 16:27, Douglas Bates wrote: I have a couple of other comments. You can write the nested grouping factors as Sundar did without having to explicitly evaluate the interaction term and drop unused levels. The lmer function, like most modeling functions in R, uses the

Re: [R] Teaching R - In front of the computer?

2005-09-20 Thread Rau, Roland
Dear all, thanks for telling me your experiences how you proceed teaching R. I am currently giving a five week course teaching R to 16 graduate-level students consisting of 12 sessions à 1.5 hours. Two weeks have passed during that course and I was just questioning myself whether my style of

[R] Running glm in batch and exporting results (AIC) to HTML

2005-09-20 Thread Jérôme Lemaître
Dear all, I'm doing univariate Poisson regressions using the glm and glm.fit functions. I have 5 independent datasets and each dataset, has one response variable and more than 20 factors to test. Currently, I run one regression at a time and manually take notes of the results in excel to have a

Re: [R] Running glm in batch and exporting results (AIC) to HTML

2005-09-20 Thread vincent
see ?write.table hih __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html

Re: [R] waveform filtering

2005-09-20 Thread tom wright
On Mon, 2005-19-09 at 10:36 -0400, tom wright wrote: I'm not an engineer so I hope I'm using the correct terminology here. I have a recorded waveform that I want to apply low and high pass filters too, are tehre already R functions existing to do this or am I going to have to program my own?

[R] Change the mirror

2005-09-20 Thread Ronaldo Reis-Jr.
Hi, Please, change the brazilian mirror http://www.termix.ufv.br/CRAN/ to http://www.insecta.ufv.br/CRAN/ in R homepage. Thanks ROnaldo -- Se dois homens no mesmo trabalho concordam o tempo todo, um deles é demais. Se discordam sem parar, então os dois são dispensáveis --Darryl F. Zanuck --

Re: [R] Change the mirror

2005-09-20 Thread Duncan Murdoch
On 9/20/2005 9:24 AM, Ronaldo Reis-Jr. [EMAIL PROTECTED] wrote: Hi, Please, change the brazilian mirror http://www.termix.ufv.br/CRAN/ to http://www.insecta.ufv.br/CRAN/ in R homepage. Requests like this should go to CRAN (to whom I've cc'd this). The admins there might not see it in

[R] Interpretation of csplit from rpart.object

2005-09-20 Thread jmoreira
Dear members of R-list, I need to reproduce the rules of a decision tree. For that I need to use the csplit information from the rpart.object. But I cannot uderstand the information because from my example I get: rpart.tree$csplit [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,]133

Re: [R] predicting residual expected survival times

2005-09-20 Thread Dieter Menne
Anne anne.piotet at gmail.com writes: Is there an implemented method to predict residual expected survival times for parametric/Cox PH models ? (I have modelled my data using the survival library) I would like to predict for a given subject (with a given profile ) having survived up to time

Re: [R] Running glm in batch and exporting results (AIC) to HTML

2005-09-20 Thread Stefano Calza
see ?html in Hmisc package HIH Ste On Tue, Sep 20, 2005 at 02:56:17PM +0200, [EMAIL PROTECTED] wrote: vincentsee vincent?write.table vincenthih vincent vincent__ vincentR-help@stat.math.ethz.ch mailing list

Re: [R] waveform filtering

2005-09-20 Thread jfontain
Quoting tom wright [EMAIL PROTECTED]: On Mon, 2005-19-09 at 10:36 -0400, tom wright wrote: I'm not an engineer so I hope I'm using the correct terminology here. I have a recorded waveform that I want to apply low and high pass filters too, are tehre already R functions existing to do this

Re: [R] waveform filtering

2005-09-20 Thread Huntsinger, Reid
fftw is a library to do FFTs (fast Fourier transform). It's excellent, but probably not necessary unless you have lots of long series and you use FFT's repeatedly (say in an iterative fitting procedure). R's fft() is essentially instantaneous for most one-shot applications. Reid Huntsinger

[R] Neat way of using R for pivoting?

2005-09-20 Thread BANNISTER, Keith
Hi, I'd like to use R to do what excel pivot tables do, and plot results. I've never used R before, and I've managed to do something, but it's quite a lot of code to do something simple. I can't help but think I'm not Doing it the R way. I could be using R for the wrong thing, in which case,

[R] Estimate predictor contribution in GAM models

2005-09-20 Thread Yves Magliulo
hi, i'm using gam() function from package mgcv. if G is my gam object, then SG=summary(G) Formula: y ~ +s(x0, k = 5) + s(x1) + s(x2, k = 3) Parametric coefficients: Estimate std. err.t ratioPr(|t|) (Intercept) 3.462e+07 1.965e+05 176.2 2.22e-16 Approximate

[R] Strange Result using weightedMedian

2005-09-20 Thread Oliver Duerr
Dear all, I found a strange result using R's weightedMedian function. Consider the following: x - c (0.2, 0.3, 0.5) w - c (1,1,2) weightedMedian(x,w) 0.3666 In cases like above, when the weights are integers, one could argue that the weighted median should be the same as the standard

[R] Course***R/Splus Programming *** 6 cities Nationwide / October 2005

2005-09-20 Thread elvis
XLSolutions Corporation (www.xlsolutions-corp.com) is proud to announce 2-day R/S-plus Fundamentals and Programming Techniques in 6 cities. www.xlsolutions-corp.com/Rfund.htm Seattle -- October 24th-25th, 2005 San Francisco ---

Re: [R] Extended Hypergeometric Distribution

2005-09-20 Thread Albatineh, Ahmed
Do you mean the generalized Hypergeometric distribution ?? * -Original Message- From: Martin Maechler [mailto:[EMAIL PROTECTED] Sent: Tuesday, September 20, 2005 5:31 AM To: Narcyz Ghinea Cc: r-help@stat.math.ethz.ch

Re: [R] Neat way of using R for pivoting?

2005-09-20 Thread Greg Snow
BANNISTER, Keith [EMAIL PROTECTED] 09/20/05 09:46AM Hi, I'd like to use R to do what excel pivot tables do, and plot results. R does not have pivot tables and I hope that it never does. My experiance with pivot tables is that they encourage poor initial design followed by

Re: [R] Transform variable number of rows per subject to column variables?

2005-09-20 Thread Kevin Bartz
Hey Bing, Reshape's the ticker -- ?reshape. For example, reshape(myFrame, idvar = ID, timevar = TEST.A) should do most of the trick. Kevin -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Bing Ho Sent: Monday, September 19, 2005 10:04 PM To:

[R] How to exclude a level from a factor

2005-09-20 Thread Qiong Yang
Hi, I could not use 'exlcude=' option in factor() to exclude a level from a existing factor. x is a factor: x [1] a b c Levels: a b c factor(x,exclude=c) [1] a b c Levels: a b c Warning message: NAs introduced by coercion However, c is not coded as NA. The following does not work either:

Re: [R] Strange Result using weightedMedian

2005-09-20 Thread Dieter Menne
Oliver Duerr Oliver.Duerr at genedata.com writes: Dear all, I found a strange result using R's weightedMedian function I could not find weightedMedian on CRAN; I know there is such a function in some Java lib, googling gave a reference to limma (Linear Models for Microarray Data by

Re: [R] Neat way of using R for pivoting?

2005-09-20 Thread hadley wickham
Hi Keith, You might want to check out my reshape package (http://had.co.nz/reshape/) which is very much pivot table inspired. I doesn't produce graphics yet, but the output is very amenable to being fed into existing R graphics function (especially lattice graphics). Hadley

Re: [R] How to exclude a level from a factor

2005-09-20 Thread Sundar Dorai-Raj
Qiong Yang wrote: Hi, I could not use 'exlcude=' option in factor() to exclude a level from a existing factor. x is a factor: x [1] a b c Levels: a b c factor(x,exclude=c) [1] a b c Levels: a b c Warning message: NAs introduced by coercion However, c is not coded as

Re: [R] Neat way of using R for pivoting?

2005-09-20 Thread Gabor Grothendieck
On 9/20/05, Greg Snow [EMAIL PROTECTED] wrote: BANNISTER, Keith [EMAIL PROTECTED] 09/20/05 09:46AM Hi, I'd like to use R to do what excel pivot tables do, and plot results. R does not have pivot tables and I hope that it never does. My experiance with pivot tables is that they

Re: [R] How to exclude a level from a factor

2005-09-20 Thread Uwe Ligges
Qiong Yang wrote: Hi, I could not use 'exlcude=' option in factor() to exclude a level from a existing factor. x is a factor: x [1] a b c Levels: a b c factor(x,exclude=c) [1] a b c Levels: a b c Warning message: NAs introduced by coercion However, c is not coded as

Re: [R] Neat way of using R for pivoting?

2005-09-20 Thread Liaw, Andy
From: Greg Snow BANNISTER, Keith [EMAIL PROTECTED] 09/20/05 09:46AM Hi, I'd like to use R to do what excel pivot tables do, and plot results. R does not have pivot tables and I hope that it never does. My experiance with pivot tables is that they encourage poor initial

Re: [R] How to exclude a level from a factor

2005-09-20 Thread Qiong Yang
Problem solved. Thanks a lot for your replies! x [1] a b c Levels: a b c factor(as.character(x),exclude=c) [1] abNA Levels: a b exclude= option may not work on factors. One has to convert the factor to character first. Qiong On Tue, 20 Sep 2005, Sundar Dorai-Raj wrote: Qiong

[R] Strange Result using weightedMedian

2005-09-20 Thread Oliver Dürr
weightedMedian is in the R.basic packag. To download see: See http://www.maths.lth.se/help/R/R.classes/#1.%20Introduction Best, Oliver __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting

Re: [R] Strange Result using weightedMedian

2005-09-20 Thread Frank E Harrell Jr
Oliver Duerr wrote: Dear all, I found a strange result using R's weightedMedian function. Consider the following: x - c (0.2, 0.3, 0.5) w - c (1,1,2) weightedMedian(x,w) 0.3666 In cases like above, when the weights are integers, one could argue that the weighted median should be

[R] help with estimating parameters with nls

2005-09-20 Thread Uri Iskin
Dear helpeRs, I have a vector containing values of incomes and I would like to estimate the three parameters of a Dagum distribution. Dagum himself recommends to use nonlinear least-squares method. I have read nls, optim (and the posting guide!) and still have not succeded. sipcf is my sorted

Re: [R] Neat way of using R for pivoting?

2005-09-20 Thread Patrick Burns
Gabor Grothendieck wrote: ... The fact that Excel has both an interactive interface and a script-based interface whereas R has only a script-based interface puts it ahead, not behind, R in at least some respects. Sorry, but I can't resist: That very much depends on if you are doing

Re: [R] Neat way of using R for pivoting?

2005-09-20 Thread Gabor Grothendieck
On 9/20/05, Patrick Burns [EMAIL PROTECTED] wrote: Gabor Grothendieck wrote: ... The fact that Excel has both an interactive interface and a script-based interface whereas R has only a script-based interface puts it ahead, not behind, R in at least some respects. Sorry, but I can't

Re: [R] Neat way of using R for pivoting?

2005-09-20 Thread Greg Snow
Gabor Grothendieck [EMAIL PROTECTED] 09/20/05 11:31AM Just one comment here lest we be arguing against a strawman. While I agree that reproducibility can be a problem with pivot tables if created interactively and this applies to just about anything you do in Excel if done interactively, it

Re: [R] Neat way of using R for pivoting?

2005-09-20 Thread Gabor Grothendieck
On 9/20/05, Greg Snow [EMAIL PROTECTED] wrote: Gabor Grothendieck [EMAIL PROTECTED] 09/20/05 11:31AM Just one comment here lest we be arguing against a strawman. While I agree that reproducibility can be a problem with pivot tables if created interactively and this applies to just

[R] annotating an axis in bwplot (lattice)

2005-09-20 Thread Sebastian Luque
Hi, I'd like to add, say, the sample size for every group in a bwplot as a parenthetical annotation to the axis. Here's a sketch: --8---cut here---start-8--- require(Hmisc) age - sample(1:100, 1000, replace = TRUE) sex - gl(2, 8, 1000, c(Male, Female)) grp -

[R] Problem with read.spss() and as.data.frame(), or: alternative to subset()?

2005-09-20 Thread Dirk Enzmann
Trying to select a subset of cases (rows of data) I encountered several problems: Firstly, because I did not read the help to read.spss() thoroughly enough, I treated the data read as a data frame. For example, dr2000 - read.spss('myfile.sav') d - subset(dr2000,RBINZ99 0) and thus received

[R] Problem with read.spss() and as.data.frame(), or: alternative to subset()? (2)

2005-09-20 Thread Dirk Enzmann
In my previous mail there is a typo. Instead of d - subset(dr2000,dr2000$RBINZ99) I used d - subset(dr2000,dr2000$RBINZ99 0) * Dr. Dirk Enzmann Institute of Criminal Sciences Dept. of Criminology Edmund-Siemers-Allee 1 D-20146 Hamburg Germany

[R] why this postscript didn't work?

2005-09-20 Thread WWei
Hi, List, I used the following codes to generate ps plots but foo.ps contains nothing. Would someone please point out what is wrong with my codes? Thanks a million! postscript('foo.ps') par(mfrow=c(2,1)) par(mfg=c(1,1)) hist(rnorm(100),col='blue') par(mfrow=c(2,2)) par(mfg=c(2,1))

Re: [R] why this postscript didn't work?

2005-09-20 Thread Marc Schwartz
On Tue, 2005-09-20 at 18:30 -0500, [EMAIL PROTECTED] wrote: Hi, List, I used the following codes to generate ps plots but foo.ps contains nothing. Would someone please point out what is wrong with my codes? Thanks a million! postscript('foo.ps') par(mfrow=c(2,1)) par(mfg=c(1,1))

Re: [R] Searchable mailing list archives -- not reachable FIXED

2005-09-20 Thread Robert King
http://tolstoy.newcastle.edu.au/~rking/R/ is now working again, following a power supply problem that kept it off the net. Unfortunately, even when it is working it is firewalled away from pings. Robert. Martin Maechler wrote: AFAIK, the correct URL --- as also used from CRAN's search page

Re: [R] help with estimating parameters with nls

2005-09-20 Thread Francisco J. Zagmutt
RSiteSearch(Dagum) Francisco From: Uri Iskin [EMAIL PROTECTED] To: r-help@stat.math.ethz.ch Subject: [R] help with estimating parameters with nls Date: Tue, 20 Sep 2005 15:44:07 -0300 Dear helpeRs, I have a vector containing values of incomes and I would like to estimate the three parameters

Re: [R] why this postscript didn't work?

2005-09-20 Thread WWei
Marc, Thank you very much for your help! Your codes worked perfectly fine here. I am using R 2.0.1 on Red Hat 9, sorry for not specifying it in my original post. I didn't know we could use layout to plot, it is much easier than par(mfg). Thanks a lot! Best, Auston Marc Schwartz [EMAIL

[R] Seu computador est? infectado com o virus [EMAIL PROTECTED]

2005-09-20 Thread security-check
__ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html