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