similar to: Multiple Lorenz curves in one diagram - populations with different "n"

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

Hi group, I would like to draw multiple Lorenz curves in a single plot using data already prepared. Here is a simple example: require("lawstat") lorenz.curve(c(1,2,3),c(4,5,4)) lorenz.curve(c(1,2,3),c(4,2,1)) This example draws two separate graphs. How can I combine them in a distinguishable way? I tried ?polygon without success... Ralf

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

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:

Dear All, I wish to calculate the Gini index (ineq from same package) and some other indices for the diameter distribution of each plot (df dgtot). dgtot: IDPlot Diameter(cm) 1 4 34.0 2 4 23.0 3 4 38.0 ... 51 5 16.0 52 5 8.0 53 5 9.0 ... 5301 140 25.0 5302 140 12.0 5303 140 7.0 I use: >

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

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 came across the package 'ineq' that computes a variety of inequality measures (e.g. gini, theil etc). I want to compute the Theil index (racial segregation) and decompose the total into sub-components (by geog levels). I think the package doesn't report the decomposition (correct me if I'm wrong). Just wonder is that available elsewhere? K.

Hi, I'm using the ineq package to calculate some concentration measures (Gini, Herfindal, ...) and I was wondering if there's around also a function to calculate standard error on these measures. If not, is anybody aware of where I can find a reference on this point? Thanks. -- ======================================================== Angelo Secchi PGP Key ID:EA280337

Ok, Thanks foe the answer, Ken: *1L, 2L etc are integers. (That is, identical to as.integer(1) , as.integer(2) etc) Using integers (instead of "numeric" type) is more efficient as here they're used as indexes and would be converted to integer anyway. Compare > is(1) ... and > is(1L) 1L:p is the sequence 1, 2, 3, ..., p (just like 1:p) Just for curiosity, what is

Hello all, I've been using binconf (package Hmisc) at a range of alpha values and noticed that using the 'Wilson' method when alpha is larger (i.e. narrow CI), results in the upper value being smaller than the lower value. The 'exact' and 'asymptotic' methods give results in the realm I'd expect. But the help file suggests: "Following Agresti and Coull, the

Dear All, I have developed a programme the anable the decomposition of the Gini index, it complets tha valuable work of Achim Zeileis, the author of the ineq package. I would like to make it to be part of all R package. How should I proceed. Must I sent it to the the Core developement team ? The proogramme is written in R. Many thanks for your advice, Best regards, Souleymane

I'm trying to call BUGS from R. But it's not working. R freezes up and BUGS gives me a strange output in the log. Just to know, BUGS is registered. The modified date on the keys file is today (Dec. 9th). It should be fully registered so that I can use it fully. And, the BUGS model is syntactically correct. Any suggestions would be very helpful. Here is my BUGS model: model {

Hello, I've come across a problem in developing a set of custom functions to calculate the number of hours of daylight at a given latitude, and the number of days a date precedes or secedes the summer solstice. I discovered an inconsistency concerning leap years between my derived values and those from the US naval databases. It seems as far as I can figure that my inconsistency arises either

Dear useRs, my last announcement is an update of the ineq package for measuring inequality, concentration and poverty. The current version is now 0.2-4. Thanks to suggestions from Rein Halbersma the Pen() function for plotting Pen's parade was improved and now allows for much more flexibility. See the help page for examples. Best wishes, Z ------------- Package: ineq Version: 0.2-4 Date:

Dear useRs, my last announcement is an update of the ineq package for measuring inequality, concentration and poverty. The current version is now 0.2-4. Thanks to suggestions from Rein Halbersma the Pen() function for plotting Pen's parade was improved and now allows for much more flexibility. See the help page for examples. Best wishes, Z ------------- Package: ineq Version: 0.2-4 Date:

Dear All, I am actually on R-devel using shell consol. When I run RStudio, it can't find libR.so in the new /lib folder where is libRblas.so libRlapack.so. At first step I configure R-devel to share library with ./configure --enable-R-shlib. when I copied libR.so from /usr/lib/R/lib/libR.so (stable version R 3.1) to /usr/local/R-devel/lib/libR.so (R-devel version), That doesn't work.

Dear All, I am actually on R-devel using shell consol. When I run RStudio, it can't find libR.so in the new /lib folder where is libRblas.so libRlapack.so. At first step I configure R-devel to share library with ./configure --enable-R-shlib. when I copied libR.so from /usr/lib/R/lib/libR.so (stable version R 3.1) to /usr/local/R-devel/lib/libR.so (R-devel version), That doesn't work.

I'm writing a paper aimed at motivating the use of GEE within the field of economics. However, after computing using the geeglm function, I noticed there's one intercept in the summary output. I assume this means the function is pooling the data. That means my code is not what I want. I want a "fixed effects" model, meaning I want the intercept to vary by cluster. Here's

I'm sure this must be trivial, but I'm a novice with R and can't work out how to handle the axes when I am constructing multiple plots on a page and try to return to a previous one to put multiple data sets it. A simple example: --- x<- 1:10 y<- (1:100)*3 par(mfcol=c(2,1)) plot(x) plot(y) par(mfg=c(1,1)) lines(x) --- The first 5 lines make two plots with a row of dots along