- R-br - DINF UFPR

Criar um VETOR de um dataframe
by Manasses Nóbrega 12 Jun '18

12 Jun '18

Olá, pessoal! Este parece ser um problema fácil de resolver, mas não estou conseguindo!! Tenho um conjunto de dados (com quase 3000 linhas) que relaciona a resposta de um voluntário sob várias condições experimentais, cada uma dela com seus respectivos níveis. Exemplo: resp V1 V2 V3 0.35 1 1 1 0.38 1 2 1 0.65 1 3 1 0.49 1 4 1 -0.25 0 1 2 -1 0 2 2 -0.73 0 3 2 -0.47 0 4 2 Para realizar uma análise, por exemplo, vou precisar de *um vetor* com as respostas apenas para o nível 1 de V2 ou, por exemplo, de* um vetor* com as respostas apenas para o nível 2 de V3. Tentei a função aggregate mas não consegui. Vocês poderiam me ajudar? Desde já agradeço a todos. Um abraço, Manasses -- *I am a brain, Watson. The rest of me is a mere appendix. * *(Sherlock Holmes)*

3 3

Re: [R-br] como usar um layout gráfico do tipo 3x2 ? (Gerson Jr)
by Rafael Tieppo 12 Jun '18

12 Jun '18

a função grid.arrange() pode ajudar: library(gridExtra)grid.arrange(plot1, plot2) mais exemplos: https://www.r-graph-gallery.com/261-multiple-graphs-on-same-page/ Rafael Tieppo State University of Mato Grosso - Department of Agricultural Engineering site: https://rafatieppo.github.io "Evite o desperdício: antes de imprimir pense na sua responsabilidade com o ambiente".

1 0

como usar um layout gráfico do tipo 3x2 ?
by Cleber N.Borges 11 Jun '18

11 Jun '18

Olá a todos, Gostaria de fazer 5 gráficos em uma mesma janela... mas não entendi como usar as funções de layout... Gostaria de algo mais ou menos como coloco abaixo: ######## # primeira linha na janela gráfica => 3 gráficos plot( rnorm(10), col=1 ); plot( rnorm(10), col=2 ); plot( rnorm(10), col=3 ); # segunda linha na janela gráfica => 2 gráficos hist( rnorm(10), col=1 ); hist( rnorm(10), col=2 ); plot( rnorm(10), col=3 ); agradeço antecipadamente pela ajuda Cleber Borges --- Este email foi escaneado pelo Avast antivírus. https://www.avast.com/antivirus

3 5

Converter character em numeric
by Edmar Caldas 06 Jun '18

06 Jun '18

Boa Tarde! estou tentando converter uma variável character em numerica e não estou conseguindo, alguem pode me ajudar por favor? > typeof(mes_ano)[1] "character"> agregar<-as.numeric(as.character(mes_ano))> is.numeric(mes_ano)[1] FALSE Edmar

2 2

linhas de densidade com legenda
by Alice Peres 05 Jun '18

05 Jun '18

Bom dia, caros membros Estou tendo dificuldade, vocês podem me ajudar a a transformar o gráfico abaixo em 4 linhas de densidade com legenda? x=c(3.42,4.01,4.04,4.26,4.10,5.61,4.90,3.67,5.30,5.05,4.44,3.50,4.09,5.05,4.40,3.99,3.86, 4.99,2.86,2.18,3.61,4.45,5.01,3.11,3.53,4.70,4.03,3.23,2.87,3.38,3.04,4.59,1.45,3.91,5.46,5.70, 3.33,3.38,5.43,3.33,4.37,3.83,3.45,5.74,3.64,3.23,4.32,2.91,4.69,4.04,3.10,3.29,4.04,3.41,4.57, 2.97,5.46,4.02,5.62,4.94,3.23,4.39,3.67,3.99,2.60,3.05,2.81,2.61,3.06,4.21,3.19,2.30,2.84,2.60, 4.32,2.55,3.23,3.66,3.26,3.90,4.60,3.99,3.42,4.53,3.20,3.46,3.19) hist(rexp(1000,1/min(x)),freq=FALSE) hist(rgamma(1000,shape=mean(x)^2,rate=mean(x)),freq=FALSE,add=TRUE) hist(rexp(1000,1/mean(x)),freq=FALSE,add=TRUE) hist(rexp(1000,1/max(x)),freq=FALSE,add=TRUE)

2 1

Re: [R-br] ggplot2_apresentar equação no gráfico
by cayanportela＠hotmail.com 05 Jun '18

05 Jun '18

Você pode gerar um modelo utilizando a função lm e extrair os coeficientes. De acordo com https://stackoverflow.com/questions/30144785/how-to-extract-intercept-and-s… , no seu caso ficaria assim: db %>% group_by(type) %>% do({ mod = lm(y ~ exp(-x), data = .) data.frame(Intercept = coef(mod)[1], Slope = coef(mod)[2]) }) Cayan _____________________________________ Sent from http://r-br.2285057.n4.nabble.com

3 3

Erro na análise de regressão
by Angélica Ricarte 04 Jun '18

04 Jun '18

