Outubro 2014 - R-br

PNAD
by Vanessa Felix 23 Out '14

23 Out '14

Pessoal, Eu tô querendo usar os dados da PNAD 2013 mas não consigo usar o pacote IBGEPesq. Fica aparecendo esse erro aqui: library(IBGEPesq) Error: package ‘IBGEPesq’ was built before R 3.0.0: please re-install it Tem alguma opção que não seja reinstalar o R? Se alguém puder ajudar... att, -- ------------------------------------------------------------------------------------------ Vanessa Felix do Nascimento Sergio Graduada em Estatística - ICE/UFJF Mestranda em Saúde Coletiva - IMS/UERJ (32) 8854-4881 / 9142-8428 (21) 97656-3448 -------------------------------------------------------------------------------------------

4 4

Comparação entre modelo linear e não linear
by Ari Clecius 23 Out '14

23 Out '14

Caros colegas, uso alguns modelos de isotermas de adsorção, para o modelo de Langmuir existem algumas linearizações, como faço para avaliar se o modelo não linear é melhor que o modelo linearizado? Att, Ari Clecius Alves de Lima Engenheiro Químico Me. Engenharia Civil (085)88412345 (085)33669042 #CRM ce=c(0.66450,0.92395,1.22510,1.53040,12.89600,32.44400,70.93000,135.54800) qe=c(0.65864,0.83815,1.09135,1.42049,11.57060,15.42198,16.41500,20.42020) #########Modelo linear#################### ce_inv=1/ce qe_inv=1/qe Langmuir_Linear=lm(qe_inv~ce_inv) summary(Langmuir_Linear) names(Langmuir_Linear) Langmuir_Linear$coefficients Langmuir_Linear$coefficients[1] qmax_L=1/Langmuir_Linear$coefficients[1] qmax_L kl_L=1/qmax*Langmuir_Linear$coefficients[2] kl_L qe_est_L=qmax_L* kl_L*ce/(1+ kl_L *ce) plot(ce,qe) lines(ce,qe_est_L,col="blue") ########################################## library(nlstools) adsorption=list(ce,qe) langmuir<-nls(qe~qmax*kl*ce/(1+kl*ce),adsorption,start=list(qmax=30,kl=0.02)) summary(langmuir) qe_Est_NL=21.40602* 0.07498 *ce/(1+ 0.07498 *ce) qe_Est_NL plot(ce,qe) lines(ce,qe_Est_NL,col="red") legend("top",legend=c("LangmuirExp","Langmuir_Est_NL"), lty=c(NA,1), col=c("black","red"),pch=c(3,NA),cex=0.6) ######################################## ############################################# plot(ce,qe) lines(ce,qe_est_L,col="blue") lines(ce,qe_Est_NL,col="red") legend("topleft",legend=c("LangmuirExp","Langmuir_Est_L","Langmuir_Est_NL"), lty=c(NA,1,1), col=c("black","blue","red"),pch=c(3,NA,NA),cex=0.6) #####################################################

1 0

Dúvida sobre análise de cluster em redes neurais
by Victor Eduardo 23 Out '14

23 Out '14

Boa noite, Ainda estou aprendendo sobre o assunto e estou com uma dúvida sobre clusterização em redes neurais utilizando o R. Eu encontrei dois pacotes que utilizam essa ferramenta, que é o nnet e kohonen, só que em nenhum deles, pelo que estava vendo, passa a formação dos clusters que tenho, ou seja, qual cluster o indivíduo pertence, como é visto no algoritmo de fuzzy c-means e no kmeans. É possível visualizar os clusters que cada indivíduo pertence? Existe algum comando que forneça essa informação dentro desses dois pacotes ou tem algum outro pacote com essa informação? Atenciosamente, Victor Eduardo

1 0

OFF- Livro Multivariada
by Andre Oliveira 22 Out '14

22 Out '14

Pessoal,qual livro de multivariada para um iniciante e em português. Achei estes .... MINGOTI, S. A. Análise de dados através de métodos de estatística multivariada: uma abordagem aplicada. Belo Horizonte: Editora UFMG, 2005. FERREIRA, Daniel Furtado.Estatística multivariada. 2. ed. Lavras : Editora UFLA, 2011 O que recomendam? André Oliveira Souza. Graduação em Matemática, mestrado em estatística aplicada.Instituto Federal de Educação, Ciência e Tecnologia do Espirito Santo. IFES

3 2

Revolution R - Paralelização
by Jefferson Ferreira-Ferreira 22 Out '14

22 Out '14

Listeiros; Alguém aí utiliza o Revolution R Enterprise? Recém o baixei. Li alguns materiais pela internet mas não estou entendendo muito bem se ele paraleliza por "default" ou se preciso usar algo como ParallelR ou doSNOW. Alguém saberia me dizer ou indicar um bom site onde posso me informar melhor sobre isso? Obrigado. Abraços -- *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*

5 6

