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 some

[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

[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

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 -

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

[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 -

[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

[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)

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. Thanks.

[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]]

[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),

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 and to Phil for

[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.

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? Eg. x -

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