- R-br - DINF UFPR

Gráfico
by Antonio Armando Ulian do Lago Albuquerque 01 Out '14

01 Out '14

Caros amigos, Sou muito novato no R, mas tenho nos últimos 6 meses lido muito e tentando realizar algumas operações muito básicas. Infelizmente ainda estou gatinhando. Utilizo o UBUNTU 14.04, 64bits, RKward. A partir do banco de dados do TSE, baixei os arquivos referentes a raça nas eleições 2014. Como não consegui abri-lo na extensão csv transformei-os em xlsx e utilizei o pacote XLConnect. Quero elaborar 2 gráficos com os dados, mas não estou conseguindo. As tabelas são as seguintes. Quais comandos devo elaborar? Summary Raça Qtde Percent 1 Amarela 120 0.46 2 Branca 14379 54.96 3 Indígena 86 0.33 4 Parda 9158 35.00 5 Preta 2422 9.26 6 Total 26165 100.00 Raça Qtde Percent Length:6 Min. : 86.0 Min. : 0.33 Class :character 1st Qu.: 695.5 1st Qu.: 2.66 Mode :character Median : 5790.0 Median : 22.13 Mean : 8721.7 Mean : 33.34 3rd Qu.:13073.8 3rd Qu.: 49.97 Max. :26165.0 Max. :100.00 Muito obrigado e desculpe pela pergunta e dificuldade tão raza a qual vos exponho. Antonio

4 4

[R-br] Parâmetros_=?UTF-8?Q?modelo_de_regress=C3?==?UTF-8?Q?=A3o_din=C3=A2mico?=
by Fatima do Nascimento Silva 30 Set '14

30 Set '14

olá caros colegas, boa noite! estou ajustando um modelo de regressão dinâmico e preciso ponderar o vetor de parâmetros beta 2 (b2) por 7, que representaria um aumento de 7° na variável temperatura. E depois faria as previsões inserindo estes novos valores de b2 no calculo dos valores ajustados. Alguém poderia me dar uma dica de como alterar isso no R? segue os comandos utilizados. desde já agradeço a atenção. ponderar o beta 2 (L(Temperatura, 1) = -4.8725) por 7. caso possa me ajudar com alguma sugestão agradeço muito. #Comandos (CMR) rm(list=ls(all=TRUE)) #pacotes utilizados require(dynlm) require(tseries) require(car) require(verification) require(MASS) require(forecast) require(TSA) require(nortest) ################################################################## (dados=read.table(file="BDAC.txt",header=T,dec=",")) attach(dados) (dados1<-ts(dados,frequency=12, start=c(2000,1),end=c(2010,12))) #modelo ajustado MOD1 <- dynlm(BEAI~trend(BEAI)+season(BEAI)+L(Temperatura,1)+L(Poluente,3),data = dados1 ) summary(MOD1) par(mfrow = c(2,2)) plot(MOD1, main="") names(MOD1) #lista de objetos do modelo #previsões p<-(forecast(MOD1$fitted.values, h=12, level=c(80,95), fan=FALSE)) plot(p) #saída do modelo > summary(MOD1) Time series regression with "ts" data: Start = 2000(4), End = 2010(12) Call: dynlm(formula = BEAI ~ trend(BEAI) + season(BEAI) + L(Temperatura, 1) + L(Poluente, 3), data = dados1) Residuals: Min 1Q Median 3Q Max -42.946 -13.816 -1.482 12.327 80.586 Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) 263.2442 45.6925 5.761 7.21e-08 *** trend(BEAI) 0.4470 0.6201 0.721 0.47241 season(BEAI)Feb -19.3407 9.5712 -2.021 0.04565 * season(BEAI)Mar 9.0604 9.6612 0.938 0.35032 season(BEAI)Apr 5.1280 9.4526 0.543 0.58853 season(BEAI)May 25.2990 9.6410 2.624 0.00988 ** season(BEAI)Jun 35.7033 11.7395 3.041 0.00292 ** season(BEAI)Jul 30.0824 11.3742 2.645 0.00933 ** season(BEAI)Aug 19.1968 11.9829 1.602 0.11192 season(BEAI)Sep 4.6108 10.3265 0.447 0.65608 season(BEAI)Oct 3.7939 10.4920 0.362 0.71832 season(BEAI)Nov -1.3031 9.4690 -0.138 0.89078 season(BEAI)Dec -7.6315 9.4334 -0.809 0.42020 L(Temperatura, 1) -4.8725 1.5284 -3.188 0.00185 ** L(Poluente, 3) 2.6678 0.8455 3.155 0.00205 ** --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 Residual standard error: 21.34 on 114 degrees of freedom Multiple R-squared: 0.618, Adjusted R-squared: 0.5711 F-statistic: 13.17 on 14 and 114 DF, p-value: < 2.2e-16 Atenciosamente, Fátima Nascimento Estatística - CONRE nº 9439 Mestre em Ciência e Engª de Petróleo-PPGCEP/UFRN Doutoranda no Programa de Pós-Graduação em Ciência e Engª de Petróleo- PPGCEP/UFRN

3 5

Números aleatórios
by Leonard de Assis 29 Set '14

29 Set '14

Existe algo pronto para gerar n sequencias de números? Por exemplo: gerar n sequências de 1 a 5, sendo que cada número tem uma probabilidade específica de cair em cada posição. Tenho mais ou menos o algoritmo na cabeça, mas gostaria de saber se já existe, para evitar retrabalho

5 5

Configuração de cluster para programação em paralelo
by Clécio da Silva Ferreira 29 Set '14

29 Set '14

Olá. Estou tentando montar no nosso laboratório, que já está em ligado em rede, uma estrutura de cluster. O objetivo é conseguir usar os pacotes snow e parallel do R. O problema inicial é na configuração, isto é, usar a função makeCluster. Digitamos, por exemplo, makeCluster("maquina1","maquina2",type="SOCK"), mas o R fica 'aguardando' possivelmente alguma resposta dos 'slaves'. Alguém já tentou e conseguiu trabalhar com este(s) pacote(s)? Atenciosamente, Prof. Clécio da Silva Ferreira Depto de Estatística - UFJF Tel: (32) 2102-3306 Página: www.ufjf.br/clecio_ferreira

3 3

desdobramento de uma analise fatorial
by helcio duarte pereira 29 Set '14

29 Set '14

Pessoal, estou tentando desdobrar uma interação pela função aov mas só dá erro. Estou trabalhando com o objeto "ex5" da função "ExpDes", cuja estrutura basica é: str(ex5) 'data.frame':160 obs. of 4 variables: $ trat : Factor w/ 4 levels "10g","15g","15t",..: 1 1 1 1 1 $ genero: Factor w/ 2 levels "F","M": 1 1 1 1 1 1 1 1 1 1 ... $ bloco : Factor w/ 20 levels "1","2","3","4",..: 1 2 3 4 5 6 7 $ sabor : int 8 9 7 9 8 8 8 5 6 7 ... o modelo de analise que usei foi: anava<-aov(formula = sabor~trat*genero + bloco, data = ex5) Response: sabor Df Sum Sq Mean Sq F value Pr(>F) trat 3 19.37 6.4562 1.9452 0.12537 genero 1 7.66 7.6563 2.3068 0.13118 bloco 19 97.82 5.1484 1.5512 0.07832 . trat:genero 3 3.92 1.3062 0.3936 0.75783 Residuals 133 441.43 3.3190 Eu sei que a interação foi não significativa, mas estou tentando desdobra-la para aprender, já que irei precisar fazer uma analise assim em breve e que provavelmente vai dar interação significativa O primeiro desdobramento(genero dentro de tratamento) rodou sem problemas: desd1<-aov(formula = ex5$sabor~ex5$trat/ex5$genero + bloco, data = ex5) summary.aov(object = desd1, split = list("ex5$trat:ex5$genero"= list("10g"=1, "15g"=2, "15t"=3, "20t"=4))) Porém o segundo: desd2<-aov(formula = ex5$sabor~ex5$genero/ex5$trat + bloco, data = ex5) summary.aov(object = desd2, split = list("ex5$genero:ex5$trat"= list("F"=1, "M"=2))) Não roda de maneira alguma. Usei o mesmo procedimento do anterior, assim, alguém sabe no que estou errando? Att. Helcio

2 1

Modelos multinível com variância complexa
by Leonardo Ferreira Fontenelle 27 Set '14

27 Set '14

