Novembro 2012 - R-br

Criar data.frame utilizando a PNAD
by Alexandre Loures 24 Nov '12

24 Nov '12

Boa noite pessoal! Sou iniciante no R e estou tentando criar um data.frame a partir dos dados da PNAD. No Windows obtive sucesso, porém, no Ubuntu 12.10 está ocorrendo o erro abaixo. library(IBGEPesq) data(dicPNAD2009) microdados <- file.path("/home/alexandre","Dados","PES2011.txt") dados<-le.pesquisa(dicionario=dicpes2009,pathname.in=microdados,pathname.out="dados.Rdata",codigos=c("V4803","V4704","V4805","V4707","V4729"), nomes=c("ANOEST","CONDATV","CONDOCP","HTRAB","PESO")) Erro em waitbar() : não foi possível encontrar a função "windows" Obrigado! -- *Alexandre Rodrigues Loures* Mestrando em Economia do Desenvolvimento (Bolsista Capes) Faculdade de Administração, Contabilidade e Economia - FACE Pontifícia Universidade Católica do Rio Grande do Sul - PUCRS Tel.: +55 32 9102-0690 (MG) ou +55 51 8241-6579 (RS) e-Mail: alexandre(a)loures.pro.br Skype: alexandre_loures Site: http://loures.pro.br

1 0

RES: Simulação de variáveis correlacionadas usando Cópulas
by Gabriel Bruno de Lemos 23 Nov '12

23 Nov '12

Mauro, estranho pq aqui em meu computar esta função está rodando normalmente. Segue o código desta função e também o output fornecido pela rotina que passei... > rMvdc function (n, mvdc) { dim <- mvdc@copula@dimension u <- rCopula(n, mvdc@copula) x <- u for (i in 1:dim) { qdf.expr <- asCall(paste0("q", mvdc@margins[i]), mvdc@paramMargins[[i]]) x[, i] <- eval(qdf.expr, list(x = u[, i])) } x } > require(copula) > > myCop.clayton <- archmCopula(family = "clayton", dim = 3, param = 2) > > myMvd <- mvdc(copula = myCop.clayton, margins = c("exp", "norm", + "chisq"), paramMargins = list(list(rate= 2), list(mean = 0, + sd = 1), list(df=5))) > > > myCop.t <- ellipCopula(family = "t", dim = 3, dispstr = "un", + param = c(0.8, 0.5,0.2), df = 8) > > myMvd <- mvdc(copula = myCop.t, margins = c("exp", "norm","chisq"), + paramMargins = list(list(rate= 2), list(mean = 0,sd = 1), list(df=5))) > > x.samp <- rMvdc(10000,myMvd) > cor(x.samp) [,1] [,2] [,3] [1,] 1.0000000 0.7194154 0.4536219 [2,] 0.7194154 1.0000000 0.1652259 [3,] 0.4536219 0.1652259 1.0000000 De: Mauro Sznelwar [mailto:sznelwar@uol.com.br] Enviada em: quarta-feira, 21 de novembro de 2012 23:08 Para: gblemos(a)usp.br Assunto: Re: [R-br] Simulação de variáveis correlacionadas usando Cópulas Tentei rodar o seu script e diz que este comando rMvdc não existe! Boa noite a todos! Já tem algum tempo que estou tentando aprender um pouco mais sobre cópulas... e confesso que não está sendo a coisa mais simples do mundo. Sinto falta de alguma aplicação (rotina no R) mais simples mostrando uma utilização aplicada deste tipo de função (principalmente para casos acima de 2 dimensões). Utilizo-a no meu trabalho para precificar portfólios... mas em um software caixa preta... Hoje consegui rodar um pequeno exemplo de uma das possíveis utilidades destas cópulas. Gerar variáveis aleatórias correlacionadas... com as marginais possuindo QUALQUER função de distribuição de probabilidade (não ficando restrito ao caso da normal multivariada). Parece que está funcionando legal... se alguém tiver mais alguns exemplos ou possíveis correções na rotina abaixo será bem vinda! Abs require(copula) ## Primeiro define-se a cópula a ser utilizada t neste caso. O comando dispstr=un especifica que desejo informar 3 correlações diferentes na matriz de correlações. myCop.t <- ellipCopula(family = "t", dim = 3, dispstr = "un", param = c(0.8, 0.5,0.2), df = 8) ## Definida a cópula, especifico quais marginais desejo que minha função possua. Neste caso defini exponencial, normal e qui-quadrado. myMvd <- mvdc(copula = myCop.t, margins = c("exp", "norm","chisq"), paramMargins = list(list(rate= 2), list(mean = 0,sd = 1), list(df=5))) x.samp <- rMvdc(10000,myMvd) cor(x.samp) par(mfrow=c(1,3)) hist(x.samp[,1]) hist(x.samp[,2]) hist(x.samp[,3]) Gabriel Bruno de Lemos Mestrando em Estatística e Experimentação Agronômica Esalq / Usp skype: gb_lemos msn: gb_lemos(a)hotmail.com (19) 8212-6999 .

