- R-br - DINF UFPR

ordenar colunas - seleção expandida
by marcelo claro de souza 18 Out '12

18 Out '12

Olá pessoal, Preciso reordenar uma coluna de um conjunto de dados, mas gostaria que todas as outras fossem reordenadas juntamente a primeira em estilo seleção expandida, assim como no excel. Tentei o comando sort , mas só consegui ordenar uma coluna de vez. Alguma dica? Muito obrigado. Abraço Marcelo

6 11

plot de uma timeseries univariada com seu modelo SARIMA fitted?
by carlos pombo sonderblohm 18 Out '12

18 Out '12

Bom dia pessoal, estou a tentar fazer um plot de um modelo SARIMA de uma serie temporal> library(stats) library(astsa) library(graphics) SL <-c(-0.289220360456518, -0.412287982738895, -0.30186163626253, -0.325463729174269, -0.950943918812091, -0.621873409193565, -0.454370839103369, -0.201717177298859, -0.0125484487133019, 0.37501627306213, 0.860628415115297, 0.573001978891618, 0.949171439258297, 0.884107931496505, 0.321528581746018, 0.197039928762375, -0.0511602427005644, -0.14797533453817, 0.0997648324025207, -0.356603216398346, -0.486741440877119, -0.796433339558246, -0.797717898751611, -1.03991357760201, -0.272391281348932, -0.280840105493716, -0.220566038983021, -0.092982026032348, -0.737941509253829, -0.25720443208744, 0.281594104947966, 0.295827677265899, 0.381173023222877, 0.410685495501984, 0.0935538548738117, 0.320870181634448, -0.176916211372866, -0.434468255568435, -0.177142353391016, -0.479723718529784, -1.18506069524595, -0.845078218877219, -0.466484365108015, -0.378310531861291, -0.568334035779147, -0.756612458676533, -0.612289956791277, -0.595067441206468, 0.227790700535617, -0.328474103547927, 0.0382893789935963, -0.107472883393398, 0.104836312989005, 0.206908642082212, 0.566811114425553, 0.226484501566576, 0.268795384938629, 0.247658319827575, 0.728293948754543, 0.704287103883026, 0.665057803685019, 1.01752383083208, 1.34104381438303, 1.10724539596799, 0.850746327229352, 0.746239914170796, 0.665531282137068, 0.469573210434806, -0.0626525410819481, -0.294702780105021, -0.160850370184094, 0.511318161080905) SL <- ts(SL, start = c (1995, 1), freq = 12) #transformação numa ts ar1 <-sarima(SL, 1,0,0,0,0,0,12) #o modelo seleccionado plot(SL, type ='o')) # o plot da serie original lines(fitted(ar1),col=2) # estou a tentar adicionar o modelo fitted a meus dados Mas não consigo fazer o grafico pretendido, quero fazer que o grafico tenha os pontos o linhas da serie original e no mesmo grafico, o modelo ajustado para ver como bate com os dados originais, alguma dica? Meu R es: R version 2.13.1 (2011-07-08) Copyright (C) 2011 The R Foundation for Statistical Computing ISBN 3-900051-07-0 Platform: x86_64-pc-mingw32/x64 (64-bit) SO, windows 7 Obrigado Carlos -- Carlos A. Pombo Sonderblohm PhD Student on Marine Science (Fisheries) Faculdade de Ciências e Tecnología Universidade do Algarve, Campus de Gambelas 8005-139 Faro Portugal Tef. 289 800 905 ext. 7605

1 0

Trabalhar com muitos dados
by Roney Fraga 17 Out '12

17 Out '12

Caros, Estou trabalhando com os dados dos censos populacionais de 2000 e 2010. A minha dificuldade é o volume de dados, pois cada censo tem aproximadamente 20 milhões de linhas. Considerando os dados estão em arquivos .csv, um único aquivo de 12 GB para o censo de 2010 e um arquivo para cada estado no censo de 2000, que totalizam 11 GB. Utilizo um notebook com processador Intel Core i5 2410M, 8 GB de ram com 3 GB de SWAP e o sistema operacional é o Arch Linux. Tenho quatro dúvidas que estão abaixo em meio da minha explicação sobre os problemas que venho enfrentando. Para ler os arquivos .csv conheço as seguintes opções: a) trabalhar com um banco de dados relacional. R.: estou com um problema com o postgresql que ainda não tenho solução, portanto, essa alternativa esta descartada no momento. b) utilizar o sqldf. R.: acredito que o sqldf não é eficiênte para trabalhar com arquivos grandes, o processo é encerrado antes de concluir a operação. Exemplo: > library(sqldf) > read.csv.sql('censo00ba1.csv', sql = "select V0102 from file", header = TRUE, sep = ",") [ reached getOption("max.print") -- omitted 1498127 rows ] Warning messages: 1: closing unused connection 4 (censo00ba1.csv) 2: closing unused connection 5 (censo00ba1.csv) 1ª dúvida: isso sempre acontece com o sqldf, ou alguém consegue utilizá-lo com para ler arquivos grandes? c) utilizar o comando cut e grep do Linux/Mac para filtrar os dados e ter um arquivo .csv menor, dessa forma é possível ler o arquivo pelo comando do R, read.csv(). R.: utilizei essa opção. Exemplo: cut -d',' -f1,4,5,6,10,79,81,84,91,92,95,97,111,113,114,149,150,188 br2010.csv > br2010b.csv Uma vez superada a dificuldade da importação dos dados deletei algumas coisas e fiquei com três arquivos com as seguintes dimensões: > dim(cs);dim(cs2000);dim(cs2010) [1] 6077327 23 # censos 2000 e 2010 [1] 2865716 23 # censo 2000 [1] 3211611 23 # censo 2010 Quanto ao processamento: d) Consigo rodar regressão simples para o censo de 2000 e para o censo de 2010 separadamente, como segue exemplo: > ma <- lm(lnrenda ~ factor(ee) + idade + idade2 + factor(tt) + factor(cor.raca) + rural + feminino + migrante + factor(UF), data=cs2000) > mb <- lm(lnrenda ~ factor(ee) + idade + idade2 + factor(tt) + factor(cor.raca) + rural + feminino + migrante + factor(UF), data=cs2010) mas não consigo rodar o mesmo modelo para os dados do arquivo 'cs', que tem 6.077.327 linhas, que é censo 2000 mais censo 2010. Pois, aparece a seguinte mensagem de erro: > mc <- lm(lnrenda ~ factor(ee) + idade + idade2 + factor(tt) + factor(cor.raca) + rural + feminino + migrante + factor(UF) + factor(ano), data=cs) Error: cannot allocate vector of size 2.0 Gb 2ª dúvida: qual procedimento posso utilizar para solucionar esse problema de rodar a regressão das 6.077.327 observações? Um amigo indicou aumentar a memória SWAP, essa é uma boa opção? 3ª dúvida: ao realizar um processo que demanda muita memória ram o R não consegue realizar outro processo pesado em seguida, pois a memória fica ocupada. Exemplo, ao rodar uma regressão em uma base de dados com 3 milhões de linhas eu não consigo realizar outra regressão pq a memória fica cheia. Dúvida, tem como 'limpar' os dados da memória para poder realizar outros cálculos? 4ª dúvida: trabalhar com paralelismo pode auxiliar nesse tipo de processamento? Grato pela atenção. Roney

