library(ineq) my_minimum <- min(x) if(my_minimum < 0) { stop('Negative values are not allowed.') } myLength <- length(x) myMaximumEntropy <- log(myLength) mySum <- sum(x) myProportion <- x/mySum myEntropy <- -sum(myProportion * log(myProportion)) myNormalizedEntropy <- myEntropy / myMaximumEntropy myDifference <- myMaximumEntropy - myEntropy myTheilEntropyIndex <- entropy(x,parameter=1,na.rm=T) myExponentialIndex <- exp(-myEntropy) myHerfindahlMeasure <- sum(myProportion^2) myHerfindahl <- conc(x,type='Herfindahl',na.rm=T) myRosenbluth <- conc(x,type='Rosenbluth',na.rm=T) myNormalizedHerfindahlMeasure <- (myHerfindahlMeasure - 1/myLength) / (1 - 1/myLength) myGini <- Gini(x,na.rm=T) myConcentrationCoefficient <- myLength/(myLength -1)*myGini myRS <- RS(x,na.rm=T) myAtkinson <- Atkinson(x,na.rm=T) myKolm <- Kolm(x,na.rm=T) myCoefficientOfVariation <- var.coeff(x,square=F,na.rm=T) mySquaredCoefficientOfVariation <- var.coeff(x,square=T,na.rm=T) bitmap(file='plot1.png') plot(Lc(x)) grid() dev.off() bitmap(file='plot2.png') plot(Lc(x),general=T) grid() dev.off() load(file='createtable') a<-table.start() a<-table.row.start(a) a<-table.element(a,'Concentration - Ungrouped Data',2,TRUE) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'Measure',header=TRUE) a<-table.element(a,'Value',header=TRUE) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'Number of Categories',header=F) a<-table.element(a,myLength,header=F) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'Maximum Entropy',header=F) a<-table.element(a,myMaximumEntropy,header=F) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'Entropy',header=F) a<-table.element(a,myEntropy,header=F) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'Normalised Entropy',header=F) a<-table.element(a,myNormalizedEntropy,header=F) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'Max. Entropy - Entropy',header=F) a<-table.element(a,myDifference,header=F) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'Theil Entropy Index',header=F) a<-table.element(a,myTheilEntropyIndex,header=F) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'Exponential Index',header=F) a<-table.element(a,myExponentialIndex,header=F) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'Herfindahl',header=F) a<-table.element(a,myHerfindahl,header=F) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'Normalised Herfindahl',header=F) a<-table.element(a,myNormalizedHerfindahlMeasure,header=F) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'Rosenbluth',header=F) a<-table.element(a,myRosenbluth,header=F) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'Gini',header=F) a<-table.element(a,myGini,header=F) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'Concentration',header=F) a<-table.element(a,myConcentrationCoefficient,header=F) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'Ricci-Schutz (Pietra)',header=F) a<-table.element(a,myRS,header=F) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'Atkinson',header=F) a<-table.element(a,myAtkinson,header=F) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'Kolm',header=F) a<-table.element(a,myKolm,header=F) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'Coefficient of Variation',header=F) a<-table.element(a,myCoefficientOfVariation,header=F) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'Squared Coefficient of Variation',header=F) a<-table.element(a,mySquaredCoefficientOfVariation,header=F) a<-table.row.end(a) a<-table.end(a) table.save(a,file='mytable.tab')
|