Dezembro 2012 - R-br

Cálculo do percentage prediction error em ZINB
by Humberto Hazin 20 Dez '12

20 Dez '12

Prezados, Estou com uma dúvida que pode parecer besteira, porém gostaria da ajuda de vocês! A dúvida é a seguinte como calcular o percentual de "Bias" (percentage prediction error) por linha de um modelo ZINB (package- pscl)? O meu objetivo é fazer um boxplot com esse percentual por ano Eu iniciei da seguinte forma: dois bancos de dados um de try.csv (10% retirados do banco original) e outro original.csv try<-read.csv(”try.csv”, header = TRUE, sep = ",", quote="\"", dec=".",fill = TRUE) https://www.dropbox.com/s/v7vruztus11nnbi/try.csv original<-read.csv(”original.csv”, header = TRUE, sep = ",", quote="\"", dec=".",fill = TRUE) https://www.dropbox.com/s/7htqc2bqi45rkku/original.csv # Criar Fatores para o banco original original$area<-rep(1,nrow(original)) original$area[original$lat<=-15&original$lat>=-45]<-2 original$A=as.factor(original$area) summary(original) original$S<-factor(original$id.group.fleet) original$Q<-factor(original$Quarter) original$Y<-factor(original$Year) original$logEf<-log(original$effort1) summary(original) # Criar Fatores para o banco try try$area<-rep(1,nrow(try)) try$area[try$lat<=-15&try$lat>=-45]<-2 try$A=as.factor(try$area) summary(try) try$S<-factor(try$id.group.fleet) try$Q<-factor(try$Quarter) try$Y<-factor(try$Year) try$logEf<-try$effort1*0 summary(try) zeroinfl.fit_original<-zeroinfl(BSH~Y+S+Q+A+offset(logEf)|Y+S+Q+A,dist="negbin",data=original) pred<-predict(zeroinfl.fit_original,newdata=try,type="response") Daqui em diante não tenho mais ideia de como gerar esse percentage prediction error Alguem poderia me ajudar! Desde já Agradeço Humberto Enviado do Email do Windows

2 1

problemas com R-br listas
by Fernando A. Souza 19 Dez '12

19 Dez '12

Caros colegas estou com um problema Gostaria de mudar a forma como eu recebo as postagem na lista. Geralmente recebo vários e-mails de todos os que postam. Recentemente um amigo fez o cadastro na lista e ele geralmente recebe um único e-mail com a listas das mensagens de um dado período. Gostaria de receber desta forma tambem. Pois evitaria a sobrecarca de minha caixa de e-mail. Algúem sabe como se faz isto? att -- Fernando Antônio de Souza Zootecnista, Msc. Nutricão Animal, UFMG lattes: http://lattes.cnpq.br/6519538815038307

3 3

Problema com strsplit()
by Augusto Ribas 19 Dez '12

19 Dez '12

Ola. Eu estou usando o pacote taxize, e estava tentando usar a função iucn_summary() deste pacote, mas ela funciona perfeitamente com os exemplos, mas quando eu vou usar um nome que digito de uma lista minha eu recebo a seguinte mensagem: ##### > iucn_summary("Molossus currentium") Erro em strsplit(pop, split = ":") : argumento modo não caractere ###### Eu olhei no site que a função busca e a espécie esta la. Não entendo o que acontece que ela não funcionou no meu caso, mas funcionou nos exemplos. Os exemplos que digo são esses: # library(taxize) iucn_summary("Lynx lynx") iucn_summary("aaaa") iucn_summary("Molossus currentium") # Eu não tenho ideia do que estou fazendo errado, mesmo quando a espécie não é encontrada no site, deveria ser retornado algo, como no exemplo do "aaaa". Eu estou usando Ubunto aqui e usando o R 2.15.2 do Emacs. -- Grato Augusto C. A. Ribas Site Pessoal: http://recologia.wordpress.com/ Lattes: http://lattes.cnpq.br/7355685961127056

2 1

Expression
by Fabio Mathias Corrêa 19 Dez '12

19 Dez '12

Pessoal, Ao utilizar o expression em um plot, o resultado não sai como o esperado. O mesmo argumento é usado no xlab e ylab, mas só funciona corretamente no xlab. CMR plot(1:10, xlab=expression(gamma[2]^"**"), ylab=expression(gamma[2]^"**")) Muito obrigado! Fábio Mathias Corrêa Universidade Estadual de Santa Cruz Departamento de Ciências Exatas e da Terra - DCET Campus Soane Nazaré de Andrade, km 16 Rodovia Ilhéus-Itabuna CEP 45662-900. Ilhéus-Bahia Tel.: 73-3680-5076

2 1

Chamada de Artigo Revista Mediações
by Ronaldo Baltar 18 Dez '12

18 Dez '12

