- R-br - DINF UFPR

lendo de csv
by Leonard de Assis 22 Fev '12

22 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 22 Fev '12

22 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

Como ajustar algumas propriedades do segundo eixo Y?
by Cleber N.Borges 20 Fev '12

20 Fev '12

Olá a todos, Gostaria de ajustar um segundo eixo Y (VAR 01), de acordo com o primeiro eixo Y (VAR 02). o que eu não consegui: 1 - ajustar o 'zero' do primeiro com o 'zero' do segundo eixo Y 2 - colocar o nome do segundo eixo (um ylab='VAR 02') Imagino que vendo o que eu quero, fica mais fácil entender, então segue um código exemplo junto. Agradeço antecipadamente por qualquer ajuda. Cleber ############# set.seed(1234) x1 <- abs(rnorm(10)) x2 <- abs(rnorm(10)) + 100 par(mar=c(5,4,4,5)) bp <- barplot( x1, xlim=c(0,12), ylim=c(0,3), xlab='', ylab='VAR 01', xaxt='n') par( new=T ) plot( bp[,1], x2, t='b', xlim=c(0,12), xlab='', ylab='', xaxt='n', yaxt='n', bty='n', ylim=c(100,103) ) axis(4, ylab='VAR 02') grid() axis(1, at=bp[,1], labels=letters[1:10])

3 3

duvidas comandos R distribuição weibull
by Vanúcia Schumacher 20 Fev '12

20 Fev '12

Pessoal, Estou começando a utilizar o R agora e estou com duvidas em alguns argumentos desse comando para plotar grafico: hist(Vento[[1]], seq(0, 18, 0.5), prob=TRUE, xlab="Velocidade (m/s)",ylab="Densidade", main="50 m") curve(dweibull(x, shape=2.223, scale=5.534), 0,18,add=T, col='red') o que significa ? seq(0, 18, 0.5), 0,18,add=T Quando insiro o comando para saber os parametros 50<-fitdistr(dados[[1]],desfun="weibull") retorna um erro Weibull values must be > 0 Alguem poderia me ajudar?

2 1

parametros weibull
by Vanúcia Schumacher 20 Fev '12

20 Fev '12

Olá pessoal, Gostaria de ajuda para calcular os parâmetros de weibull, pois as linhas de comando que estou utilizando são para calcular pelo método dos momentos, porém , está me retornando valores negativos, se alguém que entenda melhor poder ajudar agradeço.Talvez tenha uma maneira melhor de calcular dados<-read.table("mediaPOAinverno.txt",header=FALSE) med.weibull<-mean(dados) ## sample mean var.weibull<-var(dados) ## sample variance l.est<-med.weibull/var.weibull ## lambda estimate (corresponds to rate) a.est<-((med.weibull)^2)/var.weibull ## alfa estimate Valeu Vanúcia SchumacherCurso de graduação em meteorologia - UFPELBolsista do Programa de Educação Tutorial - PET

3 2

erro parametros de weibull
by Vanúcia Schumacher 19 Fev '12

19 Fev '12

Pessoal, estou tentando gerar os parametros de weibull no R, mas não estou conseguindoConforme orientaçao do Bruno obtive um erro, e não estou entendendo.alguem pode ajudar? library(MASS)> dados<-read.table("teste.txt",header=FALSE)> is.na(dados) <- dados == -999.9> vento50<-fitdistr(dados[[1]],densfun="weibull")Mensagens de aviso perdidas:In dweibull(x, shape, scale, log) : NaNs produzidos> is.na(dados) <- dados == -999> vento50<-fitdistr(dados[[1]],densfun="weibull")Mensagens de aviso perdidas:In dweibull(x, shape, scale, log) : NaNs produzidos> str(vento50)List of 5 $ estimate: Named num [1:2] 1.39 2.11 ..- attr(*, "names")= chr [1:2] "shape" "scale" $ sd : Named num [1:2] 0.191 0.292 ..- attr(*, "names")= chr [1:2] "shape" "scale" $ vcov : num [1:2, 1:2] 0.0366 0.0184 0.0184 0.0854 ..- attr(*, "dimnames")=List of 2 .. ..$ : chr [1:2] "shape" "scale" .. ..$ : chr [1:2] "shape" "scale" $ loglik : num -46.9 $ n : int 30 - attr(*, "class")= chr "fitdistr">

2 1

Re: [R-br] teste nao parametrico Mann-Kendall
by Inajara Welter Macedo 18 Fev '12

18 Fev '12

Ola, Para o tipo de resultado que quero obter, de apenas verificar se a linha de tendencia encontrada atraves de regressao lineareh significativa esse teste MAnn-Kendall eh o mais adequado, por ser muito utilizado na bibliografia para encontrar resultados muito semelhantes aos meus. Respondendo a questao da necessidade de um vetor com os dias sem chuva por ano, tenho um vetor com o numero de dias sem chuva por ano.Mas o problema nao eh esse, pois seria facil de encontrar outro tipo de vetor, ou um data.frame. Preciso mesmo entender melhor como funciona esse teste no R.Como aplica-lo a minha serie. Atte,Inajara

1 0