Abril 2020 - R-br - DINF UFPR

Re: [R-br] Busca de coordenadas de municípios usando RgoogleMaps
by sznelwar＠uol.com.br 15 Abr '20

15 Abr '20

1 0

Re: [R-br] Busca de coordenadas de municípios usando RgoogleMaps
by sznelwar＠uol.com.br 15 Abr '20

15 Abr '20

2 1

Re: [R-br] Busca de coordenadas de municípios usando RgoogleMaps
by verginio＠unesc.net 15 Abr '20

15 Abr '20

Boa tarde Flávio, Fiz um teste acrescentando o nome do país e separando por vírgula. O que ajuda a busca a entender que se trata de um município. library(RgoogleMaps) # Município, Estado, País aux_city <- c( "Vitória, ES, BR", "Serra, ES, BR", "Vila Velha, ES, BR", "Cariacica, ES, BR", "Colatina, ES, BR" ) aux_city_coord <- data.frame(Municípios = aux_city, lat = NA, long = NA) aux_city_coord <- with(aux_city_coord, data.frame(Municípios = aux_city_coord$Municípios, t( sapply(aux_city_coord$Municípios, getGeoCode) ))) aux_city_coord Atenciosamente, Max Richard Verginio _____________________________________ Sent from http://r-br.2285057.n4.nabble.com

2 1

Busca de coordenadas de municípios usando RgoogleMaps
by Flávio Fagundes 14 Abr '20

14 Abr '20

boa noite, colegas. Gostaria de uma ideia para o seguinte problema: Estou buscando as coordenadas de alguns municípios do estado de ES usando o seguinte código: library(RgoogleMaps) aux_city <- c("Vitória ES","Serra ES","Vila Velha ES","Cariacica ES","Colatina ES") aux_city_coord <- data.frame(Municípios=aux_city, lat=NA, long=NA) aux_city_coord <- with(aux_city_coord,data.frame(Municípios=aux_city_coord$Municípios, t(sapply(aux_city_coord$Municípios, getGeoCode)))) aux_city_coord Que retorna essas coordenadas. > aux_city_coord Municípios lat lon 1 Vitória ES 42.84651 -2.6724025 2 Serra ES 39.68388 -0.4288545 3 Vila Velha ES -20.32973 -40.2920284 4 Cariacica ES -20.26320 -40.4165490 5 Colatina ES -19.53467 -40.6287836 Notem que as coordenadas de Vitória e Serra estão totalmente erradas Porém, na tentativa e erro eu descobri que Vitória, não precisada do "ES" e retornou corretamente, mas Serra, sem chances de encontrar. aux_city <- c("Vitória","Serra","Serra-ES","Serra.ES") aux_city_coord <- data.frame(Municípios=aux_city, lat=NA, long=NA) aux_city_coord <- with(aux_city_coord,data.frame(Municípios=aux_city_coord$Municípios, t(sapply(aux_city_coord$Municípios, getGeoCode)))) aux_city_coord > aux_city_coord Municípios lat lon 1 Vitória -20.32009 -40.3376682 2 Serra 43.07201 1.1727739 3 Serra-ES 39.68388 -0.4288545 4 Serra.ES 39.68388 -0.4288545 Provavelmente poderá acontecer isso com outros municípios. Eu gostaria de saber se existe alguma uma forma de obter/consultar/verificar a grafia correta de, por exemplo, um ou todos os municípios de um determinado estado, sem que tenha que digitar o nome do mesmo. Mas qualquer outra forma de contornar este problema será bem vinda. Obrigado Att Flávio

1 0

Gráfico de network para o comércio internacional entre países
by Alexandre Loures 13 Abr '20

13 Abr '20