2 1

duvida no pacote ff (read.csv.ffdf)
by Vinicius Brito Rocha 23 Nov '12

23 Nov '12

Pessoal, estou tentando ler um txt com teste=read.table.ffdf(filename, first.rows=2, next.rows=10, VERBOSE=TRUE,colClasses=c("character", "character", "character","character", "character","character","character","character", "character","POSIXct", "POSIXct","character") ) o filename possui o caminho com o nome do arquivo Estou a seguinte mensagem de erro: Erro em read.table.ffdf(filename, first.rows = 2, next.rows = 10, VERBOSE = TRUE, : only ffdf objects can be used for appending (and skipping the first.row chunk) Alguém sabe o que pode ser? -- *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.

1 0

Re: [R-br] Movimento Valorização do profissional de Estatística, já!
by Richard Santos 23 Nov '12

23 Nov '12

E esclarecer que: um software de geoprocessamento não tem estatística, tem estatística básica e conglomerados. Um bioestatístico não é estatístico, é bioestatístico. Um econometrista não é estatístico, é econometrista. Um economista que sabe modelos matemáticos e econômicos, e que de modelos estatísticos só teve uma disciplina de regressão e Séries Temporais, só sabe regressão e Séries Temporais, pelo menos em que se formou. Pela grade de formação dele, não sabe mais, não sabe inferência, amostra, etc. Quem sabe amostra por cotas é porque não sabe amostra probabilística, e eu não tenho coragem de discutir a validade da aleatoriedade, sem fazer um doutorado antes sobre isso. Um estatístico necessariamente mexe com bases, e muitos com BIG DATA, outro profissional que trabalha com BIG DATA necessariamente não é um estatístico. Se as pessoas, a ciência, valoriza tanto a estatística, não vamos deixar acontecer um sequestro mal feito, e valorizá-la também. Abçs. Richard Em 23/11/12, Richard Santos<jamesrichardsantos(a)gmail.com> escreveu: > Caros, > > Devido aos comentários recentes e conclusões, que tal, quem concordar, > inclusive tendo em vista o ANO INTERNACIONAL DA ESTATÍSTICA, em 2013, > lançarmos o movimento "valorização do profissional de Estatística já"? > > Esse movimento contemplaria: > > - Só repetindo, de novo, acolher e valorizar os outros profissionais > que se especializam em alguma técnica estatística, mas não os > confundir com estatístico, e, mesmo assim, corrigi-los quando acharmos > preciso. > - Parabenizar as instituições que fazem concurso para Estatístico e > vagas privadas que contemplam o estatístico, quando o trabalho for em > sua maioria fazer estatística; > - Motivar o trabalho em equipe, economista, com estatístico, com > engenheiro, com TI, com cientista social, etc; > - Motivar os estatísticos a fazerem o concurso de Analista Geral > existentes atualmente no IPEA e no IBGE, e para nos fortalecer lá > dentro, fazendo frente aos economistas; > - Definir estatístico, na regulamentação: graduação, outro > profissional que tenha pós em estatística? Matemática com pós em > estatística? Economista com pós em estatística? Manter a definição de > estatístico e técnico em estatística; > - conscientizar nos cursos de Estatística básica, nos mini-cursos, etc; > - Fazer marketing à sociedade, à sociedade civil, aos empregadores e > intelectuais brasileiros com relação à Estatística e ao profissional > de estatística; > - Só repetindo, de novo, acolher e valorizar os outros profissionais > que se especializam em alguma técnica estatística, mas não os > confundir com estatístico, e, mesmo assim, corrigi-los quando acharmos > preciso. > > Porque se um dia no Brasil a estatística virar uma coisa que qualquer > um pode fazer ou se virar totalmente "promíscua", eu só recomendo > livros, ensino quais botões apertar, e não mudo de país não, eu vou > fazer outra coisa, aqui no Brasil mesmo, rs. > > Abraços, > Att. > Richard > Estatístico - CONRE 1 >

2 2

Concursos na Ufal
by Alejandro C. Frery 23 Nov '12

23 Nov '12

Boa tarde, Há vários editais para concursos na Ufal. Talvez nesta lista haja candidatos interessados. Ver http://www.copeve.ufal.br/index.php?opcao=concurso&idConcurso=153 Alejandro -- Alejandro C. Frery Maceió, AL - Brazil http://sites.google.com/site/acfrery/ http://www.researcherid.com/rid/A-8855-2008

1 0

R ler bases SAS
by Richard Santos 23 Nov '12

23 Nov '12

Um colega disse que "Especificamente com o formato sas7bdat tem o pacote sas7bdat. Eu já li bases do sas com esse pacote no R open source. Funciona! http://cran.r-project.org/web/packages/sas7bdat/index.html" Se ler o Censo Demográfico, que nem é um BIGDATA, eu já fico feliz. Rsrs. Lê? Rsrsrs. Abraços, Richard

1 0

Seminários dos Usuários de R no Brasil
by Ary Casado 23 Nov '12

23 Nov '12

Pelo que vejo na lista de usuários temos representantes de vários estados brasileiros e de várias instituições. Acredito que seja a hora de organizar os usuários e desenvolvedores em torno deste foco de interesse: a utilização e desenvolvimento do R tanto no meio científico como no âmbito da aplicabilidade prática. Não sei quantos interessados existem no Recife, mas acredito que já exista o interesse de iniciar encontros e seminários com o objetivo maior gerar um grupo de usuários e desenvolvedores. Arykerne Nascimento Engenheiro Químico - UFPE Especialista em Informática - UFPE Mestrando em Ciências Nucleares - UFPE

2 1

Movimento Valorização do profissional de Estatística, já!
by Richard Santos 23 Nov '12

23 Nov '12

Caros, Devido aos comentários recentes e conclusões, que tal, quem concordar, inclusive tendo em vista o ANO INTERNACIONAL DA ESTATÍSTICA, em 2013, lançarmos o movimento "valorização do profissional de Estatística já"? Esse movimento contemplaria: - Só repetindo, de novo, acolher e valorizar os outros profissionais que se especializam em alguma técnica estatística, mas não os confundir com estatístico, e, mesmo assim, corrigi-los quando acharmos preciso. - Parabenizar as instituições que fazem concurso para Estatístico e vagas privadas que contemplam o estatístico, quando o trabalho for em sua maioria fazer estatística; - Motivar o trabalho em equipe, economista, com estatístico, com engenheiro, com TI, com cientista social, etc; - Motivar os estatísticos a fazerem o concurso de Analista Geral existentes atualmente no IPEA e no IBGE, e para nos fortalecer lá dentro, fazendo frente aos economistas; - Definir estatístico, na regulamentação: graduação, outro profissional que tenha pós em estatística? Matemática com pós em estatística? Economista com pós em estatística? Manter a definição de estatístico e técnico em estatística; - conscientizar nos cursos de Estatística básica, nos mini-cursos, etc; - Fazer marketing à sociedade, à sociedade civil, aos empregadores e intelectuais brasileiros com relação à Estatística e ao profissional de estatística; - Só repetindo, de novo, acolher e valorizar os outros profissionais que se especializam em alguma técnica estatística, mas não os confundir com estatístico, e, mesmo assim, corrigi-los quando acharmos preciso. Porque se um dia no Brasil a estatística virar uma coisa que qualquer um pode fazer ou se virar totalmente "promíscua", eu só recomendo livros, ensino quais botões apertar, e não mudo de país não, eu vou fazer outra coisa, aqui no Brasil mesmo, rs. Abraços, Att. Richard Estatístico - CONRE 1

1 0

Comando install.packages(pacote)
by Richard Santos 23 Nov '12

23 Nov '12

Pessoal, Eu dei o comando install.packages('pacote'), e o R me deu as mensagens: --- Please select a CRAN mirror for use in this session --- Aviso: unable to access index for repository http://cran.fiocruz.br/src/contrib Aviso: unable to access index for repository http://www.stats.ox.ac.uk/pub/RWin/src/contrib Aviso: unable to access index for repository http://cran.fiocruz.br/bin/windows/contrib/2.15 Aviso: unable to access index for repository http://www.stats.ox.ac.uk/pub/RWin/bin/windows/contrib/2.15 Mensagens de aviso perdidas: 1: In open.connection(con, "r") : não foi possível conectar a 'cran.r-project.org' na porta 80. 2: package ‘pacote’ is not available (for R version 2.15.2) What the hell? Em 23/11/12, Richard Santos<jamesrichardsantos(a)gmail.com> escreveu: > Excelentes dicas, Flávio, obrigado pela atenção. Gostei do que foi > falado sobre a linguagem, e sobre a instalação de pacotes local. Vou > testar, realmente, para aprender o R aos poucos, também. Só uma > confirmação: essas opções do Revolution e dataframe, etc, são para ler > bases grandes também, né? > Grato. > Richard > > Em 22/11/12, Flavio Barros<flaviomargarito(a)gmail.com> escreveu: >> - sempre prefiro a programação, mas a programação do R é mais extensa >> do que o SAS e o SPSS. Por exemplo, para um comando de igualdade no >> SAS ou no SPSS se digita X4=1 e no R se digita dados$x4 != 1; >> >> Os programas em R são bem compactos na verdade (em termos de sintaxe), >> mesmo se comparados a outras linguagens como python, ruby ou mesmo o >> matlab. Vide ( >> http://blog.revolutionanalytics.com/2012/11/which-programming-language-is-t… >> ) >> >> Quando a usar X4 = 1 ou dados$x4 = 1 no R é possível das duas formas, >> basta >> 'fixar' o conjunto de dados antes com o attach(dados). Particularmente >> eu >> prefiro escrever x <- 4 (padrão do google também! >> http://google-styleguide.googlecode.com/svn/trunk/google-r-style.html#assig… >> que diga-se de passagem é um usuário ativo do R.) >> >> - para grandes bases, tem-se algumas opções como Revolution, mas mesmo >> o Revolution parece que apresenta um problema que encontrei similar ao >> Minitab, que as variáveis envolvidas tem que ter o mesmo tamanho ou >> formato. Eu não queria encontrar nenhuma limitação quanto a ler >> grandes bases; >> >> Na verdade isso não é necessário. No R existem algumas estruturas de >> dados >> como o matrix, o data.frame e o list. De fato, se você quiser ler os seus >> dados em uma matriz, todos os dados tem que ser do mesmo tipo, mas em um >> data.frame não! O formato XDF da revolution também não impõe esse tipo de >> limitação. Caso as variáveis tenham tamanhos diferentes ainda é possível >> ler como um list de forma que nem o tamanho das variáveis precisa ser >> igual. Mande um e-mail para o suporte da revolution que acredito que o >> problema foi sua implementação. Eles respondem rápido. Já recebei uma >> resposta no domingo depois de 20 minutos! >> >> - para cada procedimento que se quer fazer, pode ser necessário >> instalar um pacote diferente, eu tentei pedir para o administrador do >> meu trabalho que uma pasta minha fiique livre para eu instalar >> pacotes, mas não sei se isso vai ser feito. >> >> Acho que já foi falado aqui, mas não custa lembrar: não é necessário >> privilégio de administrador para instalar pacotes. Caso você não seja >> administrador basta dar um install.packages('pacote') que ele vai criar >> um >> biblioteca local de pacotes. Esse não é um problema. >> >> No mais, é normal para aprender qualquer software ter que gastar >> tempo, mas as questões que citei acima se quiserem podem considerar >> como gosto pessoal mesmo. Mas acho importante eu citar essas questões >> para os empregadores saberem, pelo menos. Pois creio que essas >> questões não permitem ainda dizer para todos que o R é uma opção fácil >> de programa livre, porque não é tão fácil assim. Mas ele é ótimo, não >> que não seja ótimo. Apenas acho que ainda há argumentos para que uma >> empresa compre tanto o SPSS, quando achar que precisa, quanto o SAS, >> quando for necessário. >> >> Na minha opinião existem casos e casos para migrar do SAS ou SPSS para o >> R. >> Uma questão importante é o treinamento: se a empresa já tem profissionais >> bem treinados usando SAS, muito código pronto que seria difícil migrar a >> melhor solução é continuar a usar essas ferramentas. Agora no caso de um >> consultor, uma startup ou uma empresa pública que vai optar por um pacote >> estatístico acredito que o R é sim uma melhor opção. >> >> O R é muito utilizado no Google por exemplo. Se é bom para o Google dá >> para >> usar em qualquer empresa. >> >> Não advogo todo mundo usar o R, mas a linguagem tem muitas >> características >> interessantes, é bem intuitiva para quem programa em outras linguagens >> (python, C, Java por ex) e tem custo muito menor que o SAS ou SPSS. E é >> open source. >> >> No seu caso Richard que acredito que deve ter uma boa experiência em SAS, >> trocar pode não ser a melhor opção para você. Se você tem acesso a uma >> licença do SAS ou do SPSS e tem experiência com essas ferramentas não >> vejo >> porque trocar. No entanto é interessante ir aprendendo o R em paralelo. >> >> Flavio Barros >> [image: Facebook] >> <http://www.facebook.com/flavio.barros.1650?ref=tn_tnmn> [image: >> LinkedIn] <http://www.linkedin.com/profile/view?id=61839390&trk=tab_pro> >> Contact me: [image: Google Talk] flaviomargarito(a)gmail.com >> “"Sing like no one's listening, love like you've never been hurt, dance >> like nobody's watching, and live like its heaven on earth." - Mark >> Twain<http://www.quotesdaddy.com/quote/1040448/mark-twain/sing-like-no-ones-liste…> >> ” Get this email app! >> <http://www.wisestamp.com/apps/quotes?utm_source=extension&utm_medium=email&…> >> >> Create your free signature: >> <http://r1.wisestamp.com/r/landing?promo=33&dest=http%3A%2F%2Fwww.wisestamp.…> >> CLICK >> HERE!<http://r1.wisestamp.com/r/landing?promo=33&dest=http%3A%2F%2Fwww.wisestamp.…> >> >> >> >> 2012/11/22 Richard Santos <jamesrichardsantos(a)gmail.com> >> >>> ** >>> >>> >>> Caros, >>> >>> Apenas gostaria de informar que eu, como estatístico, particularmente, >>> não vou divulgar para as empresas e instituições que o R é uma opção >>> com os mesmos níveis de facilidade do que o SAS e o SPSS. Eu não tenho >>> dúvidas de que alguém que fique expert no R consiga fazer tudo e ainda >>> mais rápido do que o SAS. >>> >>> Apenas que eu pessoalmente não vou optar pelo R agora, pelo menos por >>> enquanto, ainda que eu tenha que trazer notebook para o trabalho, por >>> alguns motivos que não tive adesão por questão de gosto pessoal: >>> >>> - sempre prefiro a programação, mas a programação do R é mais extensa >>> do que o SAS e o SPSS. Por exemplo, para um comando de igualdade no >>> SAS ou no SPSS se digita X4=1 e no R se digita dados$x4 != 1; >>> >>> - para grandes bases, tem-se algumas opções como Revolution, mas mesmo >>> o Revolution parece que apresenta um problema que encontrei similar ao >>> Minitab, que as variáveis envolvidas tem que ter o mesmo tamanho ou >>> formato. Eu não queria encontrar nenhuma limitação quanto a ler >>> grandes bases; >>> >>> - para cada procedimento que se quer fazer, pode ser necessário >>> instalar um pacote diferente, eu tentei pedir para o administrador do >>> meu trabalho que uma pasta minha fiique livre para eu instalar >>> pacotes, mas não sei se isso vai ser feito. >>> >>> No mais, é normal para aprender qualquer software ter que gastar >>> tempo, mas as questões que citei acima se quiserem podem considerar >>> como gosto pessoal mesmo. Mas acho importante eu citar essas questões >>> para os empregadores saberem, pelo menos. Pois creio que essas >>> questões não permitem ainda dizer para todos que o R é uma opção fácil >>> de programa livre, porque não é tão fácil assim. Mas ele é ótimo, não >>> que não seja ótimo. Apenas acho que ainda há argumentos para que uma >>> empresa compre tanto o SPSS, quando achar que precisa, quanto o SAS, >>> quando for necessário. >>> >>> Obrigado. >>> >>> Abraços, >>> Att. >>> Richard >>> >>> >> >

5 13

pacote fora do ar
by Rodrigo Plei 23 Nov '12

23 Nov '12

Prezados, O pacote "PBSmapping" foi retirado recentemente do CRAN. Alguém teria o arquivo zip do mesmo para compartilhar? Agradeço desde já. Rodrigo -- =8-->=8-->=8-->=8-->=8-->=8-->=8-->=8-->=8-->=8-->=8-->=8--> Rodrigo Silvestre Martins, PhD Bolsista FAPESP Pós-Doutorado Laboratório de Ecossistemas Pesqueiros (LabPesq) Universidade de São Paulo, Instituto Oceanográfico Praça do Oceanográfico, 191. Cidade Universitaria (sala 107-A/B) Butantã - São Paulo/SP, Brasil 05508-900 Tel: +55 11 3091 6549 Email: rodrigo.plei(a)gmail.com ; ocersm(a)lycos.com; rsmartins(a)usp.br CV Lattes: http://lattes.cnpq.br/5350064124902777 Google Scholar: http://scholar.google.com.br/citations?user=zX_EezEAAAAJ&hl=pt-BR

2 2