- R-br - DINF UFPR

erro na validação cruzada com pacote geoR (Éder Comunello)
by Wenceslau Teixeira 22 Out '13

22 Out '13

Message: 1 Date: Sun, 20 Oct 2013 15:01:59 -0300 From: Éder Comunello <comunello.eder(a)gmail.com> To: r-br(a)listas.c3sl.ufpr.br Subject: Re: [R-br] erro na validação cruzada com pacote geoR Message-ID: <CABmC8g=yZ1eCQXA_hNp_iOrP9v8CqNztVLBYN4dgCtHMy7Z9cQ(a)mail.gmail.com> Content-Type: text/plain; charset="iso-8859-1" Caros, Aproveitando a discussão para aprender sobre a validação de modelos em geoestatisticos. Dado que estest modelos e tecnicas são novas para mim. A krigagem é um estimador exato se o numero de pontos estimados for o mesmo dos avaliados, não ? Se eu estimar 10.000 valores para um universo de 1000 valores medidos. Os valores estimados para os 1000 pontos avaliados serão os mesmos valores ? Saudações geoestatisticas Wenceslau Teixeira Embrapa Solos _________________________________________________ Prezado Helio, boa tarde! Pensando melhor sobre a questão do resíduo, o que coloquei no email anterior não faz muito sentido, já que a krigagem é um estimador exato e a operação retornaria praticamente o próprio valor inicial (e o resíduo tenderia ser zero/nulo). O que faria sentido seria separar o efeito da tendência ('drift') e descontar esse valor dos dados originais, tendo aí sim o resíduo a ser modelado (sem o efeito da tendência). Para o produto final teria que somar o drift novamente para retornar os valores. Acredito que não seja muito difícil de fazer no R, mas só tentando fazer pra ver. Fui dar uma conferida no seu post antigo sobre o erro na validação cruzada: < http://r-br.2285057.n4.nabble.com/R-br-erro-na-validacao-cruzada-td4659526.… > Mas o link dos dados não está mais disponível... <https://www.dropbox.com/s/tfml7zjpc7hmbc3/dados250.txt> Se importaria em dispor esses dados novamente? Atte.

3 2

tapply
by geovane barbosa 22 Out '13

22 Out '13

Olá pessoal tudo bem, estou com uma dúvida. Estou sando a função tapply da seguinte forma : tapply(pressaoarterial, raca9, mean,na.rm = TRUE) # só que nesse caso ele me retorna o valor médio da pressão por raça, eu gostaria de introduzir mais uma variável que é atividade física caracterizada por 4 variáveis. Logo, seria visto a média da pressão arterial nas raças dividas por períodos de atividade física. grato #################################### #################################### Prof. Geovane Carlos Barbosa UCL - Faculdade do Centro Leste ################################### ###################################

3 2

Duvida rcpp, como remover um elemento de um vetor no código C++
by Augusto Ribas 20 Out '13

20 Out '13