Prezados, boa noite. Estou com a seguinte dúvida, se puderem me ajudar, agradeço demais: Utilizando o pacote "ExpDes.pt", fiz a seguinte análise: >dados Oleo S 1 0.0 85.71 2 0.0 85.71 3 0.0 85.71 4 0.0 85.71 5 0.0 85.71 6 0.5 85.71 7 0.5 85.71 8 0.5 85.71 9 0.5 85.71 10 0.5 85.71 11 1.0 100.00 12 1.0 100.00 13 1.0 100.00 14 1.0 100.00 15 1.0 100.00 16 1.5 71.42 17 1.5 71.42 18 1.5 71.42 19 1.5 71.42 20 1.5 71.42 > dic(dados[,1],dados[,2],quali=F,sigF=0.05) ------------------------------------------------------------------------ Quadro da analise de variancia ------------------------------------------------------------------------ GL SQ QM Fc Pr>Fc Tratamento 3 2042 680.68 4.4071e+30 1.5357e-239 Residuo 16 0 0.00 Total 19 2042 ------------------------------------------------------------------------ CV = 0 % ------------------------------------------------------------------------ Teste de normalidade dos residuos Valor-p: 5.640644e-07 ATENCAO: a 5% de significancia, os residuos nao podem ser considerados normais! ------------------------------------------------------------------------ ------------------------------------------------------------------------ Teste de homogeneidade de variancia valor-p: NaN Error in if (pvalor.hvar < 0.05) { : valor ausente onde TRUE/FALSE necessário > dic(dados[,1],dados[,2],quali=F,sigF=0.05) Porque não consigo obter um resultado para a análise de regressão? Seria porque o CV é igual a zero e porque não há variância? Se sim, poderiam me dar uma explicação sobre a relação entre a variância e a regressão? Também tentei fazer o mesmo considerando os tratamentos qualitativos (quali=T) para obter um teste de médias, porém aparece a mesma mensagem de erro. OBS: para outras variáveis (que possuem variância e CV diferente de zero) deu certo a análise. Grata, aguardo retorno. Att., Angélica Batista.

2 3

IBGE - Distorção Malha de Setores Censitários Censo 2010
by Bruno Rezende 04 Jun '18

04 Jun '18

Prezados, 1. Estou trabalhando com a malha de setores censitários fornecida pelo IBGE (CENSO 2010). 2. Inicialmente, obtive a malha em formato .shp e carreguei no R. Notei que, ao projetar essa malha no Google Earth e no Leaflet, o desenho dos polígonos não coincidia com o desenho dos mapas. 3. Pensando que a distorção acima pudesse ser causada por algum erro meu, baixei o arquivo Google Earth .kmz dos setores censitários, fornecido diretamente pelo IBGE. 4. Contudo, os problemas são os mesmos: para minha cidade, algumas partes dos mapas se sobrepõem e algumas partes destoam do mapa "real". 5. Sabemos que um setor censitário é uma subdivisão do território com a finalidade de delimitar o trabalho do agente censeador a cerca de 250 a 350 domicílios. Uma divisa de setor censitário que cruze o meio de uma quadra não permite clareza na definição do setor para o agente que está em campo. A lógica infere que utilizar os próprios contornos do terreno, ou equipamentos públicos, facilita a divisão e clareza in loco. 6. Assim, deseja-se que, na maioria das vezes, as divisas do setor censitário e dos equipamentos públicos coincidam. 7. O problema é que o mapa dos setores não "casa" com o mapa "real". Assim, coletar uma coordenada em campo e projetar no shape file de setor censitário apresentará erro, inviabilizando análises. 8. Estou analisando a cidade de Passo Fundo - RS. Em Porto Alegre, capital do estado, essas distorções aparentam não ocorrer. 9. Tentei algumas transformações do CRS, mas sem sucesso. 10. Senhores, alguém de vocês já passo por problema semelhante? Saberiam indicar como corrigí-lo?

1 0

OpenBUGS no R
by Manasses Nóbrega 04 Jun '18

04 Jun '18

Prezados, Estou iniciando meus estudos em MCMC e uso do OpenBUGS no R. Ao tentar rodar um pequeno exemplo inicial (do livro do Kinas & Andrade, Introdução à Inferência Bayesiana com R), apareceu uma mensagem de erro que não consegui resolver. Abaixo segue o código que estou usando e as mensagens de erro que apareceram: lddados <- list(n=4,x=0) params <- list("teta") inits <- list(teta=0.5) sink("ru486_mcmc.txt") cat("model { teta ~ dunif(0,1) x ~ dbin(teta,n) }") sink() library(R2OpenBUGS) library(R2WinBUGS) library(BRugs) *Error :* .onLoad falhou em loadNamespace() para 'BRugs', detalhes: chamada: if (!is.na(ob$version)) msg <- paste(msg, "version", ob$version) else msg <- paste(msg, erro: argument is of length zero In addition: Warning messages: 1: package ‘BRugs’ was built under R version 3.3.3 2: In findOpenBUGS() : OpenBUGS 3.2.2 or greater must be installed (if so, this indicates missing registry keys of OpenBUGS). Setting the environment variable 'OpenBUGS_PATH' in advance of loading 'BRugs' overwrites the path. See ?loadOpenBUGS in order to load OpenBUGS manually. 3: In is.na(ob$version) : is.na() applied to non-(list or vector) of type 'NULL' Error: package or namespace load failed for ‘BRugs’ fit <- bugs(lddados,inits,params,model.file = "ru486_mcmc.txt", n.chains = 1,n.iter = 25000,n.burnin = 20000,n.thin = 5) *Error in* bugs(lddados, inits, params, model.file = "ru486_mcmc.txt", n.chains = 1, : Cannot find the OpenBUGS program Estou usando a versão 3.5 do R, e já tentei nas formas 32 e 64 bit. Também tentei o loadOpenBUGS e não tive êxito! Agradeço a ajuda de vocês na resolução desses erros. Um abraço a todos, Manasses -- *I am a brain, Watson. The rest of me is a mere appendix. * *(Sherlock Holmes)*

2 1

Re: [R-br] Digest R-br, volume 90, assunto 3
by Cayan Portela 04 Jun '18

04 Jun '18

Você pode extrair os coeficientes utilizando a função lm. De acordo com https://stackoverflow.com/questions/30144785/how-to-extract-intercept-and-s… ,no seu caso ficaria: db %>% group_by(type) %>% do({ mod = lm(y ~ exp(-x), data = .) data.frame(Intercept = coef(mod)[1], Slope = coef(mod)[2]) }) ________________________________ De: R-br <r-br-bounces(a)listas.c3sl.ufpr.br> em nome de r-br-request(a)listas.c3sl.ufpr.br <r-br-request(a)listas.c3sl.ufpr.br> Enviado: segunda-feira, 4 de junho de 2018 15:00 Para: r-br(a)listas.c3sl.ufpr.br Assunto: Digest R-br, volume 90, assunto 3 Enviar submissões para a lista de discussão R-br para r-br(a)listas.c3sl.ufpr.br Para se cadastrar ou descadastrar via WWW, visite o endereço https://nam04.safelinks.protection.outlook.com/?url=https%3A%2F%2Flistas.in… ou, via email, envie uma mensagem com a palavra 'help' no assunto ou corpo da mensagem para r-br-request(a)listas.c3sl.ufpr.br Você poderá entrar em contato com a pessoa que gerencia a lista pelo endereço r-br-owner(a)listas.c3sl.ufpr.br Quando responder, por favor edite sua linha Assunto assim ela será mais específica que "Re: Contents of R-br digest..." Tópicos de Hoje: 1. ggplot2_apresentar equação no gráfico (Henrique Boriolo Dias) ---------------------------------------------------------------------- Message: 1 Date: Mon, 4 Jun 2018 08:50:07 +0000 (UTC) From: Henrique Boriolo Dias <henrique.bdias(a)yahoo.com.br> To: "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br> Subject: [R-br] ggplot2_apresentar equação no gráfico Message-ID: <1597288900.8454811.1528102207591(a)mail.yahoo.com> Content-Type: text/plain; charset="utf-8" Olá pessoal, Sabem se é possível plotar a equação especifica (e talvez o R2) gerada pelo geom_smooth quando se usa a opção facet_wrap?Estou interessado no valor que multiplica o 'x' da equação.Alguém poderia me dar uma luz de qual o critério é usado para o ajuste dessa equação no geom_smooth? - Dados em anexo.- Script: db <- read.csv('data.csv', sep = ';', dec = '.', header = T) tiff(filename=paste0("graph.tiff"), width = 2000, height = 1000, res=300)ggplot(db, aes(y=y, x=x)) + geom_point(aes(shape=type))+ facet_wrap(~type) + geom_smooth(method = lm, formula = y ~ exp(-x), se = FALSE)dev.off() ObrigadoHenrique -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <https://nam04.safelinks.protection.outlook.com/?url=http%3A%2F%2Flistas.inf…> -------------- Próxima Parte ---------- Um anexo não-texto foi limpo... Nome: data.csv Tipo: application/vnd.ms-excel Tamanho: 5000 bytes Descrição: não disponível URL: <https://nam04.safelinks.protection.outlook.com/?url=http%3A%2F%2Flistas.inf…> ------------------------------ Subject: Legenda do Digest _______________________________________________ R-br mailing list R-br(a)listas.c3sl.ufpr.br https://nam04.safelinks.protection.outlook.com/?url=https%3A%2F%2Flistas.in… ------------------------------ Fim da Digest R-br, volume 90, assunto 3 ****************************************

1 0

ggplot2_apresentar equação no gráfico
by Henrique Boriolo Dias 04 Jun '18

04 Jun '18

Olá pessoal, Sabem se é possível plotar a equação especifica (e talvez o R2) gerada pelo geom_smooth quando se usa a opção facet_wrap?Estou interessado no valor que multiplica o 'x' da equação.Alguém poderia me dar uma luz de qual o critério é usado para o ajuste dessa equação no geom_smooth? - Dados em anexo.- Script: db <- read.csv('data.csv', sep = ';', dec = '.', header = T) tiff(filename=paste0("graph.tiff"), width = 2000, height = 1000, res=300)ggplot(db, aes(y=y, x=x)) + geom_point(aes(shape=type))+ facet_wrap(~type) + geom_smooth(method = lm, formula = y ~ exp(-x), se = FALSE)dev.off() ObrigadoHenrique

1 0

Ajuda com amostra pequena
by Elias Carvalho 01 Jun '18

01 Jun '18

Bom dia Recebi um conjunto de dados (summary abaixo) para análise, porém vejo alguns problemas: 1. Amostra muito pequena, apenas 17 registros 2. Algumas variáveis com categorias ocorrendo apenas um ou duas vezes como sexo (SEX) por exemplo 3. Variável idade (IDA) e numero de alunos na sala (NAL) com muitas categoria 4. Variável gerenciamento (GBY) não possui variabilidade Esse é um banco de dados sobre bullying na escola, preciso analisar a capacidade do coordenador em definir bullying (DBY) reconhecer bullying (RBY) e gerenciar bullying (GBY) Alguem pode sugerir que tipo de ações eu poderia executar para conseguir fazer alguma análise sobre esses dados ou me dizer se é impossível fazer isso? -- *In Jesu et Maria* *Obrigado* *Prof. Elias Carvalho* *"Felix, qui potuit rerum cognoscere causas" (Virgil 29 BC)"Blessed is he who has been able to understand the cause of things"*

3 2

INDICADORES SOCIAIS OF TOPICS
by Edson Lira 01 Jun '18

01 Jun '18

Caros amigos, boa tarde, um amigo meu Economista me fez uma pergunta que resolvi compartilhar com vocês. "Existe alguma técnica em Estatística, específica para tratamento de Indicadores Sociais?" Eu diria que não, mas gostaria de saber se vocês conhecem algo do tipo. [ ['s.Prof. Edson Lira, Me Estatístico Manaus-Amazonas

4 4

quebra de linha do título de tabelas geradas pela função xtable
by Paulo Eduardo de Mesquita 29 Mai '18

29 Mai '18

Bom noite (dia) a todo(a)s, Paz e bem Sou um neófito em R e Latex. Estou tentando redigir um artigo relativamente longo no formato Rnw no Rstudio e esbarrei num problema que não estou sendo capaz de superar, mesmo depois de pesquisar bastante. O problema é ajustar a largura do título das tabelas geradas com a função xtable à largura da tabela propriamente dita através da inserção de quebras de linha no título, fazendo com que algo assim (que o código que escrevi á capaz de gerar): se transforme em algo assim (que eu editei como figura no MS - Paint): Um exemplo de trabalho simplificado segue anexo. Muito obrigado. -- Paulo Eduardo de Mesquita Disciplina de Infectologia - Faculdade de Medicina Universidade do Oeste Paulista Presidente Prudente - São Paulo - Brasil

4 4

Tutorial de RMarkdown
by Marcos Vital 29 Mai '18

29 Mai '18

Oi, pessoal da lista! Publiquei no meu blog um pequeno tutorial de RMarkdown, e acredito que possa ser útil para os que estão começando a aprender a usar esta ferramenta: https://cantinhodor.wordpress.com/2018/05/26/tutorial-de-rmarkdown/ Caso alguém tiver alguma sugestão, estou à disposição para corrigir, expandir ou atualizar o tutorial, ok? Abraços! -- Marcos Vinícius Carneiro Vital Universidade Federal de Alagoas Instituto de Ciências Biológicas e da Saúde Setor de Biodiversidade e Ecologia

1 0

Error: argument "ast" is missing, with no default (solve_model)
by Adriele Giaretta Biase 29 Mai '18

29 Mai '18

Boa tarde pessoal, estou trabalhando com modelos matemáticos de otimização usando o pacote ompr. A ultima linha do código (azul) retorna um erro que não estou conseguindo tratar. Copiei o erro em baixo (vermelho). Alguém já teve a experiência de encontra esse mesmo erro? rm(list=ls(all=TRUE)) #install package from github #devtools::install_github("dirkschumacher/ompr") #devtools::install_github("dirkschumacher/ompr.roi") #install.packages("ROI.plugin.glpk") library(dplyr, warn.conflicts = FALSE) library(ROI) library(ROI.plugin.glpk) library(ompr) library(ompr.roi) frig = rep(seq(1, 4,1), rep((150*40),4) ) carro = rep(rep(seq(1, 40, 1), rep(150, 40)), 4) dias = rep( seq(1:150), 4*40) lucro = (rnorm(4*40*150, 30,1)) merc = rep(1, length(150*40*4)) praca = rep(1, length(150*40*4)) conj_dados= cbind(frig, carro, dias, lucro, merc, praca) DIF = 150 Trava_param <- function(DIF){ dia <- seq(1, DIF, 1) trav <- rep(0, DIF) praca <- rep(1, DIF) Trava <- data.frame(dias = dia, praca = praca, Trava = trav) return(Trava) } Trava = Trava_param(DIF) Max_descasque <- 3 capac_cam <- 60 Carga_mercado_futuro <- 330 permissao_venda <- 1 bigM <- 5000 model <- MIPModel() %>% # venda do animal i para o frigorifico j no mercado m, # dado o peso no dia k add_variable( x[i,j,k,m], i = as.numeric(as.vector(levels( as.factor(conj_dados[,"carro"])))), #1:carro, j = as.numeric(as.vector(levels( as.factor(conj_dados[,"frig"])))), #1:frig, k = as.numeric(as.vector(levels( as.factor(conj_dados[,"dias"])))), m = as.numeric(as.vector(levels( as.factor(conj_dados[,"merc"])))),#1:merc, type = "binary" )%>% # indica a venda de animais no dia k add_variable( y[k], k = as.numeric(as.vector(levels( as.factor(conj_dados[,"dias"])))),#1:dias, type = "binary" )%>% # numero de caminhoes do tipo t indo para o frigorifico j no dia k add_variable( total[j,k], j = as.numeric(as.vector(levels( as.factor(conj_dados[,"frig"])))), #1:frig, k = as.numeric(as.vector(levels( as.factor(conj_dados[,"dias"])))), #1:dias, type = "integer", lb = 0 )%>% # indica qual mercado esta sendo usado add_variable( w[k,m], k = as.numeric(as.vector(levels( as.factor(conj_dados[,"dias"])))),#1:dias, m = as.numeric(as.vector(levels( as.factor(conj_dados[,"merc"])))), #1:merc, type = "binary" )%>% # indica se o boi i foi vendido na praca u no dia k add_variable( t[i,u,k], i = as.numeric(as.vector(levels( as.factor(conj_dados[,"carro"])))), # 1:carro, u = as.numeric(as.vector(levels( as.factor(conj_dados[,"praca"])))), #1:praca, k = as.numeric(as.vector(levels( as.factor(conj_dados[,"dias"])))), #1:dias, type = "binary" )%>% # maximizar set_objective( sum_expr( x[i,j,k,m]*lucro[i,j,k,m], i = as.numeric(as.vector(levels( as.factor(conj_dados[,"carro"])))),# 1:carro, j = as.numeric(as.vector(levels( as.factor(conj_dados[,"frig"])))),# 1: frig, k = as.numeric(as.vector(levels( as.factor(conj_dados[,"dias"])))),#1:dias, m = as.numeric(as.vector(levels( as.factor(conj_dados[,"merc"])))) ) #1:merc) - sum_expr( y[k], k = as.numeric(as.vector(levels( as.factor(conj_dados[,"dias"])))) )#1:dias) + sum_expr(130*t[i,u,k]*15, i = as.numeric(as.vector(levels( as.factor(conj_dados[,"carro"])))), #1:carro, u = as.numeric(as.vector(levels( as.factor(conj_dados[,"praca"])))), #1:praca, k = as.numeric(as.vector(levels( as.factor(conj_dados[,"dias"])))) )#1:dias) ) %>% #Um animal eh vendido uma unica vez add_constraint( sum_expr( x[i,j,k,m], i = as.numeric(as.vector(levels( as.factor(conj_dados[,"carro"])))),#1:carro, j = as.numeric(as.vector(levels( as.factor(conj_dados[,"frig"])))),#1: frig, k = as.numeric(as.vector(levels( as.factor(conj_dados[,"dias"])))), # 1:dias, m = as.numeric(as.vector(levels( as.factor(conj_dados[,"merc"])))))#1:merc) + sum_expr( t[i,u,k], i = as.numeric(as.vector(levels( as.factor(conj_dados[,"carro"])))),# 1:carro, u = as.numeric(as.vector(levels( as.factor(conj_dados[,"praca"])))),# 1: praca, k = as.numeric(as.vector(levels( as.factor(conj_dados[,"dias"])))))#1:dias, <= 1 )%>% # A capacidade do caminhao eh respeitada add_constraint( sum_expr( x[i,j,k,m], i = as.numeric(as.vector(levels( as.factor(conj_dados[,"carro"])))),#1:carro, m = as.numeric(as.vector(levels( as.factor(conj_dados[,"merc"])))))#1:merc) <= total[j,k] * capac_cam, j = as.numeric(as.vector(levels( as.factor(conj_dados[,"frig"])))),#1:frig, k = as.numeric(as.vector(levels( as.factor(conj_dados[,"dias"])))) #1:dias )%>% # Caso exista uma venda no dia k eh calculado um # descasque para mercado a vista e a prazo add_constraint( sum_expr( x[i,j,k,m], i = as.numeric(as.vector(levels( as.factor(conj_dados[,"carro"])))),#1:carro, j = as.numeric(as.vector(levels( as.factor(conj_dados[,"frig"])))),#1:frig, m = as.numeric(as.vector(levels( as.factor(conj_dados[,"merc"])))) )#1:merc) + sum_expr( t[i,u,k], i = as.numeric(as.vector(levels( as.factor(conj_dados[,"carro"])))), #1:carro, u = as.numeric(as.vector(levels( as.factor(conj_dados[,"praca"])))) )#,1:praca) <= bigM * y[k], k = as.numeric(as.vector(levels( as.factor(conj_dados[,"dias"]))))# 1:dias )%>% # Respeito ao numero maximo de descasques add_constraint( sum_expr( y[k], k = as.numeric(as.vector(levels( as.factor(conj_dados[,"dias"])))) ) #1:dias) <= Max_descasque )%>% # Desigualdade valida para ajudar na convergencia - # todos os animais devem ser vendidos add_constraint( sum_expr( x[i,j,k,m], i = as.numeric(as.vector(levels(as.factor(conj_dados[, "carro"])))), j = as.numeric(as.vector(levels(as.factor(conj_dados[, "frig"])))), k = as.numeric(as.vector(levels(as.factor(conj_dados[, "dias"])))), m = as.numeric(as.vector(levels(as.factor(conj_dados[, "merc"])))) ) == 40 )%>% # A carga vendida para mercado futuro eh de pelo menos 330 arobas add_constraint( 15*130*sum_expr( t[i,u,k], i = as.numeric(as.vector(levels( as.factor(conj_dados[,"carro"])))) ) #1:carro) <= Carga_mercado_futuro*Trava[k,u], u = as.numeric(as.vector(levels( as.factor(conj_dados[,"praca"])))), #1:praca, k = as.numeric(as.vector(levels( as.factor(conj_dados[,"dias"])))) #1:dias )%>% # Permissao para venda no mercado futuro add_constraint( sum_expr( t[i,u,k], i = as.numeric(as.vector(levels( as.factor(conj_dados[,"carro"])))), #1:carro, u = as.numeric(as.vector(levels( as.factor(conj_dados[,"praca"])))), #1:praca, k = as.numeric(as.vector(levels( as.factor(conj_dados[,"dias"])))) )#1:dias) <= bigM * permissao_venda )%>% saveRDS(model, file = "model.rds") %>% solve_model(with_ROI(solver = "glpk")) Error in on_element(push, inplace_update_ast, get_ast_value, element) : argument "ast" is missing, with no default

2 1

Erro na inserção de restrições - função add_constraint - pacote ompr
by Adriele Giaretta Biase 28 Mai '18

28 Mai '18

Bom dia a todos, pessoal estou com um erro na construção de uma restrição no modelo de pesquisa operacional. A soma da combinação de um índices de fatores deve ser igual um número. O erro é referente algum "=" ou ")" errado. Mas não consigo encontrar. i é um vetor que contem 40 índices (1:40); j é um vetor que contem 4 índices (1:4); k é um vetor que contem 150 dias (1:150); m é um numero que contem apenas um índice (1 apenas); Vejam o código: ### Script: add_constraint( sum_expr( x[i,j,k,m], i = as.numeric(as.vector(levels(as.factor(conj_dados[, "bois"])))), j = as.numeric(as.vector(levels(as.factor(conj_dados[, "frig"])))), k = as.numeric(as.vector(levels(as.factor(conj_dados[, "dias"])))), m = as.numeric(as.vector(levels(as.factor(conj_dados[, "merc"])))) ) = 40 )%>% ### Execução no Console > add_constraint( + sum_expr( + x[i,j,k,m], + i = as.numeric(as.vector(levels(as.factor(conj_dados[, "bois"])))), + j = as.numeric(as.vector(levels(as.factor(conj_dados[, "frig"])))), + k = as.numeric(as.vector(levels(as.factor(conj_dados[, "dias"])))), + m = as.numeric(as.vector(levels(as.factor(conj_dados[, "merc"])))) ) + = 40 Error: unexpected '=' in: " m = as.numeric(as.vector(levels(as.factor(conj_dados[, "merc"])))) ) =" > )%>% Error: unexpected ')' in " )" -- Adriele Giaretta Biase. Mestre em Estatística e Experimentação Agropecuária - UFLA. Doutora em Estatística e Experimentação Agronômica - ESALQ/ USP Contato: (19) 98861-0619.

1 0

GLM para variável resposta contagem e regressoras fatores.
by Alice Peres 28 Mai '18

28 Mai '18

Caros membros, bom dia! Tenho uma base de dados com 4 variáveis. Onde, a resposta é discreta, com dist. Poisson, e as demais regressoras, são fatores, onde cada um tem um número diferente de níveis. Resposta, Peso(g); Regressoras, Tratamento: C, O, S, E1, E2; Tempo: 0, 1, 2; Sexo: M, F. Para saber como as variáveis se comportam com o Peso, inicialmente, apliquei o seguinte modelo: mod<-glm(Peso~Tratamento+Tempo+Sexo+Tratamento:Sexo,family = poisson) Quando eu gero o summary(modelo) ele não tem o tratamento C e nem o sexo M. O modelo que eu fiz está errado ? Abs

2 1

Modelo linear de segundo grau
by Andre Oliveira 27 Mai '18

27 Mai '18

Sei que posso rodar dois modelos e fazer o LRT, entre outras análises, mas por algum motivo preciso da rsm! Alguém saberia o motivo deste modelo não "rodar"? x=c(1,1,1,2,2,2,3,3,3,4,4,4,5,5,5) y=c(4,8,5,9,12,9,10,13,12,12,14,13,11,11,12) plot(y~x) require(rsm) m3=rsm(y~ FO(x)+FO(x^2)) Error in dimnames(LM$B) <- list(fonm, fonm) : comprimento de 'dimnames' [1] não é igual ao tamanho do array Obg

2 2

Extraindo os raios e área de uma elipse ellipsoidhull() do pacote cluster
by ASANTOS 24 Mai '18

24 Mai '18

Caros Membros, Com fins ilustrativos e exploratório, tenho pontos correspondentes a dois tratamentos A e B espacialmente distribuídos no campo, estou delimitando o limite (zona de influência) de cada um como uma elipse, através da função ellipsoidhull() do pacote cluster conforme CRM abaixo. Gostaria de saber se alguém familiarizado com a função poderia me dizer como faço para extrair os raios e as áreas das elipses EpointsA e EpointsB? #Criação do data frame dados<-NULL #Coordenadas x dados$x<-c(0.350,0.487,0.637,0.775,0.825,0.087,0.237,0.400,0.575,0.737,0.062,0.212) #Coordenadas y dados$y<-c(0.025,0.087,0.050,0.025,0.125,0.187,0.150,0.162,0.212,0.237,0.362,0.337) #Classificação dados$Class<-c("A","A","A","A","A","A","B","B","B","B","B","B") #Visualização plot(dados$x,dados$y) text(dados$x,dados$y+0.01,labels=dados$Class,col="red") # #Usando ellipsoidhull require(cluster) PointsA<-unname(cbind(dados$x[1:6],dados$y[1:6])) EpointsA<-ellipsoidhull(PointsA) lines(predict(EpointsA), col="blue") PointsB<-unname(cbind(dados$x[7:12],dados$y[7:12])) EpointsB<-ellipsoidhull(PointsB) lines(predict(EpointsB), col="orange") # Abraços, -- ====================================================================== Alexandre dos Santos Proteção Florestal IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso Campus Cáceres Caixa Postal 244 Avenida dos Ramires, s/n Bairro: Distrito Industrial Cáceres - MT CEP: 78.200-000 Fone: (+55) 65 99686-6970 (VIVO) (+55) 65 3221-2674 (FIXO) e-mails:alexandresantosbr@yahoo.com.br alexandre.santos(a)cas.ifmt.edu.br Lattes: http://lattes.cnpq.br/1360403201088680 OrcID: orcid.org/0000-0001-8232-6722 - ResearcherID: A-5790-2016 Researchgate: www.researchgate.net/profile/Alexandre_Santos10 LinkedIn: br.linkedin.com/in/alexandre-dos-santos-87961635 Mendeley:www.mendeley.com/profiles/alexandre-dos-santos6/ ======================================================================

3 3

Transformação normal score para krigagem
by Wagner Wolff 22 Mai '18

22 Mai '18

Olá pessoal! Estou trabalhando com dados de chuva numa escala de 15min. Meu banco de dados tem muitos zeros o que deixa a distribuição dos dados assimétrica. Gostaria de interpolar espacialmente esses dados utilizando o geoR, ou seja, modelo geostatistico gaussiano. Assim, a normalização dos dados é necessaria, já tentei a boxcox e não funcionou. Pesquisando sobre o assunto encontrei uma solução na transformação normal score. O problema é como faço a retransformação pra escala original após a krigagem? Alguém tem alguma dica ou outra sugestão para esse problema? Outra característica é que os dados são zero inflacionados. Att.

3 3

TRABALHAR COM O CEP PARA CRIAR LATITUTDE E LONGITUDE EM UM BANCO
by Edson Lira 21 Mai '18

21 Mai '18

Boa tarde Listeiros, estou estudando a elaboração de mapas e me deparei com uma situação que ainda não consegui resolver. Se eu tenho um banco de dados com CEPs, como faço para a partir do CEP criar as variáveis lat lomg. No google mapas tem esse serviço, quando eu tenho um banco pequeno é fácil de resolver, mas quando o banco é grande, só para assinantes. Suponha um banco com 16000 linhas, ai é dificultoso. Eu pesquisei e encontrei um pacote chamado cepR, mas não consegui com o material que está disponível resolver o problema. Alguém poderia sugerir algo? [ ]'s. Prof. Edson Lira, Me Estatístico Manaus-Amazonas

3 3

Gerar pdf de um arquivo .rnw - Sweave
by Andre Oliveira 19 Mai '18

19 Mai '18

bom dia! ao processar meu arquivo .Rnw ocorre este erro. Pesquisei e me parece ser um erro comum, mas não achei a solução. Error in driver$finish(drobj) : the output file 'arquivo.tex' has disappeared Calls: <Anonymous> -> <Anonymous> Alguém tem alguma dica? Obrigado

1 0

Boxplot em grupo
by Marcele de Jesus Correa 08 Mai '18

08 Mai '18

Olá, boa tarde! Sou iniciante no R e preciso de muita ajuda. Tenho recorrido aos tutoriais, livros, algumas dúvidas que são postadas aqui e vídeos no youtube, que tem me ajudado bastante, porém, estou com dificuldade para fazer um gráfico boxplot envolvendo três variáveis juntas. Assim: Estou trabalhando com série temporal de três estações meteorológicas, num total de 27 variáveis, entre estas, estou analisando no momento 3 (precipitação, temperatura mínima e máxima) de cada estação. Para cada uma dessas variáveis e suas respectivas estações, elaborei um gráfico boxplot, conforme comando abaixo: boxplot(TempMin_Jagua ~ Mes, las=1, data = data_estacoes, col = "gray", notch = FALSE, ylab="Temperatura", xlab= "Mês", main= "Média Mensal da Temperatura Mínima da Estação de Jaguaruana"). Entretanto, gostaria de inserir num único gráfico a variável temperatura mínima das três estações, a fim de comparar numa mesma escala o "comportamento" dessa variável por estação, e não estou conseguindo. Li algo parecido aqui no grupo, em que foi utilizado o pacote ggplot2, tentei usar sem sucesso. Fiz-me entender? Espero que sim, porque esse negócio é bem complicado... mas estou gostando. Desde já agradeço. Cordialmente, Marcele Correa. Geographer Master Student in Climate Sciences Research Line - Bioclimatology and C limatic C hanges Advisor - Kellen Carla Lima Federal University of Rio Grande do Norte - UFRN / Brazil.

6 8

Instalação de pacote do github
by Andre Oliveira 06 Mai '18

06 Mai '18

Bom dia! faz algum tempo que este problema vem ocorrendo e dificultando as coisas para mim. Veja bem: require(devtools) Carregando pacotes exigidos: devtools install_github("tagteam/Publish") Installation failed: Problem with the SSL CA cert (path? access rights?) Pesquisando passei usar o comando! httr::set_config(httr::config(ssl_verifypeer = FALSE )) install_github("tagteam/Publish") Downloading GitHub repo tagteam/Publish@master from URL https://api.github.com/repos/tagteam/Publish/zipball/master Installation failed: Peer reports incompatible or unsupported protocol version. Não sei mais o que fazer, pois, preciso atualizar alguns pacotes e não consigo. Alguém teria uma sugestão? obg

2 1

Escalonar matrizes
by Andre Oliveira 06 Mai '18

06 Mai '18

boa tarde, alguém conhece uma função que escalona uma matriz? obg André Oliveira Souza

2 1

Emprego da PCA para seleção de variáveis mais correlacionadas no primeiro componente principal
by ASANTOS 05 Mai '18

05 Mai '18

Prezados Membros, Gostaria de selecionar as 5 variáveis resposta no primeiro componente principal mais correlacionadas usando a PCA, através da função prcomp(), conforme o post https://stats.stackexchange.com/questions/115032/how-to-find-which-variable…, sendo o meu CRM: ##Banco de dados set.seed(12345) mat <- matrix(rnorm(120,0,0.5),nrow=6,byrow=TRUE) rownames(mat) <- paste("s",1:6,sep="") colnames(mat) <- paste("g",1:20,sep="") head(mat) ## Espectros com maior correlação pca.object <- prcomp(mat,center=TRUE,scale.=FALSE) plot(pca.object) #Quero os cinco mais correlacionados topN <- 5 load.rot <- scale(pca.object$rotation) names(load.rot[,1][order(abs(load.rot[,1]),decreasing=TRUE)][1:topN]) Então cheguei as variáveis resposta, mas gostaria que alguém me explicasse o por que disso fazendo um favor, pois quando utilizo os meus dados originais, são sempre as variáveis resposta com maior valor numérico que são selecionadas, mesmo utilizando a função scale(). Obrigado, Alexandre -- ====================================================================== Alexandre dos Santos Proteção Florestal IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso Campus Cáceres Caixa Postal 244 Avenida dos Ramires, s/n Bairro: Distrito Industrial Cáceres - MT CEP: 78.200-000 Fone: (+55) 65 99686-6970 (VIVO) (+55) 65 3221-2674 (FIXO) e-mails:alexandresantosbr@yahoo.com.br alexandre.santos(a)cas.ifmt.edu.br Lattes: http://lattes.cnpq.br/1360403201088680 OrcID: orcid.org/0000-0001-8232-6722 - ResearcherID: A-5790-2016 Researchgate: www.researchgate.net/profile/Alexandre_Santos10 LinkedIn: br.linkedin.com/in/alexandre-dos-santos-87961635 Mendeley:www.mendeley.com/profiles/alexandre-dos-santos6/ ======================================================================

3 7

Escala em gráficos
by Lucas Stempkowski 05 Mai '18

05 Mai '18

Bom dia, Plotei um gráfico boxplot com o seguinte comando: boxplot(prodt~sistema, xlab="Anos sem trigo", ylab="Produtividade kg ha-¹") Como a produtividade foi variável na cultivares testadas, gostaria de padronizar a escala no eixo y para que ficasse dessa forma 1500, 2000, 2500, 3000, 3500, 4000, como faço essa alteração? Muito obrigado. -- ----------------------------------------------------- Lucas Antonio Stempkowski Eng. Agronômo - UFFS Mestrando em Produção Vegetal CAV/UDESC +55 54 99153-0124/ 98448-4633 --------------------------------------------------

2 1

Estrutura de covariância proposto por Singer, Nobre e Rocha (2017)
by Fernando Souza 30 Abr '18

30 Abr '18

No livro Singer, Nobre e Rocha (2017) disponível em: https://www.ime.usp.br/~jmsinger/MAE0610/Singer&Nobre&Rocha2017junho.pdf os autores citam como estratégia para identificação da estrutura da matriz de variãncia covariãncia entre outros a construção do seguinte gŕafico: " No caso de estudos longitudinais sem dados omissos, análise do gráfico das linhas da matriz de covariâncias intraunidades amostrais por meio de técnicas de regressão polinomial para a identificação de efeitos aleatórios." Tenho algumas suspeitas, mas não estou seguro. Como é a aparência desse gráfico? -- ========================================= Fernando Souza Zootecnista, DSc. Produção e Alimentação Animal Celular: (31)99796-8781 (Vivo) E-mail:nandodesouza@gmail.com <e-mail%3Anandodesouza(a)gmail.com> Lattes: http://lattes.cnpq.br/6519538815038307 Blog: https://producaoanimalcomr.wordpress.com/ ==========================================

1 0

Tamanho máximo de base de dados para trabalhar no R
by Charlenny Melo 26 Abr '18

26 Abr '18

Pessoal, boa tarde. Alguém sabe me informar qual a quantidade de linhas e variáveis que o R consegue ler? Estou com uma base de dados de mais de 2500 linhas, e o R está lendo apenas 1000 linhas. Alguém pode me ajudar? Obrigada.

3 2

Extrair por sorteio linhas de um data frame dentro de cada nível de um fator
by ASANTOS 26 Abr '18

26 Abr '18

Prezados Membros do R, Preciso extrair linhas de um data frame dentro de cada nível de um fator, usando a função sample_n() do pacote dplyr. Sendo meu CRM: #Data set artificial df <- data.frame(matrix(rnorm(80), nrow=40)) df$color <- rep(c("blue", "red", "yellow", "pink"), each=10) require(dplyr) #Sorteio sample_n(df, 3) #Sorteio sobre todo o data frame X1 X2 color 6 -0.2386293 -1.2611679 blue 21 0.6337288 -0.7817087 yellow 28 -0.1949911 0.6397060 yellow No entanto, não é bem isso que eu quero, na verdade, eu gostaria de realizar o sorteio de uma amostra aleatória de tamanho 3 em cada nível do fator color, ou seja, quero que seja sorteados 3 valores para cada cor. Ficando algo do tipo: X1 X2 color 3 0.22232534 1.83730750 blue 6 -0.23862925 -1.26116791 blue 7 -1.22555359 -1.60739354 blue 15 -1.12092508 -0.04395679 red 16 -0.68780179 0.51116139 red 17 -0.98583965 0.09294828 red ... Alguém poderia me ajudar a resolver isso? Obrigado, Alexandre -- ====================================================================== Alexandre dos Santos Proteção Florestal IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso Campus Cáceres Caixa Postal 244 Avenida dos Ramires, s/n Bairro: Distrito Industrial Cáceres - MT CEP: 78.200-000 Fone: (+55) 65 99686-6970 (VIVO) (+55) 65 3221-2674 (FIXO) e-mails:alexandresantosbr@yahoo.com.br alexandre.santos(a)cas.ifmt.edu.br Lattes: http://lattes.cnpq.br/1360403201088680 OrcID: orcid.org/0000-0001-8232-6722 - ResearcherID: A-5790-2016 Researchgate: www.researchgate.net/profile/Alexandre_Santos10 LinkedIn: br.linkedin.com/in/alexandre-dos-santos-87961635 Mendeley:www.mendeley.com/profiles/alexandre-dos-santos6/ ======================================================================

2 1

MAPAS
by Edson Lira 18 Abr '18

18 Abr '18

Bom dia caros amigos, nesses links tem o script e o banco que estou usando para confeccionar mapas de sarampo no município de Manaus. Eu consegui plotar o mapa por bairro. Plotei os mapas de cada região (6) separadas. Meu problema: quero montar um mapa único das 6 regiões , no script vcs verão que são 6 zonas geográficas separadas. Só não consegui montar u mapa único das 6 regiões, sem plotar os nomes dos bairros aonde o meu poligono seria as 6 regiões. SCRIPT Download mapa mao sar.R | | | | Download mapa mao sar.R | | | BANCO Download manaus sar 63 2016.xlsx | | | | Download manaus sar 63 2016.xlsx | | | [ ]'s. Prof. Edson Lira, Me Estatístico Manaus-Amazonas

1 1

Re: [R-br] [SPAM Detectado]como usar um layout gráfico do tipo 3x2 ?
by Felinto COSTA 16 Abr '18

16 Abr '18

Cleber: Tente com o multiplot: http://www.cookbook-r.com/Graphs/Multiple_graphs_on_one_page_(ggplot2)/ Caso o ggplot2 não traga a função, use este script: http://www.peterhaschke.com/r/2013/04/24/MultiPlot.html Abs. FCosta ================================================================= On Mon, 16 Apr 2018 16:03 , "Cleber N.Borges via R-br" r-br(a)listas.c3sl.ufpr.br sent: > Olá a todos, > > Gostaria de fazer 5 gráficos em uma mesma janela... > > mas não entendi como usar as funções de layout... > Gostaria de algo mais ou menos como coloco abaixo: > > ######## > > # primeira linha na janela gráfica => 3 > gráficos > plot( rnorm(10), col=1 ); plot( rnorm(10), col=2 ); plot( rnorm(10), > col=3 ); > > # segunda linha na janela gráfica => 2 > gráficos > hist( rnorm(10), col=1 ); hist( rnorm(10), col=2 ); plot( rnorm(10), > col=3 ); > > > agradeço antecipadamente pela ajuda > > Cleber Borges > > > --- > Este email foi escaneado pelo Avast antivírus. > https://www.avast.com/antivirus > _______________________________________________ > R-br mailing list > R-br > @listas.c3sl.ufpr.brhttps://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e fornea cdigo > mnimo reproduzvel. > ---- Mensagem enviada pelo sistema e-mail Onda- Powered By Onda Empresas

1 0

Plot_ly background
by Jobenil Luiz Magalhães Júnior 13 Abr '18

13 Abr '18

Prezados(as), boa tarde Gostaria de construir um gráfico usando o pacote plotly de acordo com o exemplo que aparece no link abaixo. Ao que parece devo fazer as alterações alterando "shapes". A questão que ainda não consegui resolver é como devo configurar estes parâmetros no R. Alguém pode me dar alguma dica a respeito. Agradeço desde já, Atenciosamente Jobenil Júnior https://stackoverflow.com/questions/47064886/setting- variable-yaxis-background-color-on-plot-ly-js-based-chart

2 2

arredondamento de coordenada utm
by Ana Maria Nievas 10 Abr '18

10 Abr '18

Olá! Tenho um data frame com coordenadas x e y, em utm, que estão sendo arredondadas e ficando com um número de dois dígitos, ao invés da coordenada completa. Alguma dica de resolução desse problema? -- *Ana Maria Nievas* Mestre em Ecologia Aplicada Doutoranda em Psicobiologia Laboratório de Etologia e Bioacústica - EBAC Faculdade de Filosofia, Ciências e Letras de Ribeirão Preto - FFCLRP/USP Lattes: http://lattes.cnpq.br/3529613851615012

1 0

Reproduçao dos resultados de Gabriel(1971)
by Luiz Roberto Martins Pinto 06 Abr '18

06 Abr '18

Pessoal, Alguém conseguiu reproduzir no R a Tabela 2 e a Figura 2 de Gabriel(1971)? (anexo) Abraços, Luiz Roberto. Luiz Roberto Martins Pinto Prof. Pleno/DCET/UESC Laboratório de Estatística Computacional Universidade Estadual de Santa Cruz Ilhéus-Bahia-Brasil luizroberto.uesc(a)gmail.com skype: lrmpinto http://lattes.cnpq.br/2732314327604831 "*The s**cience exists because there are patterns. * * The patterns exist because God created them*. * The statistic exists to research the patterns that God created.*"

2 5

Chi-quadrado para amostras grandes
by Elias Carvalho 06 Abr '18

06 Abr '18

Boa tarde Estou estudando correlação e em um teste para medir a intensidade de relacionamento entre duas variáveis categóricas (1 dicotômica e outra ordinal com 5 níveis). Observando o resultado abaixo percebe-se o p-value do chi-quadrado muito baixo indicando que não são independentes,. No entanto a intensidade da relação medida indica uma intensidade muito baixo: 0.03. Acredito que isso está ocorrendo pelo tamanho da amostra (63.3884 registros). Há uma solução pare resolver este problema? -- *In Jesu et Maria* *Obrigado* *Prof. Elias Carvalho* *"Felix, qui potuit rerum cognoscere causas" (Virgil 29 BC)"Blessed is he who has been able to understand the cause of things"*

2 1

github - Instalação de pacote
by Andre Oliveira 06 Abr '18

06 Abr '18

Boa tarde, estou com este problema de incompatibilidade: require(devtools) install_github("tagteam/Publish") Downloading GitHub repo tagteam/Publish@master from URL https://api.github.com/repos/tagteam/Publish/zipball/master Installation failed: Peer reports incompatible or unsupported protocol version. Alguma sugestão em apontar a solução? Obrigado

1 0

Margem de gráfico
by Andre Oliveira 05 Abr '18

05 Abr '18

Boa tarde,já mexi em parâmetros de margens, de eixos, mas sem sucesso. Preciso diminuir a distância entre as barras do IC e dos respectivos nomes de cada grupo. Segue um CMR. Alguma dica? require(Publish) data(Diabetes) x=ci.mean(bp.2s/500~frame,data=Diabetes) plot(x, leftmargin=8, rightmargin=0) Grato,André

2 1

R Day - última semana para submissão de propostas
by Fernando Mayer 03 Abr '18

03 Abr '18

O prazo para submissão de propostas de apresentação oral e tutorial para o R Day está acabando! O prazo se encerra dia 08/04 (domingo), portanto não deixe de enviar sua proposta! Já temos dois conferencistas principais confirmados: - Gabriela de Queiroz (Fundadora R Ladies) - via webconferência - Paulo Justiniano Ribeiro Jr. (LEG/DEST/UFPR) Também atualizamos o local do evento: - UFPR - Campus Jardim Botânico - Setor de Ciências da Saúde E não esqueça da data: 22 de maio de 2018. Mais informações em: http://rday.leg.ufpr.br Fernando Mayer. P/ Comissão Organizadora. --- Fernando de Pol Mayer Laboratório de Estatística e Geoinformação - LEG Departamento de Estatística - DEST Universidade Federal do Paraná - UFPR URL: http://leg.ufpr.br/~fernandomayer e-mail: fernando.mayer [@] {gmail.com, ufpr.br}

2 1

Formatar base de dados
by Andre Oliveira 01 Abr '18

01 Abr '18

bom dia, tenho vetor muito grande e preciso fazer como fiz manualmente neste CMR. notas =c(2, 3, 20, 4, 8, 9, 2, 20, 4, 4)notas[notas == 2 ] <- "A"notas[notas == 3 ] <- "B"notas[notas == 20 ] <- "C"notas[notas == 4 ] <- "D"notas[notas == 8 ] <- "E"notas[notas == 9 ] <- "F"notas Resolveu, mas é cansativo! Se alguém tiver algo melhor e puder compartilhar! gratoAndré

2 4

dificuldade em carregar pacote já instalado
by Ana Maria Nievas 31 Mar '18

31 Mar '18

Bom dia! Nunca havia enfrentado esse problema antes, mas nos últimos dias não tenho conseguido carregar pacotes já instalados. Na tentativa de nova instalação e carregamento, obtive o seguinte: > utils:::menuInstallPkgs() tentando a URL ' https://vps.fmvz.usp.br/CRAN/bin/windows/contrib/3.3/ggplot2_2.2.1.zip' Content type 'application/zip' length 2760677 bytes (2.6 MB) downloaded 2.6 MB package ‘ggplot2’ successfully unpacked and MD5 sums checked The downloaded binary packages are in C:\Users\Ana Maria\AppData\Local\Temp\RtmpYdZpFs\downloaded_packages > library(ggplot2) Error in loadNamespace(j <- i[[1L]], c(lib.loc, .libPaths()), versionCheck = vI[[j]]) : there is no package called ‘Rcpp’ Erro: package or namespace load failed for ‘ggplot2’ Acho que esse problema tem a ver com a pasta onde estão sendo salvos os arquivos de instalação, mas não sei como consertar. Alguém poderia me ajudar? Obrigada. -- *Ana Maria Nievas* Mestre em Ecologia Aplicada Doutoranda em Psicobiologia Laboratório de Etologia e Bioacústica - EBAC Faculdade de Filosofia, Ciências e Letras de Ribeirão Preto - FFCLRP/USP Lattes: http://lattes.cnpq.br/3529613851615012

3 2

Re: [R-br] Digest R-br, volume 87, assunto 22
by Alice Peres 29 Mar '18

29 Mar '18

Consegui resolver sim, obrigada. Em qui, 29 de mar de 2018 11:00, <r-br-request(a)listas.c3sl.ufpr.br> escreveu: > Enviar submissões para a lista de discussão R-br para > r-br(a)listas.c3sl.ufpr.br > > Para se cadastrar ou descadastrar via WWW, visite o endereço > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > ou, via email, envie uma mensagem com a palavra 'help' no assunto ou > corpo da mensagem para > r-br-request(a)listas.c3sl.ufpr.br > > Você poderá entrar em contato com a pessoa que gerencia a lista pelo > endereço > r-br-owner(a)listas.c3sl.ufpr.br > > Quando responder, por favor edite sua linha Assunto assim ela será > mais específica que "Re: Contents of R-br digest..." > > > Tópicos de Hoje: > > 1. IC para a diferença de médias de dois grupos independentes > (Alice Peres) > 2. Re: IC para a diferença de médias de dois grupos > independentes (Rodrigo Araújo e Castro) > 3. Re: IC para a diferença de médias de dois grupos > independentes (Rodrigo Araújo e Castro) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Wed, 28 Mar 2018 21:13:18 -0400 > From: Alice Peres <neves.peres135(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: [R-br] IC para a diferença de médias de dois grupos > independentes > Message-ID: > <CAJ+iSzAVPzA2WPfPTgtF= > ByBJL5RzaezOCiA8TTB7wSV2kJ_MA(a)mail.gmail.com> > Content-Type: text/plain; charset="utf-8" > > Olá, boa noite. Eu estou com dificuldade para gerar o intervalo de > confiança para a diferença de média dos dois seguintes grupos. Lembrando > que são independentes. Eu calculei o teste t no R, mas está dando um erro. > Preciso instalar algum pacote ? > > > grupo1=c(4.81,4.17,4.41,3.59,5.87,3.83,6.03,4.98,4.90,5.75,5.36,3.48,4.69,4.44,4.89,4.71,5.48,4.32,5.15,6.34) > > > grupo2=c(4.17,3.05,5.18,4.01,6.11,4.10,5.17,3.57,5.33,5.59,4.66,5.58,3.66,4.50,3.90,4.61,5.62,4.53,6.05,5.14) > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20180328/4ae22e59/atta… > > > > ------------------------------ > > Message: 2 > Date: Thu, 29 Mar 2018 05:34:35 +0000 (UTC) > From: Rodrigo Araújo e Castro <roacastro87(a)yahoo.com.br> > To: Alice Peres via R-br <r-br(a)listas.c3sl.ufpr.br> > Subject: Re: [R-br] IC para a diferença de médias de dois grupos > independentes > Message-ID: <1943970344.734841.1522301675527(a)mail.yahoo.com> > Content-Type: text/plain; charset="utf-8" > > Boa noite, > com base nos conjuntos de dados executei o comando t.test (grupo1,grupo2) > para executar o teste T (não necessita de pacote) e não deu erro. O > resultado foi o seguinte: > Welch Two Sample t-test > > data: grupo1 and grupo2 > t = 0.50985, df = 37.711, p-value = 0.6131 > alternative hypothesis: true difference in means is not equal to 0 > 95 percent confidence interval: > -0.3967069 0.6637069 > sample estimates: > mean of x mean of y > 4.8600 4.7265 > > > Se ainda quiser calcular o intervalo de confiança de cada grupo, pode > utilizar, pode exemplo, o pacote Publish (após instalar) da sequinte forma: > install.packages("Publish", dependencies=T) > require(Publish) > > ci.mean(grupo1) > ci.mean(grupo2) > Conseguiu resolver seu problema? Se não, favor detalhar o que faltou para > que possamos ajudá-lo. > att,Rodrigo > Em quarta-feira, 28 de março de 2018 22:13:33 BRT, Alice Peres via > R-br <r-br(a)listas.c3sl.ufpr.br> escreveu: > > Olá, boa noite. Eu estou com dificuldade para gerar o intervalo de > confiança para a diferença de média dos dois seguintes grupos. Lembrando > que são independentes. Eu calculei o teste t no R, mas está dando um erro. > Preciso instalar algum pacote ? > > grupo1=c(4.81,4.17,4.41,3.59,5.87,3.83,6.03,4.98,4.90,5.75,5.36,3.48,4.69,4.44,4.89,4.71,5.48,4.32,5.15,6.34) > > grupo2=c(4.17,3.05,5.18,4.01,6.11,4.10,5.17,3.57,5.33,5.59,4.66,5.58,3.66,4.50,3.90,4.61,5.62,4.53,6.05,5.14)_______________________________________________ > 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. > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20180329/44dadb48/atta… > > > > ------------------------------ > > Message: 3 > Date: Thu, 29 Mar 2018 05:38:22 +0000 (UTC) > From: Rodrigo Araújo e Castro <roacastro87(a)yahoo.com.br> > To: Alice Peres via R-br <r-br(a)listas.c3sl.ufpr.br> > Subject: Re: [R-br] IC para a diferença de médias de dois grupos > independentes > Message-ID: <39593106.754027.1522301902990(a)mail.yahoo.com> > Content-Type: text/plain; charset="utf-8" > > Alice, > se quiser algum material de referência, sugiro estes a seguir, dentre > outros: > > http://staff.pubhealth.ku.dk/~tag/Teaching/share/R-tutorials/ConfidenceInte… > R Handbook: Confidence Intervals > Interval Estimation of the Population Mean > > | > | > | | > Interval Estimation of the Population Mean > > Interval estimation of the population mean can be computed from the > functions of the following R packages:stats ... > | > > | > > | > > > > > > | > | > | | > R Handbook: Confidence Intervals > > Clear examples in R. Confidence intervals for means; Confidence intervals > by bootstrap; Exercises D > | > > | > > | > > > > > att,Rodrigo > Em quinta-feira, 29 de março de 2018 02:34:35 BRT, Rodrigo Araújo e > Castro <roacastro87(a)yahoo.com.br> escreveu: > > Boa noite, > com base nos conjuntos de dados executei o comando t.test (grupo1,grupo2) > para executar o teste T (não necessita de pacote) e não deu erro. O > resultado foi o seguinte: > Welch Two Sample t-test > > data: grupo1 and grupo2 > t = 0.50985, df = 37.711, p-value = 0.6131 > alternative hypothesis: true difference in means is not equal to 0 > 95 percent confidence interval: > -0.3967069 0.6637069 > sample estimates: > mean of x mean of y > 4.8600 4.7265 > > > Se ainda quiser calcular o intervalo de confiança de cada grupo, pode > utilizar, pode exemplo, o pacote Publish (após instalar) da sequinte forma: > install.packages("Publish", dependencies=T) > require(Publish) > > ci.mean(grupo1) > ci.mean(grupo2) > Conseguiu resolver seu problema? Se não, favor detalhar o que faltou para > que possamos ajudá-lo. > att,Rodrigo > Em quarta-feira, 28 de março de 2018 22:13:33 BRT, Alice Peres via > R-br <r-br(a)listas.c3sl.ufpr.br> escreveu: > > Olá, boa noite. Eu estou com dificuldade para gerar o intervalo de > confiança para a diferença de média dos dois seguintes grupos. Lembrando > que são independentes. Eu calculei o teste t no R, mas está dando um erro. > Preciso instalar algum pacote ? > > grupo1=c(4.81,4.17,4.41,3.59,5.87,3.83,6.03,4.98,4.90,5.75,5.36,3.48,4.69,4.44,4.89,4.71,5.48,4.32,5.15,6.34) > > grupo2=c(4.17,3.05,5.18,4.01,6.11,4.10,5.17,3.57,5.33,5.59,4.66,5.58,3.66,4.50,3.90,4.61,5.62,4.53,6.05,5.14)_______________________________________________ > 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. > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20180329/b879499f/atta… > > > > ------------------------------ > > Subject: Legenda do Digest > > _______________________________________________ > R-br mailing list > R-br(a)listas.c3sl.ufpr.br > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > > > ------------------------------ > > Fim da Digest R-br, volume 87, assunto 22 > ***************************************** >

1 0

IC para a diferença de médias de dois grupos independentes
by Alice Peres 29 Mar '18

29 Mar '18

Olá, boa noite. Eu estou com dificuldade para gerar o intervalo de confiança para a diferença de média dos dois seguintes grupos. Lembrando que são independentes. Eu calculei o teste t no R, mas está dando um erro. Preciso instalar algum pacote ? grupo1=c(4.81,4.17,4.41,3.59,5.87,3.83,6.03,4.98,4.90,5.75,5.36,3.48,4.69,4.44,4.89,4.71,5.48,4.32,5.15,6.34) grupo2=c(4.17,3.05,5.18,4.01,6.11,4.10,5.17,3.57,5.33,5.59,4.66,5.58,3.66,4.50,3.90,4.61,5.62,4.53,6.05,5.14)

2 2

Erro ao fazer boxplot - Arquivo SPSS
by Edmar Caldas 27 Mar '18

27 Mar '18

Bom Dia! estou trabalhando com um arquivo SPSS. Ao fazer um boxplot aparece essa mensagem. Error: 'x' and labels' must be same type. Edmar

5 8

Dúvida TRUE/FALSE
by Danielle Tamaki 23 Mar '18

23 Mar '18

Boa tarde, estou com uma dúvida na linguagem R e precisava de uma ajuda! Preciso fazer a comparação dos valores de um vetor se eles são maiores que um determinado valor. Porém, ele não faz a comparação e me retorna um erro: Error in if ((aux <= limiteinferior) == TRUE) { : valor ausente onde TRUE/FALSE necessário Como posso solucionar este problema? Obrigada pela ajuda! -- -- *Danielle Mayumi Campos Tamaki* Doutoranda em Engenharia de Produção (UNIFEI) Mestre em Engenharia de Produção (UNIFEI) Bacharel em Ciência da Computação (UNIFEI) *Celular:* +55 35 9222 2746 *Skype:* danielle.tamaki *Outro Email:* danimayu88(a)gmail.com

4 3

Erro no truncamento da query ODBC
by Pedro Emmanuel Alvarenga Americano do Brasil 23 Mar '18

23 Mar '18

Amigos de R, Começo pedindo desculpas porque não sei como postar um script que reproduza o erro que estou encarando. Estou tentando fazer uma consulta numa base SQL, e preciso consultar muitas e muitas colunas de uma vez. Muitos nomes de colunas são longos. o fim da linha de erro que encontro é... failed with -1 error code. Error on Job: ODBC statement error -1 :[Microsoft][ODBC SQL Server Driver]Truncamento Error in doTryCatch(return(expr), name, parentenv, handler) : ODBC statement error -1 :[Microsoft][ODBC SQL Server Driver][SQL Server]Unable to launch runtime for 'R' script. Please check the configuration of the 'R' runtime. Depois de muito fuçar, vi em alguns lugares que o driver do SQL usado na conexão possui uma limitação de 4096 caracteres na query. Contei aqui as minhas querys e todas possuem mais de 4700 caracteres. Antes de sair cortando as colunas da query e fatiar a tarefa, alguem saberia se há alguma possível configuração do driver do SQL que aumente essa limitação de caracteres? Abraço forte, Pedro Brasil

1 0

Oportunidade para estudantes: desenvolvimento de projeto durante o GSOC com o InterMineR
by Marcos Vital 19 Mar '18

19 Mar '18

Oi, pessoal! Recebi no Twitter, e repasso aqui. Não estou por dentro dos detalhes, então qualquer dúvida deveria ser tirada diretamente com os responsáveis, ok? A equipe do InterMineR (https://github.com/intermine/InterMineR) está procurando estudantes que queriam contribuir com o projeto durante o Google Summer of Code. Um dos objetivos é construir um aplicativo Shiny para o pacote. A chamada está aqui: http://intermine.org/gsoc/project-ideas/2018/#interminer-a-client-library-f… E nela estão os emails dos responsáveis, para o caso de alguém se interessar e quiser entrar em contato. Abraços! -- Marcos Vinícius Carneiro Vital Universidade Federal de Alagoas Instituto de Ciências Biológicas e da Saúde Setor de Biodiversidade e Ecologia

1 0

Modelo geoestatístico sobre condição, jogo SUDOKU
by Wagner Wolff 19 Mar '18

19 Mar '18

Olá pessoal da lista! Estou lidando com um problema de certa forma parecido com o jogo SUDOKU. Resumindo nesse jogo tem uma área subdividida em quadriculas ou células que devem ser preenchidas seguindo condições de que os números de 1 a 9 não se repitam nas colunas e linhas. No meu caso, eu preciso interpolar espacialmente uma variável chamada de intensidade da chuva, entretanto, eu tenho a média dessa variável entre dois pontos que ligando eles formam uma linha, figura abaixo. A questão é a seguinte, como a partir da média eu obtenho os dados pontuais sobre as linhas que satisfaça a condição de que a média permaneça até mesmo nos pontos onde as linhas se cruzam. Ou seja, onde elas se cruzam precisa ser um valor que satisfaça a média da linha 1 e da linha 2. Depois disso preciso interpolar espacialmente esses pontos na área toda. Eu obtive uma solução, mas não sei se é a mais adequada e talvez para uma banco de dados maior ela pode ser impraticável. Resumindo, eu simulei valores nas linhas com a função grf() do geoR e depois com a função optim() variei as células (pixels) até que as condições sejam satisfeitas. Segue o código reproduzível e em anexo os arquivos referente ao shapefile que representa as linhas. proj="+proj=sterea +lat_0=52.15616055555555 +lon_0=5.38763888888889 +k=0.9999079 +x_0=155000 +y_0=463000 +ellps=bessel +units=km +no_defs " ## Colocar os arquivis anexados na mesma pasta da código R Links2l <- spTransform(readOGR("Links2l.shp"), CRS(proj)) Rmean <- c(2.9424, 4.3982) ## Criando uma imagem raster vazia ras <- raster(Links2l, res=0.1, ext=extent(Links2l)); ras[] <- rep(NA, ncell(ras)) ## Celulas (pixels) interceptadas pelas links cells_links <- cellFromLine(ras, Links2l) cells <- unique(unlist(cells_links)) ## Celulas compartilhadas entre os links cells_links[[1]][which(cells_links[[1]] %in% cells_links[[2]])] ## Simulando valores a partir da media nas celulas interceptadas pelas links set.seed(234) for(i in 1:nrow(Links_ams)){ ras[cells_links[[i]]] <- apply(grf(length(cells_links[[i]]), grid=xyFromCell(ras, cells_links[[i]]), nsim=10, mean = Rmean[i],cov.pars = c(3.7, 18.7), nugget = 0.37)$data, 1, mean) } plot(ras, col = gray.colors(10, start = 0.3, end = 0.9, gamma = 2.2, alpha = NULL)) lines(Links2l) ## Funcao de otimizacao opt <- function(par, Rmean){ res <- c() for(i in 1:length(Rmean)){ res[i] <- sum(mean(par[cells %in% cells_links[[i]]])-Rmean[i])^2 } sum(res) } ## Definir espaco de iteracao lower <- list() upper <- list() for(i in 1:nrow(Links_ams)){ lower[[i]] <- rep(mean(ras[cells_links[[i]]])-3.7^0.5, ncell(cells_links[[i]])) upper[[i]] <- rep(mean(ras[cells_links[[i]]])+3.7^0.5, ncell(cells_links[[i]])) } lower <- ifelse(unlist(lower)<0, 0, unlist(lower)) upper <- ifelse(unlist(upper)<0, mean(unlist(upper))+sd(unlist(upper)), unlist(upper)) ## Otimizar e atualizar as celulas para resolver a condicao cell_upd <- optim(par=ras[cells], fn=opt,lower = lower, upper=upper, Rmean=Rmean,method="L-BFGS-B")$par ras_old <- ras ras[cells] <- cell_upd ## Comparando celulas antigas com as novas, apenas a célula compartilhada entre os links nao é linear plot(ras_old[cells_links[[1]]], ras[cells_links[[1]]], ann=T) abline(a=0, b=1) plot(ras_old[cells_links[[2]]], ras[cells_links[[2]]], ann=T) abline(a=0, b=1) ## Verificando se as condicoes foram satisfeitas mean(ras[cells_links[[1]]]) Rmean[1] mean(ras[cells_links[[2]]]) Rmean[2] ## Resultado, valores iguais "SUDOKU" resolvido! ###################################################### ## Interpolacao das celulas pts <- rasterToPoints(ras) plot(pts) boxcox(pts~1) geo_link <- as.geodata(pts, coords.col = 1:2, data.col = 3) plot(geo_link, low =T) plot(geo_link, trend='1st', low =T) plot(geo_link, trend='2nd', low =T) vario <- variog(geo_link,max.dist=10,uvec=seq(0, 10, by=0.01)) plot(vario) link_fit <- list() link_fit$cte <- likfit(geo_link, cov.model="exp", trend ="cte",ini=c(0.5,5), nugget = 0.1) link_fit$lon <- likfit(geo_link, cov.model="exp", trend =geo_link$coords[, 1],ini=c(0.5,5), nug=0.1) link_fit$st <- likfit(geo_link, cov.model="exp", trend ="1st",ini=c(0.5,5), nug=0.1) link_fit$nd <- likfit(geo_link, cov.model="exp", trend ="2nd",ini=c(0.5,20), nug=0.1) sapply(link_fit, AIC) summary(link_fit$st) ## Predição na área krigagem grid <- xyFromCell(ras, 1:ncell(ras)) link_fit$st kr.link <- krige.conv(geo_link, loc=grid, krige=krige.control(obj=link_fit$st)) ras[] <- kr.link$predict image(ras) contour(ras, add=T) lines(Links2l, col=3, lwd=2) ## Krigagem ordinaria é BLUE, consequentemente condicoes sao satisfeitas mean(ras[cells_links[[1]]]) Rmean[1] mean(ras[cells_links[[2]]]) Rmean[2] Att. Abraco -- *Wagner Wolff, **PhD* "*Luiz de Queiroz**" College of Agriculture,* University of São Paulo Pádua Dias avenue11 | 13418-900| Piracicaba-SP| Brazil Phone: +55 19 982385582 <+55%2019%2098238-5582> http://orcid.org/0000-0003-3426-308X https://github.com/wwolff7 http://buscatextual.cnpq.br/buscatextual/visualizacv.do?id=K4463141A1

1 0

comparação entre: ginv(x) com solve( t(x)%*% x ) %*% t(x)
by Cleber N.Borges 17 Mar '18

17 Mar '18

olá a todos Ainda que a dúvida seja mais conceitual do que do R propriamente.... peço ajuda aqui Minha pergunta é: Até que ponto eu posso usar essa igualdade?? Ou mesmo se posso considerar isso igual... ############################################# library( MASS ) # para usar a função ginv x <- matrix( runif(100*100), 100, 100 ) all.equal( ginv(x), solve( t(x)%*% x ) %*% t(x) ) # TRUE ############################################# O motivo de minha dúvida e vontade de usar isso nos modelos lineares com variáveis linearmente dependentes do tipo de misturas químicas, onde a soma das proporções logicamente sempre tem que ser 100%. Por hora, exponho somente essa dúvida (entre tantas :-) ) Agradeço antecipadamente por qualquer indicação nesse tema de 'ginv' Muito Obrigado Cleber Borges --- Este email foi escaneado pelo Avast antivírus. https://www.avast.com/antivirus

2 1