- R-br - DINF UFPR

OFF TOPIC: 1o Congresso de Aplicações de Linguagem R em Administração COPPEAD/UFRJ.
by Marcos Silva 14 Mar '16

14 Mar '16

Talvez possa ser de interesse para alguns: http://cara.coppead.ufrj.br/ Abs. -- Marcos F. Silva http://sites.google.com/site/marcosfs2006

1 0

OFF-TOPIC: Aniversário da lista
by Éder Comunello 13 Mar '16

13 Mar '16

Pessoal, boa tarde! Em 15 de março a r-br completará 5 anos. Meus parabéns aos mantenedores da lista e todos seus usuários por essa importante fonte de aprendizado. Abaixo vai um código que encontrei na rede e que permite tocar o "Parabéns pra você" no R. ### <code r> # browseURL(" http://stackoverflow.com/questions/31782580/how-can-i-play-birthday-music-u… ") library("dplyr") library("audio") notes <- c(A = 0, B = 2, C = 3, D = 5, E = 7, F = 8, G = 10) pitch <- "D D E D G F# D D E D A G D D D5 B G F# E C5 C5 B G A G" duration <- c(rep(c(0.75, 0.25, 1, 1, 1, 2), 2), 0.75, 0.25, 1, 1, 1, 1, 1, 0.75, 0.25, 1, 1, 1, 2) bday <- data_frame(pitch = strsplit(pitch, " ")[[1]], duration = duration) bday <- bday %>% mutate(octave = substring(pitch, nchar(pitch)) %>% {suppressWarnings(as.numeric(.))} %>% ifelse(is.na(.), 4, .), note = notes[substr(pitch, 1, 1)], note = note + grepl("#", pitch) - grepl("b", pitch) + octave * 12 + 12 * (note < 3), freq = 2 ^ ((note - 60) / 12) * 440) tempo <- 120 sample_rate <- 44100 make_sine <- function(freq, duration) { wave <- sin(seq(0, duration / tempo * 60, 1 / sample_rate) * freq * 2 * pi) fade <- seq(0, 1, 50 / sample_rate) wave * c(fade, rep(1, length(wave) - 2 * length(fade)), rev(fade)) } bday_wave <- mapply(make_sine, bday$freq, bday$duration) %>% do.call("c", .) play(bday_wave) ### </code> ================================================ Éder Comunello Agronomist (UEM), MSc in Environ. Sciences (UEM) DSc in Agricultural Systems Engineering (USP/Esalq) Brazilian Agricultural Research Corporation (Embrapa) Dourados, MS, Brazil |<O>| ================================================ GEO, -22.2752, -54.8182, 408m UTC-04:00 / DST: UTC-03:00

1 0

GlobalEnv
by Wagner Tassinari 11 Mar '16

11 Mar '16

Olá pessoal, Sei que é besteira, mas o que significa esta mensagem abaixo, pois não consigo atachar o meu banco quando ela aparece. Como devo proceder ? > attach(data) The following object is masked _by_ .GlobalEnv: cadastro Obrigado pessoal... abs -- Wagner S. Tassinari Departamento de Matemática Universidade Federal Rural do Rio de Janeiro. BR-465, Km 7 - Seropedica, RJ - Brasil CEP: 23890-000 Skype: wagner.tassinari wtassinari(a)gmail.com tassinari(a)ufrrj.br ------------------------------------------------------- "Statistical thinking will one day be as necessary for efficient citizenship as the ability to read and write." (H.G.Wellis)

2 1

Criar comparações múltiplas entre fatores
by ASANTOS 11 Mar '16

11 Mar '16

Caros listeiros, Gostaria de achar uma maneira de automatizar a criação de comparações múltiplas entre fatores, no meu caso: #Fator treat <- gl(3, 20, labels = paste("t", 1:3, sep="")) #Variables sp1_1 = rnorm(10, 5, 0.25) sp1_2 = rnorm(50, 2.5, 0.25) sp1 = c(sp1_1,sp1_2) sp2 = rnorm(60, 2.5, 0.25) sp3_1 = rnorm(10, 12, 0.25) sp3_2 = rnorm(50, 2.5, 0.25) sp3 = c(sp3_1,sp3_2) sp4 = rnorm(60, 2.5, 0.25) sp<-cbind(sp1,sp2,sp3,sp4) #Criação da matriz de contrastes Treat_Imp<-model.matrix(~ treat-1) Treat_Imp1_2<-Treat_Imp[, 1] - Treat_Imp[, 2] Treat_Imp2_3<-Treat_Imp[, 2] - Treat_Imp[, 3] Treat_Imp1_3<-Treat_Imp[, 1] - Treat_Imp[, 3] Eu acabei criando manualmente as comparações, através dos objetos Treat_Imp1_2, Treat_Imp2_3 e Treat_Imp1_3. Mas alguém poderia me dar uma ideia de como sistematizar isso, e ter a criação desses objetos de uma maneira direta e que não envolva pacotes como o multcomp? Obrigado, -- ====================================================================== 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 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO) 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 Researchgate: https://www.researchgate.net/profile/Alexandre_Santos10 LinkedIn: https://br.linkedin.com/in/alexandre-dos-santos-87961635 ======================================================================

