Re: [R] mahalanobis

2007-06-01 Thread gatemaze
On 31/05/07, Anup Nandialath <[EMAIL PROTECTED]> wrote: > oops forgot the example example > > try this line > > sqrt(mahalanobis(all, colMeans(predictors), cov(all), FALSE) Hi and thanks for the reply Anup. Unfortunately, I had a look on the example before posting but not much of a help... I did so

[R] mahalanobis

2007-05-31 Thread gatemaze
Hi, I am not sure I am using correctly the mahalanobis distnace method... Suppose I have a response variable Y and predictor variables X1 and X2 all <- cbind(Y, X1, X2) mahalanobis(all, colMeans(all), cov(all)); However, my results from this are different from the ones I am getting using another

Re: [R] normality tests [Broadcast]

2007-05-25 Thread gatemaze
Thank you all for your replies they have been more useful... well in my case I have chosen to do some parametric tests (more precisely correlation and linear regressions among some variables)... so it would be nice if I had an extra bit of support on my decisions... If I understood well from al

Re: [R] normality tests

2007-05-25 Thread gatemaze
On 25/05/07, Frank E Harrell Jr <[EMAIL PROTECTED]> wrote: > [EMAIL PROTECTED] wrote: > > Hi all, > > > > apologies for seeking advice on a general stats question. I ve run > > normality tests using 8 different methods: > > - Lilliefors > > - Shapiro-Wilk > > - Robust Jarque Bera > > - Jarque Bera

[R] normality tests

2007-05-25 Thread gatemaze
Hi all, apologies for seeking advice on a general stats question. I ve run normality tests using 8 different methods: - Lilliefors - Shapiro-Wilk - Robust Jarque Bera - Jarque Bera - Anderson-Darling - Pearson chi-square - Cramer-von Mises - Shapiro-Francia All show that the null hypothesis that

[R] part or semi-partial correlation

2007-05-22 Thread gatemaze
Is it possible to conduct part (also called semi-partial correlation) with R. The help.search produces no results and there is also nothing into the archive, well one post asking what is part correlation. Just quickly from Field [Discovering statistics using spss]: "When we do a partial correlatio

[R] partial correlation function

2007-05-22 Thread gatemaze
Hi, after reading the archives I found some methods... adopted and modified one of them to the following. I think it is correct after checking and comparing the results with other software... but if possible someone could have a look and spot any mistakes I would be grateful. Thanks pcor3 <- func

Re: [R] partial correlation significance

2007-05-18 Thread gatemaze
On 18/05/07, [EMAIL PROTECTED] <[EMAIL PROTECTED]> wrote: > > Hi, > > among the many (5) methods that I found in the list to do partial > correlation in the following two that I had a look I am getting different > t-values. Does anyone have any clues on why is that? The source code is > below. Than

[R] partial correlation significance

2007-05-18 Thread gatemaze
Hi, among the many (5) methods that I found in the list to do partial correlation in the following two that I had a look I am getting different t-values. Does anyone have any clues on why is that? The source code is below. Thanks. pcor3 <- function (x, test = T, p = 0.05) { nvar <- ncol(x) nd

[R] t value two.sided and one.sided

2007-05-10 Thread gatemaze
Hi, on a > summary(lm(y~x)) are the computed t-values for two.sided or one.sided. By looking on some tables they seem like they are for two.sided. Is it possible to have them for one.sided? If this makes sense... Thanks. [[alternative HTML version deleted]] _

Re: [R] Error in plot.new() : figure margins too large

2007-05-09 Thread gatemaze
On 09/05/07, Prof Brian Ripley <[EMAIL PROTECTED]> wrote: > > On Wed, 9 May 2007, [EMAIL PROTECTED] wrote: > > > The code is: > > > > postscript(filename, horizontal=FALSE, onefile=FALSE, paper="special", > > You have not set a width or height, so please do your homework. Thanks a lot for that an

[R] Error in plot.new() : figure margins too large

2007-05-09 Thread gatemaze
Yes, I already had a look on previous posts but nothing is really helpful to me. The code is: postscript(filename, horizontal=FALSE, onefile=FALSE, paper="special", bg="white", family="ComputerModern", pointsize=10); par(mar=c(5, 4, 0, 0) + 0.1); plot(x.nor, y.nor, xlim=c(3,6), ylim=c(20,90), pch=

Re: [R] plots - scatterplots - find index of a point in a list

2007-05-08 Thread gatemaze
On 08/05/07, John Kane <[EMAIL PROTECTED]> wrote: > > Try ?locator Thanks. Your tip also lead to another function: ?identify to add my two cents. --- [EMAIL PROTECTED] wrote: > > > Hi, > > > > is it possible to find the index of a point of a > > plot (e.g. scatterplot) in > > an easy way? > > >

[R] plots - scatterplots - find index of a point in a list

2007-05-08 Thread gatemaze
Hi, is it possible to find the index of a point of a plot (e.g. scatterplot) in an easy way? Eg. x <- c(1:5); y <- c(1:5); plot(x, y); On the plot if I move my cursor on top of a point or click on it is it possible to have its index printed or its exact value? Any clues? Thanks. [[alte

Re: [R] general question about use of list

2007-04-30 Thread gatemaze
On 30/04/07, Y G <[EMAIL PROTECTED]> wrote: > > Hi, > > is this list only related to R issues or it has a broader context > regarding questions and discussions about statistics. Is there any other > email list or forum for that? Well, just to quickly self-reply my question. From http://www.r-proj