Boa tarde, pessoal! Tive um curso introdutório aos modelos multinível, que utilizou quase exclusivamente o MLwiN. Procurei treinar em paralelo o ajuste dos modelos com lme4, o que não foi difícil para modelos com variância simples. Mas, quando o negócio chegou na variância complexa, não consegui encontrar correspondência no lme4 do que me foi ensinado no MLwiN. Nem sei se não seria necessário utilizar nlme. Um exemplo de eventual necessidade seria em algum estudo longitudinal, onde a variabilidade das medidas (por exemplo, altura; primeiro nível) de cada pessoa (segundo nível) poderia mudar ao longo do tempo. Outro exemplo seria verificar se a variabilidade da variável de reposta é diferente conforme alguma variável explanatória dicotômica ou não. Onde encontro material sobre como fazer isso no R? Obrigado, [1]Leonardo Ferreira Fontenelle References 1. http://lattes.cnpq.br/9234772336296638

2 4

fatorial com efeito aleatório
by Maurício Lordêlo 26 Set '14

26 Set '14

Caros "listeiros" do R, Estou precisando analisar os dados de um experimento já realizado, infelizmente sem a prévia consulta a um pesquisador da área de Estatística. O experimento foi realizado com explantes de bananeiras de modo que foram levados em consideração três fatores: 1. Genótipo, com dois níveis; 2. Doses ou meios, com sete diferentes concentrações; 3. Tipos de cortes do explante, com dez diferentes tipos. O objetivo é verificar qual tipo de corte induziu melhor o explante na formação do embrião considerando os setes meios e os dois genótipos. Para isso foram preparadas duas placas (caracterizando o número de repetições) com uma quantidade de explantes idênticos de acordo com o corte e o genótipo (esta quantidade de explantes em cada placa variou de 6 a 11). Cada placa com estes explantes foi introduzida nos sete diferentes meios (doses). As respostas foram obtidas registrando-se, em relação ao total de cada placa, quantos explantes oxidaram, quantos formaram calos, quantos não desenvolveram, quantos apresentaram fungos, quantos apresentaram bactéria, etc. Os dados com apenas uma variável resposta estão neste link: https://www.dropbox.com/s/23x1dr0qxoj6ghm/explantes.txt?dl=0 Comecei a esboçar um ajuste de modelo, porém surgiram dúvidas se estou fazendo de forma correta. explantes=read.table("explantes.txt", sep="", header=TRUE) attach(explantes) str(explantes) meio=as.factor(meio1) tipo=as.factor(tipo1) prop=(n_oxidaram /n_explant) tapply(prop,genotipo,mean) tapply(prop,meio,mean) tapply(prop,tipo,mean) resp=cbind(n_oxidaram, n_explant - n_oxidaram) ajuste1=glm(resp~(genotipo/meio):tipo,family=binomial()) summary(ajuste1) Um fatorial 2 x 7 x 10 reconheço que não é adequado. Posso pensar em "tipo" como efeito aleatório? Se sim, qual seria o comando para este caso e como interpretaria os resultados? Existe alguma uma outra opção para a análise? Agradeço qualquer tipo de ajuda. Maurício

2 4

criando variável idade
by Taynãna César Simões 26 Set '14

26 Set '14

Oi, pessoal?? Tenho uma variável codificadora de idade com quatro dígitos, em que: 1 - Hora 2 - Dia 3 - Mês 4 - Ano Assim, por exemplo: 3009 - nove meses, 4018 - dezoito anos, etc Quero criar uma variável que retorne apenas a idade em anos. Obrigada, Taynãna SAMeQ/FIOCRUZ

4 3

[R-br] modelo de regressão dinâmico
by Fatima do Nascimento Silva 26 Set '14

26 Set '14

olá caros colegas, boa tarde! estou ajustando um modelo de regressão dinâmico e preciso ponderar o vetor de parâmetros beta 2 (b2) por 7, que representaria um aumento de 7º na temperatura. E depois faria as previsões inserindo estes novos valores de b2. Alguém poderia me dar uma dica de como alterar isso no R? segue os comandos utilizados. desde já agradeço a atenção. MOD1 <- dynlm(BEAI~trend(BEAI)+season(BEAI)+L(Temperatura,1)+L(Poluente,3),data = dados1 ) summary(MOD1) names(MOD1) p<-(forecast(MOD1$fitted.values, h=12, level=c(80,95), fan=FALSE)) plot(p) Atenciosamente, Fátima Nascimento Estatística - CONRE nº 9439 Mestre em Ciência e Engª de Petróleo-PPGCEP/UFRN Doutoranda no Programa de Pós-Graduação em Ciência e Engª de Petróleo- PPGCEP/UFRN

2 1

[R-br] modelo de regressão dinâmico
by Fatima do Nascimento Silva 26 Set '14

26 Set '14

olá caros colegas, boa tarde! estou ajustando um modelo de regressão dinâmico e preciso ponderar o vetor de parâmetros beta zero (b0) por 7, que representaria um aumento de 7º na temperatura. E depois faria as previsões inserindo estes novos valores de b0. Alguém poderia me dar uma dica de como alterar isso no R? segue os comandos utilizados. desde já agradeço a atenção. MOD1 <- dynlm(BEAI~trend(BEAI)+season(BEAI)+L(Temperatura,1)+L(Poluente,3),data = dados1 ) summary(MOD1) names(MOD1) p<-(forecast(MOD1$fitted.values, h=12, level=c(80,95), fan=FALSE)) plot(p) Atenciosamente, Fátima Nascimento Estatística - CONRE nº 9439 Mestre em Ciência e Engª de Petróleo-PPGCEP/UFRN Doutoranda no Programa de Pós-Graduação em Ciência e Engª de Petróleo- PPGCEP/UFRN

1 0

Kruskal Wallis - identificar o tratamento diferente
by Daniela Recchia 26 Set '14

26 Set '14

Boa tarde! Já pesquisei um pouco sobre isso mas não encontrei nada, tenho 3 variáveis X,Y,Z (não normais) e estou usando o kruskal-wallis para compara-las. Recebo então que elas são estatisticamente diferentes mas não consigo saber qual delas, e a X diferente da Y, é a X diferente da Z, é a Y diferente da Z? Alguém sabe se existe alguma nova função que fornece essa informação extra na anova não paramétrica? Abraços, Daniela -- Daniela Rodrigues Recchia M.Sc. in Statistics Tel. (Mobil): +49 01578 7588382 Germany “The standard of success in life is absolutely the amount of joy you feel” Esther & Jerry Hick

3 3

Re: [R-br] Criando variável idade
by Taynãna César Simões 26 Set '14

26 Set '14

Mas como separar esta coluna unica, como nos exemplos, em quatro colunas? Obrigada ------------------------------ Em sex, 26 de set de 2014 11:00 BRT Fernando Souza escreveu: >Não é melhor você colocar cada um desses codificadores como uma coluna >em seu data frame? Fica mais simples para selecionar subgrupos através >dos argumentos subset ou comandos de navegação de data frame > >data.frame(Hora=identificador_Hora, Dia= identificador_Dia, Mes= >identificador_Mes, Ano= identificador_Ano) >Onde identificador_.... é um vetor contendo a identificação da hora, >dia, mes e ano respectivamente >On 26-09-2014 09:38, Taynãna César Simões wrote: >> Bom dia pessoal, >> >> Tenho uma variável codificadora de idade com quatro dígitos, em que: >> 1 - Hora >> 2 - Dia >> 3 - Mês >> 4 - Ano >> Assim, por exemplo: 3009 - nove meses, 4018 - dezoito anos, etc >> >> Que criar uma variável que retorne apenas a idade em anos. >> >> Obrigada, Taynãna >> CPqRR/FIOCRUZ >> >> >> _______________________________________________ >> 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. >

2 1

Criando variável idade
by Taynãna César Simões 26 Set '14

26 Set '14

Bom dia pessoal, Tenho uma variável codificadora de idade com quatro dígitos, em que: 1 - Hora 2 - Dia 3 - Mês 4 - Ano Assim, por exemplo: 3009 - nove meses, 4018 - dezoito anos, etc Que criar uma variável que retorne apenas a idade em anos. Obrigada, Taynãna CPqRR/FIOCRUZ

2 1

Pré-processamento de texto pacote 'tm'
by Fernando Gama 26 Set '14

26 Set '14

Bom dia pessoal, Estou com problemas na etapa de *preenchimento do termo *da minha análise textual, o s.o que utilizo é o mac os x 10.8.5. segue o trecho do código: tdm_corpus <- Corpus(VectorSource(corpus)) tdm_corpus<-tm_map(tdm_corpus, stemDocument, language = "portuguese") *tdm_corpus<-tm_map(tdm_corpus, stemCompletion, dictionary=corpus.copy) * As outras etapas de limpeza foram aplicadas e o stem foi realizado com sucesso! Porém ao executar a linha em negrito o seguinte erro é lançado: Warning message: In mclapply(content(x), FUN, ...) : all scheduled cores encountered errors in user code Reforço que no windows o processo de preenchimento do termo fora realizado normalmente. Além do mais, de acordo com minhas pesquisas esta maneira de preenchimento do termo é muito usual e o "warning" é muito genérico, isso acabou me motivando a buscar ajuda de vcs! Grande abraço. Fernando. -- Att, Fernando Gama da Mata Pesquisador - Instituto Tecnológico Vale