Boa noite, Listeiros! Gostaria de fazer um gráfico de rede para o comércio internacional entre um grupo de países. Nesse sentido, estou tentando replicar o seguinte exemplo: http://beyondvalence.blogspot.com/2014/05/text-mining-7-term-network-analys… Porém, no meu caso, possuo pares de países: ARE-JPN, AUS-JPN,... Ademais, gostaria que cada node fosse proporcional ao volume comercializado pelo país (variável exports). dput (df) structure(list(ido = c("ARE", "AUS", "CHN", "GBR", "GBR", "GBR", "HKG", "IDN", "IRN", "JPN", "JPN", "KOR", "KOR", "MYS", "NGA", "NLD", "NOR", "RUS", "RUS", "USA", "USA", "AUS", "AUS", "BEL", "BEL", "BEL", "CHN", "CHN", "CHN", "DEU", "DEU", "DEU", "DEU", "DEU", "DEU", "DEU", "GBR", "GBR", "GBR", "IRN", "MYS", "MYS", "NLD", "POL", "RUS", "RUS", "SGP", "SGP", "SGP", "SGP", "USA" ), idd = c("JPN", "JPN", "KOR", "CHE", "DEU", "FRA", "USA", "JPN", "JPN", "HKG", "SGP", "CHN", "JPN", "SGP", "USA", "FRA", "GBR", "DEU", "ITA", "GBR", "KOR", "CHN", "JPN", "DEU", "FRA", "NLD", "GBR", "RUS", "SGP", "AUT", "BEL", "ESP", "FRA", "NLD", "POL", "USA", "ARE", "IRL", "USA", "ITA", "SGP", "USA", "BEL", "DEU", "CHE", "TUR", "CHN", "HKG", "IDN", "MYS", "ARE"), exports = c(13112242813, 12600200044, 11292363934, 8415084770, 34788506146, 28475980395, 46984336771, 14415189665, 4968556482, 27183303751, 20820248756, 18454539579, 20466015819, 18045941449, 11499435728, 22555328417, 12375241280, 9231828795, 7254325411, 41578521231, 27901876438, 12302482046, 21703963590, 64547770665, 57659467585, 39244445982, 18976473583, 13211283269, 16632262456, 53765327000, 54292454000, 49796586000, 98354819000, 60942951000, 27729562000, 86050473000, 10567139812, 29981708162, 56706301526, 6191623821, 22078782553, 27815727284, 40914697144, 25222684658, 10781703404, 10841448818, 19757037697, 21522016786, 22102891530, 30384988869, 8476570775 )), class = "data.frame", row.names = c("1", "2", "3", "4", "5", "6", "7", "8", "9", "10", "11", "12", "13", "14", "15", "16", "17", "18", "19", "20", "21", "22", "23", "24", "25", "26", "27", "28", "29", "30", "31", "32", "33", "34", "35", "36", "37", "38", "39", "40", "41", "42", "43", "44", "45", "46", "47", "48", "49", "50", "51")) Desde já, muito obrigado pela atenção! -- *Alexandre Loures* Professor Substituto Universidade Federal da Paraíba <https://ufpb.br> Pós-graduação em Economia - PPGE <https://sigaa.ufpb.br/sigaa/public/programa/portal.jsf?lc=pt_BR&id=1875> www.rodriguesloures.com

1 0

Problemas na instalação de pacotes
by profiack＠gmail.com 11 Abr '20

11 Abr '20

Veja este link que me parece o mesmo caso. https://stackoverflow.com/questions/36346781/r-package-installation-issue <quote author='R-br mailing list'> Boa noite, Quando tento instalar pacotes : >chooseCRANmirror() Warning: failed to download mirrors file (não foi possível abrir a URL 'https://cran.r-project.org/CRAN_mirrors.csv'); using local file 'C:/PROGRA~1/R/R-36~1.3/doc/CRAN_mirrors.csv' Warning message: In download.file(url, destfile = f, quiet = TRUE) : falha em InternetOpenUrl: '' > utils:::menuInstallPkgs() Warning: unable to access index for repository https://cran-r.c3sl.ufpr.br/src/contrib: não foi possível abrir a URL 'https://cran-r.c3sl.ufpr.br/src/contrib/PACKAGES' Error in install.packages(lib = .libPaths()[1L], dependencies = NA, type = type) : argumento "pkgs" ausente, sem padrão Tentei instalar várias versões e todas estão apresentado este problema. Atenciosamente, Cleiton Fernando _______________________________________________ R-br mailing list R-br(a)listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível. </quote> Quoted from: http://r-br.2285057.n4.nabble.com/R-br-Problemas-na-instalacao-de-pacotes-t… _____________________________________ Sent from http://r-br.2285057.n4.nabble.com

1 0

Problemas na instalação de pacotes
by Cleiton Fernando Barbosa Brito 11 Abr '20

11 Abr '20