Gostaríamos de convidar os interessados e também pedir a divulgação para Chamada de artigos – Revista Mediações, com o próximo dossiê sobre "Análise quantitativa e indicadores sociais". A Comissão Editorial de Mediações, Revista de Ciências Sociais, do Departamento e Programa de Pós-Graduação em Ciências Sociais, da Universidade Estadual de Londrina, receberá artigos para o volume 18, número 1, correspondente ao primeiro semestre de 2013. Esse Dossiê se propõe a incentivar a discussão sobre: aplicação de modelos quantitativos e indicadores sociais para análise das diferenciações sociais; procedimentos estatísticos para análise multivariada; agentes e simulação social; análise de redes sociais; entre outros. Serão especialmente bem-vindos artigos que usem R como código fonte para replicação científica. Um dos objetivos deste número é ajudar a disseminação do R como ferramenta para a pesquisa em ciências sociais. Organizadores: Prof. Ronaldo Baltar e Profª Claudia Siqueira Baltar Infosoc - Laboratório de Sociologia Computacional Departamento de Ciências Sociais Universidade Estadual de Londrina Recebimento de contribuições até 28 de fevereiro de 2013. Além do dossiê, Mediações publica trabalhos sobre temas diversos das Ciências Sociais na Seção Artigos, como também Resenhas de livros relevantes em Ciências Sociais. As submissões em português podem ser enviadas através do portal da Revista http://www.uel.br/revistas/uel/index.php/mediacoes ou pelo e-mail mediacoes(a)uel.br MEDIAÇÕES – REVISTA DE CIÊNCIAS SOCIAIS Departamento de Ciências Sociais/Centro de Letras e Ciências Humanas Universidade Estadual de Londrina Campus Universitário - Caixa Postal 6001 – Londrina/PR – 86055-900 Brasil

1 0

glm - variável factor
by Vinicius Brito Rocha 18 Dez '12

18 Dez '12

Senhores, tenho um modelo glm, logistico, onde um dos coef é uma variável factor com as Unidades Federativas. modelo=glm(y ~ UF+VA1+VA2-1) summary(modelo) O problema é que ao olhar o summary dos resultados do ajustes algumas das UF's não tem significância no teste para o parâmetro. OU seja tenho cenário com RJ significativo e SP não. Como faço para excluir a classe SP do meu novo modelo? Abs Vinicius -- *Vinicius Brito Rocha.* *Estatístico e Atuário (IM / UFRJ)** Mestre em Pesquisa Operacional (COPPE / UFRJ)* www.aplicademic.blogspot.com http://twitter.com/viniciusbritor "Não se preocupe muito com as suas dificuldades em Matemática, posso assegurar-lhe que as minhas são ainda maiores." - Albert Einstein.

3 4

Dúvida_experimentação_estimação
by andrebvs＠bol.com.br 18 Dez '12

18 Dez '12

1 0

Trabalhar com pares do RecordLinkage
by Fátima Lima Paula 17 Dez '12

17 Dez '12