2 2

Mendeley gerenciador de referência
by Fernando Souza 25 Set '14

25 Set '14

Caros amigos Recentemente instalei o gerenciador de referências Mendeley em meu computador ubuntu 14.04 32 bits. Acontece que o plugin para fazer citações/refências no writer não funciona ele retorna a seguinte : "Problem formatting citation #0 Unable to render citations due to a bug in the installed version of QtWebKit" buscando na internet encontrei o seguinte: http://support.mendeley.com/customer/portal/questions/1492956-the-qtwebkit-… Alguém que é usuário do mendeley, passou por este problema? Como resolveu? Pode me ajudar? Eu não quero ter de instalar ubuntu 64 bits em meu computador. Problema também ocorre com o zoterro Standalone que nao sincroniza as referência que coloquei utilizando o zotero plugin Firefox. Agraço a atenção

1 1

remover linhas
by Leonardo Monteiro 25 Set '14

25 Set '14

Pessoal, considerando as 4 situacoes abaixo: ###1 x<- c(1,2,3,4,5) y<- c(6,7,8,9,0) ###2 x<- c(1,2,NA, 4,5) y<- c(6,7,8,9,0) ###3 x<- c(1,2,3,4,5) y<- c(6,7,NA,8,9) ####4 x<- c(1,2,NA,4,5) y<- c(6,NA,7,8,9) dados<- cbind(x,y) eu estou tentando usar o seguinte comando p remover as linhas com NA´s ## dados<- ifelse(is.na(dados[,2])==F|is.na(dados[,1])==F, dados, dados[-which(is.na(dados[,1])==T|is.na(dados[,2])==T),]) dados ### entretanto nao esta funcionando, pois em ambas situacoes, ele me retorna apenas o primeiro vetor do obj dados. alguem poderia me ajudar, pq preciso apenas das linhas q nao contem NA/ou entao manter o objeto dados, qdo nao houver NA. obrigado, -- Leonardo Amaral Monteiro - Eng° Agr° MsC Fìsica do Ambiente Agrícola - Agrometeorologia (ESALQ/USP) Doutorando em "Engenharia de Sistemas Agrícolas" - Agrometeorologia (ESALQ/USP) Piracicaba, São Paulo. Telefone - (19) 3429-4123 r.236

3 2

[R-br] Modelo de regressão dinâmico-alterar b1
by Fatima do Nascimento Silva 25 Set '14

25 Set '14

Olá caros colegas, boa noite! Bem... Estou ajustando um modelo de regressão dinâmico e, após gerar o modelo gostaria de manter os outros parâmetros fixos e ponderar o vetor b1 por 7 e com este novo vetor fazer as previsões. Alguém teria alguma dica de como mudo apenas o vetor b1 e faço as previsões? Obrigada pela atenção? Atenciosamente, Fátima Nascimento Estatística - CONRE nº 9439 Mestre em Ciência e Engª de Petróleo-PPGCEP/UFRN Doutoranda no Programa de Pós-Graduação em Ciência e Engª de Petróleo- PPGCEP/UFRN

1 0

declaração "estranha em função"
by Vinicius Brito Rocha 24 Set '14

24 Set '14

Senhores, estou estudando hadoop/map reduce no R e achei um código na pagina do RevolutionR para o Kmeans. Estou tentando entender e ei que me surgiu uma dúvida... Existe uma declaração de função para a fase de map onde o parâmetro de input é um "." Alguém sabe o que significa? Esta função é chamada dentro de outra função. O que pensei é que o "." simboliza que a sub função herda os parâmetros da "função mãe". como está declarada a função: function(., P) ### código inteiro #### # Copyright 2011 Revolution Analytics # # Licensed under the Apache License, Version 2.0 (the "License"); # you may not use this file except in compliance with the License. # You may obtain a copy of the License at # # http://www.apache.org/licenses/LICENSE-2.0 # # Unless required by applicable law or agreed to in writing, software # distributed under the License is distributed on an "AS IS" BASIS, # WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. # See the License for the specific language governing permissions and # limitations under the License. Sys.setenv(HADOOP_STREAMING="/usr/lib/hadoop-mapreduce/hadoop-streaming-2.0.0-cdh4.7.0.jar") library(rmr2) library(rhdfs) hdfs.init() # P tabela com valores ## @knitr kmeans-signature kmeans.mr = function(P, num.clusters, num.iter, combine, in.memory.combine) { ## @knitr kmeans-dist.fun # C : Centroides dist.fun = function(C, P) {apply(C,1, function(x) {colSums((t(P) - x)^2))}} ## @knitr kmeans.map *kmeans.map = function(., P) {* nearest = {if(is.null(C)) sample(1:num.clusters, nrow(P), replace = TRUE) else {D = dist.fun(C, P) nearest = max.col(-D)}} if(!(combine || in.memory.combine)) keyval(nearest, P) else {keyval(nearest, cbind(1, P))}} ## @knitr kmeans.reduce kmeans.reduce = { if (!(combine || in.memory.combine) ) function(., P) t(as.matrix(apply(P, 2, mean))) else function(k, P) keyval( k, t(as.matrix(apply(P, 2, sum))))} ## @knitr kmeans-main-1 C = NULL for(i in 1:num.iter ) { C = values( from.dfs( mapreduce( P, map = kmeans.map, reduce = kmeans.reduce))) if(combine || in.memory.combine) C = C[, -1]/C[, 1] ## @knitr end # points(C, col = i + 1, pch = 19) ## @knitr kmeans-main-2 if(nrow(C) < num.clusters) { C = rbind( C, matrix( rnorm( (num.clusters - nrow(C)) * nrow(C)), ncol = nrow(C)) %*% C) }} C} ## @knitr end ## sample runs ## out = list() for(be in c("local", "hadoop")) { rmr.options(backend = be) set.seed(0) ## @knitr kmeans-data P = do.call(rbind, rep(list(matrix(rnorm(10, sd = 10), ncol=2)),20)) + matrix(rnorm(200), ncol =2) ## @knitr end # x11() # plot(P) # points(P) out[[be]] = ## @knitr kmeans-run kmeans.mr(to.dfs(P),num.clusters = 12, num.iter = 5, combine = FALSE, in.memory.combine = FALSE) ## @knitr end } # would love to take this step but kmeans in randomized in a way that makes it hard to be completely reprodubile # stopifnot(rmr2:::cmp(out[['hadoop']], out[['local']])) -- *Vinicius Brito Rocha.*

1 0

[R-br] Modelo de regressão dinâmico-alterar b1
by Fatima do Nascimento Silva 24 Set '14

24 Set '14

Olá caros colegas, boa noite! Bem... Estou ajustando um modelo de regressão dinâmico e, após gerar o modelo gostaria de manter os outros parâmetros fixos e ponderar o vetor b1 por 7 e com este novo vetor fazer as previsões. Alguém teria alguma dica de como mudo apenas o vetor b1 e faço as previsões? Obrigada pela atenção? Atenciosamente, Fátima Nascimento Estatística - CONRE nº 9439 Mestre em Ciência e Engª de Petróleo-PPGCEP/UFRN Doutoranda no Programa de Pós-Graduação em Ciência e Engª de Petróleo- PPGCEP/UFRN

1 0

[R-br] Modelo de regressão dinâmico-alterar b1
by Fatima do Nascimento Silva 24 Set '14

24 Set '14

Olá caros colegas, boa noite! Bem... Estou ajustando um modelo de regressão dinâmico e, após gerar o modelo gostaria de manter os outros parâmetros fixos e ponderar o vetor b1 por 7 e com este novo vetor fazer as previsões. Alguém teria alguma dica de como mudo apenas o vetor b1 e faço as previsões? Obrigada pela atenção? Atenciosamente, Fátima Nascimento Estatística - CONRE nº 9439 Mestre em Ciência e Engª de Petróleo-PPGCEP/UFRN Doutoranda no Programa de Pós-Graduação em Ciência e Engª de Petróleo- PPGCEP/UFRN

1 0

Alocação de memória
by Ari Clecius 23 Set '14

23 Set '14

