Fevereiro 2012 - R-br

Re: [R-br] Criação de função fortran em pacote (base)
by Junior Beleti 24 Fev '12

24 Fev '12

Olá Paulo, o problemas está realmente na chamada a função mvrnorm, pois essa invoca o código objeto eigen.o (função eigen.f). Meu objetivo é dividir o processamento e ao invés de realizar uma chamada, fazer diversas chamadas em paralelo por meio da função mclapply. O problema é que mesmo com a mclapply, quando o código objeto é carregado na memória, apenas uma instancia de tal código pode ser processada por vez. Estaria correta essa minha ultima colocação Benilton? Sobre a chamada mclapply(1:1, function(i) { ... }, mc.cores=4) o elemento (1:1) diz respeito ao número de voltas que a função vai executar e não a quantidade de cores que estou disponibilizando, estou certo? Entendi sua colocação e vou efetuar testes. Assim que tiver respostas lhes comunico. Obrigado pelas dicas e pela atenção. Carlos Beleti.

3 2

Converter arquivo shape file em poligono file
by Antônio Silva 24 Fev '12

24 Fev '12

Prezados! Como faço para converter um arquivo shape file e arquivo um arquivo que o pacote maps possa ler. Eu creio que o nome do arquivo que o pacote maps ler é poligono file. Abs, Antônio Alcirley

1 0

abrir arquivo
by Luiz André Zardo 23 Fev '12

23 Fev '12

Prezados, estou com dificuldades para abrir um arquivo Rfile tento utilizar o comando load("...\\Pesquisa.Rfile") porém a seguinte mensagem aparece: Erro: número mágico do arquivo de restauração inválido (arquivo pode estar corrompido) -- nenhum dado carregadoAlém disso: Warning message:file 'Pesquisa.Rfile' has magic number 'pesqu' Use of save versions prior to 2 is deprecated O que pode estar ocorrendo???? Obrigado Luiz Zardo.

2 2

WebScraping - Como passar informações ao site antes do readLines
by Paulo Nogueira 23 Fev '12

23 Fev '12

