- R-br - DINF UFPR

Nome de Função
by David Feitosa 23 Out '14

23 Out '14

Caríssimos, boa noite. Existe modo de se obter o nome de uma função em tempo de execução? Explicando melhor: eu tenho uma função da forma minhaFuncao <-function(){ #codigo } e uma variável que recebe essa função: f <- minhaFuncao Eu gostaria de obter o nome dela para ser usado em um gráfico Em Java, eu sei que poderia usar a API de Reflection para, em tempo de execução, analisar a classe, seus métodos etc e obter o nome para uso posterior. Há como fazer o equivalente em R? No meu caso, bastaria saber qual o nome da função que a variável "f" carrega. Grato. Atenciosamente, David Feitosa (\_(\ (=°;°) (("")("")

4 8

Teste de estimadores de regressão
by Adriano Borges Costa 23 Out '14

23 Out '14

Caros, Estou fazendo um exercício teórico de regressão para estimar o efeito de um programa a partir da base de dados disponível no código baixo. library(foreign) lalonde <- read.dta(" http://dl.dropboxusercontent.com/u/7453595/BD_LAB/lalonde.dta") #carrega a base de dados str(lalonde) # variáveis # As variáveis re74, re75 e re78 contêm informações sobre o rendimento para os anos de 1974, 1975 e 1978, respectivamente. # A variável de tratamento t indica se o indivíduo participou do programa (t=1) ou se não participou do programa (t=0). # As demais variáveis contém características dos trabalhadores Trata-se de um experimento aleatório e quero mostrar que há aleatoriedade no experimento. Assim, estou rodando a seguinte regressão: regressao1 <- lm(re78 ~ t, data=lalonde) summary(regressao1) E em seguida incluindo os controles em uma outra regressão e gostaria de mostrar que o estimador de t (tratamento) varia pouco (cerca de U$ 100,00). regressao2 <- lm(re78 ~ t+age+educ+nodegree+black+hisp+married, data=lalonde) summary(regressao2) Assim, gostaria de fazer um teste para verificar se o estimador de t na regressão 1 é estatisticamente diferente da mesmo estimador na regressão 2. Está certo eu fazer da seguinte forma? install.packages("car") library(car) linearHypothesis(regressao2,"t=1794.3", test = "F") Desde já obrigado pelas ajudas. At. -- Adriano Borges Costa

3 3

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 22 Out '14

22 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 22 Out '14

22 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 21 Out '14

21 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