- R-br - DINF UFPR

Ajuda com pacote "seas"
by Inajara Welter Macedo 06 Jan '12

06 Jan '12

Ola, Estou trabalhando com dados de precipitacao, de 10 anos, e quero obter resultados de periodos de seca. Preciso contar quantos dias consecutivos sem chuva houve, ao todo, e por ano e depois por mes. Meus dados eram uma matriz de 4 colunas, chamado "mod1", cada coluna corresponde ao ano, mes, dia, e precipitacao, mas transformei para data.frame, para poder utilizar o pacote "seas": as.data.frame(mod1) Tentei usar a funcao "interarrival": CCS<-interarrival (mod1, var = "Precip", p.cut = 0, inv = TRUE) # CCS= casos consecutivos secos Mas, ocorre um erro: Erro em seas.df.check(x, orig, var) : a ‘date’ column must exist in ‘mod1’ Nao sei como colocar essa coluna de data que pede, porque existem as colunas de ano , mes, e dia. Como consigo converter essas colunas em uma unica coluna de data?? Obrigada, Inajara - UFPEL/RS

5 4

Ajuda para nomear variáveis para uma função repetida
by Maria Papa 05 Jan '12

05 Jan '12

Olá pessoal, boa tarde, Gostaria da ajuda de vcs para nomear variáveis. Eu tenho este "for" que é aplicado para as três coordenadas de um ponto (XYZ). O procedimento é idêntico para os três, por isso, eu gostaria de saber se tem alguma forma de usar este mesmo conjunto mundando apenas o nome das coordenadas, no exemplo abaixo, eu faço com X e tenho os resutlados para X, depois e preciso para Y e Z. Eu fiz algumas tentativas com "for," mas não deu certo. for (j in 1:n) { X_e = rep(X[j], length(m))+ (runif(m, -U,U)) funX = function(X_e){mean(X_e)} X_boot = one.boot(X_e,funX,R) X_boot_ic= boot.ci(X_boot,conf = c(0.95),type = c("norm")) ic_inf_X [j]= X_boot_ic$normal[2] ic_sup_X [j] = X_boot_ic$normal[3] X_b[j]= X_boot$t0 } os resultados que eu preciso no final e portanto devem ser trocados o nome das coordenadas são os intervalos de confiança ic_inf_X, ic_sup_X e X_b Desde já agradeço a todos. Célia

3 5

Re: [R-br] Ajuda com pacote "seas"
by Ivan Bezerra Allaman 05 Jan '12

05 Jan '12

A mensagem de erro é clara. Uma coluna com datas deve existir em mod1. Para uma ajuda mais efetiva, nos envie parte dos seus dados e um código reproduzível. Um ?interarrival ajuda a entender melhor o uso da função. (S,f,P) Allaman \begin{signature} <<>>= Prof. Dr. Ivan Bezerra Allaman Universidade Estadual de Santa Cruz Departamento de Ciências Exatas e Tecnológicas Ilhéus/BA - Brasil Fone: +55 73 3680-5076 E-mail: ivanalaman@yahoo.com.br/ivanalaman@gmail.com @ \end{signature}

1 0

SOMAR VARIAVEIS
by Luiz André Zardo 05 Jan '12

05 Jan '12

PREZADOS Preciso criar uma variavel Y que é a soma de outras 3 variaveis (X1+X2+X3), porém quando uma delas é algum NA não consigo fazer a soma das outras duas. como resolvo isso, como nos casos 2 e 3 abaixo. caso *X1* *X2* *X3* *Y* 1 1 2 1 4 2 4 2 NA NA 3 NA 2 1 NA 4 2 0 3 5

2 1

Leitura de arquivos
by Emerson Cotta Bodevan 04 Jan '12

04 Jan '12

Prezados, tenho um arquivo .txt com uma sequência enorme de números. Por ex.: 0120304405560085673947684658495604932570932847509836536793586593420574092398564230562439065423095624359842350234567230945 7340127459042759834265093248560147237561092436509864230543609586723096756759867243670934675.... . . . ...87350274523475982734596239467358927609650972345792034759237459623767235967526 chego a ter milhares de números! Veja que não há espaço entre os dígitos da sequência! Gostaria de ler esse arquivo para um vetor mas, ao ler, quero que o R entenda que cada dígito deva ser armazenado numa posição diferente no vetor criado. Por ex. (no arquivo acima): 0 1 2 0 3 0 4 4 0 5 5 6 0 0 8... cada um pertencendo a uma posição no vetor criado. Agradeço qualquer ajuda, -- *Emerson*

4 7

Re: [R-br] Teste de homocedasticidade e ANOVA
by Ivan Bezerra Allaman 04 Jan '12

04 Jan '12

