> require(lpSolve) # carrega pacote
lpSolve previamente instalado Carregando pacotes exigidos:
lpSolve > require(xlsx) # carrega pacote
xlsReadWrite previamente instalado Carregando pacotes exigidos: xlsx > setwd('c:/Users/Usuário/UFG/MESTRADO
2017/MESTRADO 2019/R') # define o diretório de trabalho >
library(readxl) >
library(openxlsx) Attaching package: ‘openxlsx’ The following objects are masked from ‘package:xlsx’: createWorkbook,
loadWorkbook, read.xlsx, saveWorkbook, write.xlsx > library(xlsx) > dados <-
read_excel("C:/Users/Usuário/UFG/MESTRADO 2017/MESTRADO
2019/R/dados.xlsx", col_names=TRUE) > View(dados) > dados # A tibble: 18 x 5
UTILITY OPEX NETWORK MWH CUSTOMERS <chr> <dbl> <dbl> <dbl> <dbl> 1
ELETROPAULO 1249143613 45213 39922710 5987873 2 CEMIG 1682334644 460219 37476802 6832546 3 CPFL - PAULISTA 497290782 89879 25267579 3502793 4 COPEL 1018866491 224817 23525040 3628209 5 LIGHT 557206112 58074 22902552 3640182 6 CELESC 721455274 144896 18105811 2237127 7 COELBA 436436014 215001 14286757 4622046 8 ELEKTRO 414602018 107116 13398558 2123670 9 CPFL - PIRATININGA 195789961 22236 13013378 1367488 10 BANDEIRANTE 286832273 27496 12536237 1482518 11 CELPE 350651684 120428 10001560 2994259 12 AMPLA 436532756 51050 9506961 2365558 13 CELG 691472253 199494 9344291 2213198 14 RGE 186357415 84997 7993103 1226079 15 COELGE 316166876 120300 7929212 2744830 16 ESCELSA 241433335 56960 7897969 1185432 17 AES SUL 206122962 76133 7616460 1150518 18 CEEE 385990997 71892 7277929 1438072 > inputs<- data.frame(dados[,2])
# seleciona coluna 2 do objeto data > outputs <-
data.frame(dados[,c(3,4,5)]) # seleciona colunas 3, 4 e 5 do objeto data > N<- dim(dados)[1] # número de
DMUs > s <- dim(inputs)[2] # número
de inputs > m <- dim(outputs)[2] # número
de outputs >
f.rhs <- c(rep(0,1,N),1) # RHS > f.dir <-
c(rep("<=",1,N),"=") # direção das restrições > aux <-
cbind(-1*inputs,outputs) # matriz com os coeficientes tecnológicos em (6) > for (i in 1:N) { + f.obj <-
c(0*rep(1,s),as.numeric(outputs[i,])) # coeficientes da função objetivo + f.con <- rbind(aux
,c(as.numeric(inputs[i,]), + rep(0,1,m))) # adiciona a restrição bTz=1 + results <-
lp("max",as.numeric(f.obj),f.con,f.dir,f.rhs,scale=1,compute.sens=TRUE)
# resolve PPL +
multipliers <- results$solution #
multiplicadores +
efficiency <- results$objval #
indices de eficiência +
duals <- results$duals # ??? +
if (i==1) { +
weights <- multipliers +
effcrs <- efficiency +
lambdas <- duals [seq(1,N)] +
} else { +
weights <-
rbind(weights,multipliers) +
effcrs <- rbind(effcrs ,
efficiency) +
lambdas <-
rbind(lambdas,duals[seq(1,N)]) + } + } > ################################################################### > planilha <-
cbind(effcrs,weights) > #ESSA PARTE NÃO ESTÁ
FUNCIONANDO################################### > rownames(planilha) <- dados
[,1] # nomeas das linhas (empresas) Error in dimnames(x) <- dn : comprimento
de 'dimnames' [1] não é igual ao tamanho do array
|