Tabela de frequência com amplitudes desproporcionais
by Flávio Fagundes 21 Out '14

21 Out '14

Prezados(as), Por um acaso existe alguma função semi-pronta no R que produz uma tabela de frequência para intervalo de classes com amplitudes de tamanhos desproporcionais? Imagino ser simples a programação (vou começar uma agora), mas pode ser que já exista. Por exemplo: Tenho os seguintes dados: dados<-c(0, 7 ,7 ,8 ,8 , 9, 9, 10, 10, 10, 11, 11, 12, 13, 14, 14, 16, 17, 18, 18, 19, 19, 20, 22, 25) A função fdt me retorna: > tb<-fdt(dados, k=4) > tb Class limits f rf rf(%) cf cf(%) [0,6.31) 1 0.04 4 1 4 [6.31,12.6) 12 0.48 48 13 52 [12.6,18.9) 7 0.28 28 20 80 [18.9,25.2) 5 0.20 20 25 100 Gostaria de obter a mesma tabela acima, mas para as seguintes classes 0 |---- 10 10 |---- 13 13 |---- 20 20 |---- 25 utilizando, algo como o parâmetro breaks=c(0,10,13,20) Obrigado Flávio

3 3

Re: [R-br] skater e deldir
by Rafael Tieppo 21 Out '14

21 Out '14

Elias, Pelo o que eu entendi, vc usou a funcao deldir(), no caso Voronoi para fazer os polígono da área. Se for isso acho que passei pelo mesmo processo. Vou passar algumas etapas e uma função que improvisei e no meu caso deu certo. Vou explicar os procedimentos: 1) com a funcao deldir() vc faz seu polígonos na área sss <- deldir( .....) entao vc tem o objeto sss 2) Para gerar uma lista do objeto "sss" use a funcao tile.list() sss_list <- tile.list(sss) 3) Agora desse objeto que é uma lista, vc tem que extrair o conjunto de polígonos para isso improvisei uma função, grave ela no seu ambiente: FUN_VOR_MAT <- function (A,TT) #TT é o arquivo, N o número do POL { AA <- 1 BB <- 2 NROW <- length (A[[TT]]$x) + 1 #add primeiro ponto VORONOI_MATRIX <- matrix(0,NROW,2) VORONOI_MATRIX[1:NROW,AA] <- A[[TT]]$x[1:NROW] VORONOI_MATRIX[1:NROW,BB] <- A[[TT]]$y[1:NROW] VORONOI_MATRIX[NROW,AA] <- A[[TT]]$x[1] #+ 1 add primeiro ponto VORONOI_MATRIX[NROW,BB] <- A[[TT]]$y[1] #+ 1 add primeiro ponto if(TT == NROW){ print("17/10/2014, rafaeltieppo(a)yahoo.com.br ")} return(VORONOI_MATRIX) } 4) Agora vc vai usar a funcao: sss_areas <- lapply(1:length(sss_list), function (x) ( FUN_VOR_MAT(sss_list,x)) ) 5) #Formando Polygons sss_pols <- lapply(1:length(sss_areas), function(x) Polygons(list(Polygon(sss_areas[[x]])), ID = paste("reg", x, sep = ""))) 6) Para converter em SpatialPolygins sss_SPATIAL_POLS <- as.SpatialPolygons.PolygonsList(sss_pols, proj4string=CRS(as.character(NA))) Se tudo correu bem é só comecar a usar o poly2nb(sss_SPATIAL_POLS) Boa sorte Rafael Tieppo

2 1

Juntar colunas no R
by Giselle Davi 20 Out '14

20 Out '14