Boa noite, Quando tento instalar pacotes : >chooseCRANmirror() Warning: failed to download mirrors file (não foi possível abrir a URL 'https://cran.r-project.org/CRAN_mirrors.csv'); using local file 'C:/PROGRA~1/R/R-36~1.3/doc/CRAN_mirrors.csv' Warning message: In download.file(url, destfile = f, quiet = TRUE) : falha em InternetOpenUrl: '' > utils:::menuInstallPkgs() Warning: unable to access index for repository https://cran-r.c3sl.ufpr.br/src/contrib: não foi possível abrir a URL 'https://cran-r.c3sl.ufpr.br/src/contrib/PACKAGES' Error in install.packages(lib = .libPaths()[1L], dependencies = NA, type = type) : argumento "pkgs" ausente, sem padrão Tentei instalar várias versões e todas estão apresentado este problema. Atenciosamente, Cleiton Fernando

2 1

Re: [R-br] Subscribe - Erro ao exportar Dados para o Excel
by verginio＠unesc.net 04 Abr '20

04 Abr '20

Bom dia Márcio Boareto, Seu código funcionará, se alterar a linha: rownames(planilha) <- dados [,1] # nomeas das linhas (empresas) #por rownames(planilha) <- dados[[1]] # nomeas das linhas (empresas) ### O diferença está na library usada para leitura da planilha dados.xlsx. Foram carregadas três `xlsx`, `openxlsx`e `readxl`. Mas foi utilizada apenas esta última, que carrega como tibble. dados <- readxl::read_excel(path = "dados.xlsx", col_names = T) Atenciosamente, Max Richard Verginio _____________________________________ Sent from http://r-br.2285057.n4.nabble.com

2 1

Re: [R-br] Subscribe - Erro ao exportar Dados para o Excel
by Daniel Guimarães Tiezzi 04 Abr '20

04 Abr '20