Antes de sermos mecanicistas, vamos nos ater no problema como um todo. Se você tomou medidas em mesmo animal durante intervalos de tempo, essas medidas não são independentes, infrigindo um dos primeiros pressupostos da ANAVA. Se você utilizar as médias dessas medidas, obtendo um valor para cada animal, aí tudo bem, pois a unidade experimental será o animal. No entanto, você irá ter perda de informações. (S,f,P) Allaman \begin{signature} <<>>= Prof. Dr. Ivan Bezerra Allaman Universidade Estadual de Santa Cruz Departamento de Ciências Exatas e Tecnológicas Ilhéus/BA - Brasil Fone: +55 73 3680-5076 E-mail: ivanalaman@yahoo.com.br/ivanalaman@gmail.com @ \end{signature}

3 2

Problema para reinstalar pacote (EspDes) no Linux Ubuntu 10.04
by Andre Oliveira 04 Jan '12

04 Jan '12

Pessoal estou com problemas para reinstalar o EspDes !! Alguém saberia dizer o que pode ser? install.packages("EspDes", dep=T) Installing package(s) into ‘/home/andre/R/i486-pc-linux-gnu-library/2.14’ (as ‘lib’ is unspecified) --- Please select a CRAN mirror for use in this session --- Carregando a interface Tcl/Tk ... pronto Mensagens de aviso perdidas: In getDependencies(pkgs, dependencies, available, lib) : package ‘EspDes’ is not available (for R version 2.14.1) André Oliveira Souza obrigado

2 2

Teste de homocedasticidade e ANOVA
by Flávia Machado 04 Jan '12

04 Jan '12

Pessoal, trabalho com dieta de muriquis (primata da Mata Atlântica) e quero, com os dados que obtive, verificar qual foi o alimento mais consumido em um determinado período de tempo (todo período de estudo, estação seca e e. chuvosa). São quatro categorias alimentícias (Fruto, Folha, Broto e Flor) e calculei, para cada uma delas, as frequências média de consumo para os três períodos de tempo. Então, cada linha do meu arquivo corresponde a frequência de consumo de cada um desses em um determinado dia, no determinado período. Fiz uma ANOVA, mas não soube interpretar o valor do teste e o bloxplot. De qualquer forma, pode ter dado errado porque eu não sabia que tinha que fazer o teste de homocedasticidade. Não entendi o que preciso colocar para fazer o "levene's test". Na ajuda do R diz: levene.test(y, group, option = c("mean", "median","trim.mean"), trim.alpha = 1) Arguments your data groupfactor of the data optionThe option must be "mean" corresponding to classical Levene's test (default), "median" corresponding to the robust Brown-Forsythe Levene-type test or "trim.mean" corresponding to the robust Levene-type test using the trimmed mean. trim.alphathe fraction (0 to 0.5) of observations to be trimmed from each end of 'x' before the mean is computed. Porém, continuo sem entender... Alguém pode me dar um help? Abraços! Att. -- Flávia Machado

3 3

Escala de y
by Fernanda Balbino 04 Jan '12

04 Jan '12

Olá, boa tarde, Gostaria de saber como alterar em um gráfico a escala em y. Por exemplo, no caso gostaria que ela fosse de 10 em 10 unidades. Obrigada. Fernanda Balbino.

2 2

Alocação de Memória
by Edson Lira 04 Jan '12

04 Jan '12

Amigos, tenho um arquivo com milhões de registro. Uso a rotina abaixo. require(RODBC) channel<-odbcConnect("saddw",uid="saddw",case="tolower") #SOROLOGIA GERAL channel<-odbcConnect("saddw",uid="saddw",case="tolower") > soro<- sqlQuery(channel, "select doador_dw.sexo,doador_dw.idade,doador_dw.qtde_doacoes, + doador_dw.numerogente,doacao_dw.numerogente,sorologias_dw.numerogente, + sorologias_dw.exame,sorologias_dw.patologia,sorologias_dw.resultado,doacao_dw.tipodoacao, + sorologias_dw.metodo,sorologias_dw.basicasorologia,sorologias_dw.dataexamesorologia + from saddw.doador_dw, saddw.doacao_dw,saddw.sorologias_dw + where doador_dw.numerogente=sorologias_dw.numerogente + and doacao_dw.numerogente=doador_dw.numerogente + and to_char(sorologias_dw.dataexamesorologia,'DD/MM/YYYY') >= to_date('01/01/2010','DD/MM/YYYY') + and to_char(sorologias_dw.dataexamesorologia,'DD/MM/YYYY') <= to_date('31/12/2011','DD/MM/YYYY')") tenho esse erro, como resolver? Erro: não é possível alocar vetor de tamanho 26.8 Mb Meu computador tem 2 giga de ram. Uso o windows 7 32 bits. A versão do R é a 2.14 Edson Lira Estatístico Manaus-Amazonas

3 2