2 3

I Congresso de Aplicações de R em Administração
by Celina Rebello 11 Mar '16

11 Mar '16

O I Congresso de Aplicações de R em Administração - CARA - visa disseminar o universo de estudos na área de Administração que utilizam a linguagem R como forma de alavancar a eficiência e qualidade dos trabalhos elaborados. Haverá espaço para exposição de trabalhos das diversas áreas de Administração, que tenham sido desenvolvidos utilizando a linguagem R. O foco não deverá ser a discussão do código desenvolvido, mas sim a análise proposta, e do poder estatístico oferecido pela linguagem. Interessados devem submeter até três artigos pelo sistema disponível neste site, observando o calendário do evento. O CARA é destinado tanto a profissionais que nunca escreveram um código em R, bem como a profissionais experientes que desejarem conhecer trabalhos que têm sido desenvolvidos utilizando-se desta poderosa linguagem. Serão também oferecidas duas atividades: Workshop Análise de Business Intelligence, para profissionais que possuem um conhecimento básico de R; Minicurso Apresentação do R, voltado àqueles que desejarem ter o primeiro contato com a linguagem. Submeta seus trabalhos em: http://cara.coppead.ufrj.br/submissao-de-trabalhos Saiba Mais em :http://cara.coppead.ufrj.br/home

1 0

Bons valores iniciais e não convergência - regressão não linear
by Adriele Pereira 08 Mar '16

08 Mar '16

Olá, bom dia a todos. Estou tentando ajustar o modelo não linear Richards ao seguinte conjunto de dados: y<-c(32.2,45.1,63.3,87.7,107.6,120.4,131.8,148.9,159.3,166.8,173.4,187.4,195.5,198.8,211.0,221.4,226.6,231.2,237.7,243.1,246.6,250.3) x<-c(90,180,270,360,450,540,630,720,810,900,990,1080,1170,1260,1350,1440,1530,1620,1710,1800,1890,1980) Ao utilizar o shiny, obtive os seguintes valores iniciais: alpha=270 , beta= 402, k=0.0016 , d=0.8. O problema é que o gráfico fica muito bom com esses valores, mas ao colocar esse valores iniciais na nls aparece a mensagem de erro: ajuste=nls(y~(alpha*(1 + (d - 1)*exp(-k*(x-beta)))^(1/(1-d))), start=c(alpha= 270, beta= 402 , k= 0.0016, d=0.8)) Error in numericDeriv(form[[3L]], names(ind), env) : Obtido valor faltante ou infinito quando avaliando o modelo Já tentei mudar a parametrização do modelo e aumentar o número de iterações e mesmo assim não ajustou. Alguém já se deparou com um problema parecido? Desde já agradeço.

2 1

Soma com intervalo definido dentro de um loop
by Yury Duarte 07 Mar '16

07 Mar '16

Bom dia colegas de R! Preciso somar uma serie de dados com tamanho definido igual a "a" repetições. Tentei usar o comando for() para realizar essa tarefa e somar os valores de x1 até xa; x2 até x(a+1) e assim sucessivamente até xn. Usando a estrutura básica da função soma fora do comando for não tive problemas: sum(Nasa$PP.Liquida[1:140]) Quando tentei aplicar esse formato no comando for não consegui ter sucesso: PP = for(i in Nasa$PP.Liquida){ sum(Nasa$PP.Liquida[i:(i+139)]) } Acredito que possa ser um problema de entendimento conceitual do funcionamento desse comando. Alguém poderia me ajudar nessa questão? Desde já, agradeço a todos pelo apoio! Att Yury Duarte Engenheiro Agrônomo - ESALQ/USP

3 3

Aumentar o tamanho de fontes em gráfico ggplot2
by Marcelo Laia 05 Mar '16

05 Mar '16

Olá, Rodei o exemplo de gráfico disponível em http://www.r-bloggers.com/how-to-plot-three-categorical-variables-and-one-c… Tive dois problemas: 1. não consegui incluir o Título do gráfico. É informado que a função opts não existe. 2. as fontes ficaram muito pequenas. Se eu insiro a função geom_text(size = 12, ....aes(....)) dá erro. Alguém saberia em que local do código eu posso alterar as fontes globalmente (de todo o gráfico)? Obrigado. -- Marcelo

2 1

Re: [R-br] problema ao conveter data:hora armazenado em texto (função strptime)
by Éder Comunello 04 Mar '16

04 Mar '16