Prezados colegas da lista, estava tentando rodar este exemplo quando tive problema de alocação de memória, alguém poderia me ajudar a resolver? Att, Ari Clecius Alves de Lima Engenheiro Químico Me. Engenharia Civil (085)88412345 (085)33669042 #CMR processador: Intel(R) Core(TM) i3 CPU M380 @2.53GHz Memória instalada(RAM): 4,00 GB (utilizável:2,93) Sistema oeracional de 32 Bits R version 3.1.1 (2014-07-10) -- "Sock it to Me" Copyright (C) 2014 The R Foundation for Statistical Computing Platform: i386-w64-mingw32/i386 (32-bit) N <- 100000000 x <- runif(N) y <- runif(N) ok <- sum(x^2 + y^2 < 1) mc.pi <- 4 * ok / N sprintf("pi = %.6lf\n", mc.pi) #1ª tentativa >N <- 100000000 >x <- runif(N) > y <- runif(N) Error: cannot allocate vector of size 762.9 Mb In addition: Warning messages: 1: In runif(N) : Reached total allocation of 1535Mb: see help(memory.size) 2: In runif(N) : Reached total allocation of 1535Mb: see help(memory.size) 3: In runif(N) : Reached total allocation of 1535Mb: see help(memory.size) 4: In runif(N) : Reached total allocation of 1535Mb: see help(memory.size) > ok <- sum(x^2 + y^2 < 1) Error: cannot allocate vector of size 762.9 Mb In addition: Warning messages: 1: Reached total allocation of 1535Mb: see help(memory.size) 2: Reached total allocation of 1535Mb: see help(memory.size) 3: Reached total allocation of 1535Mb: see help(memory.size) 4: Reached total allocation of 1535Mb: see help(memory.size) > mc.pi <- 4 * ok / N Error: object 'ok' not found > sprintf("pi = %.6lf\n", mc.pi) Error in sprintf("pi = %.6lf\n", mc.pi) : object 'mc.pi' not found memory.size(2900) 2ª tentativa > N <- 100000000 > x <- runif(N) Error: cannot allocate vector of size 762.9 Mb > y <- runif(N) Error: cannot allocate vector of size 762.9 Mb > ok <- sum(x^2 + y^2 < 1) Error: cannot allocate vector of size 762.9 Mb > mc.pi <- 4 * ok / N Error: object 'ok' not found > sprintf("pi = %.6lf\n", mc.pi) Error in sprintf("pi = %.6lf\n", mc.pi) : object 'mc.pi' not found

2 1

alternativa à função xcorr do Matlab no R
by Paulo Nogueira Starzynski 23 Set '14

23 Set '14

Olá colegas de lista. Tenho em mãos um programa escrito para Matlab e estou tentando fazer uma versão para R. O projeto faz parte de uma migração de longo prazo que talvez façamos no longo prazo. Como muitos conhecem "de vista" os registros de eletrocardiograma vou contextualizar meu problema, que consiste em detectar o instante em que ocorrem os batimentos do coração. Tenho em mãos um script em Matlab que detecta os "picos do ECG" (que na verdade se chama complexo QRS). Ele faz isso baseado na correlação cruzada entre o vetor contendo o eletrocardiograma completo e um vetor menor contendo apenas um complexo QRS. A função xcorr do Matlab faz "correr" o vetor menor através do maior, apontando os picos de correlação no instante da coincidência (ocorrência do complexo QRS). Estou tentando fazer o mesmo no R com através da *stats:: ccf*, mas sem sucesso, porque ela exige séries de tamanhos iguais. Alguém já resolveu esse problema? Procurei bastante e não achei alternativa. Aceito sugestões e se for preciso posso enviar um CRM. Obrigado! *==================================* *Paulo Nogueira Starzynski, B.Sc.*Estudante de Mestrado em Fisiologia Geral Laboratório de Fisiologia Teórica (LFT) Instituto de Biociências Universidade de São Paulo

3 2

alinhamento de texto em grafico
by Leonardo Monteiro 23 Set '14

23 Set '14

Prezados, eu gostaria de saber se é possivel alinhar o texto no grafico, por exemplo, alinhando pelo valor do eixo 'x': x<- 1:10 y<- 1/x plot(x, y) text(6, 0.8, 'exemplo1') text(6, 0.9, 'exemplo 50') os textos 'exemplo1' e 'exemplo50' ficaram desalinhados.. é possivel alinhar "a esquerda", ou a primeira letra do texto por um valor do eixo x? obrigado,

2 1

