[R] a R function for sort a data frame.

2005-04-01 Thread Mario Morales
Is there a R function for sort a data frame by a variable ? I know sort a vector, but I don't know sort a data frame by a column. Can you help me ? the sort() function don't work with data frame. __ R-help@stat.math.ethz.ch mailing list

Re: [R] programming conventions

2005-04-01 Thread Brian D Ripley
On Thu, 31 Mar 2005, Sebastian Luque wrote: I'm trying to make my R scripts more readable by others, so I searched for some R programming conventions and found the following two: - http://www.ci.tuwien.ac.at/Conferences/useR-2004/Keynotes/Maechler.pdf - http://www.maths.lth.se/help/R/RCC

Re: [R] about get stdout from other program

2005-04-01 Thread Brian D Ripley
On Fri, 1 Apr 2005, Michael S wrote: if I want to use stdout from other language as my R program input ,which is the best way for design the API,using Pipe function or produce a temporary file,using scan function to read the file ? A file for portability, since pipes do not work well on

[R] error in kmeans

2005-04-01 Thread Asha Jayanthi
I am trying to generate kmean of 10 clusters for a 165 x 165 matrix. i do not see any errors known to me. But I get this error on running the script Error: empty cluster: try a better set of initial centers the commands are M -matrix(scan(R_mutual,n = 165 * 165),165,165,byrow = T) cl -

Re: [R] error messages on R CMD check

2005-04-01 Thread Uwe Ligges
David Firth wrote: Dear Johannes I have noticed the same complaint, and you are right that it can be unconnected with faulty S3 methods, etc., in your package. For example, in my relimp package (current version 0.9-1, new on CRAN today) the DESCRIPTION file has Suggests: tcltk and that works

Re: [R] French Curve

2005-04-01 Thread Martin Maechler
dream == dream home [EMAIL PROTECTED] on Wed, 30 Mar 2005 12:27:08 -0800 writes: dream Dear R experts, Did someone implemented French Curve dream yet? Or can anyone point me some papers that I can dream follow to implement it? Are you talking about splines ? I vaguely

[R] Classification of an image

2005-04-01 Thread Poizot Emmanuel
Dear all, I need to do a automatic classification of a raster file (image) using training samples. I would like to know if there is a library able to do such a work. Thanks Emmanuel Poizot __ R-help@stat.math.ethz.ch mailing list

[R] glm with poisson errors

2005-04-01 Thread Christoph Scherber
Dear R users, I have data on percentage leaf area damaged (in classes, e.g. 1%, 5%, 10%) in plants. My two questions are: (1) Could I use a glm with poisson errors on these data? (2) Could I still use this glm with poisson errors after arcsine transformation of the data? Thank you very much

Re: [R] about get stdout from other program

2005-04-01 Thread Jan T. Kim
On Fri, Apr 01, 2005 at 01:29:01AM +, Michael S wrote: if I want to use stdout from other language as my R program input ,which is the best way for design the API,using Pipe function or produce a temporary file,using scan function to read the file ? For the scan function, it makes no

[R] optim problem, nls regression

2005-04-01 Thread Kervahu Anne
Hi, I try to fit a non linear regression by minimising the sum of the sum of squares. The model is number[2]-(x/number[1])^number[3] Number [2] and number [1] change as the data changes but for all the set of data number[3] must be identical. I have 3 set of data (x1,y1), (x2,y2), (x3,y3).

[R] anova and post hoc

2005-04-01 Thread Paolo Ariano
Hi all i'm trying to do some statistic on my data from neuronal migration and i've to use anova and post-hoc test, does someone have any tips'trick about a pdf file to read about ? I've a control and five differents conditions and for every condition i've a different number of data, am i on the

RE: [R] Re:Plotting to A4 and replacing x-axis with actual years.

2005-04-01 Thread Marshall Mdoka
1)I have figures in a 3 rows and 3 columns. I tried to plot with: postscript(trend_jfm.eps, horizontal=FALSE, onefile=TRUE) I am getting the output cut. I am in linux (slackware) version 9 I think. I have tried to add in some of the extra conditions like for A4 as in the manual but still