8 11

Re: [R-br] Informação sobre efeito marginal para Modelos Logit Ordenado Generalizado
by Ivan Bezerra Allaman 17 Out '12

17 Out '12

fitted(mod) e terás o que quer! (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-5596 E-mail: ivanalaman@yahoo.com.br/ivanalaman@gmail.com @ \end{signature}

4 4

RecordLinkage
by Fátima Lima Paula 16 Out '12

16 Out '12

Prezados, alguém que seja íntimo do pacote RecordLinkage para me dar uma luz? Tenho um banco de dados de internados e quero buscar reinternação. Quais os tratamentos que preciso fazer no banco? Como vou relacionar um banco com ele mesmo, vou usar o comando compare.dedup? Não sei quando usar fonética. Desculpe, mas não manjo nada disso. Obrigada Fátima

2 2

Função dist()
by Alisson Lucrecio 16 Out '12

16 Out '12

Caro colegas da lista r-br, Boa tarde. Eu preciso aplicar a função dist() para calcular a distancia euclidiana de uma planilha. Quando aplico a função ela considera a rownames, porem eu nao quero que ela considere a rownames e considere uma coluna com os tratamento que estão na primeira coluna da minha tabela. Como faço isso? Obrigado. Alisson Lucrécio da Costa

2 1

Informação sobre efeito marginal para Modelos Logit Ordenado Generalizado
by Abel Brasil Ramos da Silva 16 Out '12

16 Out '12

Olá a todos, Estou atualmente trabalhando com modelos logit ordenados generalizados, estou trabalhando no R atráves do pacote VGAM (vector generalized linear and additive model) na função vglm(). Bom, eu consigo estimar o modelo com ou sem hipótese de regressão paralela (proportional odds model). Meu objetivo é estimar o efeito marginal da probabilidade, para isso, utilizo a função margeff() contida no pacote. Porém, a função executa o efeito para cada indivíduo, sem erro-padrão e estatística de teste. Então, gostaria de saber se existe outro modo (ou pacote) de estimar tais efeitos marginais de forma completa. P.S: já utilizei as funções clm (pacote ordinal) e polr() (pacote MASS) e não obtive sucesso. Agradeço a ajuda. Abel Brasil Estatístico - UFC

4 6

Base Geográfica de Coordenadas
by Listeiro 037 16 Out '12

16 Out '12

Olá a todos. Há algum tempo tinha perguntado sobre uma base geográfica de coordenadas. Lá no site do IBGE não consegui muita coisa, então continuei procurando. Daí encontrei este site: http://www.maxmind.com/en/worldcities E este link para download: http://www.maxmind.com/download/worldcities/worldcitiespop.txt.gz Possui coordenadas de localidades do mundo todo. No site é dito que constam estes dados: Country Code ASCII City Name City Name Region Population Latitude Longitude Dos quais os mais importantes seriam as coordenadas. Como ele possui 3173959 linhas, 40 MB em arquivo de texto compactado e 144 MB quando descompactado, fica difícil verificar a correção de todas essas coordenadas, mas digo que vi alguma linha duplicada. Creio haver nesta relação cerca de 235 códigos de países diferentes. O texto pode ser trabalhado e modificado com algum script tipo Perl, AWK ou no próprio R, até para ser transformado em SQL. Se alguém souber algo sobre a qualidade destes dados, eu agradeceria.

1 0

Família APPLY...
by andrebvs 15 Out '12

15 Out '12

6 7

get.best.arima <- function (Error in solve.default(res$hessian * n.used) : Lapack routine dgesv: system is exactly singular)
by carlos pombo sonderblohm 15 Out '12

15 Out '12

Caros amigos, estou a correr uma função para obter o *best fitting model* baseado no AIC. Esta função esta no livro de Introductory time series with R, de Metcalfe and Cowperwait na pagina 144, dispoinivel no site: http://tur-www1.massey.ac.nz/~pscowper/ts/scripts.R a serie a modelar es structure(c(-0.289220360456516, -0.412287982738894, -0.301861636262529, -0.325463729174269, -0.95094391881209, -0.621873409193565, -0.454370839103367, -0.201717177298859, -0.0125484487133019, 0.375016273062132, 0.860628415115299, 0.573001978891618, 0.949171439258297, 0.884107931496505, 0.321528581746019, 0.197039928762376, -0.0511602427005631, -0.14797533453817, 0.0997648324025207, -0.356603216398344, -0.486741440877118, -0.796433339558245, -0.79771789875161, -1.03991357760201, -0.272391281348931, -0.280840105493716, -0.220566038983019, -0.0929820260323471, -0.737941509253828, -0.257204432087439, 0.281594104947967, 0.295827677265899, 0.381173023222876, 0.410685495501985, 0.0935538548738126, 0.320870181634449, -0.176916211372864, -0.434468255568435, -0.177142353391016, -0.479723718529783, -1.18506069524595, -0.845078218877218, -0.466484365108014, -0.37831053186129, -0.568334035779147, -0.756612458676532, -0.612289956791276, -0.595067441206467, 0.227790700535618, -0.328474103547926, 0.0382893789935963, -0.107472883393397, 0.104836312989006, 0.206908642082212, 0.566811114425555, 0.226484501566576, 0.268795384938631, 0.247658319827575, 0.728293948754544, 0.704287103883027, 0.66505780368502, 1.01752383083208, 1.34104381438303, 1.10724539596799, 0.850746327229352, 0.746239914170797, 0.665531282137069, 0.469573210434808, -0.0626525410819472, -0.29470278010502, -0.160850370184093, 0.511318161080906), .Tsp = c(1995, 2000.91666666667, 12), class = "ts") a função es # Chapter 7, p144 get.best.arima <- function(x.ts, maxord = c(1,1,1,1,1,1)) { best.aic <- 1e8 n <- length(x.ts) for (p in 0:maxord[1]) for(d in 0:maxord[2]) for(q in 0:maxord[3]) for (P in 0:maxord[4]) for(D in 0:maxord[5]) for(Q in 0:maxord[6]) { fit <- arima(x.ts, order = c(p,d,q), seas = list(order = c(P,D,Q), frequency(x.ts)), method = "CSS") fit.aic <- -2 * fit$loglik + (log(n) + 1) * length(fit$coef) if (fit.aic < best.aic) { best.aic <- fit.aic best.fit <- fit best.model <- c(p,d,q,P,D,Q) } } list(best.aic, best.fit, best.model) } O erro esta aqui: > best.arima.SLCENT <- get.best.arima(SLCENT, maxord = c(2,2,2,2,2,2)) Error in solve.default(res$hessian * n.used) : Lapack routine dgesv: system is exactly singular Muito obrigado Carlos P.D.: meu R es R version 2.13.1 (2011-07-08) Copyright (C) 2011 The R Foundation for Statistical Computing ISBN 3-900051-07-0 Platform: x86_64-pc-mingw32/x64 (64-bit) -- Carlos A. Pombo Sonderblohm PhD Student on Marine Science (Fisheries) Faculdade de Ciências e Tecnología Universidade do Algarve, Campus de Gambelas 8005-139 Faro Portugal Tef. 289 800 905 ext. 7605

1 0

como transformar uma matriz em uma lista de submatrizes?
by Cleber N.Borges 15 Out '12

15 Out '12

Olá a todos, Como eu faço para transformar uma matrix com n linhas e 4 colunas em uma lista de comprimento n onde cada objeto dessa lista é uma submatriz 2x2 (formada a partir de cada linha da matriz original)? Agradeço antecipadamente por qualquer ajuda. Cleber N.Borges # exemplo set.seed(12345) m <- matrix( runif(2*4), ncol=4, nrow=2 ) m list( matrix(m[1,],nc=2,nr=2), matrix(m[2,], nc=2,nr=2))

3 2

Linkage
by Fátima Lima Paula 15 Out '12

15 Out '12

Prezados, preciso fazer um linkage em banco de dados para identificar pacientes que foram reinternados. Alguém sabe me dizer se posso fazer isso no R? Obrigada Fátima

2 2

Re: [R-br] como transformar uma matriz em uma lista de submatrizes?
by Ivan Bezerra Allaman 14 Out '12

14 Out '12

lista <- list() for(i in 1:dim(m)[1]){ lista[[i]] <- matrix(m[i,],nc=2,nr=2) } lista Estou com preguiça de usar o lapply! \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-5596 E-mail: ivanalaman@yahoo.com.br/ivanalaman@gmail.com @ \end{signature}

1 0

Filtrar por data
by Gabriel Bruno de Lemos 14 Out '12

14 Out '12

Boa tarde! Tenho um banco de dados com observações diárias durante um período de 50 anos. Desejo filtrar alguns períodos específicos, tais como todos os dias 10/04 de cada ano. Alguém saberia como realizar tal busca por data? A forma que utilizei para resolver no jeitinho brasileiro foi usar o Excel para criar uma coluna com o dia e outra com o mês de cada data. Em seguida, no R, faço uma busca usando o comando ## Datas ini=as.Date("1961-01-01",origin="1899-12-30",format="%Y-%m-%d") fim=as.Date("2010-12-31",origin="1899-12-30",format="%Y-%m-%d") serie=seq.Date(from=ini,to=fim,by="day") ## Criei o dados_excel no Excel usando as funções Dia() e Mês. Em seguida criei a tabela no R adicionando estes dois vetores à série Data. tabela <- data.frame(serie,dados_excel) head(tabela) # dados.Data Dia Mes Ano #1 1961-01-01 1 1 1961 #2 1961-01-02 2 1 1961 #3 1961-01-03 3 1 1961 #4 1961-01-04 4 1 1961 #5 1961-01-05 5 1 1961 #6 1961-01-06 6 1 1961 head(tabela[tabela$Dia==10 & tabela$Mes==4,]) # dados.Data Dia Mes Ano #100 1961-04-10 10 4 1961 #465 1962-04-10 10 4 1962 #830 1963-04-10 10 4 1963 #1196 1964-04-10 10 4 1964 #1561 1965-04-10 10 4 1965 #1926 1966-04-10 10 4 1966 Este mecanismo funcionou mas gostaria de saber se existe algo mais automatizado. Obs: O código acima só será reproduzível se for criado os vetores Dia e Mes (que eu não fiz no R) Att,

2 6

Gráfico ARIMA
by kaue veras 14 Out '12

14 Out '12

Boa tarde Pessoal, Então estou fazendo um trabalho em que faço previsões em Holt-Winters e ARIMA, para Holt-Winters consegui gerar um grafico em que ele me da a linha de observado e esperado pelo HW, mas não consigo fazer o mesmo pelo ARIMA. Alguém teria alguma idéia do que fazer pra plotar as linhas de observado e esperado pelo modelo ARIMA? Desde já agradeço. Att, Kauê P. Veras da Cunha Estatística - UERJ MSN: kaueveras(a)hotmail.com E-mail: kaueveras(a)gmail.com Cel: (21) 8254-9601

1 0

Re: [R-br] plug-in para modelos não-lineares em r commander
by Gustavo Marcatti 13 Out '12

13 Out '12

Ola Walmes, obrigado pela resposta. Eu também prefiro trabalhar com o R atraves de scripts, é bem mais flexivel. O interesse no r comander foi um pedido de um colega que precisará fazer algumas análises no R e não tem familiaridade alguma com o software. Quanto a interface para modelos não lineares eu ja vi em um tutorial na internet, é bem parecida com a de modelos lineares. Como nao encontrei nada ainda pode ser que seja pessoal de algum usuário. Gustavo

1 0

Dendrograma com níveis defidos
by Alisson Lucrecio 13 Out '12

13 Out '12

Caros colegas da lista r-br, Como construir um dendrograma para classificar solos com valores com valores de order, suborder, greatgroup, subgroup, class e series? O problemas é que esses valores não são numéricos e eles segue um ordem, exemplo: o dendrograma começa em order e termina em series, os níveis já estão definidos, apenas quero colocar-los juntos e ligados por um dendrograma. Obrigado. Alisson Lucrécio da Costa

1 0

plug-in para modelos não-lineares em r commander
by Gustavo Marcatti 13 Out '12

13 Out '12

Olá pessoal, tenho um colega que está precisando ajustar uns modelos não-lineares no r commander, já vi em alguns tutoriais que o r commander tem uma interface para isso (não sei se tal funcionalidade é do r commander ou é pessoal de algum usuário), porém não consegui encontrar na instalação padrão. Se possível gostaria de saber se existe essa interface para modelos mistos também. Desde já agradeço. Gustavo Marcatti Eng Florestal

2 1

Re: [R-br] Digest R-br, volume 20, assunto 16
by João Augusto Batista 13 Out '12

13 Out '12

Olá Elias, muito obrigado pela ajuda! EU baixei tanto da pasta de municípios quanto pras áreas de ponderação. Estou começando agora esse novo projeto e ainda vou determinar o que é melhor pro trabalho! Assim que possível vou testar a sua solução, e pelo visto via dar certo. Mais uma vez, muito obrigado pela atenção! Att, João Augusto

1 0

maptools
by João Augusto P Batista 13 Out '12

13 Out '12

Olá! Esta é minha primeira mensagem, vou tentar fazê-la o mais inteligível possível. Estou tentando utilizar os mapas do IBGE, através da função readShapePoly(), mas sempre dá erro. Os arquivos na pasta dos mapas tem extensão .dbf, .shp, .shx, .prj. Segue a rotina e o output: > library(descr) > library(Hmisc) > library(foreign) > library(maptools) > mapa <- readShapePoints("Mapas IBGE/Municípios/33MEE250GC_SIR") Shapefile type: PolygonM (25), # of Shapes: 6 Error in read.shape(filen = fn, verbose = verbose, repair = repair) : Shapefile type not (yet) handled by this function Já procurei nas mensagens anteriores, mas nenhuma respondeu o meu problema, e já tentei readOGR() e readShapePoints(), sem sucesso. O que fazer? Muito obrigado, João Augusto.

2 1

Valores de probabilidade de Dunnett
by tiago souza marçal 11 Out '12

11 Out '12

Alguem sabe me enformar se existe no R uma função que retorne os valores de probabilidade de Dunnett? Pois eu só consegui encontrar a de Tukey qtukey() e a de Student qt(). Att. Tiago de Souza.

4 4

VALORES TABELADOS DE DUNNETT
by tiago souza marçal 11 Out '12

11 Out '12

Bom dia pessoal, gostaria de saber se alguém sabe como calcular os valores de probabilidade tabelados de Dunnett no R? Att. Tiago.

3 2

Re: [R-br] ordenar colunas em data.frame por data
by Ivan Bezerra Allaman 11 Out '12

11 Out '12

CMR <- data.frame(dates = format(seq(ISOdate(2003,1,1), by='month', length=36), format='%d-%m-%Y'), var=rpois(36,10)) summary(CMR) str(CMR) with(CMR,CMR[order(as.Date(CMR$dates,format='%d-%m-%Y'),decreasing=TRUE),]) with(CMR,CMR[order(as.Date(CMR$dates,format='%d-%m-%Y'),decreasing=FALSE),]) (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-5596 E-mail: ivanalaman@yahoo.com.br/ivanalaman@gmail.com @ \end{signature}

2 1

shapefile novamente
by Heloíse Pavanato 11 Out '12

11 Out '12

Olá pessoal, Repetindo um pouco as dúvidas sobre a leitura de arquivo .shp, não consigo ler o arquivo com nenhuma das funções readShapePoly e readOGR. Ocorrem os seguintes erros: > ae <- readShapePoly('area_de_estudo_2011.shp') Erro em read.dbf(filen1) : unable to open DBF file > ae <- readOGR('area_de_estudo_2011.shp') Erro em readOGR("area_de_estudo_2011.shp") : missing layer Alguém sabe como solucionar? Obrigada, Heloise

4 5

[OFF] Ajuda na obtenção de artigo
by Andre Oliveira 10 Out '12

10 Out '12

Alguém tem este artigo em pdf para compartilhar? Fontelles MJ, Simões MG, Almeida JC, Fontelles RGS. Metodologia da pesquisa: diretrizes para o cálculo do tamanho da amostra. Rev Paran Med. 2010;24:57-64. André Oliveira Souza

5 6

Pacote em R para trabalhar com power laws
by Felipe Buchbinder 10 Out '12

10 Out '12

Prezados, Alguém conhece um pacote em R que seja voltado para trabalhar com *power laws*? Sei que o pacote igraph faz um ajuste de *power law*, mas só para distribuilções discretas e não me consta que ele tenha qualquer função para avaliar a qualidade do ajuste. obrigado, Felipe

2 1

ordenar colunas em data.frame por data
by Adevair Henrique Da Fonseca 10 Out '12

10 Out '12

Prezados, estou com dificuldades em ordenar os dados por data, os dados que estou usando foram filtrados de um banco maior e nesta alterou a sequencia de datas. Tentei varias alternativas inclusive a postada hoje. > amb_l1[with(amb_l1, order(d_coleta)),] ponto ano d_coleta amb_l 1 6 3 3-05-11 0 2 6 4 4-01-12 0 3 6 1 4-12-08 0 4 6 4 7-08-12 2 5 6 1 8-05-09 74 6 6 3 9-08-11 0 7 6 3 10-11-10 145 8 6 2 13-05-10 154 9 6 2 13-07-10 1352 10 6 2 14-10-09 0 11 6 3 15-02-11 0 12 6 1 19-01-09 0 13 6 4 19-04-12 0 14 6 4 26-11-11 0 15 6 2 28-01-10 0 16 6 1 30-07-09 226 str(amb_l1) 'data.frame': 16 obs. of 4 variables: $ ponto : Factor w/ 35 levels "1","2","3","4",..: 6 6 6 6 6 6 6 6 6 6 ... $ ano : Factor w/ 4 levels "1","2","3","4": 3 4 1 4 1 3 3 2 2 2 ... $ d_coleta: Date, format: "3-05-11" "4-01-12" ... $ amb_l : int 0 0 0 2 74 0 145 154 1352 0 ... sou curuja nesta lista! obrigado pela atencao. Abracos, Adevair

1 0

Re: [R-br] arumento add=T
by Ivan Bezerra Allaman 10 Out '12

10 Out '12

O erro ocorre porque "add" não é um parâmetro da função "plot". Você deve estar se confundindo com a função "curve" ou outra no qual o parâmetro "add" é utilizado para acrescentar linhas ou qualquer coisa no gráfico "plot". Para juntar dois "plots" vc tem que usar um "par" após o primeiro plot. Por exemplo: plot(seq(1:10),seq(1:10)) par(new=TRUE) plot(rnorm(100),col='blue') E não esqueça de um CMR da próxima vez para que de fato a ajuda seja eficiente. (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-5596 E-mail: ivanalaman@yahoo.com.br/ivanalaman@gmail.com @ \end{signature}

4 4

Evitar gravação
by Carlos Mendonça 10 Out '12

10 Out '12

Meus amigos, o meu problema é como fazer com que um arquivo não seja gravado se na etapa anterior do programa acontece um erro. No exemplo, eu não gostaria que o arquivo .csv fosse gravado pois aconteceu um erro na etapa anterior do programa. Exemplo: > arquivo1= merge(arquivo1, arquivo99, by = "code6", all.x = TRUE) > > arquivo1l$2004t1 = ifelse(is.na(arquivo1$i2004t1), 100, arquivo1$i2004t1) Erro: constante numérica inesperado em "serie_bmovel$2004" > > write.csv2(arquivo1, "\\arquivo1.csv", row.names = FALSE)

4 4

Rotina para classificação de imagens
by Mauro Sznelwar 10 Out '12

10 Out '12

Queria saber se já existe alguma bibioteca no R que trata sobre classificação de imagens. Mesmo que não tenha coisa pronta. se existe em algum lugar estas rotinas.

2 1

Erro em read.spss
by Pedro Emmanuel Alvarenga Americano do Brasil 10 Out '12

10 Out '12

Amigos de R, É frequente comigo quando importo dados no formato utilizado pelo SPSS (*.sav) a função read.spss do pacote foreign retornar um erro e retornar o banco tambem. O Erro é sempre algo assim: Mensagens de aviso perdidas: 1: In read.spss(filename) : C-toc modificado_diagrama.sav: Unrecognized record type 7, subtype 8 encountered in system file 2: In read.spss(filename) : C-toc modificado_diagrama.sav: Unrecognized record type 7, subtype 18 encountered in system file 3: In read.spss(filename, to.data.frame = TRUE, trim.factor.names = TRUE) : C-toc modificado_diagrama.sav: Unrecognized record type 7, subtype 8 encountered in system file 4: In read.spss(filename, to.data.frame = TRUE, trim.factor.names = TRUE) : C-toc modificado_diagrama.sav: Unrecognized record type 7, subtype 18 encountered in system file Neste exemplo, o banco veio com o número de linhas e de colunas esperadas. A questão é que eu não faço ideia o que esse erro significa e por isso não sei se há algum risco de esses erros terem modificado o banco na importação. Alguém tem algum comentário que me passe alguma sensação de frescor e felicidade? Abraço forte, Dr. Pedro Emmanuel A. A. do Brasil Curriculum Lattes: http://lattes.cnpq.br/6597654894290806 Instituto de Pesquisa Clínica Evandro Chagas Fundação Oswaldo Cruz Rio de Janeiro - Brasil Av. Brasil 4365, CEP 21040-360, Tel 55 21 3865-9648 email: pedro.brasil(a)ipec.fiocruz.br email: emmanuel.brasil(a)gmail.com ---Apoio aos softwares livres www.zotero.org - gerenciamento de referências bibliográficas. www.broffice.org ou www.libreoffice.org - textos, planilhas ou apresentações. www.epidata.dk - entrada de dados. www.r-project.org - análise de dados. www.ubuntu.com - sistema operacional

4 4

mensagem de aviso perdidas
by Fátima Lima Paula 09 Out '12

09 Out '12

Prezados, quero trabalhar com um banco que é formado por outros menores. Ao todo são 34 bancos com aproximadamente 30.000 linhas cada um que foram unidos formando 3 bancos grandes. Como não é necessário trabalhar com todas as linhas, fui eliminando o que não precisava e fiquei com 3 bancos com as seguintes dimensões: dim(banco2008idosonovo) [1] 76855 43 > dim(banco2009idosonovo) [1] 77870 43 > dim(banco2010idosonovo) [1] 63195 43 Quando vou juntar os bancos para poder fazer minhas análises, aparece o seguinte: bancocompleto=rbind(banco2008idosonovo,banco2009idosonovo,banco2010idosonovo) Mensagens de aviso perdidas: 1: In unique.default(c(all.levs[[jj]], lij)) : Reached total allocation of 1535Mb: see help(memory.size) 2: In unique.default(c(all.levs[[jj]], lij)) : Reached total allocation of 1535Mb: see help(memory.size) 3: In rbind(deparse.level, ...) : Reached total allocation of 1535Mb: see help(memory.size) 4: In rbind(deparse.level, ...) : Reached total allocation of 1535Mb: see help(memory.size) 5: In as.vector(levels(x)[x], mode) : Reached total allocation of 1535Mb: see help(memory.size) Alguém poderia me ajudar? Obrigada Fátima

2 2

Razão entre normais
by Eric Ferreira 09 Out '12

09 Out '12

Prezados colegas, Alguém conhece algum pacote facilita a obtenção da densidade da razão entre duas normais? (Digo, se X~N(m1,s1) e Y~N(m2,s2), então estou interessado na densidade de Z = X/Y) Grato. -- Dr Eric B Ferreira Exact Sciences Department Federal University of Alfenas Brazil

3 3

Re: [R-br] Razão entre normais
by Ivan Bezerra Allaman 09 Out '12

09 Out '12

Apr,oveitando o gancho do Walmes, se o seu interesse for na razão entre duas normais padrão que dará um cauchy, você pode tentar: ?dcauchy ?qcauchy que é da "stats" (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-5596 E-mail: ivanalaman@yahoo.com.br/ivanalaman@gmail.com @ \end{signature}

1 0

arumento add=T
by Stephan M. CARVALHO 09 Out '12

09 Out '12

Prezados, Encontrei um erro com o argumento add=T na função plot. Até então eu utilizava esse argumento p/ juntar duas funções plot, só que agora sempre ocorre erro. Procurando no help não encontrei nada. Por acaso alguém sabe o motivo? att, Stephan -- *************************************************** STEPHAN M. CARVALHO "Honeybee protection and ecotoxicology" Universidade Estadual Paulista Centro de Estudo de Insetos Sociais Av. 24A, 1515 Bela Vista 13.506-900 Rio Claro/SP Brazil smalfitano(a)uol.com.br msn: pefasmc(a)hotmail.com / Skype: smalfitano +55 19 3526-4284 (university) +55 19 3534-8523 (fax) +55 35 9103-9822 (mobile) +55 35 9949-0707 (mobile) Antes de imprimir, pense no ambiente! Before printing, think in the environment! Avant d'imprimer, pensez à l'environnement! ***************************************************

2 1

Erro em multiplicação - estouro de inteiros
by Pedro Emmanuel Alvarenga Americano do Brasil 08 Out '12

08 Out '12

Amigos de R, estou surpreso com o erro que encontrei. Alguem consegue propor alguma solução? > lX [1] 250000 > lX/2 [1] 125000 > lX * lX/2 [1] NA Mensagens de aviso perdidas: In lX * lX : NAs produzidos por estouro de inteiro Abraço forte, Dr. Pedro Emmanuel A. A. do Brasil Curriculum Lattes: http://lattes.cnpq.br/6597654894290806 Instituto de Pesquisa Clínica Evandro Chagas Fundação Oswaldo Cruz Rio de Janeiro - Brasil Av. Brasil 4365, CEP 21040-360, Tel 55 21 3865-9648 email: pedro.brasil(a)ipec.fiocruz.br email: emmanuel.brasil(a)gmail.com ---Apoio aos softwares livres www.zotero.org - gerenciamento de referências bibliográficas. www.broffice.org ou www.libreoffice.org - textos, planilhas ou apresentações. www.epidata.dk - entrada de dados. www.r-project.org - análise de dados. www.ubuntu.com - sistema operacional

6 7

Fwd: Erro em multiplicação - estouro de inteiros
by Jodavid Ferreira 08 Out '12

08 Out '12

Olá, Pedro, fiz aqui no meu computador e obtive o resultado normalmente > lX[1] 250000> lX/2[1] 125000> lX*lX/2 [1] 3.125e+10 Mas vc pode tentar fazer o inverso > lX/2 * lX[1] 3.125e+10 para que primeiro o R faça a divisão e depois multiplique! Em 8 de outubro de 2012 11:43, Pedro Emmanuel Alvarenga Americano do Brasil [via R-br] <ml-node+s2285057n4656458h33(a)n4.nabble.com> escreveu: Amigos de R, > > estou surpreso com o erro que encontrei. Alguem consegue propor alguma > solução? > > > lX > [1] 250000 > > > lX/2 > [1] 125000 > > > lX * lX/2 > [1] NA > Mensagens de aviso perdidas: > In lX * lX : NAs produzidos por estouro de inteiro > > Abraço forte, > > Dr. Pedro Emmanuel A. A. do Brasil > Curriculum Lattes: http://lattes.cnpq.br/6597654894290806 > Instituto de Pesquisa Clínica Evandro Chagas > Fundação Oswaldo Cruz > Rio de Janeiro - Brasil > Av. Brasil 4365, > CEP 21040-360, > Tel 55 21 3865-9648 > email: [hidden email]<http://user/SendEmail.jtp?type=node&node=4656458&i=0> > email: [hidden email]<http://user/SendEmail.jtp?type=node&node=4656458&i=1> > > ---Apoio aos softwares livres > www.zotero.org - gerenciamento de referências bibliográficas. > www.broffice.org ou www.libreoffice.org - textos, planilhas ou > apresentações. > www.epidata.dk - entrada de dados. > www.r-project.org - análise de dados. > www.ubuntu.com - sistema operacional > > > _______________________________________________ > R-br mailing list > [hidden email] <http://user/SendEmail.jtp?type=node&node=4656458&i=2> > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça > código mínimo reproduzível. > > ------------------------------ > If you reply to this email, your message will be added to the discussion > below: > > http://r-br.2285057.n4.nabble.com/R-br-Erro-em-multiplicacao-estouro-de-int… > To unsubscribe from R-br, click here<http://r-br.2285057.n4.nabble.com/template/NamlServlet.jtp?macro=unsubscrib…> > . > NAML<http://r-br.2285057.n4.nabble.com/template/NamlServlet.jtp?macro=macro_view…> > -- Jodavid Ferreira Responsável pela comissão de Comunicação do MJ da RCCPB - http://jovem.rccpb.com.br Membro da equipe Estadual do Ministério de Comunicação da RCCPB - http://rccpb.com.br Responsável pelo Ministério de Comunicação do Grupo de Oração Boa Esperança - http://grupoboaesperanca.com Bacharelando em Estatística - UFPB Estagiário no Laboratório de Estatística Aplicada ao Processamento de Imagens e Geoprocessamento - http://de.ufpb.br/~leapig Site Pessoal - http://jodavidferreira.com Twitter - @JodavidFerreira <http://twitter.com/JodavidFerreira>

1 0

[off-topic] Divulgação de Bolsa Pós Doutorado - Embrapa Arroz e Feijão
by Walmes Zeviani 08 Out '12

08 Out '12

Saudações membros da lista, Fui solicitado a divulgar na lista detalhes sobre uma bolsa de pós Doutorado cedida pela Embrapa Arroz e Feijão. Para não enviar o documento colei na integra o texto do mesmo. Contato no final da mensagem. --------------------------------------------------------------------------------------------- Oportunidade Pós Graduação O subprojeto “A methodological development of an online tool for the identification of Target Population of Environments: improving the predictions of agricultural production using crop models”, coordenado pela Embrapa Arroz e Feijão, e parte do projeto CCAFS (CGIAR Program on Climate Change, Agriculture and Food Security) disponibiliza uma bolsa de pós graduação para atuar na área de melhoramento vegetal/estatística. O objetivo desse subprojeto é definir estratégias para agrupar a região produtora de arroz de terras altas em megaambientes e posteriormente determinar ambientes alvo para populações (“Target Population of Environment – TPE”) com a finalidade de explorar adaptações genotípicas específicas. Descrição das atividades: a) Reanalisar a série histórica de dados dos ensaios de valor de cultivo e uso (VCU) da Embrapa Arroz e Feijão considerando co-variáveis climáticas; b) Desenvolver algoritmo para classificação dos ensaios (VCUs); c) Determinar genótipos que se destacaram em ambiente não ótimo. Conhecimentos requeridos: Conhecimento em programação no ambiente R; Conhecimento na aplicação de modelos mistos e modelo AMMI (Additive Main effects and Multiplicative Interaction); O interessado irá desenvolver o trabalho na Embrapa Arroz e Feijão, Goiás. Início em Janeiro de 2013. Valor da bolsa R$ 3.700,00 (Pós Doutorado). Duração de 1 ano e 10 meses. Interessados devem enviar CV para aos cuidados de alexbh <at> cnpaf.embrapa.br --------------------------------------------------------------------------------------------- À disposição. Walmes. ========================================================================== Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 VoIP: (3361 3600) 1053 1173 e-mail: walmes(a)ufpr.br skype: walmeszeviani twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================

1 0

Distância entre texto em uma legenda!
by Ivan Bezerra Allaman 08 Out '12

08 Out '12

Boa noite senhores! Já quebrei a cabeça para descobrir como eu faço para manipular a distância entre texto em uma legenda. Vejam, a distância entre o texto e um quadradinho (como no CMR abaixo) é manipulável com o argumento "inset", mais não é isso que quero. Gostaria de manipular a distância entre dois quadradinhos como no exemplo abaixo. Segue um CMR: plot(1:20,1:20,type='n',axes=F,xlab='',ylab='') legend(1,10,c("Red", "UFLA", "Commercial", "Chitralada"), col=rep(1,4), density=c(0,60,20,100), angle=c(45,45,45,0),cex=1,horiz=TRUE,xpd=TRUE) Desde já, grato! \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-5596 E-mail: ivanalaman@yahoo.com.br/ivanalaman@gmail.com @ \end{signature}

2 2

Re: [R-br] Dúvida função Anova pacote car - Medidas repetidas
by Ivan Bezerra Allaman 07 Out '12

07 Out '12

Acho uma grande bobagem querer usar modelos mistos para tudo! Que sentido há usar modelo aleatório se o objetivo não é estimar componentes da variância? Me parace que usar modelo misto é "xique", e não lembramos a essência de se utilizar tal método. \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-5596 E-mail: ivanalaman@yahoo.com.br/ivanalaman@gmail.com @ \end{signature}

3 3

Problema com download de imagens
by Lucas Cusinato 07 Out '12

07 Out '12

Boa tarde, Estou tentando montar uma base de imagens de câmeras de trânsito usando o site da prefeitura de Salvador. Encontrei uns endereços que são atualizados frequentemente com imagens de diferentes câmeras. Eu consigo baixar essas imagens sem problemas usando o navegador (firefox e chrome); mas, quando faço download usando R, o arquivo vem corrompido ou danificado (apesar de que o tamanho é sempre coerente - uns 22kb). Já testei em dois computadores com sistemas operacionais e ip's diferentes, mas dá na mesma sempre. O que está acontecendo? #endereço da imagem: http://www.transalvador.salvador.ba.gov.br/Cam/c02/snap_c1.jpg download.file(' http://www.transalvador.salvador.ba.gov.br/Cam/c02/snap_c1.jpg', 'teste.jpg') #mais uma: http://www.transalvador.salvador.ba.gov.br/Cam/c10/snap_c1.jpg download.file(' http://www.transalvador.salvador.ba.gov.br/Cam/c10/snap_c1.jpg', 'teste2.jpg') # o código inicial, que tem o mesmo problema: n <- 0 while (TRUE){ download.file(' http://www.transalvador.salvador.ba.gov.br/Cam/c06/snap_c1.jpg', paste(as.numeric(Sys.time()), '.jpg', sep='')) n <- n + 1 cat("\n", n, "\n\n") Sys.sleep(30) }

1 1

Re: [R-br] Diferença Pacotes multicore e parallel
by Junior Beleti 06 Out '12

06 Out '12

Olá Benilton, entendi. É que eu sempre utilizei o pacote multicore, mas em experimentos recentes com o pacote parallel obtive melhores desempenhos. Daí surgiu a dúvida. Obrigado, Junior.

1 0

Dúvida função Anova pacote car - Medidas repetidas
by Diego Pujoni 06 Out '12

06 Out '12

Olá pessoal, estou realizando uma ANOVA com medidas repetidas e estou utilizando a função "Anova" do pacote "car". Medi o biovolume de algas a cada dois dias durante 10 dias (no banco de dados abaixo só coloquei até o 4° dia). Tenho 2 tratamentos ("c","t") e o experimento foi realizado em tréplicas ("A","B","C"). > Pa2 Day Type Replicate logbiovolume 1 0 c A 19.34 2 0 c B 18.27 3 0 c C 18.56 4 0 t A 18.41 5 0 t B 18.68 6 0 t C 18.86 7 2 c A 18.81 8 2 c B 18.84 9 2 c C 18.52 10 2 t A 18.29 11 2 t B 17.91 12 2 t C 17.67 13 4 c A 19.16 14 4 c B 18.85 15 4 c C 19.36 16 4 t A 19.05 17 4 t B 19.09 18 4 t C 18.26 Pa2.teste = within(Pa2,{group = factor(Type) time = factor(Day) id = factor(Replicate)}) matrix = with(Pa2.teste,cbind(Pa2[,VAR][group=="c"],Pa2[,VAR][group=="t"])) matrix [,1] [,2] [1,] 19.34 18.41 [2,] 18.27 18.68 [3,] 18.56 18.86 [4,] 18.81 18.29 [5,] 18.84 17.91 [6,] 18.52 17.67 [7,] 19.16 19.05 [8,] 18.85 19.09 [9,] 19.36 18.26 [10,] 19.63 18.96 [11,] 19.94 18.06 [12,] 19.54 18.37 [13,] 19.98 17.96 [14,] 20.99 17.93 [15,] 20.45 17.74 [16,] 21.12 17.60 [17,] 21.66 17.33 [18,] 21.51 18.12 model <- lm(matrix ~ 1) design <- factor(c("c","t")) options(contrasts=c("contr.sum", "contr.poly")) aov <- Anova(model, idata=data.frame(design), idesign=~design, type="III") summary(aov, multivariate=F) Univariate Type III Repeated-Measures ANOVA Assuming Sphericity SS num Df Error SS den Df F Pr(>F) (Intercept) 12951.2 1 6.3312 17 34775.336 < 2.2e-16 *** design 19.1 1 17.3901 17 18.697 0.0004606 *** --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 O problema é que eu acho que esta função não está levando em consideração os dias, nem as réplicas. Como faço para introduzir isto na função. Vocês conhecem alguma função correspondente não paramétrica para este teste? Tipo um teste de Friedman com dois grupos (tratamento e réplica) e um bloco (tempo)? Pelo que entendi, a função friedman.test só faz com um grupo apenas. Muito Obrigado Diego PJ

4 3

Diferença Pacotes multicore e parallel
by Junior Beleti 05 Out '12

05 Out '12

Boa tarde a todos, gostaria de saber qual seria ou se existe ligação entre os pacotes "multicore" e "parallel". Por exemplo, posso utilizar os dois com a função mclapply, apenas instanciando-os: require(multicore) OU require(parallel), que consigo executar a função mclapply normalmente. Na documentação do multicore, vi que o pacote parallel tem a descrição de "Avaliar uma expressão de forma assíncrona em processos separados", e já o multicore possui a descrição de "pacote para processamento paralelo de código R". Qual seria a diferença dos dois pacotes no uso da função mclapply, por exemplo? Att, Junior.

2 1

Ordered probit binary
by Daniel Marcelino 05 Out '12

05 Out '12

Olá caros, li un texto esta semana em que se analizou dados binários como ordered probit e não simples probit. Para mim isto é novo, dei uma googlada mas não achei muita coisa interessante no R, agluem sabe alguma coisa sobre isso? Daniel -- -------------------------------------------------- View my research on these pages: http://www.crcee.umontreal.ca http://danielmarcelino.com http://ssrn.com/author=1640522 -------------------------------------------------- My contacts: Land Phone 1+514 343 6111 #3799 Skype: d.marcelino 3200 Jean Brillant, Office C5071 Montreal, QC; H3T 1N8 Canada ----------------------------------------------------- "Small steps toward a much better world"

4 4

Eclipse e R
by Emerson Cotta Bodevan 05 Out '12

05 Out '12

Olá alguém da lista usa o Eclipse como "editor" do R? Se sim, quais as vantagens de usá-lo, por exemplo, em relação ao EMACS? Att., Emerson.

1 0

Filtrar observações
by Carlos Mendonça 04 Out '12

04 Out '12

Prezados, tenho um arquivo tradutor com várias observações para um mesmo código. Preciso ficar só com a primeira observação de cada código desse arquivo tradutor. Como posso fazer isto? *Arquivo tradutor* codigo produto 1 valor3 > a 3 25 > a 5 10 > a 2 15 > b 15 85 > b 5 15 > c 2 5 > c 4 7 > c 2 5 > c 2 3 > *Resultado desejado * codigo produto 1 valor3 a 3 25 b 15 85 c 2 5 > > >

3 5

ccf (cross correlation function) entre duas time series com 'internal NAs'
by carlos pombo sonderblohm 04 Out '12

04 Out '12

Prezados colegas, estou a trabalhar com varias series temporais >library(astsa) e library(TSA); -no meu R version 2.13.1 (2011-07-08) no SO Windows7 64bit- tenho a seguinte dúvida, quando aplico a função ccf entre duas series, que devo fazer com os Na´s nesta função: 1.- na.action = na.pass ou 2.- na.action = na.contiguous estive a revisar a informação de ajuda da função ccf e encontrei isto : 'By default, no missing values are allowed. If the na.action function passes through missing values (as na.pass does), the covariances are computed from the complete cases. This means that the estimate computed may well not be a valid autocorrelation sequence, and may contain missing values. Missing values are not allowed when computing the PACF of a multivariate time series.' se comprendí bem a explicação, não devo usar na.pass devido a que estaria a estimar uma sequencia de correlações NÂO CONTINUA? com possiveis 'missing values' ou na.pass daria jeito porque estaria so a estimar aquelas que apresentan valores? e os lags vão a corresponder a cada intervalo de tempo na que se encontra minhas series (depasso seja dito são valores mensais)? Bom, adjunto aqui un ficheiro com varias duas series, OCTOPUS e PLOBO (ficheiro -> 'duvida.csv' , como podem ver, PLOBO apresenta internals NAs, eu fiz a ccf com cada na.action para ver o correlograma, e apresenta diferenças enormes! alguma dica? #scripts para visualizar o problema (duvida) OCT<- ts(duvida$OCTOPUS, start = c(1990,1), freq=12) PLOBO <- ts(duvida$PLOBO, start = c(1990,1), freq=12) ccf(PLOBO, OCTOPUS) ccf(PLOBO, OCTOPUS, na.action = na.pass) ccf(PLOBO, OCTOPUS, na.action = na.contiguous) muito obrigado Carlos -- Carlos A. Pombo Sonderblohm PhD Student on Marine Science (Fisheries) Faculdade de Ciências e Tecnología Universidade do Algarve, Campus de Gambelas 8005-139 Faro Portugal Tef. 289 800 905 ext. 7605

3 4

ggplot2 renomear strip elementos
by Alisson Lucrecio 04 Out '12

04 Out '12

Olá pessoal da r-br, Eu usava as funções ggplotGrob, grid.ls e getGrob para extrair os elementos do grid e renomear-los da forma que queria, mas o pacote ggplot2 teve uma alteração. Alguém sabe como fazer as alterações do strip.elementes com o ggplot2? Quero trocar os nomes LVAd, LVdf e PVAd para YRL, DRL e YRA, respectivamente. O link do arquivo segue abaixo: http://dl.dropbox.com/u/102629803/r%20br%20ajudar.rar Obrigado. Alisson Lucrécio da Costa

3 3

Literatura sobre análises de sobrevivência
by Cassiano 03 Out '12

03 Out '12

Olá, Gostaria de saber se alguém conhece algum livro texto, que seja bom e completo, para análises de sobrevivência usando o R. Att, Cassiano -- ======================================= Cassiano S. Rosa FaEnge - UEMG Av. Brasilia, nº 1304 - Bairro Baú, João Monlevade/MG CEP 35930-314 Tel:. 31-38593200 www.faenge.uemg.br Laboratório de Orthoptera Universidade Federal de Viçosa -UFV Av. Ph Rolfs s/n, Viçosa-MG 36570-000 Tel: 31-38992548 Curriculo lattes: http://lattes.cnpq.br/3030825329622014 =======================================

3 2