Prezados, preciso trabalhar com os pares resultantes do RecordLinkage e não sei como fazer. Enviei email para o autor do pacote, mas ele disse que terei que programar, pois não existe uma maneira automática de fazer isso. Bem, aí vai, meu desafio: Tenho dois bancos ( um retirado do outro) e vou fazer o linkage para obter as readmissões. No resultado sempre haverá um par de cada um, com mesma data de internação, já que o banco menor é retirado do maior (o menor é de internados por fratura de fêmur e o maior é de internados por todas as causas). Aí vai um exemplo criado para você verem como é a saída dos pares: nome=c("Ricardo","Maria","Tereza","Pedro","José", "Germano") sobrenome=c("Cunha","Andrade","Silva","Soares","Silva","Lima") idade=c(67,89,78,65,68,67) ano=c(1945,1923,1934,1947,1944,1945) dataint=c(20120907,20120703,20120301,20120805,20121004,20121209) banco1=cbind(nome,sobrenome,idade,ano,dataint) banco1=data.frame(banco1) banco1 nome=c("Maria","Lúcia","Paulo","Marcos", "Ricardo", "Germanio") sobrenome=c("Andrade","Silva","Soares","Pereira","Cunha","Lima") idade=c(67,88,78,60,68,80) ano=c(1945,1924,1934,1952,1944,1932) dataint=c(20121208,20121103,20120302,20120105,20121004,20121209) banco2=cbind(nome,sobrenome,idade,ano,dataint) banco2=data.frame(banco2) banco2 str(banco1) banco1$nome=as.character(banco1$nome) banco1$sobrenome=as.character(banco1$sobrenome) banco1$idade=as.character(banco1$idade) banco1$ano=as.character(banco1$ano) banco1$dataint=as.character(banco1$dataint) banco2$nome=as.character(banco2$nome) banco2$sobrenome=as.character(banco2$sobrenome) banco2$idade=as.character(banco2$idade) banco2$ano=as.character(banco2$ano) banco2$dataint=as.character(banco2$dataint) library(RecordLinkage) soundexBR<-function(termo){ termo<-toupper(termo) # 1. Retire toda pontuação da palavra; termo<-gsub("[ÃÁÀÂÄ]","A",termo) termo<-gsub("[ÉÈÊ?Ë]","E",termo) termo<-gsub("[ÍÌÏÎI]","I",termo) termo<-gsub("[ÓÕÒÔÖ]","O",termo) termo<-gsub("[ÚÙÛÜU]","U",termo) termo<-gsub("Ç","C",termo) #1.a Substituir todas as letras duplas por uma única letra #http://www.archives.gov/genealogy/census/soundex.html termo<-gsub("([A-Z])\\1", "\\1", termo) # 2. Fique com a primeira letra; N<-nchar(termo) # Modificação PT-BR (Coeli e Camargo) termo<-ifelse(substr(termo,1,2)=="WA",sub("W","V",termo),termo) termo<-ifelse(substr(termo,1,1)=="H",substr(termo,2,N),termo) termo<-ifelse(substr(termo,1,2)=="KA"|substr(termo,1,2)=="KO"|substr(termo,1,2)=="KU",sub("K","C",termo),termo) termo<-ifelse(substr(termo,1,1)=="Y",sub("Y","I",termo),termo) termo<-ifelse(substr(termo,1,2)=="CE"|substr(termo,1,2)=="CI",sub("C","S",termo),termo) termo<-ifelse(substr(termo,1,2)=="GE"|substr(termo,1,2)=="GI",sub("G","J",termo),termo) ############## Segunda Parte ####################### termo.1<-substr(termo,1,1) termo<-substr(termo,2,N) # 3. Mude todas as ocorrências das letras a seguir por '0' (zero): 'A','E','I','O','U','H','W','Y'; termo<-gsub("[A,E,I,O,U,H,W,Y]",0,termo) # 4. Mude as letras restantes de acordo com a tabela abaixo: # Número Letra # 1 'B','F','P','V' termo<-gsub("[B,F,P,V]",1,termo) # 2 'C','G','J','K', 'Q','S','X','Z' termo<-gsub("[C,G,J,K,Q,S,X,Z]",2,termo) # 3 'D','T' termo<-gsub("[D,T]",3,termo) # 4 'L' termo<-gsub("L",4,termo) # 5 'M','N' termo<-gsub("[M,N]",5,termo) # 6 'R' termo<-gsub("R",6,termo) # 5. Remova todos os zeros da string resultante; termo<-gsub(0,"",termo) # Eliminar dois numeros iguais consecutivos por um único numero #http://www.archives.gov/genealogy/census/soundex.html termo<-gsub("([0-9])\\1", "\\1", termo) # Remontar termo<-paste(termo.1,termo,sep="") # 6. Preencha a string resultante com zeros à direita e retorne desta forma: #<letra maiuscula><digito><digito><digito> termo<-paste(termo,"0000",sep="") termo<-substr(termo,1,4) return(termo) } soundexBR("termo") soundexBR("terno") pares=compare.linkage(banco1,banco2,blockfld=list(c(1,2,4),c(1,2)), phonetic=c(1,2),phonfun=soundexBR,strcmp=FALSE, strcmpfun=jarowinkler, exclude=FALSE,identity1=NA,identity2=NA, n_match=NA, n_non_match=NA) editMatch(pares) pares pesos=epiWeights(pares, e=0.01,f=pares$frequencies) hist(pesos$Wdata,plot=F) getPairs(pesos,max.weight=Inf, min.weight=-Inf) Agora, como tenho as datas da internação de todos, preciso verificar, par a par, quem tem data de internação diferente. Estes serão as readmissões. Existe uma forma de eu criar uma coluna para identificar os pares? E depois, testar se as datas de internação são diferentes? Obrigada Fátima

2 1

ifelse composto
by Edson Lira 17 Dez '12

17 Dez '12

Bom dia caros r-listeiros, estou com um probleminha. Tenho um banco com a variável Naturalidade que tem as seguintes respostas: MANAUS, BAHIA, RIO DE JANEIRO, ANAMÃ, BORBA,SALVADOR,CAREIRO... Quero criar uma variável chamda de prodedência na qual terei os seguintes níveis: MANAUS=CAPITAL ANAMÃ,BORBA,CAREIRO=INTERIOR BAHIA,RIO DE JANEIRO,SALVADOR=OUTRO ESTADO. COMO FAZER ISSO COM O ifelse? [ ]'s. Edson Lira Estatístico Manaus-Amazonas

4 5

média acumulada
by Leonard de Assis 17 Dez '12

17 Dez '12

Existe pronto um comando estilo cumsum que retorne ao invés de soma, a média? exemplo: set.seed(1) op <- par(mfrow = c(2, 2)) x <-rnorm(500) y <- cumsum(x) ts.plot(x) ts.plot(y) # para eu fazer o que quero, tenho que fazer: y.cmean <-cumsum(x)/(1:length(x)) ts.plot(y.cmean) par(op) existe algum jeito mais fácil? -- []s Leonard de Assis http://about.me/ldeassis

2 1