Bom dia Marcio Segue o código e a planilha final. Eu só converti a sua table dados para o formato CSV pois não tenho instalado o pacote xlsx. Eu fiz isso no Excel mesmo. Veja se ficou OK. Daniel Daniel Tiezzi, MD, PhD Oncologia / Mastologia Professor Associado - Livre Docente Departamento de Ginecologia e Obstetrícia Setor de Mastologia e Oncologia Ginecológica Faculdade de Medicina de Ribeirão Preto - USP Tel.: 16 3602-2488 e-mail: <mailto:dtiezzi@usp.br>dtiezzi(a)usp.br <mailto:dtiezzi@usp.br> > On 4 Apr 2020, at 10:31, Marcio Boareto <marcioboareto(a)gmail.com> wrote: > > Prezado Prof Daniel > > > Ficaria muito agradecido se puder me ajudar. > Segue a planilha "dados", os resultados, o artigo e o script do R studio, se tiver como reproduzir. > > O artigo onde estou tentando reproduzir o modelo (está dando erro, pois o pacote de leitura da planilha foi descontinuado (xlsReadWrite) e estou tentando utilizar outro pacote para exportar o modelo para uma planilha. > > > > > Desde já agradeço vossa atenção. > > > > > Em sáb., 4 de abr. de 2020 às 10:22, Daniel Guimarães Tiezzi <dtiezzi(a)usp.br <mailto:dtiezzi@usp.br>> escreveu: > Ok, seria melhor VC enviar parte da sua planilha para podemos reproduzir o erro e resolver para VC. > > > On Sat, 4 Apr 2020, 10:17 Marcio Boareto, <marcioboareto(a)gmail.com <mailto:marcioboareto@gmail.com>> wrote: > Bom dia > > Prezado Prof Daniel > > Não consegui resolver e o erro ainda persiste. Estou aprendo recentemente a usar o R para cálculo de eficiência em análise envoltória de dados, estou com dificuldades em exportar os resultados para uma planilha. > > Muito obrigado pela ajuda! > > Em sáb., 4 de abr. de 2020 às 09:49, Daniel Guimarães Tiezzi <dtiezzi(a)usp.br <mailto:dtiezzi@usp.br>> escreveu: > Bom dia > > O erro gerado diz que a dimensão do objeto que você está tentando inserir como rownames do seu dataframe são diferentes. > > Verifique isso: > > > dim(planilha)[1] > > length(dados [ ,1]) > > Esses valoras devem ser iguais. > > Daniel > > > > Daniel Tiezzi, MD, PhD > Oncologia / Mastologia > Professor Associado - Livre Docente > Departamento de Ginecologia e Obstetrícia > Setor de Mastologia e Oncologia Ginecológica > Faculdade de Medicina de Ribeirão Preto - USP > Tel.: 16 3602-2488 > e-mail: dtiezzi(a)usp.br <mailto:dtiezzi@usp.br> > > > >> On 4 Apr 2020, at 09:02, Marcio Boareto por (R-br) <r-br(a)listas.c3sl.ufpr.br <mailto:r-br@listas.c3sl.ufpr.br>> wrote: >> >> Bom dia >> >> Preciso gerar uma planilha em excel dos dados gerados no modelo, porém está dando erro ao tentar colocar o nome das empresas na coluna 1 da planilha, segue o script: está dando erro no comando rownames >> >> #modelo retirado do artigo disponível em: https://www.aedb.br/seget/arquivos/artigos13/44218525.pdf <https://www.aedb.br/seget/arquivos/artigos13/44218525.pdf> >> > 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 >> >> Segue a planilha gerada, na coluna A preciso constar os nomes da coluna A (UTILITY) >> efficiency OPEX NETWORK MWH CUSTOMERS >> effcrs 0,58270216 8,00548E-10 0 6,07396E-09 5,6817E-08 >> efficiency 0,579460003 5,94412E-10 9,73842E-07 3,50295E-09 0 >> efficiency.1 0,885425355 2,0109E-09 0 1,52572E-08 1,42719E-07 >> efficiency.2 0,518585846 9,81483E-10 7,95216E-07 1,11913E-08 2,10933E-08 >> efficiency.3 0,775513183 1,79467E-09 0 1,36166E-08 1,27372E-07 >> efficiency.4 0,523933705 1,38609E-09 1,28427E-06 1,86596E-08 0 >> efficiency.5 1 2,29129E-09 1,85644E-06 2,61264E-08 4,92426E-08 >> efficiency.6 0,687899797 2,41195E-09 1,95421E-06 2,75023E-08 5,18359E-08 >> efficiency.7 1 5,10751E-09 0 7,6844E-08 0 >> efficiency.8 0,698435209 3,48636E-09 0 2,64519E-08 2,47436E-07 >> efficiency.9 0,82245379 2,85183E-09 0 2,16376E-08 2,02402E-07 >> efficiency.10 0,549836367 2,29078E-09 0 1,73807E-08 1,62583E-07 >> efficiency.11 0,585646 1,44619E-09 2,93566E-06 0 0 >> efficiency.12 1 5,36603E-09 8,79133E-06 3,16227E-08 0 >> efficiency.13 0,819757504 3,16289E-09 0 0 2,98655E-07 >> efficiency.14 0,669679648 4,14193E-09 3,35587E-06 4,72283E-08 8,90153E-08 >> efficiency.15 0,840551996 4,85147E-09 3,93075E-06 5,53189E-08 1,04264E-07 >> efficiency.16 0,44597084 2,59073E-09 2,09906E-06 2,95408E-08 5,56781E-08 >> >> Se puderem ajudar, fico agradecido. >> Atenciosamente >> >> Márcio Boareto >> > >> >> _______________________________________________ >> R-br mailing list >> R-br(a)listas.c3sl.ufpr.br <mailto:R-br@listas.c3sl.ufpr.br> >> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br <https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br> >> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia <http://www.leg.ufpr.br/r-br-guia>) e forneça código mínimo reproduzível. > > > > -- > Márcio Boareto > Administrador > CRA-GO 15993 > Mestre em Gestão Organizacional pela UFG Catalão. > Especialista em Gestão Financeira e Controladoria - UFG Catalão > Graduado em Administração de Empresas - UFG Catalão > 64 98137-9049 whatsapp > Lattes:http://lattes.cnpq.br/8524063708922235 <http://lattes.cnpq.br/8524063708922235> > ORCID: https://orcid.org/0000-0003-0398-8774 <https://orcid.org/0000-0003-0398-8774> > > > > > -- > Márcio Boareto > Administrador > CRA-GO 15993 > Mestre em Gestão Organizacional pela UFG Catalão. > Especialista em Gestão Financeira e Controladoria - UFG Catalão > Graduado em Administração de Empresas - UFG Catalão > 64 98137-9049 whatsapp > Lattes:http://lattes.cnpq.br/8524063708922235 <http://lattes.cnpq.br/8524063708922235> > ORCID: https://orcid.org/0000-0003-0398-8774 <https://orcid.org/0000-0003-0398-8774> > > > <resultscrs.xlsx><script -seget.R><dados.xlsx><44218525.pdf><script.docx>

1 0

Subscribe - Erro ao exportar Dados para o Excel
by Marcio Boareto 04 Abr '20

04 Abr '20

Bom dia Preciso gerar uma planilha em excel dos dados gerados no modelo, porém está dando erro ao tentar colocar o nome das empresas na coluna 1 da planilha, segue o script: está dando erro no comando rownames #modelo retirado do artigo disponível em: https://www.aedb.br/seget/arquivos/artigos13/44218525.pdf > 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 *1*249*143*613 *45*213 39*922*710 5*987*873 2 CEMIG *1*682*334*644 *460*219 37*476*802 6*832*546 3 CPFL - PAULISTA 497*290*782 *89*879 25*267*579 3*502*793 4 COPEL *1*018*866*491 *224*817 23*525*040 3*628*209 5 LIGHT 557*206*112 *58*074 22*902*552 3*640*182 6 CELESC 721*455*274 *144*896 18*105*811 2*237*127 7 COELBA 436*436*014 *215*001 14*286*757 4*622*046 8 ELEKTRO 414*602*018 *107*116 13*398*558 2*123*670 9 CPFL - PIRATININGA 195*789*961 *22*236 13*013*378 1*367*488 10 BANDEIRANTE 286*832*273 *27*496 12*536*237 1*482*518 11 CELPE 350*651*684 *120*428 10*001*560 2*994*259 12 AMPLA 436*532*756 *51*050 9*506*961 2*365*558 13 CELG 691*472*253 *199*494 9*344*291 2*213*198 14 RGE 186*357*415 *84*997 7*993*103 1*226*079 15 COELGE 316*166*876 *120*300 7*929*212 2*744*830 16 ESCELSA 241*433*335 *56*960 7*897*969 1*185*432 17 AES SUL 206*122*962 *76*133 7*616*460 1*150*518 18 CEEE 385*990*997 *71*892 7*277*929 1*438*072 > 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 Segue a planilha gerada, na coluna A preciso constar os nomes da coluna A (UTILITY) efficiency OPEX NETWORK MWH CUSTOMERS effcrs 0,58270216 8,00548E-10 0 6,07396E-09 5,6817E-08 efficiency 0,579460003 5,94412E-10 9,73842E-07 3,50295E-09 0 efficiency.1 0,885425355 2,0109E-09 0 1,52572E-08 1,42719E-07 efficiency.2 0,518585846 9,81483E-10 7,95216E-07 1,11913E-08 2,10933E-08 efficiency.3 0,775513183 1,79467E-09 0 1,36166E-08 1,27372E-07 efficiency.4 0,523933705 1,38609E-09 1,28427E-06 1,86596E-08 0 efficiency.5 1 2,29129E-09 1,85644E-06 2,61264E-08 4,92426E-08 efficiency.6 0,687899797 2,41195E-09 1,95421E-06 2,75023E-08 5,18359E-08 efficiency.7 1 5,10751E-09 0 7,6844E-08 0 efficiency.8 0,698435209 3,48636E-09 0 2,64519E-08 2,47436E-07 efficiency.9 0,82245379 2,85183E-09 0 2,16376E-08 2,02402E-07 efficiency.10 0,549836367 2,29078E-09 0 1,73807E-08 1,62583E-07 efficiency.11 0,585646 1,44619E-09 2,93566E-06 0 0 efficiency.12 1 5,36603E-09 8,79133E-06 3,16227E-08 0 efficiency.13 0,819757504 3,16289E-09 0 0 2,98655E-07 efficiency.14 0,669679648 4,14193E-09 3,35587E-06 4,72283E-08 8,90153E-08 efficiency.15 0,840551996 4,85147E-09 3,93075E-06 5,53189E-08 1,04264E-07 efficiency.16 0,44597084 2,59073E-09 2,09906E-06 2,95408E-08 5,56781E-08 Se puderem ajudar, fico agradecido. Atenciosamente Márcio Boareto >

2 2