[R] new package - paraNormal distribution

2005-04-01 Thread Barry Rowlingson
I've just completed a package for generating numbers from the paraNormal distribution. This is the distribution of numbers that you are currently thinking of. Sample usage: rParaNorm(1) [1] 3 - which is what I was thinking of when I pressed 'Return'. for generating multiple numbers, a

Re: [R] hexbin and grid - input data values as coordinates

2005-04-01 Thread Martin Maechler
Paul == Paul Murrell [EMAIL PROTECTED] on Fri, 01 Apr 2005 10:45:16 +1200 writes: Paul Hi Adaikalavan Ramasamy wrote: Dear all, I am trying to use hexbin and read the very interesting article on grid (

[R] (no answer)

2005-04-01 Thread Frank E Harrell Jr
I wish to perform brain surgery this afternoon at 4pm and don't know where to start. My background is the history of great statistician sports legends but I am willing to learn. I know there are courses and numerous books on brain surgery but I don't have the time for those. Please direct me

[R] R FAQ 2.8

2005-04-01 Thread Robin Hankin
Hi [thanks to everyone for advice on Recall() and sapply()] The R FAQ section 2.8 discusses how to cite R in publications, but does not (AFAICS) tell me how to describe R in a sentence. To wit, in my latest paper (destined for Rnews) one sentence reads: The R programming language (3) has been

[R] [R-pkgs] perturb package for evaluating collinearity

2005-04-01 Thread John Hendrickx
I've uploaded the R package perturb to CRAN. Perturb contains two programs for evaluating collinearity. Colldiag calculates condition indexes and variance decomposition proportions to detect and track down collinear sets of variables. Perturb takes a different approach. It re-estimates the model

RE: [R] a R function for sort a data frame.

2005-04-01 Thread Doran, Harold
You could do something like this dataframe[order(dataframe$variable),] -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Mario Morales Sent: Thursday, March 31, 2005 10:23 PM To: r-help@stat.math.ethz.ch Subject: [R] a R function for sort a data frame.

RE: [R] a R function for sort a data frame.

2005-04-01 Thread bogdan romocea
dfr - data.frame(sample(1:50,10),sample(1:50,10)) colnames(dfr) - c(a,b) dfr - dfr[order(dfr$a),] dfr - dfr[order(-dfr$a),] -Original Message- From: Mario Morales [mailto:[EMAIL PROTECTED] Sent: Thursday, March 31, 2005 10:23 PM To: r-help@stat.math.ethz.ch Subject: [R] a R function for

Re: [R] a R function for sort a data frame.

2005-04-01 Thread Dimitris Rizopoulos
look at ?order(), e.g., dat[order(dat[,1]),] will sort your data.frame wrt the first column. I hope it helps. Best, Dimitris Dimitris Rizopoulos Ph.D. Student Biostatistical Centre School of Public Health Catholic University of Leuven Address: Kapucijnenvoer 35, Leuven, Belgium Tel:

Re: [R] a R function for sort a data frame.

2005-04-01 Thread TEMPL Matthias
Hi, x - data.frame(A=rnorm(10), B=round(runif(10,0,10)), C=rnorm(10)) #one way: #increasing order x.sort.in - x[order(x$B),] x.sort.in #decreasing order: x.sort.de - x[rev(order(x$B)),] x.sort.de Best, Matthias -Ursprüngliche Nachricht- Von: [EMAIL PROTECTED] [mailto:[EMAIL

Re: [R] error in kmeans

2005-04-01 Thread Ingmar Visser
Hi Asha, kmeans is a non-deterministic routine. The help page says the following about the centers argument: centers: Either the number of clusters or a set of initial cluster centers. If the first, a random set of rows in 'x' are chosen as the initial centers. Hence,

Re: [R] error in kmeans

2005-04-01 Thread Uwe Ligges
Asha Jayanthi wrote: I am trying to generate kmean of 10 clusters for a 165 x 165 matrix. i do not see any errors known to me. But I get this error on running the script Error: empty cluster: try a better set of initial centers the commands are M -matrix(scan(R_mutual,n = 165 *

Re: [R] a R function for sort a data frame.

2005-04-01 Thread Duncan Murdoch
On Thu, 31 Mar 2005 22:22:32 -0500, Mario Morales [EMAIL PROTECTED] wrote : Is there a R function for sort a data frame by a variable ? I know sort a vector, but I don't know sort a data frame by a column. Can you help me ? the sort() function don't work with data frame. This is a FAQ, but

Re: [R] French Curve

2005-04-01 Thread Martyn Plummer
On Fri, 2005-04-01 at 09:30 +0200, Martin Maechler wrote: dream == dream home [EMAIL PROTECTED] on Wed, 30 Mar 2005 12:27:08 -0800 writes: dream Dear R experts, Did someone implemented French Curve dream yet? Or can anyone point me some papers that I can dream follow to

Re: [R] a R function for sort a data frame.

2005-04-01 Thread Romain Francois
Le 01.04.2005 05:22, Mario Morales a écrit : Is there a R function for sort a data frame by a variable ? I know sort a vector, but I don't know sort a data frame by a column. Can you help me ? the sort() function don't work with data frame. take a look at ?order as in : A -

RE: [R] a R function for sort a data frame.

2005-04-01 Thread Liaw, Andy
http://finzi.psych.upenn.edu/R/Rhelp02a/archive/39835.html (First hit using RSiteSearch(sort.data.frame) in R-2.1.0 alpha.) Andy From: Mario Morales Is there a R function for sort a data frame by a variable ? I know sort a vector, but I don't know sort a data frame by a column. Can you

[R] French Curve

2005-04-01 Thread Ken Knoblauch
I remember that my father had a French curve: it was a plastic template used for drawing which had several smooth edges of varying curvature. You could use it to draw a wide variety of curved shapes. No doubt the French called it something else. Nobody, up and down the corridor here, of age to

[R] plot axis appearance problem

2005-04-01 Thread Yogesh K. Tiwari
Hello, When I plot any data with simple plot command in R, for example :- plot(time,co2,ylim=c(350,380),xlim=c(1993,2003),xlab=NA,ylab=NA,type=p,col=5) Then the first value of x-axis(350) and y-axis(1993) never starts from origin, always they sifted from the origin. Is there any command that

RE: [R] new package - paraNormal distribution

2005-04-01 Thread abunn
What a great package. I imagined the answer before hand, so I guess it's right... R rParaNorm(1) [1] 0+42i R This will really help me through a Bayesian analysis I have to do today. Thanks! __ R-help@stat.math.ethz.ch mailing list

RE: [R] Classification of an image

2005-04-01 Thread Berton Gunter
Search CRAN! -- Bert Gunter -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Poizot Emmanuel Sent: Thursday, March 31, 2005 11:59 PM To: r-help@stat.math.ethz.ch Subject: [R] Classification of an image Dear all, I need to do a automatic

RE: [R] R FAQ 2.8

2005-04-01 Thread Berton Gunter
Inline below. -- Bert Gunter Genentech Non-Clinical Statistics South San Francisco, CA The R FAQ section 2.8 discusses how to cite R in publications, but does not (AFAICS) tell me how to describe R in a sentence. To wit, in my latest paper (destined for Rnews) one sentence reads:

[R] fatal error unused tempdir

2005-04-01 Thread Chris Ward
I am running R on an XP machine. Lately I have been unable to start R in any mode without getting the 'Fatal Error: cannnot find unused tempdir name' message. After reading previous postings I went in and deleted all my Rtmp* directories in /documents and settings/USER/local/tmp which didn't

RE: [R] R FAQ 2.8

2005-04-01 Thread Wiener, Matthew
I have used the R statistical computing environment. No-one has ever asked me to change it, but maybe someone else has something better. Hope this helps, Matt Wiener -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Robin Hankin Sent: Friday, April 01,

Re: [R] error in kmeans

2005-04-01 Thread Prof Brian Ripley
On Fri, 1 Apr 2005, Asha Jayanthi wrote: I am trying to generate kmean of 10 clusters for a 165 x 165 matrix. i do not see any errors known to me. But I get this error on running the script Error: empty cluster: try a better set of initial centers the commands are M -matrix(scan(R_mutual,n = 165

Re: [R] R FAQ 2.8

2005-04-01 Thread Uwe Ligges
Robin Hankin wrote: Hi [thanks to everyone for advice on Recall() and sapply()] The R FAQ section 2.8 discusses how to cite R in publications, but does not (AFAICS) tell me how to describe R in a sentence. To wit, in my latest paper (destined for Rnews) one sentence reads: The R programming

Re: [R] (no answer)

2005-04-01 Thread David Forrest
On Fri, 1 Apr 2005, Frank E Harrell Jr wrote: I wish to perform brain surgery this afternoon at 4pm and don't know where to start. My background is the history of great statistician sports legends but I am willing to learn. I know there are courses and numerous books on brain surgery but I

Re: [R] error in kmeans

2005-04-01 Thread Gavin Simpson
Asha Jayanthi wrote: I am trying to generate kmean of 10 clusters for a 165 x 165 matrix. i do not see any errors known to me. But I get this error on running the script Error: empty cluster: try a better set of initial centers the commands are M -matrix(scan(R_mutual,n = 165 *

Re: [R] Re:Plotting to A4 and replacing x-axis with actual years.

2005-04-01 Thread Uwe Ligges
Marshall Mdoka wrote: 1)I have figures in a 3 rows and 3 columns. I tried to plot with: postscript(trend_jfm.eps, horizontal=FALSE, onefile=TRUE) I am getting the output cut. I am in linux (slackware) version 9 I think. I have tried to add in some of the extra conditions like for A4 as in the

Re: [R] optim problem, nls regression

2005-04-01 Thread Uwe Ligges
Kervahu Anne wrote: Hi, I try to fit a non linear regression by minimising the sum of the sum of squares. The model is number[2]-(x/number[1])^number[3] Number [2] and number [1] change as the data changes but for all the set of data number[3] must be identical. I have 3 set of data (x1,y1),

RE: [R] a R function for sort a data frame.

2005-04-01 Thread Francisco J. Zagmutt
dat-data.frame(ID=seq(1,10),x=runif(10,2,50)) dat-dat[sort.list(dat$ID),] Francisco From: bogdan romocea [EMAIL PROTECTED] To: [EMAIL PROTECTED] CC: r-help@stat.math.ethz.ch Subject: RE: [R] a R function for sort a data frame. Date: Fri, 1 Apr 2005 06:07:05 -0800 (PST) dfr -

Re: [R] French Curve

2005-04-01 Thread Marc Schwartz
On Fri, 2005-04-01 at 17:00 +0200, Martyn Plummer wrote: On Fri, 2005-04-01 at 09:30 +0200, Martin Maechler wrote: dream == dream home [EMAIL PROTECTED] on Wed, 30 Mar 2005 12:27:08 -0800 writes: dream Dear R experts, Did someone implemented French Curve dream yet? Or

[R] boot function

2005-04-01 Thread Christina D Smith
The boot function documentation says, The data as a vector, matrix or data frame. If it is a matrix or data frame then each row is considered as one multivariate observation. From trying several things, I've concluded that this means it resamples the rows but calculates the statistic on the

[R] CI for Ratios of Variance components in lme?

2005-04-01 Thread Berton Gunter
My apologies if this is obvious: Is there a simple way (other than simulation or bootstrapping) to obtain a (approximate)confidence interval for the ratio of 2 variance components in a fitted lme model? -- In particular, if there are only 2 components (1 grouping factor). I'm using nlme but lme4

Re: [R] R FAQ 2.8

2005-04-01 Thread Dave Braze
In a paper for which I recently received reviews, I have the sentence Statistical analyses were carried out with the R statistical system (R Development Core Team, 2004). Reviewers did not complain (about that). The particular journal falls in the domain of applied psychology. -Dave Braze On

[R] Error in colMeans ... what's wrong with my data?

2005-04-01 Thread Hartmut Weinrebe
Having searched and searched I still haven't found what's the problem with my data (I've attached the relevant file). Every time I tried to use the CANCOR-Function I got error messages. So I turned to check my two sets of variables separately by using the colMeans-Function. With one set no

Re: [R] plot axis appearance problem

2005-04-01 Thread Prof Brian Ripley
?par look at xaxs. Or just read `An Introduction to R' more carefully! On Fri, 1 Apr 2005, Yogesh K. Tiwari wrote: Hello, When I plot any data with simple plot command in R, for example :- plot(time,co2,ylim=c(350,380),xlim=c(1993,2003),xlab=NA,ylab=NA,type=p,col=5) Then the first value of

[R] Ordering scales in xYplot.Hmisc

2005-04-01 Thread Jose A. Hernandez
Dear R community, I am using xYplot() from the Hmisc package. The package works great to plot means + CI. But I am having issues handling the scales. I am plotting Soil Clay content vs Soil depth by land use. Usually in this type of graphs it is better to place the variable soil depth in the

RE: [R] French Curve

2005-04-01 Thread Liaw, Andy
From: Ken Knoblauch I remember that my father had a French curve: it was a plastic template used for drawing which had several smooth edges of varying curvature. You could use it to draw a wide variety of curved shapes. No doubt the French called it something else. Nobody, up and

RE: [R] plot axis appearance problem

2005-04-01 Thread Liaw, Andy
From: Yogesh K. Tiwari Hello, When I plot any data with simple plot command in R, for example :- plot(time,co2,ylim=c(350,380),xlim=c(1993,2003),xlab=NA,ylab=N A,type=p,col=5) ^ R can only do what you _ask_ it to do, not what you want it to do. You have

Re: [R] Stratified Bootstrap question

2005-04-01 Thread Qian An
Dear Tim, Thank you so much for addressing all my questions. This really makes sense. I talked with my advisor yesterday about how to do bootstrapping for my scenario: random clinic + random subject within clinic. She suggested that only clinic are independent units, so I can only resample

[R] step error

2005-04-01 Thread vasilis pappas
Could anyone tell me what am I doing wrong? pro-function(indep,dep){ + d-data.frame(indep) + form-formula(lm(dep~.,data=d)) + forward-step(lm(dep~X1,data=d),scope=form,trace=0,direction='f') + return(forward) + } pro(m,q) Error in inherits(x, data.frame) : Object d not found Where q is a

Re: [R] plot axis appearance problem

2005-04-01 Thread Marc Schwartz
On Fri, 2005-04-01 at 18:13 +0200, Yogesh K. Tiwari wrote: Hello, When I plot any data with simple plot command in R, for example :- plot(time,co2,ylim=c(350,380),xlim=c(1993,2003),xlab=NA,ylab=NA,type=p,col=5) Then the first value of x-axis(350) and y-axis(1993) never starts from

Re: [R] hexbin and grid - input data values as coordinates

2005-04-01 Thread Adaikalavan Ramasamy
Thank you to Paul Murrell and Martin Maechler for their help. pushHexport() and the rest of the codes have done the trick. I spent the afternoon trying to code up something that might be used as grid.abline() and grid.grid() before I read Martin's suggestion. Sigh. But here it is anyway in case

Re: [R] French Curve

2005-04-01 Thread Charles C. Berry
Ken Knoblauch knoblauch at lyon.inserm.fr writes: I remember that my father had a French curve: it was a plastic template used for drawing which had several smooth edges of varying curvature. Yes. These were used by graphic artists, among others. Some pictures of these are posted at

Re: [R] French Curve

2005-04-01 Thread Ted Harding
On 01-Apr-05 Martyn Plummer wrote: On Fri, 2005-04-01 at 09:30 +0200, Martin Maechler wrote: dream == dream home [EMAIL PROTECTED] on Wed, 30 Mar 2005 12:27:08 -0800 writes: dream Dear R experts, Did someone implemented French Curve dream yet? Or can anyone point me some

[R] [R-pkgs] Version 0.95-1 of the lme4 package

2005-04-01 Thread Douglas Bates
This may sound like an April Fool's joke but it isn't really. I just uploaded version 0.95-1 of the lme4 package for R. This is the package with the new versions of the lme (linear mixed-effects) function for R and the GLMM (generalized linear mixed models) function. One distinguishing feature

RE: [R] French Curve

2005-04-01 Thread Ken Knoblauch
Here is what my colleague dug up and his reaction to it, afterwards, and then one of mine: Pistolet, nom masculin DESSIN. Instrument de tracé permettant de dessiner les lignes courbes dans les tracés géométriques. Synon. curvigraphe, virgule. Les outils de l'écrivain plumiste sont: la plume, le

Re: [R] Error in colMeans ... what's wrong with my data?

2005-04-01 Thread Federico Calboli
On Fri, 2005-04-01 at 20:27 +0200, Hartmut Weinrebe wrote: Having searched and searched I still haven't found what's the problem with my data (I've attached the relevant file). Every time I tried to use the CANCOR-Function I got error messages. So I turned to check my two sets of variables

[R] Vernier Caliper function vernier()

2005-04-01 Thread Ted Harding
Hi Folks, I don't think I'm up to implementing a function for drawing French Curves, as requested by dream. However, the discussion about it took me back so vividly to the old days that, paranormally, I felt once again the urge to ascertain magnitudes as it used to be, and should be, done. I

Re: [R] Ordering scales in xYplot.Hmisc

2005-04-01 Thread Deepayan Sarkar
On Friday 01 April 2005 12:41, Jose A. Hernandez wrote: Dear R community, I am using xYplot() from the Hmisc package. The package works great to plot means + CI. But I am having issues handling the scales. I am plotting Soil Clay content vs Soil depth by land use. Usually in this type of

Re: [R] French Curve

2005-04-01 Thread Marc Schwartz
On Fri, 2005-04-01 at 20:07 +0100, [EMAIL PROTECTED] wrote: snip I still have some, from the 1950s ... The curves in the edges are supposed to be segments of logarithmic spirals (which ensures a kind of self-similarity on different scales). A nice picture is at

Re:[R] (no answer)

2005-04-01 Thread vernel3
Ben wrote: I wish to perform brain surgery this afternoon at 4pm and don't know where to start. My background is the history of great statistician sports legends but I am willing to learn. I know there are courses and numerous books on brain surgery but I don't have the time for those. ... It

Re: [R] Vernier Caliper function vernier()

2005-04-01 Thread Rajarshi Guha
On Fri, 2005-04-01 at 21:31 +0100, [EMAIL PROTECTED] wrote: Its definition is below. Test it by entering vernier(pi+sqrt(.182)) By following the instructions you should find that the result is x = 3 Units + 5 Tenths + 7 Hundredths which is good enough for anyone. It generalises

Re: [R] French Curve

2005-04-01 Thread Ted Harding
On 01-Apr-05 Marc Schwartz wrote: On Fri, 2005-04-01 at 20:07 +0100, [EMAIL PROTECTED] wrote: snip [...] Splines, in the drawing-office sense, were long narrow (about 1/4 inch wide) strips of thin springy metal with, along their length, little flanges at right-angles to the plane of the

Re: [R] French Curve

2005-04-01 Thread Marc Schwartz
On Fri, 2005-04-01 at 22:56 +0100, [EMAIL PROTECTED] wrote: On 01-Apr-05 Marc Schwartz wrote: On Fri, 2005-04-01 at 20:07 +0100, [EMAIL PROTECTED] wrote: snip [...] Splines, in the drawing-office sense, were long narrow (about 1/4 inch wide) strips of thin springy metal with,

[R] R mailing list archive difficulty

2005-04-01 Thread François Pinard
Hi, people! This is my first babble on this list, please be kind! :-) Last Tuesday, I wrote to the (likely) Webmaster of the R site to report a little problem, but also to ask for advice about how to get a bulk copy of the mailing list archives, from 2002 to now. While I quite understand that

Re: [R] Stratified Bootstrap question

2005-04-01 Thread Tim Hesterberg
Qian wrote: I talked with my advisor yesterday about how to do bootstrapping for my scenario: random clinic + random subject within clinic. She suggested that only clinic are independent units, so I can only resample clinic. But I think that since subjects are also independent within clinic, shall

Re: [R] R mailing list archive difficulty

2005-04-01 Thread Jonathan Baron
On 04/01/05 18:15, Franois Pinard wrote: I would like to get hold on a copy of R mailing lists archives, for local, off-line, progressive perusal (I find Web-based browsing of email extremely inefficient). So, I recursively got archives from `ftp://ftp.stat.math.ethz.ch/Mail-archives/'. The

[R] Importing surface data

2005-04-01 Thread Wayne Walker
Hello, I'm attempting to create a surface plot similar to that produced with the volcano data. Athough I understand persp() and can plot the volcano data just fine, I'm struggling with how best to format my own data file for subsequent importing. I've read through the available materials, but

RE: [R] Survey of moving window statistical functions - still looking f or fast mad function

2005-04-01 Thread Vadim Ogranovich
Hi, First, let me thank Jaroslaw for making this survey. I find it quite illuminating. Now the questions: * the #1 solution below (based on cumsum) is numerically unstable. Specifically if you do the runmean on a positive vector you can easily get negative numbers due to rounding errors.

Re: [R] Survey of moving window statistical functions - still looking f or fast mad function

2005-04-01 Thread Gabor Grothendieck
Jaroslaw's article was great. In fact, it was used as the basis for rapply and some optimized special cases that will be included in the R 2.1.0 version of zoo (which have been coded but not yet released). Regarding numerically stable summation, check out the idea behind the following which I

Re: [R] R mailing list archive difficulty

2005-04-01 Thread Franois Pinard
[Franois Pinard] I would like to get hold on a copy of R mailing lists archives [...] [Jonathan Baron] the archives themselves are at [...] (and probably listed in other places). Thanks a lot for bringing me to this URL. I've been there before, but did not realise on first visit that the

[R] using GAM to assess the linearity in logistic regression

2005-04-01 Thread ronggui
as agresti(2002) points out that we had better to screen the data to see if the the logit(pi) and the predictor has linear realtionship in logistic regressin.and i find some materials in MASS and the refernce of s-plus.but it is a bit simple and i can not exactly master the means to assess

Re: [R] using GAM to assess the linearity in logistic regression

2005-04-01 Thread Wensui Liu
I am a little confused about what you asked. If you want to assess the linearity in logistic regression, why do you want to use GAM instead of GLM? As far as I understand, GAM is used to capture nonlinearity rather linearity. Am I right here? On Apr 1, 2005 10:19 PM, ronggui [EMAIL

Re: [R] using GAM to assess the linearity in logistic regression

2005-04-01 Thread ronggui
maybe the idea is simle,but the details is beyond me.you are right,gam can capture the non-linearity.but if the results from gam shows little evidence on on-linearity,then we can assume linearity exists. am i right? from agresti(2002): ... Before fitting the model and making such

Re: [R] Error in colMeans ... what's wrong with my data?

2005-04-01 Thread Prof Brian Ripley
On Fri, 1 Apr 2005, Hartmut Weinrebe wrote: Having searched and searched I still haven't found what's the problem with my data (I've attached the relevant file). Every time I tried to use the CANCOR-Function I got error messages. So I turned to check my two sets of variables separately by using