Prezados,Preciso juntar colunas no R !!! Tenho um arquivo assim: D01 D02 D03 D04 D05 D06 D07 D08 D09 D10 F01 F02 F03 F04 F05 F06 F07 F08 F091 G G G G G G . G G G G G G G G G G . .2 C C C C C T C C C C C C C C C C C C C3 T T T T T T T T T T T T T T T T T T T4 G G C G G G C G G G G G C G G C G G G5 C C C C C C C C C C C C C C C C C C C6 T T T T T T T T T T T T T T T C T T T Quando peço para juntar as colunas o arquivo ( Parte do arquivo na verdade) passa a ser : D01 D01 D01 D02 D01 D03 D01 D04 D01 D05 D01 D06 D01 D07 D01 D08 D01 D091 1 1 1 1 1 1 0 1 12 2 2 2 2 2 3 2 2 23 2 2 2 2 2 2 2 2 24 1 1 3 1 1 1 3 1 15 2 2 2 2 2 2 2 2 26 2 2 2 2 2 2 2 2 2 D01 D10 D01 F01 D01 F02 D01 F03 D01 F04 D01 F05 D01 F06 D01 F07 D01 F081 1 1 1 1 1 1 1 1 02 2 2 2 2 2 2 2 2 23 2 2 2 2 2 2 2 2 24 1 1 1 3 1 1 3 1 15 2 2 2 2 2 2 2 2 26 2 2 2 2 2 2 3 2 2 MAS, eu preciso juntar as colunas dessa forma: | D01 | com | F01 | , | D02 | com | F01 | , | D03 | com | F01 | , | D04 | com | F01 | , | D05 | com | F01 | , | D06 | com | F01 | , | D07 | com | F01 | , | D08 | com | F01 | , | D09 | com | F01 | , | D10 | com | F01 | , | | D01 | com | F02 | , | D02 | com | F02 | , | D03 | com | F02 | , | D04 | com | F02 | , | D05 | com | F02 | , | D06 | com | F02 | , | D07 | com | F02 | , | D08 | com | F02 | , | D09 | com | F02 | , | D10 | com | F02 | , | | D01 | com | F03 | , | D02 | com | F03 | , | D03 | com | F03 | , | D04 | com | F03 | , | D05 | com | F03 | , | D06 | com | F03 | , | D07 | com | F03 | , | D08 | com | F03 | , | D09 | com | F03 | , | D10 | com | F03 | , | | D01 | com | F04 | , | D02 | com | F04 | , | D03 | com | F04 | , | D04 | com | F04 | , | D05 | com | F04 | , | D06 | com | F04 | , | D07 | com | F04 | , | D08 | com | F04 | , | D09 | com | F04 | , | D10 | com | F04 | , | | D01 | com | F05 | , | D02 | com | F05 | , | D03 | com | F05 | , | D04 | com | F05 | , | D05 | com | F05 | , | D06 | com | F05 | , | D07 | com | F05 | , | D08 | com | F05 | , | D09 | com | F05 | , | D10 | com | F05 | , | | D01 | com | F06 | , | D02 | com | F06 | , | D03 | com | F06 | , | D04 | com | F06 | , | D05 | com | F06 | , | D06 | com | F06 | , | D07 | com | F06 | , | D08 | com | F06 | , | D09 | com | F06 | , | D10 | com | F06 | , | | D01 | com | F07 | , | D02 | com | F07 | , | D03 | com | F07 | , | D04 | com | F07 | , | D05 | com | F07 | , | D06 | com | F07 | , | D07 | com | F07 | , | D08 | com | F07 | , | D09 | com | F07 | , | D10 | com | F07 | , | | D01 | com | F08 | , | D02 | com | F08 | , | D03 | com | F08 | , | D04 | com | F08 | , | D05 | com | F08 | , | D06 | com | F08 | , | D07 | com | F08 | , | D08 | com | F08 | , | D09 | com | F08 | , | D10 | com | F08 | , | | D01 | com | F09 | , | D02 | com | F09 | , | D03 | com | F09 | , | D04 | com | F09 | , | D05 | com | F09 | , | D06 | com | F09 | , | D07 | com | F09 | , | D08 | com | F09 | , | D09 | com | F09 | , | D10 | com | F09 | , | Alguém poderia me ajudar ??? Desde já agradeço a atenção !!!! Att, Giselle

5 6

Krigagem: Obter valores pontuais a partir da sup. interpolada
by Helder Gramacho 20 Out '14

20 Out '14

Bom dia pessoal, Para código reproduzível desta dúvida pode-se utilizar o exemplo disponível na página dos tutoriais do geoR: http://leg.ufpr.br/geoR/tutorials/ctc.R Ao final deste tutorial chegamos ao mapa da distribuição espacial obtido por krigagem. [image: Imagem inline 1] Gostaria de saber se é possível, obter o valor da variável para coordenadas específicas por exemplo na figura acima: x=5; y=8 x=25; y=10 x=15; y=15 Se é possível, como posso implementar no R. Desde já agradeço toda colaboração, *Hélder Gramacho * Recife-PE / *agrohelder(a)gmail.com <agrohelder(a)hotmail.com>*

2 4

skater e deldir
by Hélio Gallo Rocha 20 Out '14

20 Out '14

Caros da lista, Estou tentando implementar o uso de SKATER do pacote spdep. Gostaria de contar com a boa vontade dos listeiros para resolver esta dificuldade. Para iniciar o skater tenho de dividir minha área em poligonos. Para tal estou usando o pacote deldir. Abaixo um exemplo da própria página CRAN. require(deldir) set.seed(42) x <- runif(20) y <- runif(20) d <- deldir(x,y) (ttt <- triang.list(d) ) plot(ttt,border="red",showrect=TRUE) (sss <- tile.list(d)) plot(sss) plot(ttt,add=TRUE,border="blue") O objeto ttt lista os conjuntos de três pontos que formam o poligono Teria uma forma automatizada para gerar estes poligonos no formato shape, pois o CRM do skater que tenho utiliza o comando "readShapePoly" do spdep. Tentei exportar para o formato XLS, mas não foi possível pois o objeto ttt não é um data.frame Agradeço desde já -- Hélio Gallo Rocha IFSULDEMINAS - Câmpus Muzambinho

2 1