duvidas ( rotacionar a escala no eixo y para que fique os numeros fiquem na horizontal
by Eric Leandro 22 Set '14

22 Set '14

Pessoal boa noite, Como faço para rotacionar a escala do eixo y para que os numeros aparecem horizontalmente em relacão ao eixo y. E como faço para colocar o eixo 0 na vertical do gráfico att. plot(declividademediam,TSS, xlab="Declividade média (m)",ylab="Conc SST (mg/L)",pch=16,xlim=range(0:1200),ylim=range(0:400)) Att Eric Pereira. Bsc. Engenharia Ambiental - UfV / Bsc Enviromental Engineer UFV (University of Viçosa - Brazil) Mestre em Clima e Ambiente - INPA / Msc Climate and Environment (INPA) Programa de Grande Escala da Biosfera-Atmosfera na Amazônia (LBA project)

2 1

duvidas ( rotacionar a escala no eixo y para que fique os numeros fiquem na horizontal
by Eric Leandro 22 Set '14

22 Set '14

Pessoal boa noite, Como faço para rotacionar a escala do eixo y para que os numeros aparecem horizontalmente em relacão ao eixo y. E como faço para colocar o eixo 0 na vertical do gráfico att. plot(declividademediam,TSS, xlab="Declividade média (m)",ylab="Conc SST (mg/L)",pch=16,xlim=range(0:1200),ylim=range(0:400)) Att Eric Pereira. Bsc. Engenharia Ambiental - UfV / Bsc Enviromental Engineer UFV (University of Viçosa - Brazil) Mestre em Clima e Ambiente - INPA / Msc Climate and Environment (INPA) Programa de Grande Escala da Biosfera-Atmosfera na Amazônia (LBA project)

1 0

Agrupar linhas pela primeira letra
by Fátima Lima Paula 19 Set '14

19 Set '14

Olá pessoal. Tenho um df assim: id=c("ana","pedro","maria","teresa","monica","jorge","mario","raquel","marcus") idade=c("67","68","76","84","60","62","79","81","80") diag=c("A110","A112","A134","A560","B123","B165","C123","D123","D345") df1=data.frame(id,idade,diag) Quero criar uma nova coluna chamada categ de forma que fique assim: os diag que começam com A ficam "categA" os diag que começam com B ficam "categB" E assim por diante. Alguém poderia me ajudar? Obrigada Fátima -- "Minha felicidade depende da qualidade dos meus pensamentos"

2 5

geracao de numeros aleatorios representativos de estacoes virtuais
by Leonardo Monteiro 18 Set '14

18 Set '14

Pessoal, tenho a ideia de gerar pares de numeros, como representativos de localizacao geografica (lat, long) para alguns locais do Brasil...de preferencia uma grade regular para um estado especifico do Brasil.. existe algum pacote que faça isso? obrigado, Leonardo

2 1

Plotar perfil de umidade no solo ao longo do ano
by Wenceslau Teixeira 18 Set '14

18 Set '14

Caros colegas do R, muito bacana o script disponibilizado para traçar os perfis de umidade do solo. Será útil para nós. Só gostaria de trocar as cores, com valores mais alto (umido) para azul e valores mais baixos (solo seco) para vermelho, assim o gráfico fica mais cognitivo. Alguém sabe como mudar as cores ? Cordialmente, Wenceslau Teixeira Embrapa Solos library(repmis) # reads text data directly from dropbox - no need to download any file # read data url <- 'https://dl.dropboxusercontent.com/u/27700634/precip.txt' tmp <- repmis::source_data(url, sep = '', header = TRUE) # convert julian day to date date <- as.Date(tmp$julian, origin='2011-12-31') data <- cbind(date, tmp) head(data) # define vector with depth of soil layers depths <- c(0.05,0.10,0.20,0.30, 0.40,0.60,0.80,1.00, 1.50,2.00,2.50) # teste com lattice cols <- ncol(data):4 levelplot(as.matrix(data[, cols]), aspect = "fill", xlab = "date", ylab = "soil Moist.", scales = list(x = list(at = which(date %in% pretty(date)), labels = pretty(date)), y = list(at = length(cols):1, labels = depths)), contour = T) # Plot soil moisture profile # como são profundidades (abaixo da superfície) use a convenção de z negativo # acrescentei a chuva, precisa ajustar ylabels, mtext() resolve x11() op <- par(mar = c(3,4,12,3)) plot <- filled.contour(date, sort(-depths), as.matrix(data[, cols]), ylab = expression(theta(m^3~~m^-3))) par(new=T,mar=c(23,4,1,8)) with(data, plot(date, precip, type = "h", xaxt = "n", xpd = T, xaxs="i", ylab = "Precip (mm)", xlab = "")) par(op) -- Wenceslau Teixeira wgt007(a)gmail.com skype: wgt007

4 4

Re: [R-br] Variáveis contínuas explicando Variáveis categóricas
by Jefferson Ferreira-Ferreira 18 Set '14

18 Set '14

Daniel e demais, Aparentemente não resolve, pois recebo NAs para diversos slopes. (ver summary do modelo abaixo). E acho que isso não pode estar certo. Segue aqui um exemplo dos meus dados. O que eu gostaria é de ver o quanto as variáveis categóricas explicam minha media_CBT e minha media_CCS. Alguma ideia??? agua_quente detergentes agua_trat Filtr_leite Higie_resf pre_dipping toalha pos_dipping vacuo pulsador Aliment_pos media_CBT media_CCS 1 0 0 1 0 1 1 1 1 1 1 931 623.3333333 1 0 0 1 1 0 1 1 1 1 1 163.6666667 259 1 0 0 1 1 1 1 1 0 0 0 690.3333333 306.3333333 1 0 0 1 1 1 1 1 0 0 1 690.3333333 306.3333333 summary(m) Call: lm(formula = dados$media_CBT ~ dados$agua_quente + dados$detergentes + dados$agua_trat + dados$Filtr_leite + dados$Higie_resf + dados$pre_dipping + dados$toalha + dados$pos_dipping + dados$vacuo) Residuals: 1 2 3 4 0.000e+00 -3.155e-30 4.019e-14 -4.019e-14 Coefficients: (7 not defined because of singularities) Estimate Std. Error t value Pr(>|t|) (Intercept) 4.043e+02 8.988e-14 4.499e+15 <2e-16 *** dados$agua_quente NA NA NA NA dados$detergentes NA NA NA NA dados$agua_trat NA NA NA NA dados$Filtr_leite NA NA NA NA dados$Higie_resf -2.407e+02 6.962e-14 -3.457e+15 <2e-16 *** dados$pre_dipping 5.267e+02 6.962e-14 7.565e+15 <2e-16 *** dados$toalha NA NA NA NA dados$pos_dipping NA NA NA NA dados$vacuo NA NA NA NA --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 Residual standard error: 5.684e-14 on 1 degrees of freedom Multiple R-squared: 1, Adjusted R-squared: 1 F-statistic: 4.872e+31 on 2 and 1 DF, p-value: < 2.2e-16 Em 18 de setembro de 2014 15:36, <dtiezzi(a)usp.br> escreveu: > Jefferson > > Montando o modele assim > > model <- lm(test$x ~ test$c1 + test$c2 + test$c3 + test$c4) > > não resolve seu problema? > > Daniel > > > ------------------------------ > > *De: *"Jefferson Ferreira-Ferreira" <jecogeo(a)gmail.com> > *Para: *r-br(a)listas.c3sl.ufpr.br > *Enviadas: *Quinta-feira, 18 de Setembro de 2014 13:33:57 > *Assunto: *Re: [R-br] Variáveis contínuas explicando Variáveis categóricas > > > > Professor daniel, > Obrigado pela resposta. Ainda não testei essa alternativa. Mas a > desvantagem dessa abordagem seria o grande número de modelos a serem > ajustados, visto que tenho 11 variáveis explicativas categóricas para duas > variáveis resposta numéricas contínuas. Isso me faria ajustar 22 modelos > que teriam como premissa a independencia entre as variáveis explicativas. > > O que quero dizer é: o quanto a combinação das variáveis categóricas c1 c2 > c3 c4 c5 c6... etc (biárias =0 ou 1) explicam meu x (numérico contínuo). > > Será que existe um modo de eu ajustar dois modelos? Tipo com a variável > resposta X em função de todas as varíaveis explicativas categóricas e outro > modelo com a variável resposta Y em função de todas as variáveis > explicativas categóricas? > > > > c1 c2 c3 c4 x y > 0 1 1 1 931 623 > 0 1 0 0 163 259 > 1 0 1 0 690 306 > 1 0 1 0 690 306 > > > -- > > *Jefferson Ferreira-Ferreira* > > Geógrafo – GEOPROCESSAMENTO IDSM | Coordenadoria de TI > > > Jefferson.ferreira(a)mamiraua.org.br > > *Instituto de Desenvolvimento Sustentável Mamirauá* > > Ministério da Ciência, Tecnologia e Inovação > > Telefone: +55 97 3343-9710 > > *Google Maps* - Mapas deste e-mail: > > Exibir mapa ampliado > <https://maps.google.com.br/maps?q=-3.355557,-64.731151&ll=-3.355471,-64.731…> > > > *Contatos particulares:* > *(55) 9615-0100* > > > > > Em 18 de setembro de 2014 12:56, Daniel Tiezzi <dtiezzi(a)usp.br> escreveu: > >> Você precisa fazer uma regressão. >> >> Segue um modelo >> >> # Regression analyses, standardized >> model1.z <- lm(scale(PE$endurance) ~ scale(PE$age)) >> summary(model1.z) >> confint(model1.z) >> >> model2.z <- lm(scale(PE$endurance) ~ scale(PE$activeyears)) >> summary(model2.z) >> confint(model2.z) >> >> model3.z <- lm(scale(PE$endurance) ~ scale(PE$age) + >> scale(PE$activeyears)) >> summary(model3.z) >> confint(model3.z) >> >> # Conduct a model comparison NHST to compare the fit of model2.z to the >> fit of model3.z >> anova(model2.z, model3.z) >> >> >> Acho que seria assim >> >> >> Daniel >> >> >> >> >> Daniel Tiezzi, MD, PhD >> Professor Associado >> 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)fmrp.usp.br >> >> >> >> >> >> On Sep 18, 2014, at 12:47 PM, Jefferson Ferreira-Ferreira < >> jecogeo(a)gmail.com> wrote: >> >> >> Prezados; >> >> Uma dúvida básica de um iniciante em análises estatísticas. >> Tenho uma série de variáveis categóricas binárias e duas variáveis >> contínuas. Eu gostaria de saber o quando minhas variáveis categóricas >> explicam minhas duas variáveis contínuas. Por exemplo: >> >> c1 c2 c3 c4 x y >> 0 1 1 1 931 623 >> 0 1 0 0 163 259 >> 1 0 1 0 690 306 >> 1 0 1 0 690 306 >> >> >> A pergunta é: o quanto a combinação das variáveis c explicam os valores >> de x e y? Ou, qual a correlação entre as variáveis c e as variáveis x e y? >> >> Podem me dar alguma ideia de análises possíveis? >> Obrigado. >> >> >> -- >> >> *Jefferson Ferreira-Ferreira* >> Geógrafo – GEOPROCESSAMENTO IDSM | Coordenadoria de TI >> >> Jefferson.ferreira(a)mamiraua.org.br >> *Instituto de Desenvolvimento Sustentável Mamirauá* >> >> Ministério da Ciência, Tecnologia e Inovação >> Telefone: +55 97 3343-9710 >> >> *Google Maps* - Mapas deste e-mail: >> >> Exibir mapa ampliado >> <https://maps.google.com.br/maps?q=-3.355557,-64.731151&ll=-3.355471,-64.731…> >> >> >> >> *Contatos particulares:* >> *(55) 9615-0100* >> _______________________________________________ >> 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. >> >> >> >> _______________________________________________ >> 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. >> > > > > > _______________________________________________ > 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. > > >

3 6

Retornar NA nos que deram errado
by Rodrigo Coster 18 Set '14

18 Set '14

Caros, Estou usando o pacote survey para calcular a mediana de alguns dados vindos de amostras complexas. O problema é que dependendo das variáveis que peço, algumas combinações apresentam menos de 2 observações válidas, impossibilitando o cálculo (apesar deu achar que 1 observação já seria suficiente, já que ela tem peso maior que 1, mas enfim...) Alguém saberia como fazer o comando não parar, mas sim retornar NA nesses valores? Segue um exemplo simplificado (estou usando os dados da amostra do censo demográfico) require(survey) exemplo <- data.frame(ID = 1:10, var = rnorm(10), grupo = factor(c('A', 'A', 'A', 'A', 'B', 'B', 'B', 'C', 'D', 'D')), peso = rchisq(10, 5)) exemplo[10, 2] <- NA amostra <- svydesign(ids = ~ ID, data = exemplo, weights = ~ peso) svyby(~ var, ~ grupo, amostra, svyquantile, quantiles = .5)

2 3

[OT]Neural Network
by Tropidurus Torquatus 18 Set '14

18 Set '14

Caros R-users, Alguém trabalha com redes neurais no R? Caso positivo gostaria de trocar algumas idéias e experiências. Se alguém estiver interessado me envie e-mail pvt para não encher a caixa de mensagem do grupo com um assunto off topic. Atenciosamente Thales

3 2

posicao de graficos em figura com multiplos graficos
by Leonardo Monteiro 17 Set '14

17 Set '14

Pessoal, eu tenho que plotar 10 graficos em um unico figura... e nao estou usando nenhum pacote especifico. para isso estou fazendo: par(mfrow = c(3,4)) plot... assim, as posicoes sao completadas, até que na terceira linha (ultima) fica 'sobrando' as duas ultimas colunas.. alguem poderia me ajudar a centralizar as 2 ultimos 'graficos' nas colunas 2 e 3, respectivamente? Obrigado, Leonardo.

4 3

Plotar perfil de umidade no solo ao longo do ano
by Thiago V. dos Santos 17 Set '14

17 Set '14

Olá pessoal, Estou tentando fazer um gráfico de contorno no R mas ainda não encontrei a função certa para isso. Basicamente, eu tenho uma simulação de umidade em 11 camadas do solo para cada dia do ano, e gostaria de mostrar isso na forma de um perfil de solo. Por favor leia os comentários ao reproduzir os meus dados com o código abaixo: library(repmis) # reads text data directly from dropbox - no need to download any file # read data url <- 'https://dl.dropboxusercontent.com/u/27700634/precip.txt' tmp <- repmis::source_data(url, sep = '', header = TRUE) # convert julian day to date date <- as.Date(tmp$julian, origin='2011-12-31') data <- cbind(date, tmp) head(data) # now, convert soil layers to matrix and transpose it mat <- t(as.matrix(data[, 4:14])) # essa é a matriz que eu quero plotar. Reparem que ela já está organizada como um "perfil", # com as linhas representando as camadas de solo e as colunas representando os dias do ano A minha primeira tentativa foi usar a função filled.contour do pacote graphics. Eu defino um vetor com o nome das camadas do solo e então eu tento plotar, mas a função diz que a matriz está com as dimensoes invertidas: # define vector with depth of soil layers depths <- c(0.05,0.10,0.20,0.30, 0.40,0.60,0.80,1.00, 1.50,2.00,2.50) # Plot soil moisture profile plot <- filled.contour(data$julian, depths, mat) #Error in .filled.contour(x, y, z, levels, col) : dimension mismatch Obviamente eu consigo transpor a matrix para forçar o gráfico, mas a figura final não é o que eu preciso - o perfil fica de "cabeça para baixo". plot <- filled.contour(data$julian, depths, t(mat)) Como essa função só aceita os eixos em ordem crescente, eu não consigo inverter o eixo y para mostrar a primeira camada do solo no topo do gráfico. Poranto, eu aprecio qualquer sugestão de funções e/ou pacotes para construir o gráfico que eu preciso. Saudações, -- Thiago V. dos Santos PhD student Land and Atmospheric Science University of Minnesota http://www.laas.umn.edu/CurrentStudents/MeettheStudents/ThiagodosSantos/ind… Phone: (612) 323 9898

4 5

Abrindo duas Shapes no R
by Paulo Pimenta 15 Set '14

15 Set '14

Prezados, Estou tentando abrir dois shapefiles ao mesmo tempo, ou seja, sobrepor um sobre o outro e não consigo. Abaixo esta parte do meu código: ******************************************************************************************************** library(maptools) #dados <- readShapeSpatial(file.choose()) dados <- readShapeSpatial("F:/qgis_1/exercicio_buffer_L1/estrada_2000.shp") dados_1 <- readShapePoints("F:/qgis_1/exercicio_buffer_L1/L1_SAD69UTMZ23.shp") summary(dados)summary (dados_1) plot(dados)points (dados_1) ******************************************************************************************************** Tentei usar o comando merge e juntar as duas shapes, mas ainda assim não dá certo. Por outro lado não sei como plotar os dois shapes. Agradeço desde já a ajuda! Abraços.

3 2

set locale no Linux
by Manoel Galdino 15 Set '14

15 Set '14

Olá, Eu tentei mudar o local do R (time) pra Brasil, mas não consegui. Alguém consegue me ajudar? Segue abaixo o que tentei fazer e o máximo de info possível para o resultado ser reproduzível. Eu tentei: Sys.setlocale("LC_TIME", "Portuguese_Brazil.UTF-8")[1] ""Warning message:In Sys.setlocale("LC_TIME", "Portuguese_Brazil.UTF-8") : OS reports request to set locale to "Portuguese_Brazil.UTF-8" cannot be honored Tentei outras variantes de Brasil (como "br", "BR", mas não funcionou). Aí, procurando no StackOverflow, descobri como recuperar os localles possíveis (no linux, pelo que entendi): locales <- system("locale -a", intern = TRUE)> library(stringr)> unique(str_split_fixed(locales, "\\.", 2)[, 1]) [1] "C" "POSIX" "en_AG" "en_AU" "en_BW" "en_CA" "en_DK" "en_GB" "en_HK" "en_IE" "en_IN" "en_NG" "en_NZ" "en_PH" "en_SG" "en_US" "en_ZA" "en_ZM" "en_ZW" Como se vê pela lista acima, não parece haver opção de pt ou br. Alguma ideia de como resolver isso? Session info do R: > sessionInfo()R version 3.0.2 (2013-09-25) Platform: x86_64-pc-linux-gnu (64-bit) locale: [1] LC_CTYPE=en_US.UTF-8 LC_NUMERIC=C LC_TIME=C LC_COLLATE=C LC_MONETARY=C LC_MESSAGES=C LC_PAPER=C [8] LC_NAME=C LC_ADDRESS=C LC_TELEPHONE=C LC_MEASUREMENT=C LC_IDENTIFICATION=C attached base packages: [1] grid stats graphics grDevices utils datasets methods base other attached packages: [1] stringr_0.6.2 RGoogleAnalytics_0.1.1 httr_0.5 lubridate_1.3.2 extrafont_0.16 scales_0.2.3 [7] httpuv_1.3.0 reshape_0.8.5 gridExtra_0.9.1 plyr_1.8 rjson_0.2.13 RCurl_1.95-4.1 [13] bitops_1.0-6 ggplot2_0.9.3.1 data.table_1.9.3 RJDBC_0.2-1 rJava_0.9-4 DBI_0.2-7 loaded via a namespace (and not attached): [1] MASS_7.3-29 RColorBrewer_1.0-5 Rcpp_0.11.1 Rttf2pt1_1.3.2 colorspace_1.2-4 dichromat_2.0-0 digest_0.6.3 extrafontdb_1.0 [9] gtable_0.1.2 jsonlite_0.9.11 labeling_0.2 memoise_0.1 munsell_0.4.2 proto_0.3-10 reshape2_1.2.2 tools_3.0.2 Meu OS é Linux GNU Ubuntu 12.04 Grato, Manoel -- Manoel Galdino https://sites.google.com/site/galdinomcz/

2 3

Erro com FactoMiner e Rcmdr
by Ederson Civardi 15 Set '14

15 Set '14

Consegui solucionar o problema atualizando os pacotes e reinstalando o R. Em 15/09/2014 12:00, <r-br-request(a)listas.c3sl.ufpr.br> escreveu:

1 0

Criação de data.frame com valores médios de variaveis
by Gilenio Borges Fernandes 15 Set '14

15 Set '14

Prezados(as) Tenho um banco de dados com as seguintes variáveis (colunas): Tipo Trat Tempo Triplicata acidezKOH e acidezOLEIC Segue o CMR com a estrutura do meu banco de dados: d1=structure(list(Tipo = structure(c(1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L), .Label = c("TI1", "TI2", "TI3"), class = "factor"), Trat = structure(c(1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L), .Label = c("TR1", "TR2", "TR3"), class = "factor"), Tempo = c(0L, 0L, 0L, 3L, 3L, 3L, 6L, 6L, 6L, 0L, 0L, 0L, 3L, 3L, 3L, 6L, 6L, 6L, 0L, 0L, 0L, 3L, 3L, 3L, 6L, 6L, 6L, 0L, 0L, 0L, 3L, 3L, 3L, 6L, 6L, 6L, 0L, 0L, 0L, 3L, 3L, 3L, 6L, 6L, 6L, 0L, 0L, 0L, 3L, 3L, 3L, 6L, 6L, 6L, 0L, 0L, 0L, 3L, 3L, 3L, 6L, 6L, 6L, 0L, 0L, 0L, 3L, 3L, 3L, 6L, 6L, 6L, 0L, 0L, 0L, 3L, 3L, 3L, 6L, 6L, 6L), Triplicata = c(1L, 2L, 3L, 1L, 2L, 3L, 1L, 2L, 3L, 1L, 2L, 3L, 1L, 2L, 3L, 1L, 2L, 3L, 1L, 2L, 3L, 1L, 2L, 3L, 1L, 2L, 3L, 1L, 2L, 3L, 1L, 2L, 3L, 1L, 2L, 3L, 1L, 2L, 3L, 1L, 2L, 3L, 1L, 2L, 3L, 1L, 2L, 3L, 1L, 2L, 3L, 1L, 2L, 3L, 1L, 2L, 3L, 1L, 2L, 3L, 1L, 2L, 3L, 1L, 2L, 3L, 1L, 2L, 3L, 1L, 2L, 3L, 1L, 2L, 3L, 1L, 2L, 3L, 1L, 2L, 3L ), acidezKOH = c(2.29, 2.29, 2.38, 6.16, 6.16, 5.92, 4.82, 5.04, 5.04, 3.37, 3.54, 3.69, 5.92, 5.92, 5.71, 5.8, 5.75, 5.72, 3.16, 3.02, 3.22, 5.33, 5.33, 5.3, 5.66, 5.75, 5.7, 0.23, 0.23, 0.23, 0.3, 0.3, 0.3, 0.34, 0.31, 0.31, 0.26, 0.29, 0.26, 0.33, 0.36, 0.33, 0.43, 0.4, 0.43, 0.26, 0.26, 0.29, 0.3, 0.33, 0.3, 0.34, 0.31, 0.31, 0.09, 0.09, 0.09, 0.33, 0.33, 0.33, 0.37, 0.37, 0.34, 0.14, 0.14, 0.14, 0.3, 0.3, 0.3, 0.31, 0.34, 0.34, 0.57, 0.57, 0.57, 0.33, 0.33, 0.33, 0.8, 0.77, 0.8), acidezOLEIC = c(1.15, 1.15, 1.2, 3.1, 3.1, 2.97, 2.42, 2.54, 2.54, 1.69, 1.78, 1.85, 2.98, 2.98, 2.87, 2.92, 2.89, 2.87, 1.59, 1.52, 1.62, 2.68, 2.68, 2.66, 2.85, 2.89, 2.86, 0.11, 0.11, 0.11, 0.15, 0.15, 0.15, 0.17, 0.16, 0.16, 0.13, 0.14, 0.13, 0.17, 0.18, 0.17, 0.22, 0.2, 0.21, 0.13, 0.13, 0.14, 0.15, 0.17, 0.15, 0.17, 0.16, 0.16, 0.04, 0.04, 0.04, 0.17, 0.17, 0.17, 0.19, 0.19, 0.17, 0.07, 0.07, 0.07, 0.15, 0.15, 0.15, 0.17, 0.16, 0.16, 0.29, 0.29, 0.29, 0.17, 0.17, 0.17, 0.4, 0.39, 0.4)), .Names = c("Tipo", "Trat", "Tempo", "Triplicata", "acidezKOH", "acidezOLEIC"), class = "data.frame", row.names = c(NA, -81L)) Preciso criar um banco de dados (data.frame) com as colunas: Tipo Trat Tempo e mais as colunas acidezKOH e acidezOLEIC com a média das três Triplicatas. Tentei fazer com a função tapply, mas, não resolve porque não cria uma estrutura de data.frame. Gostaria de ter auxilio da lista para este problema. Desde já agradeço. Cordialmente, -- Gilenio Borges Fernandes Professor Associado IV (Aposentado) Universidade Federal da Bahia Instituto de Matemática Departamento de Estatística Av. Adhemar de Barros, s/n – Ondina. 40.170-110 - Salvador - BA, Brasil Tel.: (071)3283-6340/6341/6337 Fax: (071)3283-6336 URL: http://lattes.cnpq.br/6764860618464860

5 4

Erros com FactoMiner e Rcmdr
by Ederson Civardi 14 Set '14

14 Set '14

Prezados, consegui resolver os problemas reistalando o R e pacotes. Em 14/09/2014 12: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. Erros com FactoMineR e Rcmdr (Ederson Civardi) > 2. Erro com FactoMineR e Rcmdr (Ederson Civardi) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Sat, 13 Sep 2014 15:34:17 -0300 > From: Ederson Civardi <civardi.eac(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: [R-br] Erros com FactoMineR e Rcmdr > Message-ID: > < > CACA7RZec6DoHgaMTAhU8qQo8+HGBS0_DtcWgsOzNf2WHgeLK6w(a)mail.gmail.com> > Content-Type: text/plain; charset="utf-8" > > Prezados, > > Estou encontrando dificuldades de abrir o Rcmdr e o FactoMineR, onde os > utilizo para análise multivariada. > Meu sistema operacional é o Win 8.1 e > R version 3.0.3 (2014-03-06) -- "Warm Puppy" > Copyright (C) 2014 The R Foundation for Statistical Computing > Platform: x86_64-w64-mingw32/x64 (64-bit) > > O erro que aparece é o seguinte: > > The downloaded binary packages are in > C:\Users\Ederson\AppData\Local\Temp\Rtmpo7YjZQ\downloaded_packages > > library(FactoMineR) > Error in get(Info[i, 1], envir = env) : lazy-load database 'P' is corrupt > Além disso: Mensagens de aviso perdidas: > In get(Info[i, 1], envir = env) : internal error -3 in R_decompress1 > Erro: package or namespace load failed for ?FactoMineR? > > library(Rcmdr) > Carregando pacotes exigidos: splines > Carregando pacotes exigidos: RcmdrMisc > Carregando pacotes exigidos: car > Error in get(Info[i, 1], envir = env) : lazy-load database 'P' is corrupt > Além disso: Mensagens de aviso perdidas: > In get(Info[i, 1], envir = env) : internal error -3 in R_decompress1 > Erro: package ?car? could not be loaded > > Alguém poderia me ajudar ou dizer o que está acontecendo. > Desde já agradeço. > > Ederson > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20140913/813d0643/atta… > > > > ------------------------------ > > Message: 2 > Date: Sat, 13 Sep 2014 15:42:47 -0300 > From: Ederson Civardi <civardi.eac(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: [R-br] Erro com FactoMineR e Rcmdr > Message-ID: > <CACA7RZeqzSL-= > y5dK4J57wLgPRW+h7iFX2Bk37z6VKeGuVm2ug(a)mail.gmail.com> > Content-Type: text/plain; charset="utf-8" > > Já com a versão mais atual do R: > R version 3.1.1 (2014-07-10) -- "Sock it to Me" > Copyright (C) 2014 The R Foundation for Statistical Computing > Platform: x86_64-w64-mingw32/x64 (64-bit) > > Acontece o mesmo problema: > > > library(FactoMineR) > Error in get(Info[i, 1], envir = env) : > não foi possível abrir o arquivo > 'C:/Users/Ederson/Documents/R/win-library/3.1/leaps/R/leaps.rdb': No such > file or directory > Erro: package or namespace load failed for ?FactoMineR? > > library(Rcmdr) > Carregando pacotes exigidos: splines > Carregando pacotes exigidos: RcmdrMisc > Carregando pacotes exigidos: car > Carregando pacotes exigidos: sandwich > Error in get(Info[i, 1], envir = env) : > não foi possível abrir o arquivo > 'C:/Users/Ederson/Documents/R/win-library/3.1/abind/R/abind.rdb': No such > file or directory > Erro: package ?RcmdrMisc? could not be loaded > > Como faço para solucionar esse problema? Alguém pode me ajudar? > > Desde já agradeço. > Ederson > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20140913/823ce1aa/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 45, assunto 21 > ***************************************** >

1 0

Erro com FactoMineR e Rcmdr
by Ederson Civardi 13 Set '14

13 Set '14

Já com a versão mais atual do R: R version 3.1.1 (2014-07-10) -- "Sock it to Me" Copyright (C) 2014 The R Foundation for Statistical Computing Platform: x86_64-w64-mingw32/x64 (64-bit) Acontece o mesmo problema: > library(FactoMineR) Error in get(Info[i, 1], envir = env) : não foi possível abrir o arquivo 'C:/Users/Ederson/Documents/R/win-library/3.1/leaps/R/leaps.rdb': No such file or directory Erro: package or namespace load failed for ‘FactoMineR’ > library(Rcmdr) Carregando pacotes exigidos: splines Carregando pacotes exigidos: RcmdrMisc Carregando pacotes exigidos: car Carregando pacotes exigidos: sandwich Error in get(Info[i, 1], envir = env) : não foi possível abrir o arquivo 'C:/Users/Ederson/Documents/R/win-library/3.1/abind/R/abind.rdb': No such file or directory Erro: package ‘RcmdrMisc’ could not be loaded Como faço para solucionar esse problema? Alguém pode me ajudar? Desde já agradeço. Ederson

1 0

Erros com FactoMineR e Rcmdr
by Ederson Civardi 13 Set '14

13 Set '14

Prezados, Estou encontrando dificuldades de abrir o Rcmdr e o FactoMineR, onde os utilizo para análise multivariada. Meu sistema operacional é o Win 8.1 e R version 3.0.3 (2014-03-06) -- "Warm Puppy" Copyright (C) 2014 The R Foundation for Statistical Computing Platform: x86_64-w64-mingw32/x64 (64-bit) O erro que aparece é o seguinte: The downloaded binary packages are in C:\Users\Ederson\AppData\Local\Temp\Rtmpo7YjZQ\downloaded_packages > library(FactoMineR) Error in get(Info[i, 1], envir = env) : lazy-load database 'P' is corrupt Além disso: Mensagens de aviso perdidas: In get(Info[i, 1], envir = env) : internal error -3 in R_decompress1 Erro: package or namespace load failed for ‘FactoMineR’ > library(Rcmdr) Carregando pacotes exigidos: splines Carregando pacotes exigidos: RcmdrMisc Carregando pacotes exigidos: car Error in get(Info[i, 1], envir = env) : lazy-load database 'P' is corrupt Além disso: Mensagens de aviso perdidas: In get(Info[i, 1], envir = env) : internal error -3 in R_decompress1 Erro: package ‘car’ could not be loaded Alguém poderia me ajudar ou dizer o que está acontecendo. Desde já agradeço. Ederson

1 0

Pacote Animation, uso de memória
by David Feitosa 12 Set '14

12 Set '14

Caríssimos, boa noite. Segue a resposta do Robert Iquiapaza e fui usar o "animation" para gerar uma animação. O resultado ficou interessante com o saveHTML(...), pois ele gera uma página com um player e a animação etc. Entretanto, quando levo o código gerado para outra máquina, ocorre lentidão ou, até mesmo, travamento. Analisei os recursos da máquina e vi que o código que é gerado, carrega todas as imagens para a memória de forma incremental. Ao executar na minha máquina em casa com 8GB de RAM, fica tudo OK, mas, em outra, que uso no laboratório e que possui apenas 2GB "livre" ocorre este problema... Olhei o código fonte da página para tentar uma solução de contorno, mas não obtive uma "luz". Alguém já passou por situação similar e sabe de algum parâmetro, ou ajuste posterior, que permitisse carregar as imagens novas e descartar as antigas de forma a otimizar o uso da memória? Atenciosamente, David Feitosa (\_(\ (=°;°) (("")("")

1 0

Veracidade das Loterias
by andrebvs＠bol.com.br 12 Set '14

12 Set '14

6 5

Seleção de modelos
by Cassiano 12 Set '14

12 Set '14

Prezados, Preciso selecionar um modelo com 5 variáveis preditoras e uma resposta que não atende a pressupostos de normalidade. Estou usando poisson. Qual o melhor (ou mais curto) caminho para fazer a seleção de modelos no R. MASS e stepAIC não dá certo. Obrigado Cassiano -- ======================================= Cassiano S. Rosa Professor Visitante Universidade do Estado do Mato Grosso - UNEMAT Programa de Pós-graduação em Ecologia e Conservação Br 158, Km 655 - Caixa Postal 08 Antiga Sede da FAB UNEMAT – *Campus* Universitário de Nova Xavantina Nova Xavantina - MT CEP: 78690-000 Telefone: (66)3438-1224 ou 3438-2670 - FAX: (66) 3438-1224 Curriculo lattes: http://lattes.cnpq.br/3030825329622014 Scholar Google: http://scholar.google.com/citations?user=JmELzOMAAAAJ =======================================

3 2

Duplicar coluna
by Giselle Davi 12 Set '14

12 Set '14

Prezados, Poderiam me dizer como eu conseguira duplicar 20 colunas . E além disso, posicioná-las ao lado da coluna de origem ? Desde já agradeço. Att, Giselle

2 3

Re: [R-br] Veracidade das Loterias
by andrebvs＠bol.com.br 12 Set '14

12 Set '14

4 3

Mudanca da classe de dados e modificacao de valor desses
by Leonardo Monteiro 12 Set '14

12 Set '14

Prezados, segue o codigo para o exemplo: setwd('C:/Users/Leonardo/Desktop/') ##Baixando os dados da internet setInternet2(FALSE) #desativa setInternet2(NA) #mostra o status base<- ' http://power.larc.nasa.gov/cgi-bin/cgiwrap/solar/agro.cgi?email=agroclim%40… ' url<- paste0(base,sep=-22.7, '&lon=', sep = -48, '&ms=1&ds=1&ys=', sep = 1983, '&me=12&de=31&ye=', sep = 1986, '&p=toa_dwn&p=swv_dwn&p=lwv_dwn&p=T2M&p=T2MN&p=T2MX&p=RH2M&p=DFP2M&p=RAIN&p=WS10M&submit=Submit') download.file(url, 'teste.txt', method = "auto", quiet = FALSE, mode = "w", cacheOK = TRUE) ## dados<- read.table('teste.txt', sep= '', dec = '.', skip=24); head(dados) dados<- dados[,-c(3,5, 10)]; head(dados) colnames(dados)<- c('ano', 'dj', 'Qg', 'tmed', 'tmin', 'tmax', 'ur', 'chuva', 'u10'); head(dados) dados[180:190,] # dados onde comeca a haver dados de Qg (os dados estao corretos) dados$Qg[dados$Qg == '-']<- NA; dados[180:190,] ##apresentando os dados de Qg ja com NA class(dados$Qg) # dados como factor dados$Qg<- as.numeric(dados$Qg); class(dados$Qg) ##Ok, eles estao como 'numeric' dados[180:190,] #na transformacao, os valores de Qg mudaram ##### Alguem poderia me ajudar a inserir os valores NA onde estao apresentados '-', entretanto sem que modifique o valor original dos dados (os que nao sao NA)?? obrigado, Leonardo

2 1

Re: [R-br] Veracidade das Loterias
by andrebvs＠bol.com.br 11 Set '14

11 Set '14

1 0

erro em write.table
by Leonardo Monteiro 11 Set '14

11 Set '14

Pessoal, to tentando gerar um arquivo com: getwd() 'C:/Users/Leonardo/dados' a<- 1:10 write.table(a, 'a.txt', dir = 'C:/Users/Leonardo/Desktop') entretanto ta dando erro.. alguem sugere algo para salvar o arquivo em outro dir? obrigado Leonardo.

4 3

Re: [R-br] erro em write.table
by Robert Iquiapaza 11 Set '14

11 Set '14

Use write.table(a, ''C:/Users/Leonardo/Desktop/a.txt') # se Desktop existir vai funcionar me parece que não existe o parâmetro ‘dir’ nessa função. Sds From: David Feitosa Sent: Thursday, September 11, 2014 5:44 PM To: r-br(a)listas.c3sl.ufpr.br Subject: Re: [R-br] erro em write.table Qual a MSG de erro? Em 11/09/2014 17:40, "Leonardo Monteiro" <monteiroleonardo6(a)gmail.com> escreveu: Pessoal, to tentando gerar um arquivo com: getwd() 'C:/Users/Leonardo/dados' a<- 1:10 write.table(a, 'a.txt', dir = 'C:/Users/Leonardo/Desktop') entretanto ta dando erro.. alguem sugere algo para salvar o arquivo em outro dir? obrigado Leonardo. _______________________________________________ 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. -------------------------------------------------------------------------------- _______________________________________________ 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

Oportunidade para cientísta de dados
by Fabrício Barth 11 Set '14

11 Set '14

Caros(as), Peço desculpas pelo "spam". Mas, acredito que pode ser relevante para algumas pessoas da lista. A VAGAS.com está com uma oportunidade para cientísta de dados. Segue a descrição da oportunidade: *Principais responsabilidades: * - Desenvolver sistemas utilizando técnicas de aprendizagem de máquina (modelos descritivos e preditivos) a serem aplicados na construção de softwares para a área de Recrutamento & Seleção, por exemplo: sistemas de recomendação de vagas, entre outros. *Perfil desejado:* - Experiência profissional ou acadêmica no desenvolvimento de modelos descritivos (por exemplo: clustering e regras de associação), preditivos (por exemplo: modelos classificação ou regressão) ou sistemas de recomendação. - Experiência profissional ou acadêmica no desenvolvimento de sistemas capazes de tratar um grande volume de dados heterogêneos (dados estruturados ou não). *Diferenciais:* - Mestrado ou Doutorado em Exatas tais como: Computação, Engenharia, Matemática, Física. - Experiência com mais de uma das seguintes linguagens de programação: R, Python, Ruby ou Java. - Já ter utilizado banco de dados não relacionais, por exemplo: Neo4J e MongoDB. Mais informações podem ser obtidas em http://site.vagas.com.br/PagVagaDirSS.asp?v=1031491 Abraços, -- Fabrício J. Barth http://fbarth.net.br http://blog.fbarth.net.br

1 0