Paulo, bom dia! Não consegui reproduzir o erro. Qual a saída de str(lista.exames)? ### <code r> lista.exames <- structure(list(id = 1:3, `Data de Liberação:` = c("10/06/2007 01:03:34", "10/06/2007 01:03:34", "10/06/2007 01:03:34")), .Names = c("id", "Data de Liberação:"), row.names = c(NA, -3L), class = "data.frame") str(lista.exames) # 'data.frame': 3 obs. of 2 variables: # $ id : int 1 2 3 # $ Data de Liberação:: chr "10/06/2007 01:03:34" "10/06/2007 01:03:34" "10/06/2007 01:03:34" class(lista.exames[,"Data de Liberação:"]) # [1] "character" dt.hr <- lista.exames[,"Data de Liberação:"][1:3]; dt.hr # [1] "10/06/2007 01:03:34" "10/06/2007 01:03:34" "10/06/2007 01:03:34" strptime(dt.hr, format = "%d/%m/%Y %H:%M:%S", tz="America/Sao_Paulo") # [1] "2007-06-10 01:03:34 BRT" "2007-06-10 01:03:34 BRT" "2007-06-10 01:03:34 BRT" ### </code> ================================================ Éder Comunello Agronomist (UEM), MSc in Environ. Sciences (UEM) DSc in Agricultural Systems Engineering (USP/Esalq) Brazilian Agricultural Research Corporation (Embrapa) Dourados, MS, Brazil |<O>| ================================================ GEO, -22.2752, -54.8182, 408m UTC-04:00 / DST: UTC-03:00 Em 3 de março de 2016 17:03, Paulo Nogueira Starzynski <paulons(a)gmail.com> escreveu: > Pessoal, > estou tentando fazer uma conversão simples de data e hora armazena em um > vetor de caracters para formato "POSIXlt" "POSIXt". > Porém, está acontecendo algo que não entendo ao usar a função strptime. > > Se eu passo o vetor com os dados diretamete, retorna NA. > Se eu copio e colo os valores na função, dá certo... ou o mesmo usando o > dput, como abaixo... vejam... > > > > class(lista.exames[,"Data de Liberação:"])[1] "character"> dt.hr <- lista.exames[,"Data de Liberação:"][1:10]> strptime(dt.hr, format = "%d/%m/%Y %H:%M:%S") [1] NA NA NA NA NA NA NA NA NA NA > > > # não deu certo > > > dput(lista.exames[,"Data de Liberação:"][1:10])c("10/06/2007 01:03:34", "10/06/2007 01:03:34", "10/06/2007 01:03:34", > "10/06/2007 01:03:34", "10/06/2007 01:03:34", "10/06/2007 01:03:34", > "10/06/2007 01:03:34", "10/06/2007 01:03:34", "10/06/2007 01:03:34", > "10/06/2007 01:03:34")> dt.hr <- c("10/06/2007 01:03:34", "10/06/2007 01:03:34", "10/06/2007 01:03:34",+ "10/06/2007 01:03:34", "10/06/2007 01:03:34", "10/06/2007 01:03:34",+ "10/06/2007 01:03:34", "10/06/2007 01:03:34", "10/06/2007 01:03:34",+ "10/06/2007 01:03:34")> strptime(dt.hr, format = "%d/%m/%Y %H:%M:%S") [1] "2007-06-10 01:03:34 BRT" "2007-06-10 01:03:34 BRT" "2007-06-10 01:03:34 BRT" > [4] "2007-06-10 01:03:34 BRT" "2007-06-10 01:03:34 BRT" "2007-06-10 01:03:34 BRT" > [7] "2007-06-10 01:03:34 BRT" "2007-06-10 01:03:34 BRT" "2007-06-10 01:03:34 BRT" > [10] "2007-06-10 01:03:34 BRT" > > > > # deu certo!! > > O que pode estar ocorrendo? Não basta o meu vetor original ser do tipo > character? > > > Abraços, > Paulo > > _______________________________________________ > 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. >

1 0

Verificar duplicados na integra
by Wagner Tassinari 04 Mar '16

04 Mar '16

Olá pessoal, Gostaria de verificar as duplicidade na integra, o comando "duplicated" mostra os duplicados a partir do segundo registro. Mas eu gostaria de identificar desde o primeiro registro as duplicidade. Por exemplo, o que ocorre fazendo o comando "duplicated": df <- read.table(header=TRUE, text=' label value A 4 B 3 C 6 B 3 B 1 A 2 A 4 A 4 ') df = data.table(df, key = "value") df$duplic = duplicated(df, by="label") df label value duplic 1: B 1 FALSE 2: A 2 FALSE 3: B 3 TRUE 4: B 3 TRUE 5: A 4 TRUE 6: A 4 TRUE 7: A 4 TRUE 8: C 6 FALSE Mas na realidade gostaria que aparecesse da seguinte forma: df label value duplic 1: B 1 TRUE 2: A 2 TRUE 3: B 3 TRUE 4: B 3 TRUE 5: A 4 TRUE 6: A 4 TRUE 7: A 4 TRUE 8: C 6 FALSE Obrigado pela ajuda de vocÊs galera. ABS -- Wagner S. Tassinari Departamento de Matemática Universidade Federal Rural do Rio de Janeiro. BR-465, Km 7 - Seropedica, RJ - Brasil CEP: 23890-000 Skype: wagner.tassinari wtassinari(a)gmail.com tassinari(a)ufrrj.br ------------------------------------------------------- "Statistical thinking will one day be as necessary for efficient citizenship as the ability to read and write." (H.G.Wellis)

3 2