similar to: R example of Mahalanobis, Kolmogrov, ROC, Gini, Delta

Dear helper. I am a student in SFSU. I had a class that introduce a new open-source "R" in this semester. It's, however, hard to understand how to make "R" command. Especially, at this time we( including my friend) have to explain some concepts (e.g. mahalanobis distance, kolmogorov-smirnow statistics, ROC curve, Gini coefficient, and the delta approach) We have to

Hello, Thank you very much for your help. How can I draw a Lorenz curve with several replications ? Here is an example with 4 replications: hosts=c(23,31,19,10,7,7,3, 39,40,8,3,6,2,2, 47,17,8,10,6,11,1, 30,30,10,0,15,15,0) parasites=rep(seq(from=0,to=6,by=1),4) replications=c(rep(1,7),rep(2,7),rep(3,7),rep(4,7)) test <- cbind(parasites,hosts,replications) Should I

Hello dear R-help members, I would appreciate any help in understanding how the rpart function computes the "improve" (which is given in fit$split) when using the split='information' parameter. Thanks to Professor Atkinson help, I was able to find how this is done in the case that split='gini'. By following the explanation here:

http://www.nabble.com/file/p24569511/mahalanobis.txt mahalanobis.txt http://www.nabble.com/file/p24569511/concentrations.txt concentrations.txt Dear Forum members, I have a problem calculating mahalanobis distances. My data file mahalanobis.txt and categories file concentrations.txt are attached. I do the following steps: x <- as.matrix(read.table("mahalanobis.txt", header=TRUE))

Dear R-experts, Searching the help archives I found a recommendation to do multivariate outlier identification by mahalanobis distances based on a robustly estimated covariance matrix and compare the resulting distances to a chi^2-distribution with p (number of your variables) degrees of freedom. I understand that compared to euclidean distances this has the advantage of being scale-invariant.

Dear all Why isn'it possible to calculate Mahalanobis distances with R for a matrix with 1 row (observations) more than the number of columns (variables)? > mydata <- matrix(runif(12,-5,5), 4, 3) > mahalanobis(x=mydata, center=apply(mydata,2,mean), cov=var(mydata)) [1] 2.25 2.25 2.25 2.25 > mydata <- matrix(runif(420,-5,5), 21, 20) > mahalanobis(x=mydata,

Hello R helpers, I'm trying to use Mahalanobis distance to calculate distance of two time series, to make some comparations with euclidean distance, DTW, etc, but I'm having some dificults. I have, for example, two objects: s.1 <- c( 5.6324702, 1.3994353, -3.2572327, -3.8311846, -1.2248719, 0.9894694, -2.2835332, -5.1969285, -5.2823988, -3.1499400, -1.7307950, 2.8221209,

I have a data set that contains income data and a group identifier. Sort of like: DATA Group,Income A,2300 B,6776 A,6668 A,6768 B,9879 C,5577 A,7867 (etc),(etc) I am trying to compute the gini coefficient for each group. I have tried the following and none seem to do the trick: 1) attach(DATA) by(DATA, group, function(x) gini(income)) 2) attach(data) tapply(income, group,

Hi there, I am trying to compute Gini coefficients for vectors containing income classes. The data I possess look loke this: yit <- c(135, 164, 234, 369) piit <- c(367, 884, 341, 74 ) where yit is the vector of income classes, and fit is the vector of associated frequencies.(This data is from Rustichini, Ichino and Checci (Journal of Public Economics, 1999) ). In ineq pacakge, Gini( )

If I use the Ineq library and the Gini function in this way: >Gini(c(100,0,0,0)) I obtain the result 0.75 instead of 1 (that is the perfect inequality). I think Gini's formula in Ineq is based on a formula as reported here: http://mathworld.wolfram.com/GiniCoefficient.html but in the case of perfect inequality: x_1=.......=x_n-1 =0 x_n>0 these formula are equal to 1 - 1/n, not to

I am a new R user. i have a question about Mahalanobis distance.actually i have 300 rows and 7 columns. columns are different measurements, 300 rows are genes. since genes can classify into 4 categories. i used dist() with euclidean distance and cmdscale to do MDS plot. but find out Mahalanobis distance may be better. how do i use Mahalanobis() to generate similar dist object which i can use

Hi listers, Does it necessary to install any package in order to use the GINI or INEQ functions. If I use the following command the R tells me that didn't find the GINI function. x<-c(541, 1463, 2445, 3438, 4437, 5401, 6392, 8304, 11904, 22261) G<-gini(x) Thanks in advance, Marcio -- View this message in context:

Hi, help on 'mahalanobis' (in the stats package in Rv2.2.0) now says: "Description: Returns the Mahalanobis distance of all rows in 'x' and the vector mu='center' with respect to Sigma='cov'. This is (for vector 'x') defined as D^2 = (x - mu)' Sigma^{-1} (x - mu)" It does return D^2 as written. However,

Dear R ExpeRts, I'm having memory difficulties using mahalanobis distance to trying to cluster in R. I was wondering if anyone has done it with a matrix of 6525x17 (or something similar to that size). I have a matrix of 6525 genes and 17 samples. I have my R memory increased to the max and am still getting "cannot allocate vector of size" errors. My matrix "x" is

I want to calculate the probability that a group will include a particular point using the squared Mahalanobis distance to the centroid. I understand that the squared Mahalanobis distance is distributed as chi-squared but that for a small number of random samples from a multivariate normal population the Hotellings T2 (T squared) distribution should be used. I cannot find a function for

Hello, I would like to build a Lorenz curve and calculate a Gini coefficient in order to find how much parasites does the top 20% most infected hosts support. Here is my data set: Number of parasites per host: parasites = c(0,1,2,3,4,5,6,7,8,9,10) Number of hosts associated with each number of parasites given above: hosts = c(18,20,28,19,16,10,3,1,0,0,0) To represent the Lorenz curve: I

Dear all, I am working with randomForest package and I am interested in examining the "Gini importance" measures that are used as a general indicator of feature relevance. Is there a possibility of getting the Gini measure that is being estimated in each tree by the output of the getTree() function? Thanks a lot, Chrysanthi [[alternative HTML version deleted]]

Hello all, I am a 2 month newbie to R and am stumped. I have a data set that I've run multivariate stats on using the manova function (I included the data set). Now it comes time for a table of effect sizes with significance. The univariate tests are easy. Where I run into trouble filling in the table of effect sizes is the Mahalanobis D as an effect size. I've included the table so

Dear List, I am unsure if this is specifically a R question or a stats question? I thought i would ask here and if i get no replies it will answer that! I am trying to calculate Gini coefficients in R, based on a slight modification of the typical equation that i have seen in a paper. -------------- next part -------------- A non-text attachment was scrubbed... Name: PastedGraphic-2.pdf Type:

Hello all, I am a 2 month newbie to R and am stumped. I have a data set that I've run multivariate stats on using the manova function (I included the data set). Now it comes time for a table of effect sizes with significance. The univariate tests are easy. Where I run into trouble filling in the table of effect sizes is the Mahalanobis D as an effect size. I've included the table so