- R-br - DINF UFPR

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

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

21 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

instalar versão antiga do R no Windows
by Rodrigo Plei 18 Fev '12

18 Fev '12

Prezados, Preciso instalar a versão 2.9.2 do R no meu computador (Windows &) para rodar um pacote enviado por um colega. Como faço para instalar essa versão antiga (que eu já baixei do CRAN). Agradeço desde logo a ajuda. 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

2 4

correção
by Fernando Antonio de souza 17 Fev '12

17 Fev '12

Caros amigos apenas uma correção. no e-mail anterior coloquei os meus dados na função deput() para torná-lo reproduzível entretanto a função é dput() por isto estou enviando novamente. Peço desculpas pela distração estou tentando rodar um curva de crescimento (Modelo gompertz) para predizer o conteúdo de cálcio (UtGest1f$UTGEST) no útero gestante em função da idade gestacional (UtGest1f$Gest). Entretanto ao rodar a função nls ela anda me retornando a seguinte mensagem de erro: " Erro em assign(i, temp, envir = env) : tentativa de usar um nome de variável com comprimento zero Não sei onde corrigir este problema. Alguém pode me orientar. Um grande abraço a todos #conteúdo de calcio no útero gestante gestação simples UtGest1f<-subset(tabcalcio, Manejo==2 & Fetos==1,na.rm=T) str(UtGest1f) MedUtGest1NR<-tapply(UtGest1f$UTGEST,UtGest1f$Gest,mean) MedUtGest1NR attach(UtGest1f) modelo<-nls(UTGEST~a-b*exp(-c * Gest),start=list(a=52,b=48,c=-0,009143)) plot(MedUtGest1NR,main="relacionamento calcio utero gestante e idade gestacional,gestação multipla",xlab="idade gestacional",ylab="conteúdo de cálcio (g)",na.rm=T) lines(MedUtGest1NR) modelo<-nls(UTGEST~a-b*exp(-c * Gest),start=list(a=52,b=48,c=-0,009143)) dput(UtGest1f)

1 0

dúvida função nls
by Fernando Antonio de souza 17 Fev '12

17 Fev '12

Caros amigos estou tentando rodar um curva de crescimento (Modelo gompertz) para predizer o conteúdo de cálcio (UtGest1f$UTGEST) no útero gestante em função da idade gestacional (UtGest1f$Gest). Entretanto ao rodar a função nls ela anda me retornando a seguinte mensagem de erro: " Erro em assign(i, temp, envir = env) : tentativa de usar um nome de variável com comprimento zero Não sei onde corrigir este problema. Alguém pode me orientar. Um grande abraço a todos #conteúdo de calcio no útero gestante gestação simples UtGest1f<-subset(tabcalcio,Manejo==2 & Fetos==1,na.rm=T) str(UtGest1f) MedUtGest1NR<-tapply(UtGest1f$UTGEST,UtGest1f$Gest,mean) MedUtGest1NR modelo<-nls(UTGEST~a-b*exp(-c * Gest),start=list(a=52,b=48,c=-0,009143)) plot(MedUtGest1NR,main="relacionamento calcio utero gestante e idade gestacional,gestação multipla",xlab="idade gestacional",ylab="conteúdo de cálcio (g)",na.rm=T) lines(MedUtGest1NR) modelo<-nls(UTGEST~a-b*exp(-c * Gest),start=list(a=52,b=48,c=-0,009143)) deput(UtGest1f)

1 0