Fevereiro 2012 - R-br

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

teste nao parametrico Mann-Kendall
by Inajara Welter Macedo 17 Fev '12

17 Fev '12

Ola, Preciso aplicar o teste Mann-Kendall em meus dados para verificar se ha ou nao tendecia significativa de aumento. Estou trabalhado com uma serie de dados de eventos de seca.Fiz a regressao linear da serie, e plotei a reta de tendencia linear dos casos de seca por ano, agora quero verificar se essa tendencia eh realmente significativa, ouseja, se ha um aumento dos casos de seca nos ultimos anos. Minha serie:ano CCS *CCS = casos consecutivos de seca, i.e, numero de dias sem chuva por ano2000 02001 12002 22003 62004 3... ...2009 5 Preciso de um pacotes especifico para utilizar o teste Mann-Kendall?? como faco para aplicar esse teste na minha serie? Obrigada, Inajara

4 4

Re: [R-br] Erro: 'Sigma' is not positive definite
by Junior Beleti 16 Fev '12

16 Fev '12

Obrigado Benilton, vou verificar o determinante de tais matrizes. Paulo, estarei verificando também a construção de tais matrizes. Sobre a decomposição por autovalores/autovetores, sabe me indicar uma função que realize tais procedimentos. Grato pelas orientações.. Carlos Beleti.

2 1