Boa tarde R-Br. Mais uma vez peço a vocês uma orientação de caminho para descobrir como solucionar problemas relacionados a web scraping. Hoje consigo ler códigos HTML de ṕáginas atráves do readLines e depois extrair os dados que me interessam utilizando funções como grep, gsub e regexpr. Evolui nisso e tenho conseguido resultados bastante interessantes. O problema é que nem sempre o que é exibido no navegador, está disponível no código fonte da página. Por exemplo, a página inicial do site Webmotors ( http://www.webmotors.com.br/index.html) apresenta um box de pesquisa. Escolhendo, na página, os campos marca = Chevrolet, modelo = Astra, estado = AM e cidade = Manaus e clicando em BUSCAR, encontro esse resultado: http://www.webmotors.com.br/Webmotors/Compra/carrosResultado/carros-resulta… E assim consigo extrair com o R as informações dos veículos Astra a venda em Manaus. Porém, como posso fazer para o R informar ao site que desejo fazer essa busca? Gostaria de automatizar essa parte. Sei que na URL resultado basta trocar os campos como &modelo=1042 ou marca=2, etc... mas também não sei quais são os possíveis valores para esses campos. Ressalto que esse problema é um exemplo. Gostaria de aprender a fazer chamada às funções (acredito que são as .js) da página para construir as páginas com os resultados e então importar com o readLines para guardar os dados de meu interesse. Espero ter sido claro no questionamento. Agradeço qualquer ajuda. Abraços, Paulo Nogueira Starzynski

3 8

FW: [R] HELP ERROR Weibull values must be > 0
by Vanúcia Schumacher 23 Fev '12

23 Fev '12

Vanúcia SchumacherCurso de graduação em meteorologia - UFPELBolsista do Programa de Educação Tutorial - PET From: vanucia-schumacher(a)hotmail.com To: tura(a)centroin.com.br; r-help(a)r-project.org Date: Tue, 21 Feb 2012 22:03:09 +0000 Subject: Re: [R] HELP ERROR Weibull values must be > 0 Tura, meus dados possuem mais de 700 elementos, fica muito inviavel passa-los aqui por summarya questão e que nesses dados eu tenho valores ausentes e valores 0, li algumas sugestões em que deve-se mensurar os valores em 0 para 0.00001,sendo assim ele roda tranquilo , mas retorna um aviso de Nas, mas me retorna os parametros. Não sei se estou fazendo a coisa certa, mas como estou querendo ajustar a distribuição m dados de vento paraprodução de energia eólica, conforme varios artigos que li, essa seria a melhor distribuição. Se puder ajudar fico muito agradecida. library(MASS)> dados<-read.table("inverno.txt",header=FALSE)> vento50<-fitdistr(dados[[1]],densfun="weibull")Mensagens de aviso perdidas:1: In dweibull(x, shape, scale, log) : NaNs produzidos2: In dweibull(x, shape, scale, log) : NaNs produzidos3: In dweibull(x, shape, scale, log) : NaNs produzidos4: In dweibull(x, shape, scale, log) : NaNs produzidos5: In dweibull(x, shape, scale, log) : NaNs produzidos6: In dweibull(x, shape, scale, log) : NaNs produzidos> vento50 shape scale 0.61473436 2.43536064 (0.01895222) (0.13202750) Vanúcia SchumacherCurso de graduação em meteorologia - UFPELBolsista do Programa de Educação Tutorial - PET > Date: Tue, 21 Feb 2012 18:40:47 -0200 > From: tura(a)centroin.com.br > To: vanucia-schumacher(a)hotmail.com > Subject: Re: [R] HELP ERROR Weibull values must be > 0 > > On 02/21/2012 01:48 PM, Vanúcia Schumacher wrote: > > > > GUYS, > > > > I NEED HELP WITH ERROR: > > > > library(MASS) > >> dados<-read.table("mediaRGinverno.txt",header=FALSE) > >> vento50<-fitdistr(dados[[1]],densfun="weibull") > > Erro em fitdistr(dados[[1]], densfun = "weibull") : > > Weibull values must be> 0 > > WHY RETURN THIS ERROR? WHAT CAN I DO? > > Vanúcia > > Vou lhe ajudar dentro do possível mas antes de tudo se seus dados tem > valores menores ou iguais a zero eles não pode ser ajustados numa > distribuição weibull pois esta só permite valores >= 0. > > Que tal você me mandar o resultado do summary(dados) > > Abracos Tura [[alternative HTML version deleted]] ______________________________________________ R-help(a)r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.

3 2

lendo de csv
by Leonard de Assis 23 Fev '12

23 Fev '12

Possuo um csv (+- 350 000 registros) com a seguinte 'xara': "137831","4/12/2006 17:30:47","0","ENTRE 50 e 100","","","-1","1" "137549","4/12/2006 11:13:26","0","ENTRE 50 e 100","","","-1","1" "137936","4/12/2006 19:58:57","0","ENTRE 50 e 100","","","-1","1" "137661","4/12/2006 13:51:43","0","ENTRE 50 e 100","","","-1","1" "138046","4/12/2006 23:35:00","0","ENTRE 100 e 200","","","-1","1" "140275","8/12/2006 17:15:22","0","ENTRE 50 e 100","","","-1","1" "137837","4/12/2006 17:36:06","0","ENTRE 50 e 100","","","-1","1" "138210","5/12/2006 12:02:20","0","ENTRE 50 e 100","","","-1","1" depois de alguns milhares de regístros (105296 regístros depois), ele adquire o seguinte formato: "2144767","6/3/2008","0","ENTRE 1600 e 3200","LJ","SITE","1","1" "2144768","6/3/2008","0","ENTRE 200 e 400","LJ","SITE","1","1" "2144769","6/3/2008","0","ENTRE 100 e 200","LJ","SITE","2","1" "2144770","6/3/2008","0","ENTRE 200 e 400","LJ","SITE","1","1" "2144771","6/3/2008","0","ENTRE 400 e 800","LJ","TVEN","1","1" "2144772","6/3/2008","0","ENTRE 200 e 400","LJ","SITE","1","1" "2144773","6/3/2008","0","ENTRE 100 e 200","LJ","SITE","1","1" "2144774","6/3/2008","0","ENTRE 400 e 800","LJ","SITE","1","1" mas o problema não chegou a ser essa mudança repentina no formato da data (creio que existirá um sério problema referente a isto ainda), o problema é um pouco antes: rodei o seguinte código: reg1 <- read.table( file="CSV/REG0001.csv", header=FALSE, colClasses= c('integer', 'POSIXct', 'integer', 'factor', 'factor', 'factor', 'integer', 'integer'), sep=',', quote='\"', dec=',' ) aí, ao executar, sou surpreendido com o seguinte erro: Error in scan(file, what, nmax, sep, dec, quote, skip, nlines, na.strings, : scan() expected 'an integer', got '"155638"' "155638" é a 1ª linha, 1ª variável do CSV aí vem a pergunta: Como fazer o R entender que precisa ignorar estas '"' ao ler o csv? -- []s Leonard de Assis assis <dot> leonard <at> gmail <dot> com

3 5

porque meu grafico sempre gera o mesmo para arquivos diferentes?
by Vanúcia Schumacher 23 Fev '12

23 Fev '12

Pessoal, Estou com problema num script do R, que quando quero gerar um grafico em png ou mesmo plotar em tela, sempre plota o mesmo grafico, mesmo trocando os parametros e dados, porque está acontecendo isso?? ja tentei trocar de diretório, usar outro terminal e nada, plota sempre o mesmo grafico!!!!! Se alguem puder me orientar, agradeço

2 1

rotulo z da função persp
by Gustavo Marcatti 22 Fev '12

22 Fev '12

Boa noite a todos! Como distanciar o título do eixo z (zlab) dos números do eixo para não ficar sobreposto? Assim como no exemplo abaixo: x <- seq(-10, 10, length= 30) y <- x f <- function(x,y) { r <- sqrt(x^2+y^2); 10 * sin(r)/r } z <- outer(x, y, f) z[is.na(z)] <- 1 z <- z * 100 par(mar=c(0.1,2.2,0.1,1.4),cex=1.6) persp(x, y, z, theta = 30, ticktype = "detailed", xlab="eixo X",ylab="eixo Y", zlab = "eixo Z", col = "gray") Outra dúvida: Tem como fazer um gráfico assim utilizando a função scatterplot3d do pacote scatterplot3d. Se alguém tiver outra sugestão de funções para fazer gráficos 3d, pode se manifestar. Já fiz algumas tentativas com wireframe do lattice, porém tive problemas, principalmente com os títulos dos eixos que não ficam na angulação correta. Desde já agradeço! Gustavo Marcatti Eng. Florestal - UFV

2 2

CRAN Task Views
by Benilton Carvalho 21 Fev '12

21 Fev '12

eye candy para alguns: http://www.maths.lancs.ac.uk/~rowlings/R/TaskViews/ b

1 0

ERRO COM VALORES 0 E NAs
by Vanúcia Schumacher 21 Fev '12

21 Fev '12

PESSOAL, ESTOU TENTANDO GERAR A DISTRIBUIÇÃO DE WEIBULL PARA PLOTAR GRAFICO, MAS ESTOU TENDO UM ERRO QUANDO USO COMANDO: library(MASS)> dados<-read.table("inverno.txt",header=FALSE)> vento50<-fitdistr(dados[[1]],densfun="weibull")Mensagens de aviso perdidas:1: In dweibull(x, shape, scale, log) : NaNs produzidos2: In dweibull(x, shape, scale, log) : NaNs produzidos3: In dweibull(x, shape, scale, log) : NaNs produzidos4: In dweibull(x, shape, scale, log) : NaNs produzidos5: In dweibull(x, shape, scale, log) : NaNs produzidos6: In dweibull(x, shape, scale, log) : NaNs produzidos> vento50 shape scale 0.61473436 2.43536064 (0.01895222) (0.13202750) ALGUEM SABE ME DIZER O QUE SIGNIFICA OU O QUE ESTOU FAZENDO ERRADO? QUANDO TENTO RODAR MEU ARQUIVO DE DADOS COM VALORES AUSENTES OU 0, DA OUTRO ERRO WEIBULL RECISA SER >0 E QUANDO SUBSTITUO ESSES VALORES POR 0.000001 ELE DEVOLVE O ERRO ACIMA. PRECISO DE UM HELP NISSO!!!!!!! BRIGADO

2 1