Ola pessoal. Eu estou com uma dificuldade em como fazer uma chamada recursiva usando o Rcpp. Por curiosidade, eu estava fazendo alguns códigos aqui usando o Rcpp, só por experiencia. Ai eu estava implementando o código para o insertion sort ( http://pt.wikipedia.org/wiki/Insertion_sort) Minha duvida é a seguinte: #Suponha que eu tenho um vetor qualquer: vetor<-sample(100) vetor #Iterativamente, eu fiz um código para usar com o Rcpp assim: library(Rcpp) cppFunction(" NumericVector insertionsortC(NumericVector vetor) { int n = vetor.size(); double aux; int i , j; for(i=1;i<n;i++) { aux=vetor[i]; j=i-1; while(j>=0 && vetor[j]>aux) { vetor[j+1]=vetor[j]; j=j-1; } vetor[j+1]=aux; } return vetor; } ") #Funcionou beleza. Ai eu fui tentar repetir a tentativa, recursivamente cppFunction(" NumericVector insertionsortRC(NumericVector vetor, int n) { double aux; int i; if(n>1) { insertionsortRC(vetor,n-1); aux=vetor[n-1]; i=n-1; while(vetor[i-1]>aux && i>=0 ) { vetor[i]=vetor[i-1]; i--; } vetor[i]=aux; } return vetor; } ") #E assim funcionou blz, mas eu tive que começar a mandar o argumento n, #que é o tamanho do vetor #o que eu queria é fazer algo aqui, na chamada recursiva #insertionsortRC(vetor,n-1); #sei la, se eu conseguir tirar o ultimo elemento, eu posso suprimir esse n-1 #Eu não entendo de Rcpp e C++ muito bem. #Mas tem uma forma simples de mudar o tamanho do vetor, ou arrancar o ultimo elemento dele #algo tipo vetor.pop() #pra fazer a chamada ficar tipo #vetor.pop(); #insertionsortRC(vetor); #que ai no inicio da chamada eu uso o int n = vetor.size(); igual ao formato iterativo. #Bem, espero não estar sendo muito confuso na pergunta e que a função funcione. #Mas imagino que deve ser super simples isso. Agradeço a atenção desde ja :) -- Grato Augusto C. A. Ribas Site Pessoal: http://recologia.com.br/ <http://augustoribas.heliohost.org> Github: https://github.com/Squiercg Lattes: http://lattes.cnpq.br/7355685961127056

3 5

(sem assunto)
by Fátima Lima Paula 20 Out '13

20 Out '13

Pessoal, o Fernando me ajudou a agregar com 2 variáveis. Agora preciso agregar estratificando. Ou seja, idade=c("1","3","2","2","1","2","1","3","1","3") sexo <- c("1","1","2","2","1","2","1","1","1","1") tcont <- c(100,320,24,256,134,290,18,34,15,121) read <- c("s","n","n","n","s","s","s","n","n","s") dados <- data.frame(sexo,tcont,read) dados dados$read <- as.numeric(ifelse(read == 's', 1, 0)) dados Agora quero que o resultado fique assim: sexo idade tcont read 1 1 1 2 1 2 3 1 3 4 2 1 5 2 2 6 2 3 O tcont será a soma e a coluna de read terá a quantidade de readmissões estratificada. Alguém pode me ajudar? Obrigada Fátima -- "Minha felicidade depende da qualidade dos meus pensamentos"

7 11

mensagem de erro em simulacao
by Maurício Lordêlo 19 Out '13

19 Out '13

#Caros, #Estou fazendo uso da função "bild" pertencente ao pacote com mesmo nome. #Nos dados que são usados como exemplo do pacote (locust) temos: install.packages("bild") require(bild) locust2 <- bild(move ~ time+feed, data = locust,start = NULL, aggregate = feed, dependence = "MC1R") summary(locust2) #O programa retorna os resultados sem problemas, como mostrado abaixo: Call: bild(formula = move ~ time + feed, data = locust, aggregate = feed,start = NULL, dependence = "MC1R") Number of profiles in the dataset: 24 Number of profiles used in the fit: 24 Log likelihood: -1537.923 AIC: 3085.845 Coefficients: Label Value Std. Error t value p-value (Intercept) 1 -1.048305 0.30730336 -3.411 0.000647 time 2 1.572124 0.14600334 10.768 0.000000 feed1 3 -3.161945 0.43072021 -7.341 0.000000 log.psi1 4 1.070874 0.09833972 10.890 0.000000 Random effect (omega): Value Std. Error -0.03532051 0.34195846 Message: 0 #Porém, se faço warnings() #Ele me retorna Mensagens de aviso: 1: In sqrt(prob * (1 - prob)) : NaNs produzidos 2: In sqrt(prob * (1 - prob)) : NaNs produzidos 3: In sqrt(prob * (1 - prob)) : NaNs produzidos 4: In sqrt(prob * (1 - prob)) : NaNs produzidos #Mesmo com estas mensagens consigo obter os resultados. #O problema surge quando peço para simular muitos ajustes deste modelo. #Por exemplo, n.time=4 #numero de medidas no tempo n.trat=3 #numero de tratamentos n.rep=10 #numero de repeticoes ns=15 #numero de simulacoes n.ind <- n.trat * n.rep #numero de individuos n <- n.trat * n.rep * n.time #numero total de registros id <- rep(1:n.ind, each = n.time) #identificacao do individuo time <- rep(1:n.time, times = n.ind) #tempo trat <- rep(factor(LETTERS[1:n.trat]), each = n.rep*n.time) #tratamentos dados=data.frame(id,time,trat) #estrutura do banco de dados #funcao para gerar resposta binaria z <- function() { y <- integer(n) for(i in 1:n) y[i] <- rbinom(1,1,.5) return(y) } repet <- replicate(ns, z(), simplify = FALSE) #funcao para obter as estimativas coef_bild <- function(y) { dat <- data.frame(dados, y=z()) coefs=bild(y ~ trat+time, data = dat, start=NULL,dependence="MC1R")@coefficients coefs_bild=coefs[1:4] } #lista dos coeficientes ajustados coefs_ajust <- do.call(rbind, lapply(repet, coef_bild)) #obtencao das medias dos coeficientes estimados media_coefs <- apply(coefs_ajust, 2, mean) #Com apenas ns=15 ele me dá as mensagens de aviso sem eu pedir (as mesmas do exemplo para o conjunto de dados "locust"), mas me retorna o que preciso. #Quando coloco ns=2000 (que é a quantidade que preciso) demora bastante tempo, não retorna o que preciso e emite a mensagem: Error in solve.default(Info) : system is computationally singular: reciprocal condition number = 1.69179e-16 In addition: There were 50 or more warnings (use warnings() to see the first 50) #As 50 ou mais mensagens são as mesmas que aparecem quando rodo com uma quantidade menor . Alguém tem alguma dica ou forma de conseguir estes resultados para o número de simulações que desejo? # Caso nao encontre uma solução, verifiquei que com ns=100 ele faz. Tentei fazer um "for de 1 a 20, com ns=100" das funções que defini como "z" e "coef_bild" e depois concatenar, mas não consegui. Em último caso, peço ajuda a quem saiba fazer para este caso específico. #Agradeço antecipadamente, #Maurício

1 0

Materiais de programacao gratuitos
by Benilton Carvalho 19 Out '13

19 Out '13

https://github.com/vhf/free-programming-books/blob/master/free-programming-… divirtam-se

3 2

Notação potências de dez (Gráficos)
by Eric Leandro 18 Out '13

18 Out '13

Boa tarde pessoal Tenho uma dúvida quanto a intervalo com potências de dez no R. Quando dou este comando : plot(areadrenagem,vazao, xlab="Área drenagem (m²)",ylab="Q(m³/s)",pch=16,xlim=range(0:5e+12),ylim=range(0:150000)), No eixo x aparece 1e+12..... 5e+12, como eu faria para aparecer neste formato : 2.1012 ..... 5.1012. ATT Eric Leandro Silva Pereira Mestrando em Clima e Ambiente Instituto Nacional de Pesquisas da Amazônia Programa de Grande Escala da Biosfera-Atmosfera na Amazônia Av. André Araújo, 2936 - Aleixo, CEP 69060-001, Manaus - AM.

4 14

Substituir NA em data.frame 1 por dados data.frame 2
by Alisson Lucrecio 18 Out '13

18 Out '13

Caro Colegas da r-br, Boa noite. Eu tenho 2 data.frame, uma com dados originais publicado e outra como dados gerados a partir de dados de estação climatológicas do Brasil. A primeira data.frame tem alguns NA, os quais preciso substituir pelos valores calculado, assim preciso fazer uma função que encontre o valor e substitua pelo valor na data.frame 2, entretanto essas duas data.frame não tem e numero de linha e ordem. Obrigado. Alisson Lucrécio da Costa

4 7

função 'relógio'
by Teodoro Calvo 17 Out '13

17 Out '13

Olá, mais uma vez recorro a vocês. hehe Gostaria de saber se existe alguma função "relógio", onde eu possa criar uma estrutura de repetição que seja executada a cada 'x' minutos por exemplo. Muito obrigado. Att Téo Calvo.

3 3

Duvida sobre rlm()
by Luiz Henrique da Conceicao Leal 17 Out '13

17 Out '13

Pessoal, estou com dúvidas para interpretar os resultados abaixo. O intercepto estimado (em cada caso) é uma média ponderada? Quais seriam os pesos? Desde já agradeço a ajuda. summary(rlm(y~1,method="M",scale.est="Huber",k2=1.345))$coefficients summary(rlm(y~1,method="M",weights=1/uncer^2,scale.est="Huber"))$coefficients summary(rlm(y~1,method="MM",weights=1/uncer^2))$coefficients Luiz Henrique.

1 0

Aggregate
by Fátima Lima Paula 17 Out '13

17 Out '13

Pessoal, estou querendo agregar umas variáveis, mas não sei. Imaginemos o seguinte data frame: sexo=c("1","1","2","2","1","2","1","1","1","1") tcont=c("100","320","24","256","134","290","18","34","15","121") read=c("s","n","n","n","s","s","s","n","n","s") df=data.frame(sexo,idade,read) df Quero agregar essas variáveis de forma que fique a primeira coluna sexo, a segunda coluna a soma do tcont e a terceira coluna a contagem dos read = sim. Alguém pode me ajudar, por favor? Obrigada Fátima -- "Minha felicidade depende da qualidade dos meus pensamentos"

5 5

auxílio: conjunto de dados e servidor
by Maurício Lordêlo 17 Out '13

17 Out '13

Caro usuário do R, Estou desenvolvendo meu trabalho de doutorado e pretendo obter resultados através de simulações. Consegui fazer todo o procedimento no R através do auxílio dos membros desta lista. No entanto, gostaria de também fazer análises através de dados reais. Como estou com dificuldade em obter estes dados e sei que existem muitos aqui que são de várias áreas diferentes, resolvi solicitar a ajuda nesta lista. O trabalho consiste em um delineamento experimental inteiramente casualizado com resposta binária medida no tempo onde se tem interesse em avaliar o efeito de três ou mais tratamentos. Em último caso, mesmo que a resposta não seja binária, usarei alguma recurso que possa transformá-la. Caso tenha ou conheça alguém que possa ter este conjunto de dados (e lógico, que tenha interesse em disponibilizá-lo) ou até mesmo uma fonte onde possa encontrar ficarei muito agradecido. Aproveitando a oportunidade, como tenho muitos cenários para simular, alguém conhece alguma instituição que disponibilize o servidor para que estas simulações sejam rodadas? Desde já agradeço a quem possa contribuir. Abraço, Maurício Lordêlo

3 3

Re: [R-br] Distribuição de Bernoulli
by Ivan Bezerra Allaman 14 Out '13

14 Out '13

Se foi lhe proposto trabalhar com a função dbern, então se trata de outro software e não o R. Caso exista esta função em outro software, ela é extramente restrita, uma vez que um variável aleatória binomial são "n" realizações de uma variável bernoulli. Outra coisa, sua regressão não é multivariada pelo amor de Deus!!!!!!!! O que caracteriza uma análise uni ou multivariada não são as variáveis preditoras!! Caso ainda queira que alguém lhe ajude, nos informe qual a mensagem de erro ou ainda melhor, nos envie um CMR. Cordialmente, (S,f,P) Allaman

3 3

Correlação policórica
by Felipe Buchbinder 14 Out '13

14 Out '13

Professores, bom dia, Existe alguma função para calcular uma matriz de correlação policórica a partir de uma matriz de respostas a um questionário de escala Likert? A função polychor, do pacote polycor, me retorna apenas um número... Obrigado, Felipe

4 9

Remover NULL de um output da função lapply
by Tiago Souza Marçal 14 Out '13

14 Out '13

Bom dia pessoal, gostaria de saber como remover os NULL do CMR abaixo. x=c("b","a","c") u=lapply(1:3,function(i){ if(x[i]=="a"){ 1} } ) att. Tiago. ################################################################# Tiago de Souza Marçal - Graduando em Agronomia pelo CCA-UFES Bolsista de Iniciação Científica da área de Genética e Melhoramento de Plantas #################################################################

2 1

[OFF-Topic] Dados climáticos de Paris durante maratona de 20 km
by walmes . 12 Out '13

12 Out '13

Saudações à todos, Nos meus cursos para motivar os alunos estou considerando dados reais e maiores que os dados disponíveis nos livros. Dessa vez pedi para que trabalhassem com os dados da maratona de paris cujas planilhas dos últimos anos estão disponíveis em http://www.20kmparis.com/web/resultats_precedents.asp Muito conhecimento se adquire da análise desses dados. É muito gostoso trabalhar esses dados em sala de aula. São infinitas possibilidades de gráficos. Para ver alguns gráficos baixe descompacte o arquivo e rode a aula10.R http://dl.dropboxusercontent.com/u/48140237/ce083.zip No entanto eu gostaria de verificar o efeito do clima sobre os resultados das corridas (no tempo médio, no tempo recorde, nos homens e mulheres, etc). Para isso eu preciso de séries climáticas (ou valores pontuais) de Paris no dia da prova. As datas das provas estão nos panfletos de divulgação da mesma em http://www.20kmparis.com/web/affiches.asp Alguém saberia indicar um site de onde eu pudesse pegar dados climáticos (temperatura e umidade relativa) nessas datas? A idéia e fazer merge de bases de dados, e inclusive testes de hipótese simples envolvendo às variáveis. Muitos maratonistas participam vários anos e tem como ver o desempenho individual. Fiz uma busca no google mas nada que me desse séries históricas. Talvez seja do conhecimento dos que trabalham com séries temporais a origem de dados assim. Quem souber fontes de dados assim, por gentileza me indique. A maratona acontece novamente agora dia 13 de outubro. Quero que os alunos façam uma previsão do tempo recorde por meio da temperatura prevista para Paris ou coisas do tipo. É uma forma de estimulá-los a conhecer os dados. Grato à todos. Walmes. ========================================================================== Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 skype: walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================

6 8

Script de Bernouilli
by Gisele Correia 12 Out '13

12 Out '13

Olá!! Alguém, por gentileza, poderia me dizer por que o modelo está dando erro? Segue: # Gerando dados tabela<-read.csv("morfo.csv",sep=";",head=T) tabela$SX<-as.factor(tabela$SX) #criar a lista de dados para o jags lddados <- list(y=tabela$SX, N=length(tabela$SX), CT=tabela$CT, CF=tabela$CF, CP=tabela$CP, CC=tabela$CC, FO=tabela$FO, AC=tabela$AC, PD=tabela$PD, CD=tabela$CD, BD=tabela$BD, PP=tabela$PP, CPT=tabela$CPT, BP=tabela$BP, PV=tabela$PV, BV=tabela$BV, LC=tabela$LC, CV=tabela$CV, PA=tabela$PA, BA=tabela$BA, CA=tabela$CA, AP=tabela$AP, LP=tabela$LP, CM=tabela$CM, DOH=tabela$DOH, DOV=tabela$DOV) #criar modelo sink('morfo.txt') cat(' model{ for(i in 1:N){ # loop sobre as observações y[i] ~ dbern(p[i]) # saída binária logit(p[i]) <- beta0 # regressão multivariada + beta1*CT[i] + beta2*CF[i] + beta3*CP[i] + beta4*CC[i] + beta5*FO[i] + beta6*AC[i] + beta7*PD[i] + beta8*CD[i] + beta9*BD[i] + beta10*PP[i] + beta11*CPT[i] + beta12*BP[i] + beta13*PV[i] + beta14*BV[i] + beta15*LC[i] + beta16*CV[i] + beta17*PA[i] + beta18*BA[i] + beta19*CA[i] + beta20*AP[i] + beta21*LP[i] + beta22*CM[i] + beta23*DOH[i] + beta24*DOV[i] } # Prioris beta0 ~ dnorm(0,1.0E-6) beta1 ~ dnorm(0,1.0E-6) beta2 ~ dnorm(0,1.0E-6) beta3 ~ dnorm(0,1.0E-6) beta4 ~ dnorm(0,1.0E-6) beta5 ~ dnorm(0,1.0E-6) beta6 ~ dnorm(0,1.0E-6) beta7 ~ dnorm(0,1.0E-6) beta8 ~ dnorm(0,1.0E-6) beta9 ~ dnorm(0,1.0E-6) beta10 ~ dnorm(0,1.0E-6) beta11 ~ dnorm(0,1.0E-6) beta12 ~ dnorm(0,1.0E-6) beta13 ~ dnorm(0,1.0E-6) beta14 ~ dnorm(0,1.0E-6) beta15 ~ dnorm(0,1.0E-6) beta16 ~ dnorm(0,1.0E-6) beta17 ~ dnorm(0,1.0E-6) beta18 ~ dnorm(0,1.0E-6) beta19 ~ dnorm(0,1.0E-6) beta20 ~ dnorm(0,1.0E-6) beta21 ~ dnorm(0,1.0E-6) beta22 ~ dnorm(0,1.0E-6) beta23 ~ dnorm(0,1.0E-6) beta24 ~ dnorm(0,1.0E-6) } ') sink() #carregar a biblioteca jags library(rjags) #definir os parâmetros para as posteriores params <- c("beta0", "beta1", "beta2", "beta3", "beta4", "beta5", "beta6", "beta7", "beta8", "beta9", "beta10", "beta11", "beta12", "beta13", "beta14", "beta15", "beta16", "beta17", "beta18", "beta19", "beta20", "beta21", "beta22", "beta23", "beta24") #inicializar cadeia inicio <- list(beta0=0, beta1=0, beta2=0, beta3=0, beta4=0, beta5=0, beta6=0, beta7=0, beta8=0, beta9=0, beta10=0, beta11=0, beta12=0, beta13=0, beta14=0, beta15=0, beta16=0, beta17=0, beta18=0, beta19=0, beta20=0, beta21=0, beta22=0, beta23=0, beta24=0) #chamar o modelo no jags morfo.m <- jags.model("morfo.txt",data=lddados,inicio,n.chains=1,n.adapt=1000) Abraço, Gisele

1 0

Distribuição de Bernoulli
by Gisele Correia 12 Out '13

12 Out '13

Olá!! Alguém teria um exemplo de distribuição de Bernoulli no R? Até agora não encontrei nenhum trabalho/problema utilizando a função dbern(p). Pretendo resolver o problema utilizando cadeias de Markov. Desde já agradeço, Gisele

2 3

Função-escalonar
by Maria de Fatima Ferreira Almeida 11 Out '13

11 Out '13

Gostaria de saber se existe alguma função do R que faz escalonamento de matrizes. Alguém poderia me ajudar?

4 6

Fw: [STAT-MATH] Oportunidade Profissional - Consultor Estatístico
by pjobenil＠gmail.com 11 Out '13

11 Out '13

Repassando aos colegas. Enviado do Email do Windows De: nati.0401 Enviado: ‎quinta-feira‎, ‎10‎ de ‎outubro‎ de ‎2013 ‎15‎:‎29 Para: STAT-MATH(a)yahoogrupos.com.br Consultor Estatístico – Foco no software R Será responsável por atuar com pesquisas baseadas em dados estatísticos, mapeamento de concorrência, mensuração de análises de mercado, entre outras atividades. Profissional com experiência e profundos conhecimentos no software estatístico R. Ensino Superior completo. Local: zona Oeste de São Paulo Remuneração: a combinar Interessados, encaminhar CV atualizado para nmarques(a)fitrh.com.br ou então, retornar no tel.: (11) 3149-9128. __._,_.___ | através de email | Responder através da web | Adicionar um novo tópico Mensagens neste tópico (1) Atividade nos últimos dias: Novos usuários 13 Visite seu Grupo Yahoo! Grupos Trocar para: SÃ³ Texto, Resenha DiÃ¡ria • Sair do grupo • Termos de uso . __,_._,___

1 0

eyefit
by Hélio Gallo Rocha 08 Out '13

08 Out '13

Caros da Lista Instalei o R numa máquina nova e rodando o pacote geoR dáerro na função eyefit: Erro em eyefit(vario_esc) : não foi possível encontrar a função "tclVar" Procurei esta função para instalar, achei tcltk2, mas não é o caso. Em que outro pacote pode ter esta função? Grato a todos -- Hélio Gallo Rocha IFSULDEMINAS - Câmpus Muzambinho

1 1

Regras de Associação Indiretas
by Flavio Barros 08 Out '13

08 Out '13

Alguém sabe se existe alguma implementação em R para o seguinte algortimo: http://link.springer.com/chapter/10.1007%2F3-540-32392-9_9#page-1 qualquer implementação para regras de associação indiretas serviria. Eu sei que existe o pacote *arules* mas ele não encontra regras de associação indiretas. Flavio Barros [image: Facebook]<http://s.wisestamp.com/links?url=http%3A%2F%2Fwww.facebook.com%2Fflavio.bar…> [image: LinkedIn]<http://s.wisestamp.com/links?url=http%3A%2F%2Fwww.linkedin.com%2Fprofile%2F…> Contact me: [image: Google Talk] flaviomargarito(a)gmail.com “"Only as high as I reach can I grow, only as far as I seek can I go, only as deep as I look can I see, only as much as I dream can I be." - Karen Ravn<http://s.wisestamp.com/links?url=http%3A%2F%2Fwww.quotesdaddy.com%2Fquote%2…> ” Get this email app! <http://s.wisestamp.com/links?url=http%3A%2F%2Fwww.wisestamp.com%2Fapps%2Fqu…> Create your free signature: <http://s.wisestamp.com/links?url=http%3A%2F%2Fr1.wisestamp.com%2Fr%2Flandin…> CLICK HERE<http://s.wisestamp.com/links?url=http%3A%2F%2Fr1.wisestamp.com%2Fr%2Flandin…>

3 2

Re: [R-br] Usar R p/ importar informações da web
by Eder Comunello 05 Out '13

05 Out '13

Teodoro, bom dia! Procure pelos pacotes httr, XML e RCurl. Se o site não exigir login/senha, será bem mais simples e você poderá operar com os dois primeiros. Caso seja necessário passar o login/senha vale o tópico indicado anteriormente (RCurl). < http://r-br.2285057.n4.nabble.com/R-br-r-baixando-dados-inmet-td4660459.html > Abaixo tem uns exemplos com os dois primeiros pacotes que peguei na net: ### ======================================= library(XML) library(httr) theurl <- "http://www.rout.gr/index.php?name=Rout&file=results&year=2011" doc <- htmlParse(GET(theurl, user_agent("Mozilla"))) results <- xpathSApply(doc, "//*/table[@id='table_results_r_1']") results <- readHTMLTable(results[[1]]) #rm(doc) htmlTreeParse(txt, asText = TRUE) ### ======================================= library(XML) res = htmlTreeParse("http://www.omegahat.org/RSXML/", useInternalNodes=TRUE) xpathApply(res, "//h1", xmlValue) xpathApply(res, "//pre", xmlValue) #[[1]] #[1] "An XML package for the S language" ### ======================================= -- Éder Comunello <e <comunello.eder(a)gmail.com>comunel(a)gmail.com> Dourados, MS - [22 16.5'S, 54 49'W] Em 2 de outubro de 2013 19:43, Teodoro Calvo <teocalvo2(a)gmail.com> escreveu: > Olá, boa noite. > > Como posso extrair uma parte de um texto de determinado site, utilizando o > R ? > É possível ? Existe algum material ? > > Obrigado desde já. > > Att Téo Calvo. > ______________________________**_________________ > R-br mailing list > R-br(a)listas.c3sl.ufpr.br > https://listas.inf.ufpr.br/**cgi-bin/mailman/listinfo/r-br<https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br> > Leia o guia de postagem (http://www.leg.ufpr.br/r-br-**guia<http://www.leg.ufpr.br/r-br-guia>) > e forneça código mínimo reproduzível. >

3 2

Posição das Barras - ggplot2
by Sérgio Henrique almeida da silva ju 05 Out '13

05 Out '13

Prezados Estou fazendo um gráfico usando o ggplot e não estou conseguindo diminuir os espaços entre as barras. Só consigo fazer isso aumentando o tamanho das barras, mas gostaria que elas ficassem finas e não tão distantes. Já usei o position_dodge(width=0.5), mas não me ajudou. Alguém sabe como eu posso fazer isso? Meu CMR: seguro <- matrix(0,4,3) seguro <- as.data.frame(seguro) seguro[,1] <- c("2010","2011","2012","2013") seguro[,2] <-c(89,86,87,88) seguro[,3] <-c("89%","86%","87%","88%") names(seguro)[c(1)]<-c("Ano") ggplot(seguro, aes(Ano, V2, fill=Ano))+geom_bar(width=0.3,stat="identity", position="identity", aes(fill=Ano)) + scale_y_continuous(labels=percent) + geom_text(data=seguro,aes(x=Ano,label=V3),vjust=0) Abraços -- Sérgio Henrique Almeida da Silva Junior Doutorando em Epidemiologia em Saúde Pública Escola Nacional de Saúde Pública Sérgio Arouca - ENSP/FIOCRUZ http://lattes.cnpq.br/1611345552843383 Tel: (21) 68463637 http://www.linkedin.com/profile/view?id=250437145&trk=tab_pro

1 0

Re: [R-br] correção de tendência direcional
by Éder Comunello 05 Out '13

05 Out '13

Aline, bom dia! Conseguiu resolver o problema com a modelagem dos dados geoestatísticos? Caso ainda esteja precisando de ajuda, haveria algum dado que você poderia disponibilizar (para reproduzir o problema)? -- Éder Comunello <c <comunello.eder(a)gmail.com>omunello.eder(a)gmail.com> Dourados, MS - [22 16.5'S, 54 49'W]

1 0

tentando compilar o R no Windows (com Rtools)
by Cleber N.Borges 04 Out '13

04 Out '13

olá a todos, estou tentando compilar o R no Win 7 e compila uma parte pequena pq qdo tenta entrar para compilar as bibliotecas da base o script não consegue criar diretorios temporarios. tomei o cuidado de instalar o Rtools como administrador e também chamei o prompt (terminal) como administrador se alguém tiver algo pra ajudar já agradeço antecipadamente. segue algum detalhe abaixo. obrigado cleber ############################################# Arquivo de LOG http://klebyn.ploud.com/arquivo_log/log C:\Rsrc> C:\Rsrc> C:\Rsrc>tar -xf R-3.0.2.tar.gz C:\Rsrc>where basename cat cmp comm cp cut date diff du echo expr gzip ls makeinfo C:\Rtools\bin\basename.exe C:\Rtools\bin\cat.exe C:\Rtools\bin\cmp.exe C:\Rtools\bin\comm.exe C:\Rtools\bin\cp.exe C:\Rtools\bin\cut.exe C:\Rtools\bin\date.exe C:\Rtools\bin\diff.exe C:\Rtools\bin\du.exe C:\Rtools\bin\echo.exe C:\Rtools\bin\expr.exe C:\Rtools\bin\gzip.exe C:\Rtools\bin\ls.exe C:\Rtools\bin\makeinfo.exe C:\Program Files (x86)\MiKTeX 2.9\miktex\bin\makeinfo.exe C:\Rsrc>where mkdir mv rm rsync sed sort texindex touch uniq C:\Rtools\bin\mkdir.exe C:\Rtools\bin\mv.exe C:\Rtools\bin\rm.exe C:\Rtools\bin\rsync.exe C:\Rtools\bin\sed.exe C:\Rtools\bin\sort.exe C:\Windows\System32\sort.exe C:\Rtools\bin\texindex.exe C:\Program Files (x86)\MiKTeX 2.9\miktex\bin\texindex.exe C:\Rtools\bin\touch.exe C:\Rtools\bin\uniq.exe C:\Rsrc>sort --version sort (GNU coreutils) 8.15 Packaged by Cygwin (8.15-1) Copyright (C) 2012 Free Software Foundation, Inc. License GPLv3+: GNU GPL version 3 or later <http://gnu.org/licenses/gpl.html>. This is free software: you are free to change and redistribute it. There is NO WARRANTY, to the extent permitted by law. Written by Mike Haertel and Paul Eggert. C:\Rsrc>cd R-3.0.2\src\gnuwin32 C:\Rsrc\R-3.0.2\src\gnuwin32>make all recommended > compilaR.log ######################################################################### ######################################################################### http://klebyn.ploud.com/arquivo_log/log ######################################################################### ######################################################################### connections.c: In function 'do_readbin': connections.c:3759:8: warning: dereferencing type-punned pointer will break strict-aliasing rules [-Wstrict-aliasing] connections.c:3761:8: warning: dereferencing type-punned pointer will break strict-aliasing rules [-Wstrict-aliasing] connections.c:3769:4: warning: dereferencing type-punned pointer will break strict-aliasing rules [-Wstrict-aliasing] connections.c:3784:4: warning: dereferencing type-punned pointer will break strict-aliasing rules [-Wstrict-aliasing] connections.c:3788:4: warning: dereferencing type-punned pointer will break strict-aliasing rules [-Wstrict-aliasing] pcre_exec.c: In function 'pcre_exec': pcre_exec.c:7190:20: warning: 'match_partial' may be used uninitialized in this function [-Wuninitialized] localtime.c: In function 'timesub.isra.2': localtime.c:1407:5: warning: assuming signed overflow does not occur when assuming that (X + c) < X is always false [-Wstrict-overflow] localtime.c:1411:8: warning: assuming signed overflow does not occur when assuming that (X - c) > X is always false [-Wstrict-overflow] localtime.c: In function 'time2sub.constprop.10': localtime.c:1566:8: warning: assuming signed overflow does not occur when assuming that (X + c) < X is always false [-Wstrict-overflow] localtime.c:1581:5: warning: assuming signed overflow does not occur when assuming that (X + c) < X is always false [-Wstrict-overflow] localtime.c:1593:9: warning: assuming signed overflow does not occur when assuming that (X + c) < X is always false [-Wstrict-overflow] localtime.c:1599:8: warning: assuming signed overflow does not occur when assuming that (X - c) > X is always false [-Wstrict-overflow] localtime.c:1619:5: warning: assuming signed overflow does not occur when assuming that (X - c) > X is always false [-Wstrict-overflow] cannot create /tmp/R4428: directory nonexistent mv: cannot stat `/tmp/R4428': No such file or directory make[3]: *** [mkR1] Error 1 make[2]: *** [all] Error 2 make[1]: *** [R] Error 1 make: *** [all] Error 2 C:\Rsrc\R-3.0.2\src\gnuwin32>

1 0

[OT] III CONBRATRI - Congresso Brasileiro de Teoria de Resposta ao Item - Dez2013 - Belem/PA
by Leandro Marino 04 Out '13

04 Out '13

Repassando: Informamos a todos da realização do 3º Congresso Brasileiro de Teoria da Resposta ao Item (CONBRATRI), a ocorrer em Belém do Pará, de 04 a 06 de dezembro de 2013, no Hotel Hilton. O evento é uma realização da Associação Brasileira de Avaliação Educacional (ABAVE) com apoio da Associação Brasileira de Estatística (ABE), e visa reunir pesquisadores e estudantes da área. O foco principal será Testes Adaptativos Computadorizados/Informatizados (Computerized Adaptive Testing - CAT), apresentando o que há de mais recente nesta área. Considerada uma das áreas da Estatística em maior evidência na atualidade, a TRI vem ganhando cada vez mais espaço em diversas áreas, e tem sido apontada como a tendência das avaliações em larga escala, também apontando para concursos públicos, associada à informatização dos testes. O evento contará com um minicurso, conferencistas internacionais e nacionais, além de seções técnicas orais e pôster. Também haverá concursos em nível de Iniciação Científica e Dissertação de Mestrado na área de TRI. Os conferencistas convidados são: * Wim van der Linden, CTB/McGraw-Hill, EUA * Michal Kosinski, University of Cambridge, U.K. * Mariana Curi (ICMC/USP) * Fernando Moura (UFRJ) * Ricardo Primi (Univ. São Francisco) Mais informações em:http://www.conbratri.org. Os emails para contato são conbratri(a)conbratri.org A submissão de trabalhos inicia em 16/09 e encerra em 20/10. As regras de submissão estarão na página do evento. Haverá desconto para os associados, em dia com a sua anuidade, nas seguintes associações: ABAVE, ABE e Instituto Brasileiro de Avaliação Psicológica (IBAP). Gostaríamos de convidar os membros da nossa comunidade para participar do 3º CONBRATRI em Belém do Pará, bem como na divulgação do evento. A seguir, alguns links sobre a cidade. http://pt.wikipedia.org/wiki/Belém http://www.paraturismo.pa.gov.br/ http://www.belemdopara.tur.br/ http://www.marajo.tur.br/ Atenciosamente, Heliton Tavares (Coordenador, UFPA) Ruben Klein (Presidente da ABAVE) Dalton Andrade (UFSC) <http://www.showmetech.com.br> *LEANDRO MARINO* | Showmetech Estatístico | Fotógrafo Mobile: 55 21 9845-7707 Email: contato(a)leandromarino.com.br Profiles SMT: facebook <http://www.facebook.com/showmetech> | twitter<http://www.twitter.com/showmetech> | Google+ <https://plus.google.com/115440786884851522659> <http://www.leandromarino.com.br>

1 0

Eixo x gráfico
by Fátima Lima Paula 03 Out '13

03 Out '13

Pessoal, fiz um gráfico usando o comando plot. O eixo x deveria ser "60-64 65-69 70-74 75-79 80+" Só que se coloco idade=c("60-64","65-69","70-74","75-79","80+"), no momento de plotar ele me retorna Erro em plot.window(...) : need finite 'xlim' values Além disso: Mensagens de aviso perdidas: 1: In xy.coords(x, y, xlabel, ylabel, log) : NAs introduzidos por coerção 2: In min(x) : nenhum argumento não faltante para min; retornando Inf 3: In max(x) : nenhum argumento não faltante para max; retornando -Inf Só consigo plotar se coloco idade=c("1","2","3","4","5") Como resolvo esse problema? Obrigada Fátima -- "Minha felicidade depende da qualidade dos meus pensamentos"

3 3

Integração entre R e Calc
by Andre Oliveira 03 Out '13

03 Out '13

Pessoas boa noite, gostaria de saber se existe alguma pacote ou forma de integrar R + calc? André Oliveira Souza

3 2

Fortran 90 + R: como carregar a dll
by Cleber N.Borges 03 Out '13

03 Out '13

olá a todos, estou tentando aprender o básico de como construir um código em Fortran e usá-lo no R. Para tanto, criei o segundo código em Fortran: ------------------------------------------------------------------------------------------------ subroutine contador ! Compilar como DLL para uso no R implicit none integer n, i integer, allocatable :: r(:) read *, n allocate( r(n) ) !!!!!!!!!!!!!!!!!!!!!!!!!!!!! do i = 1, n r(i) = i enddo !!!!!!!!!!!!!!!!!!!!!!!!!!!!! end subroutine ------------------------------------------------------------------------------------------------ No prompt do DOS, eu fui tentar usar o comando: R CMD SHLIB contador.f90 e obtive o seguinte erro: ------------------------------------------------------------------------------------------------ C:\Users\CLEBER\Documents\code\contador>R CMD SHLIB contador.f90 C:/R/etc/x64/Makeconf:196: warning: overriding commands for target `.m.o' C:/R/etc/x64/Makeconf:189: warning: ignoring old commands for target `.m.o' gfortran -m64 -O2 -mtune=core2 -c contador.f90 -o contador.o gfortran -m64 -shared -s -static-libgcc -o contador.dll tmp.def contador.o -Ld:/RCompile/CRANpkg/extralibs64/local/lib/x64 -Ld:/RCompile/CRANpkg/extralibs64/local/lib -LC:/R/bin/x64 -lR /usr/bin/sh: line 8: gfortran -m64: command not found ------------------------------------------------------------------------------------------------ não entendi o motivo do erro... Gostaria de alguma ajuda nesse Ponto (1) <- *** Contornei o problema "Copiando & Colando" os comandos no prompt. O arquivo * tmp.def * não existe... mas então eu o criei por que havia compilado anteriormente com a IDE Code Blocks, e sabia da existência de um fulano *libcontador.def* +++++++++++++++++++++++++++++ o conteúdo de def é apenas: EXPORTS contador_ @1 +++++++++++++++++++++++++++++ ------------------------------------------------------------------------------------------------ gfortran -m64 -O2 -mtune=core2 -c contador.f90 -o contador.o gfortran -m64 -shared -s -static-libgcc -o contador.dll tmp.def contador.o -Ld:/RCompile/CRANpkANpkg/extralibs64/local/lib -LC:/R/bin/x64 -lR ------------------------------------------------------------------------------------------------ Após os comandos acima é gerado o arquivo contador.dll No R, eu fui tentar carregar a dll, assim: > > dyn.load('contador.dll') > > dlls <- getLoadedDLLs() > dlls Filename Dynamic.Lookup base base FALSE utils C:/R/library/utils/libs/x64/utils.dll FALSE methods C:/R/library/methods/libs/x64/methods.dll FALSE grDevices C:/R/library/grDevices/libs/x64/grDevices.dll FALSE graphics C:/R/library/graphics/libs/x64/graphics.dll FALSE stats C:/R/library/stats/libs/x64/stats.dll FALSE contador C:/Users/CLEBER/Documents/code/contador/contador.dll TRUE > > getDLLRegisteredRoutines('contador') data frame with 0 columns and 0 rows > > Porém não há registro da subrotina 'contador'... :-( Usei o comando: .Fortran('contador_', n=as.integer( 10 ) ) mas daí trava o R. Se alguém puder me ajudar a entender isso, agradeço desde já. cleber

4 7

correção de tendência direcional
by Aline Kieskoski 01 Out '13

01 Out '13

Olá,estou trabalhando em um conjunto de dados e verifiquei a presença de tendência direcional linear dos dados em relação ao eixo Y, quando vou fazer o ajuste de modelos estou utilizando a função likfit da seguinte forma: dexp.ml<-likfit(d,ini=c(0.003,50),lambda=1, lik.method= "ML", cov.model="exp", trend=~d$coords[,2]) utilizei o comando trend sendo que d$coords[,2] está referindo-se as coordenadas do eixo Y. O R até ajusta um modelo e responde com os parâmetros estimados beta0 e beta1, porém, para a escolha do melhor modelo ajustado estou utilizando a validação cruzada, quando vou rodar a validação o R me responde com: xvalid: performing cross-validation at location ... 1, Erro em solve.default(ttivtt, crossprod(ivtt, as.vector(data))) : sistema é computacionalmente singular: condição recíproca número = 5.7415e-23 Testei o ajuste para os modelos exponencial, gaussiano e matérn, todos utilizando o método de máxima verossimilhança. Pensei que o problema pudesse ser nos valores dos chutes iniciais, porém, fiz o teste trocando esses valores e mesmo assim a validação não funcionou. Att, Aline.

1 0

Produção de Gráficos: PCA
by Vinícius Lionel Mateus 01 Out '13

01 Out '13

Prezados, Estou começando a aplicar PCA usando R, mas anteriormente o fazia no Statistica. Tenho utilizado os pacotes 'psych' e 'GPArotation' pois prefiro utilizar rotações em Varimax, porém os gráficos não ficaram ao meu gosto. Depois de algum esforço, eu consegui trocar as labels do meu gráfico, mas além disto gostaria de retirar o histograma, e plotar gráficos da dispersão das amostras(labels) e das variáveis entres os componentes, em 2D e 3D. Segue meu CMR. Desde já, obrigado. # instalando os pacotes install.packages("psych") install.packages("GPArotation") #carregando require (psych) install.packages (GPArotation) #entrada de dados (https://www.dropbox.com/s/l8v9g6qub27oig0/enqalog.xls) read.table(file = "clipboard", header = TRUE, sep = "\t", dec = ",")-> d #entrada de dados1 read.table(file = "clipboard", header = TRUE, sep = "\t", dec = ",")-> enqa enqa -> d attach(d) names(d) summary(d) #PCA(Varimax)/psych dpca<-principal(d[,2:16], nfactors=3, rotate="varimax", scores=T) #Scree plot plot (dpca$values, type="b", ylab="Eigenvalue", xlab="Component", lab=c(5,5,5)) #biplot biplot(dpca, labels=d[,1], cex=c(0.75,1)) Vinícius -- -- Atenciosamente, VINÍCIUS LIONEL MATEUS, M.Sc (http://lattes.cnpq.br/6501001637020665) Bacharel em Química - Doutorando em Química Analítica Laboratório de Química Atmosférica - Departamento de Química Pontifícia Universidade Católica - Rio de Janeiro (PUC - Rio) Rua Marquês de São Vicente, 225, Gávea - Rio de Janeiro, RJ - Brasil CEP.: 22451-900 Telefone: (+55) (21) 3527-1327 (+55) (21) 9358-8051 www.puc-rio.br

2 2

Subscrito em textos
by Tiago Souza Marçal 01 Out '13

01 Out '13

Boa noite, gostaria de saber se ao escrever um texto em R é possível representar 2^ns na forma 2ns. CMR cat("\n 2^ns \n") Att. Tiago. ################################################################# Tiago de Souza Marçal - Graduando em Agronomia pelo CCA-UFES Bolsista de Iniciação Científica da área de Genética e Melhoramento de Plantas #################################################################

2 1

Desdobramento de interações em dados com distribuição de Poisson inflacionada por zeros
by ASANTOS 30 Set '13

30 Set '13

Bom dia Pessoal, Estou quebrando a cabeça para tentar desdobrar interações em dados com distribuição de Poisson inflacionada por zeros, na verdade é uma situação muito recorrente em entomologia, onde em algumas situações temos zero insetos e em outras muitos, considerando-se dados de contagens. Fiz uma pequena simulação abaixo para ilustrar o problema (Variáveis observadas Trat e Epoca e variável resposta y), se acho um modelo significativo que tem uma parte com distribuição de Poisson e outra com Binomial, como faço o desdobramento? Uma para cada distribuição ou tem alguma abordagem diferente? Segue CRM abaixo: ## Desdobramento de dados com distribuição de Poisson inflacionado por zeros --- require(pscl) y1<- c(mapply(rpois, lambda=c(5,20,45), MoreArgs=list(n=400)))##Criação da variável resposta Poisson y2<- c(mapply(rbinom, size=c(1,0,0), prob=c(0.5,1,1), MoreArgs=list(n=200)))##Criação da variável resposta Binomial y<-c(y1,y2) trat <- as.factor(gl(3,600)) ##Criação dos tratamentos tempo<- as.factor(rep(gl(6,100),3)) ### Criação da variável tempo dados<-as.data.frame(cbind(trat,tempo, y)) #------------------------------------------------------------------------------- # Análise de variância do dados inflacionados summary(m1 <- zeroinfl(y ~ trat | epoca, data = dados)) ## Modelo completo mnull <- update(m1, . ~ 1) ### Modelo nulo pchisq(2 * (logLik(m1) - logLik(mnull)), df = 2, lower.tail = FALSE) ## Teste de Chi o modelo completo foi significativo ## Comparando com GLM Poisson sem inflação por zeros --------------------------- summary(p1 <- glm(y ~ trat + epoca, family = poisson, data = dados)) ##Teste de Vuong---------------------------------------------------------------- vuong(p1, m1) ## O GLM Poisson estava tão mal ajustado que o Poisson iflacionado é a única opção # ## E agora, como desdobrar isto? Faço um desdobramento para parte de poisson e outra para parte binomial? 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 ======================================================================

4 13

[OFF] Gerenciador de referencias
by Andre Oliveira 30 Set '13

30 Set '13

Alguém me indica um bom gerenciador de referencias. Sem ser o JabRef, pois, ele não abre aqui no meu ubuntu? André Oliveira Souza

9 15

[OFF-TOPIC] regex
by Roney Fraga Souza 30 Set '13

30 Set '13

Caros, Tenho um arquivo em formato txt com abstracts de artigos científicos, segue exemplo: We examine the productive efficiency of 70 Indian commercial banks during the early stages (1986-1991) of the ongoing period of liberalization. We use data envelopment analysis to calculate radial technical efficiency scores. We then use stochastic frontier analysis to attribute variation in the calculated efficiency scores to three sources: a temporal component, an ownership component, and a random noise component. We find publicly-owned Indian banks to have been the most efficient, followed by foreign-owned banks and privately-owned Indian banks. We also find a temporal improvement in the performance of foreign-owned banks, virtually no trend in the performance of privately-owned Indian banks, and a temporal decline in the performance of publicly-owned Indian banks. We attempt to explain these patterns in terms of the government's evolving regulatory policies. (C) 1997 Elsevier Science B.V Recently, considerable attention has been focused on the performance of various airlines and air carriers in terms of efficiency. Although it is obvious that air carriers use airports, few studies have focused on airport operational efficiency. This empirical study evaluates the operational efficiencies of 44 major U.S. airports using data envelopment analysis and some of its recent developments. Various airport characteristics are evaluated to determine their relationship to an airport's efficiency. Efficiency measures are based on four resource input measures including airport operational costs, number of airport employees, gates and runways, and five output measures including operational revenue, passenger flow, commercial and general aviation movement, and total cargo transportation. The results of this study have operational as well as public policy implications. (C) 2000 Elsevier Science B.V. All rights reserved considerando que cada parágrafo esta em uma linha, ou seja, a quebra de linha existe apenas no final de cada abstract, preciso apagar tudo que estiver a partir do '(C)' até a quebra de linha. No primeiro parágrafo preciso eliminar '(C) 1997 Elsevier Science B.V', no segundo parágrafo '(C) 2000 Elsevier Science B.V. All rights reserved'. Uso um macbook, qualquer sugestão via R ou terminal é bem vinda. Atenciosamente Roney

2 2

Re: [R-br] ajuda normatização de grid
by Eder Comunello 27 Set '13

27 Set '13

Prezado Elias, Fornecendo dados junto com seu código, permitindo torná-lo reproduzível (CMR), você obtém ajuda mais facilmente. Suponho que o 17 seja em função da variável que você quer representar no grid que está exportando (predição). Estou supondo que 17 seja o índice da última coluna adicionada (cbind). Um objeto SpatialPixelsDataFrame permite multivariáveis, mas pode ser que o formato escolhido exija a seleção de uma única variável pra compor um grid. ================================================ Éder Comunello Em 27 de setembro de 2013 07:43, ELIAS MENDES COSTA MENDES escreveu: > > Obrigado Éder, era isso mesmo, agora que rodei um comando para tirar todos > os NAs deu certo. > Já aproveitando o email talvez pudesse me ajudar mais uma vez. No script > que estou usando o comando que uso para salvar o mapa predito em um formato > que > abra no SagaGis é: writeGDAL(grid_V[,17], "V_rna37.sdat", > drivername="SAGA", mvFlag=-99999) Você sabe explicar o porque do 17? > > Os comandos anteriores a esse foram. > > # fazer o mapa para V%, ou aplicar o modelo para a area de estudo > grid_V <- cbind(grid1, predito) # usar "predito" ou "test1" > names(grid_V) > gridded(grid_V) <- ~X+Y # transformar "grid_V" de DataFrame para > SpatialPixelsDataFrame > proj4string(grid_V) <- CRS(utm) # definir projecao > names(grid_V) > spplot(grid_V, "predito") # plotar no R apenas a coluna "predito" > >

1 0

Re: [R-br] ajuda normatização de grid
by Eder Comunello 27 Set '13

27 Set '13

Prezado Elias, boa noite! Pra ajudar de fato seria necessário ter acesso aos dados que você está manipulando (ou uma amostra). Mesmo assim vou arriscar um palpite. Você removeu NA's com base somente na coluna 14 (gride[,14]) e por fim, normaliza de 4 a 18 (grid1[,4:18]). É provável que outras colunas além da 14 tenham valores de NA e estejam causando o erro. Para verificar o número de ocorrências de NA's em grid3: length(grid3[is.na(grid3)]) Espero que ajude, -- Éder Comunello <e <comunello.eder(a)gmail.com>comunel(a)gmail.com> Dourados, MS - [22 16.5'S, 54 49'W]

1 0

[OFF] Formatar classe de uma distribuição
by Andre Oliveira 26 Set '13

26 Set '13

Pessoal afim de ficar bem formatado eu uso o símbolo \vdash no LaTeXpara montar as classes de uma distribuição de dados contínuos, porém estou precisando de construir um intervalo incluindo as duas extremidades para a última classe. Alguém conhece algum símbolo em LaTeX que faça isto? Eu usei \vdash \dashv mas isto é um gambiarra feia! obrigado André Oliveira Souza

1 0

Análise de Série Temporal - Modelo
by Edson Lira 26 Set '13

26 Set '13

Caros amigos, estou com uma dúvida, estou analisando uma série de dados de doadores aptos de 2000 a 2011. Observando a série e fazendo o teste das raizes, conclui que a mesma é estacionária sem haver a necessidade da diferenciação. Usando a função auto.arima do pacote forecast, ele me indica um modelo com sazonalidade, porém observando os gráficos de correlação e correlação parcial não identifico o modelo sugerido pelo pacote, alguém poderia me dar uma sugestão? No link abaixo está a saída com os gráficos. http://www.datafilehost.com/d/88912c37 [ ]'s. Edson Lira Estatístico Manaus-Amazonas

3 3

Banco de dados muito grande
by Edson Lira 26 Set '13

26 Set '13

Caros amigos, boa noite, tenho um banco de dados com 10 mb e 176 mil linhas com 11 colunas. Estou usando: tra<-read.xlsx("/Users/edsonlira/Minuza/Dados_geral.xlsx",1,header=T,encoding= "UTF-8") Me mostra o seguinte erro: Erro em .jcall("RJavaTools", "Ljava/lang/Object;", "invokeMethod", cl, : java.lang.OutOfMemoryError: Java heap space O que devo fazer? [ ]'s. Edson Lira Estatístico Manaus-Amazonas

8 15

Divida em qual analise usar para resolver um problema de pesquisa
by Alisson Lucrecio 25 Set '13

25 Set '13

Caro Colegas da r-br, Boa tarde. Estou tentando ajudar um colega de universidade, mas não sei qual a ferramenta mais apropria para resolver esse problema, assim gostaria da ajuda de vocês. Um amigo fez análise 14 elementos de terra rara em 20 fertilizante por dois métodos de abertura, assim fatorial 20 (Fertilizantes) x 2 (Métodos de abertura) com 3 repetições. Ele quer responder 2 peguntas: os fertilizantes são diferentes entre si considerando todos os elementos de terra rara? Os métodos de abertura apresentam o mesmos valor de elemento de terra rara? Se não, em quais fertilizantes, os métodos de abertura superestima ou subestima o valor quando comparado ao outro método. Obrigado Alisson Lucrécio da Costa

1 0

combinar variáveis - PNAD
by Luciane Maria Pilotto 25 Set '13

25 Set '13

Olá pessoal, gostaria de ajuda para combinar variáveis. Estou trabalhando com dados da PNAD 2008 e preciso criar uma nova variável com a combinação de outras 3 variáveis. As variáveis são: 1321 (tem direito a plano de saúde) - com respostas 1,3 e 5 1342 (dá direito à assistencia odontologica) - respostas 2 e 4 1344 (plano exclusivo para odontologia) - respostas 2 e 4 Variável a ser criada (tipo de plano): Sem plano: se 1321=5 Plano com odontologia: se 1321=1 ou 1321=3 e 1342=2 e 1344=2; 1321=1 ou 1321=3 e 1342=2 e 1344=4 Plano sem odontologia: 1321=1 ou 1321=3 e 1342=4 e 1344=4 Atenciosamente, Luciane Maria Pilotto

1 0

Combinar arrays de diferentes comprimentos em data.frame
by Thiago V. dos Santos 24 Set '13

24 Set '13

Olá pessoal, Como resultado da leitura de algumas images de satélite, tenho o seguinte array no R > dim(ncdata) [1] 18 15 16 2555 > length(ncdata) [1] 11037600 onde as dimensões representam 18 é a quantidade de longitudes, 15 é a quantidade de latitudes, 16 é a quantidade de tipos de vegetação (floresta tropical, cerrado etc), 2555 é a quantidade de datas em que as imagens foram feitas. Por outro lado, tenho os seguintes arrays: > length(nclon) [1] 18 > length(nclat) [1] 15 > length(pft) [1] 16 > length(nctime) [1] 2555 E, em suma, o que eu quero fazer é criar um data.frame com a seguinte estrutura: nclon1 nclat1 pft1 ncdate1 (primeiro valor do ncdata) nclon1 nclat1 pft1 ncddate2 (segundo valor do ncdata) . . . nclon18 nclat15 pft16 nctime2555 (último valor do ncdata) Como a quantidade de linhas é muito grande, o ideal seria até usar alguma técnica que permita paralelização usando os cores do pc. Alguém conhece uma maneira (realmente) eficiente de fazer isso? Saudações, -- Thiago V. dos Santos PhD student Land and Aé qtmospheric Science University of Minnesota http://www.laas.umn.edu/CurrentStudents/MeettheStudents/ThiagodosSantos/ind… Phone: (612) 323 9898

2 2

Re: [R-br] r baixando dados inmet
by Éder Comunello 24 Set '13

24 Set '13

Bom dia! Também tenho interesse no tópico, pois costumava baixar os arquivos um a um, diretamente no site! :D Examinando o arquivo da lista, tem uma discussão sobre o RCurl em: <http://r-br.2285057.n4.nabble.com/R-br-RCurl-td4659610.html> Seguindo a discussão de Daniel Marcelino e Raphael Saldanha, fiz umas adaptações no código postado por eles para o caso em questão: ### <begin> library(RCurl) # Links para login e dados myURL1 <- "http://www.inmet.gov.br/projetos/rede/pesquisa/inicio.php" myURL2 <- " http://www.inmet.gov.br/projetos/rede/pesquisa/gera_serie_txt.php?&mRelEsta…,,,,, " ### como é só um teste, encurtei a série pesquisada (mRelDtInicio=01/01/2012) # Dados para acesso myParams=list( mCod="seuemail", ### alterar! mSenha="suasenha", ### alterar! btnProcesso = " Acessar ") myCurl <- getCurlHandle() curlSetOpt(cookiejar="cookies.txt", useragent="Mozilla/5.0", followlocation=TRUE, curl=myCurl) login <- postForm(myURL1, .params=myParams, curl=myCurl) dados <- getURLContent(myURL2, curl=myCurl) head(strwrap(dados, 70), 30) ### <end> O código acima faz o login e baixa os dados de um único exemplo (myURL2). Tá testado e funcionando!!! Logicamente requer cadastro no BDMEP... Agora falta entrar com a limpeza e importação do código armazenado na variável dados e automatizar a geração de links para dados (basicamente mudar mRelEstacao, mRelDtInicio e .mRelDtFim em cada link). Mas acredito que já seja um começo! :D Éder Comunello ### Uma saída parcial do conteúdo armazenado em dados: > head(strwrap(dados, 70), 30) [1] "<html><meta http-equiv=\"Content-Type\" content=\"text/html;" [2] "charset=iso-8859-1\"> <head> <title>Estação - 83595</title> \r <script src=\"facebox/js/jquery.js\"" [4] "type=\"text/javascript\"></script>\r <link" [5] "href=\"facebox/styles/facebox.css\" media=\"screen\" rel=\"stylesheet\"" [6] "type=\"text/css\" >\r \r <script src=\"facebox/js/facebox.js\"" [7] "type=\"text/javascript\"></script>\r <script type=\"text/javascript\">\r" [8] "jQuery(document).ready(function($) {\r $('a[rel*=facebox]').facebox({\r" [9] "loading_image : 'loading.gif',\r close_image : 'closelabel.gif'})\r" [10] "$('area[rel*=facebox]').facebox({\r loading_image : 'loading.gif',\r" [11] "close_image : 'closelabel.gif'})\r })\r </script>\r <script" [12] "type=\"text/javascript\">\r function showFoto(src) {\r //criando um link" [13] "falso, para o LightBox captar as informações\r var a_false =" [14] "document.createElement(\"A\");\r a_false.rel = \"lightbox\";\r a_false.href" [15] "= src;\r //chamando a função do lightbox responsável pelo início das" [16] "ações\r myLightbox.start(a_false);\r //anulando o link\r a_false =" [17] "null;\r }\r </script>\r </head> <pre>-------------------- BDMEP - INMET" [18] "-------------------- Estação : AIMORES - MG (OMM: 83595) Latitude" [19] "(graus) : -19.49 Longitude (graus) : -41.07 Altitude (metros): 82.74" [20] "Estação Operante Inicio de operação: 01/06/1972 Periodo solicitado" [21] "dos dados: 01/01/2012 a 01/01/2013 Os dados listados abaixo são os" [22] "que encontram-se digitados no BDMEP -------------------- Obs.: Os" [23] "dados aparecem separados por ; (ponto e vírgula) no formato txt." [24] "Para o formato planilha XLS, <A HREF=\"instrucao.html\" TARGET=\"_top\"" [25] "rel=\"facebox\">siga as instruções</a> --------------------" [26] "Estacao;Data;Hora;TempBulboSeco;UmidadeRelativa;PressaoAtmEstacao;DirecaoVento;VelocidadeVentoNebulosidade;" [27] "83595;01/01/2012;0000;25.2;83;1004.8;14;1;8;" [28] "83595;01/01/2012;1200;26;84;1004.9;9;1;10;" [29] "83595;01/01/2012;1800;27;91;1000;9;1;10;" [30] "83595;02/01/2012;0000;25;92;1003;14;1;10;"

2 2

r baixando dados inmet
by Alisson Lucrecio 24 Set '13

24 Set '13

Caro Colegas, Boa noite. Eu estou tentando baixar uma série de dados meteorológico do inmet, assim preciso automatizar o processo, alguém sabe como o r poderia ler um arquivo dessa forma e salva-lo em .txt? O dado é proveniente do site do inmet http://www.inmet.gov.br/projetos/rede/pesquisa/gera_serie_txt.php &mRelEstacao=83595&btnProcesso=serie&mRelDtInicio=01/01/1961&mRelDtFim=01/01/2013&mAtributos=1,,,,1,1,,1,1,,,1,,,,, -------------------- BDMEP - INMET -------------------- Estação : AIMORES - MG (OMM: 83595) Latitude (graus) : -19.49 Longitude (graus) : -41.07 Altitude (metros): 82.74 Estação Operante Inicio de operação: 01/06/1972 Periodo solicitado dos dados: 01/01/1961 a 01/01/2013 Os dados listados abaixo são os que encontram-se digitados no BDMEP -------------------- Obs.: Os dados aparecem separados por ; (ponto e vírgula) no formato txt. Para o formato planilha XLS, siga as instruções -------------------- Estacao;Data;Hora;TempBulboSeco;UmidadeRelativa;PressaoAtmEstacao;DirecaoVento;VelocidadeVentoNebulosidade; 83595;01/06/1972;1200;19.8;84;1011.6;36;1;8; 83595;01/06/1972;1800;27.2;47;1007.8;9;1;9; 83595;02/06/1972;0000;23;67;1009.5;14;3;10; 83595;02/06/1972;1200;22;76;1011;32;1;4; 83595;02/06/1972;1800;29;45;1006.4;32;1;4; 83595;03/06/1972;0000;19;82;1009;14;1;10; 83595;03/06/1972;1200;22.2;78;1011.4;32;1;5; 83595;03/06/1972;1800;30.8;47;1005.6;32;1;8; 83595;04/06/1972;0000;20.2;85;1008.6;23;1;10; 83595;04/06/1972;1200;23.8;78;1010.7;36;1;8; 83595;04/06/1972;1800;30.8;47;1007.4;32;1;6; 83595;05/06/1972;0000;25;81;1008.8;14;5;10; 83595;05/06/1972;1200;22;83;1010;32;1;4; 83595;05/06/1972;1800;29;57;1007.2;14;1;4; 83595;06/06/1972;0000;24;75;1008.6;14;3;10; 83595;06/06/1972;1200;21;81;1010.2;32;1;3; 83595;06/06/1972;1800;29;54;1005.8;36;1;8; 83595;07/06/1972;0000;24.5;73;1008.4;14;1;10; 83595;07/06/1972;1200;25.4;71;1009.6;9;1;7; 83595;07/06/1972;1800;25.8;71;1007;5;3;9; 83595;08/06/1972;0000;23;76;1008.7;14;1;10; 83595;08/06/1972;1200;23.5;80;1012.3;32;1;5; 83595;08/06/1972;1800;30;52;1007.2;32;1;6; 83595;09/06/1972;0000;23;76;1009.2;14;1;10; 83595;09/06/1972;1200;24.5;77;1010.3;5;1;5; 83595;09/06/1972;1800;29;53;1006.3;5;1;6; 83595;10/06/1972;0000;24;80;1008.7;14;3;10; 83595;10/06/1972;1200;23;82;1010.8;0;0;3; 83595;10/06/1972;1800;29.8;50;1005.2;5;1;5; 83595;11/06/1972;0000;24;77;1008.4;14;1;10; 83595;11/06/1972;1200;23;80;1010;5;1;3; 83595;11/06/1972;1800;30;50;1007.1;23;1;4; 83595;12/06/1972;0000;21;77;1009;18;1;10; 83595;12/06/1972;1200;20;87;1010.4;23;1;10; Obrigado. Alisson Lucrécio da Costa

3 3

Re: [R-br] Vetor de caracter como objeto em um pacote
by Ivan Bezerra Allaman 24 Set '13

24 Set '13

Complexo!! Já tentou system(paste(....)) ou cat(..., file = Sys.getenv('APPDATA')) ?? O problema de lhe ajudar, é que não tem como gerar um CRM (eu acho) para testar aqui!! Aí fica mais difícil só na imaginção!! Abraços!!

3 2

Re: [R-br] Notação potências de dez (Gráficos)
by Ivan Bezerra Allaman 22 Set '13

22 Set '13

Em 2006 um sabido respondeu no r-help, sciNotation <- function(x, digits = 1) { if (length(x) > 1) { return(append(sciNotation(x[1]), sciNotation(x[-1]))) } if (!x) return(0) exponent <- floor(log10(x)) base <- round(x / 10^exponent, digits) as.expression(substitute(base %.% 10^exponent, list(base = base, exponent = exponent))) } plot(1:1000, axes = FALSE, type = "l", frame.plot = TRUE) axis(1, at = axTicks(1), label = sciNotation(axTicks(1), 1)) Deve ser o que procuras!!! (s,f,p) Allaman

1 0

Re: [R-br] Vetor de caracter como objeto em um pacote
by Ivan Bezerra Allaman 22 Set '13

22 Set '13

Não tive problema com as dependências! De fato agora não apareceu mesmo, ou seja, obtive o seguinte: trPaths [1] "\\Tinn-R\\tmp\\" "\\Tinn-R\\tmp\\search.txt" [3] "\\Tinn-R\\tmp\\objects.txt" "\\Tinn-R\\tmp\\file.r" [5] "\\Tinn-R\\tmp\\selection.r" "\\Tinn-R\\tmp\\block.r" [7] "\\Tinn-R\\tmp\\lines.r" "\\Tinn-R\\tmp\\reformat-input.r" [9] "\\Tinn-R\\tmp\\reformat-output.r" Abraços!

1 0

Re: [R-br] Vetor de caracter como objeto em um pacote
by Ivan Bezerra Allaman 22 Set '13

22 Set '13

Ok! Agora que entendi que Sys lalala é uma função. Pois bem!! Estou no windows agora e reproduzi o que você postou, e não tive NENHUM PROBLEMA, ou seja, ao rodar trPaths <- paste(paste(as.vector(Sys.getenv('APPDATA')), '\\Tinn-R\2013/09/22 - 12:20:51mp\\', sep=''), c('', 'search.txt', 'objects.txt', 'file.r', 'selection.r', 'block.r', 'lines.r', 'reformat-input.r', 'reformat-output.r'), sep='') obtive trPaths [1] "C:\\Users\\ivan\\AppData\\Roaming\\Tinn-R3/09/22 - 12:20:51mp\\" [2] "C:\\Users\\ivan\\AppData\\Roaming\\Tinn-R3/09/22 - 12:20:51mp\\search.txt" [3] "C:\\Users\\ivan\\AppData\\Roaming\\Tinn-R3/09/22 - 12:20:51mp\\objects.txt" [4] "C:\\Users\\ivan\\AppData\\Roaming\\Tinn-R3/09/22 - 12:20:51mp\\file.r" [5] "C:\\Users\\ivan\\AppData\\Roaming\\Tinn-R3/09/22 - 12:20:51mp\\selection.r" [6] "C:\\Users\\ivan\\AppData\\Roaming\\Tinn-R3/09/22 - 12:20:51mp\\block.r" [7] "C:\\Users\\ivan\\AppData\\Roaming\\Tinn-R3/09/22 - 12:20:51mp\\lines.r" [8] "C:\\Users\\ivan\\AppData\\Roaming\\Tinn-R3/09/22 - 12:20:51mp\\reformat-input.r" [9] "C:\\Users\\ivan\\AppData\\Roaming\\Tinn-R3/09/22 - 12:20:51mp\\reformat-output.r" Logo, Qual é o problema?????????????????????????????????????????

2 2