- R-br - DINF UFPR

gráfico de áreas
by D'Apolito 29 Mar '12

29 Mar '12

Olá, Alguém saberia informar uma função ou pacote que desenhe gráficos de áreas, como estes exemplos: a parte verde http://www.pnas.org/content/102/31/10904/F3.large.jpg a parte azul http://www.ncdc.noaa.gov/paleo/drought/images/grimpollen.jpg grato Carlos

3 3

Tabela Personalizada
by Sérgio Henrique almeida da silva ju 28 Mar '12

28 Mar '12

Olá Estou aqui mais uma vez para pedir ajuda aos meus amigos. Gostaria de gerar uma tabela como abaixo: Ano do óbito por uma doença X 2007 2008 2009 UF Rio de Janeiro Município Nova Iguaçu 10 9 11 Rio de janeiro 20 21 19 Petrópolis 15 16 14 São Paulo São Paulo 50 55 53 Santos 5 7 10 Guarulhos 12 10 9 Minas Gerais Belo Horizonte 25 22 20 Juiz de Fora 10 12 15 Ipatinga 9 10 11 Os números representam a contagem dos óbitos em cada ano Estou usando o comando: ftable(dados[c("uf", "municipio", "anoobt")]) Existe um comando melhor para isso? Como faço para exportar essa tabela? Abraços -- Sérgio Henrique Almeida da Silva Junior Doutorando em Epidemiologia em Saúde Pública Escola Nacional de Saúde Pública Sérgio Arouca - ENSP/FIOCRUZ http://lattes.cnpq.br/1611345552843383 Tel: (21) 94429486/78101651 id: 123*20942

3 3

duvida formatação de arquivo .csv em R
by Vitor Aguiar 28 Mar '12

28 Mar '12

Boa tarde, senhores, alguém pode por favor ajudar na dúvida abaixo? Tenho uma planilha de excel (.csv) com 30 colunas e dezenas de milhares de linhas. Estou tentando edita-la no excel mas como é muito grande o programa sempre trava, talvez culpa do meu computador com 2Gb de RAM. Como fazer isso usando R? Nas minhas 30 colunas, eu quero trabalhar com 2 a 2, ou seja, tenho 15 grupos de 2 colunas. Nessas duas colunas tenho valores numéricos como nesse exemplo: coluna 1 coluna 2 13 15 11 10 9 18 14 9 Quando você compara cada linha, as vezes eu tenho valores maiores na 2ª coluna, e as vezes o valor maior está na 1ª coluna e o menor na 2º. Eu queria transferir todos os valores menores para a primeira coluna e todos os maiores para a 2ª. E assim para todas as 30 colunas (15 grupos de 2 colunas). No excel dá pra fazer com: coluna 1 coluna 2 MIN(A1:B1) MAX(A1:B1) Desde já agradeço, Vitor Rezende da Costa Aguiar --------------------------------------------------- PhD student in Biotechnology Rede Nordeste de Biotecnologia Universidade Federal do Espírito Santo --------------------------------------------------- Current Adress: Department of Integrative Biology University of California, Berkeley 2033 Valley Life Sciences Building, office 4134 Berkeley, CA - USA 94720 Phone: 1 (510)-643-0060

3 9

Colocar Linha no Meio de um Gráfico no R
by Jodavid Ferreira 28 Mar '12

28 Mar '12

Estou precisando de algum comando que coloque uma linha reta na horizontal no meio de um gráfico. Tem esse exemplo aqui em baixo. O Comando lines, coloca uma linha, mas ela não fica na horizontal! x=seq(-5,5,0.5) y=x^2 plot(x,y,type='l',main="Teste",xlab='X', ylab='Y',col='red') lines(x,x+10,col='blue') -- Jodavid Ferreira Responsável pela comissão de Comunicação do MJ da RCCPB - http://jovem.rccpb.com.br Responsável pelo Ministério Jovem e 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

3 3

Cálculo de Volume de Sólido Irregular
by Tiago Vieira 27 Mar '12

27 Mar '12

Olá a todos, Você sabem me informar se há alguma função, já implementado no R, que calcule o volume de um sólido irregular. As dados que tenho é um grid (X, Y e Z) e uma variável indicadora que informa se o ponto pertence ou não ao sólido, resultados de uma krigagem indicativa. Att, Tiago Vieira.

3 6

Criar um banco selecionando casos
by Sérgio Henrique almeida da silva ju 27 Mar '12

27 Mar '12

Olá Tenho um banco com um conjunto de variáveis e gostaria de criar um outro banco atraves do primeiro. Exemplo: Id Sexo x y 1 M 1 0 2 M 0 1 3 F 0 1 4 F 1 0 5 M 0 0 6 F 0 0 7 M 0 0 8 F 0 1 9 F 0 0 10 M 1 1 Quero criar um outro banco se as variáveis do primeiro banco forem x ou y iguais a 1 No SPSS eu uso o seguinte comando, como faço para rodar no R? SELECT IF (x = 1 OR y = 1). EXECUTE. Obrigado mais uma vez!!! -- Sérgio Henrique Almeida da Silva Junior Doutorando em Epidemiologia em Saúde Pública Escola Nacional de Saúde Pública Sérgio Arouca - ENSP/FIOCRUZ http://lattes.cnpq.br/1611345552843383 Tel: (21) 94429486/78101651 id: 123*20942

5 6

Ajuda definir efeitos e contraste polinomial
by andre neves 27 Mar '12

27 Mar '12

Prezados, Ajustei o seguinte modelo: ResLme=lme (Differ~Medication,random=~1|Person/Period,data=Ch30pr11,method="REML") e gostaria de pedir ajuda para incluir na mesma tabela de análise da REML, os efeitos (linear, quadrático, etc) com os contrastes polinomiais como descrito no link: _https://stat.ethz.ch/pipermail/r-help/2006-January/085724.html _ Um abraço. André =========================== Aviso de confidencialidade: =========================== - Esta mensagem da Empresa Brasileira de Pesquisa Agropecuaria (Embrapa), empresa publica federal regida pelo disposto na Lei Federal n. 5.851, de 7 de dezembro de 1972, eh enviada exclusivamente a seu destinatario e pode conter informacoes confidenciais, protegidas por sigilo profissional. Sua utilizacao desautorizada eh ilegal e sujeita o infrator as penas da lei. Se voce a recebeu indevidamente, queira, por gentileza, reenvia-la ao emitente, esclarecendo o equivoco. ===================== Confidentiality note: ===================== - "This message from Empresa Brasileira de Pesquisa Agropecuaria (Embrapa), a government company established under Brazilian law (5.851/72), is directed exclusively to its addressee and may contain confidential data, protected under professional secrecy rules. Its unauthorized use is illegal and may subject the transgressor to the law's penalties. If you are not the addressee, please send it back, elucidating the failure".

2 1

Problemas FixedBreaks
by Julianna Trindade 27 Mar '12

27 Mar '12

Bom dia, Estou com problemas na hora de implementar meu intervalo de classes pra mapas coropleticos. O código utilizado é: brks <- classIntervals(dados$ano_2007, n = 7,style="fixed", fixedBreaks=c(-1, 0, 20, 40, 60, 80, 100)) e meu objetivo é que os números sejam agrupados assim: "Sem informação", "0", "]0,20]", "]20,40]", "]40,60]", "]60,80]","Acima de 80" Já utilizei IntervalClosure = right mas não dá certo. Alguém pode me ajudar a entender como o R determina os intervalos?? Se arredondando ou não? Alguém já implementou algo parecido e pode postar o exemplo?? -- Julianna Trindade "Exalta-te, Senhor, na Tua força! Nós cantaremos e louvaremos o Teu poder." Salmos 21; 13

1 0

Ajuda definir efeitos e contraste polinomial
by andre.neves＠cpatc.embrapa.br 27 Mar '12

27 Mar '12

Prezados, Ajustei o seguinte modelo: ResLme=lme (Differ~Medication,random=~1|Person/Period,data=Ch30pr11,method="REML") e gostaria de pedir ajuda para incluir na mesma tabela de análise da REML, os efeitos (linear, quadrático, etc) com os contrastes polinomiais como descrito no link: https://stat.ethz.ch/pipermail/r-help/2006-January/085724.html Um abraço. André =========================== Aviso de confidencialidade: =========================== - Esta mensagem da Empresa Brasileira de Pesquisa Agropecuaria (Embrapa), empresa publica federal regida pelo disposto na Lei Federal n. 5.851, de 7 de dezembro de 1972, eh enviada exclusivamente a seu destinatario e pode conter informacoes confidenciais, protegidas por sigilo profissional. Sua utilizacao desautorizada eh ilegal e sujeita o infrator as penas da lei. Se voce a recebeu indevidamente, queira, por gentileza, reenvia-la ao emitente, esclarecendo o equivoco. ===================== Confidentiality note: ===================== - "This message from Empresa Brasileira de Pesquisa Agropecuaria (Embrapa), a government company established under Brazilian law (5.851/72), is directed exclusively to its addressee and may contain confidential data, protected under professional secrecy rules. Its unauthorized use is illegal and may subject the transgressor to the law's penalties. If you are not the addressee, please send it back, elucidating the failure".

1 0

linas repetidas após unir duas data.frames com comando merge
by Lucas Fernandes de Souza 27 Mar '12

27 Mar '12

A questão da memória, resolvi dividindo a planilha por tratamentos. Porém, quando faço o merge, reparo que ocorre repetição de várias linhas. Abaixo um exemplo das planilhas que quero unir: p1 - dim: 3326 10 TRAT RUN year DOY das dap laid yield CWAD DTTD 14144 5 125 1979 32 216 0 0 0 0 14.78 14145 5 125 1979 33 217 1 0 0 0 15.94 14146 5 125 1979 34 218 2 0 0 0 16.24 14147 5 125 1979 35 219 3 0 0 0 17.57 14148 5 125 1979 36 220 4 0 0 0 15.94 14149 5 125 1979 37 221 5 0 0 0 17.47 p2 - dim: 9999 6 RUN PRED SRAD PARD TMXD TMND 55840 125 0.0 16.2 32.4 26.6 6.7 55841 125 0.0 15.7 31.4 25.9 7.2 55842 125 0.2 14.6 29.2 25.9 9.9 55843 125 0.0 15.2 30.4 27.2 9.9 55844 125 0.0 15.5 31.0 27.4 9.6 55845 125 0.0 15.3 30.6 27.1 9.7 Como resultado deveria vir uma nova planilha com dim: 3326 15, porém, o resultado é uma planilha gigantesca, com várias linhas repetidas: p3 - dim: 1073748 15 RUN TRAT year DOY das dap laid yield CWAD DTTD PRED SRAD PARD TMXD TMND 1 125 5 1979 32 216 0 0 0 0 14.78 0.0 16.2 32.4 26.6 6.7 2 125 5 1979 32 216 0 0 0 0 14.78 0.0 15.7 31.4 25.9 7.2 3 125 5 1979 32 216 0 0 0 0 14.78 0.2 14.6 29.2 25.9 9.9 4 125 5 1979 32 216 0 0 0 0 14.78 0.0 15.2 30.4 27.2 9.9 5 125 5 1979 32 216 0 0 0 0 14.78 0.0 15.5 31.0 27.4 9.6 6 125 5 1979 32 216 0 0 0 0 14.78 0.0 15.3 30.6 27.1 9.7 Alguma sugestão? obrigado Lucas -- Lucas F. de Souza Eng. Agrícola - CONAB M.Sc. Produção Vegetal - UENF Doutorando Engenharia de Sistemas Agrícolas - ESALQ/USP tel: (19) 8160-7068

2 1

Ajuda com criação de loop para transforma lista em matrix de presença e ausencia
by Augusto Ribas 26 Mar '12

26 Mar '12

Ola pessoal. Estou aqui denovo pra pedir uma ajudinha. Eu estou com uns dataframes com 2 colunas, uma são nome de especies de hospedeiro e outra nome de especie de parasita. Ai a lista é na verdade parasita tal ocorreu em hospedeiro fulano Até aqui tudo bem, agora eu empaquei pq queria transformar isso em uma matriz de presensa e ausencia, com colunas sendo parasitos e linhas hospedeiro, ai onde tem interação entre os 2 fica preenchido com 1 e onde nao tem com 0. Eu fiz um exemplo com uma matriz pequenininha e como eu estou tentando resolver. Pensei em usar um teste de logica pra ir preenchendo a matriz mas eu não sei como fazer isso de modo eficiente. Pq o maximo que cheguei foi testar cada celula da matriz com toda a lista, e como a lista é meio grande é meio perdido fazer assim. Mas segue o exemplo que acho que ta facil entender, se alguem tiver uma solução facil na manga faria minha semana muito feliz :) #------------------------------------------------------------------------------------------------- #exemplo de dados dados<-data.frame(parasito=c("a","b","b","c","d","d"), hospedeiro=c("1","1","1","2","2","3")) #como eles estão agora dados #Eu preciso de uma matriz de presença e ausencia #dai faço uma matriz matriz<-matrix(data=NA,nrow=length(unique(dados$parasito)), ncol=length(unique(dados$hospedeiro)), dimnames=list(c(as.character(unique(dados$parasito))), c(as.character(unique(dados$hospedeiro))))) #blz agora preciso que a lista em dados seja representada #atraves dessa matriz matriz #Eu imaginei um teste assim #Linha 1 coluna 1 matriz[1,1]<-ifelse(colnames(matriz)[1]==dados$hospedeiro[1] & rownames(matriz)[1]==dados$parasito[1] ,1,0) #aqui eu testei se a primeira linha do data.frame corresponde a essa interação, como corresponde ele coloca 1 matriz #linha 1 Coluna 2 matriz[1,2]<-ifelse(colnames(matriz)[2]==dados$hospedeiro[1] & rownames(matriz)[1]==dados$parasito[1] ,1,0) #o mesmo pra linha 1 coluna 2, ele coloca 0 pq a a primeira linha do data.frame não corresponde a essa interação matriz #Blz, mas eu não consegui pensar em um jeito de faze ele rodar pra #toda a matriz que eu tenho de forma direta, eu queria monta um loop #ou algo do tipo que eu so fala-se o tamanho da matriz e ele testa-se #logo tudo, ou melhor acho que seria faze algo que olhe a lista #e teste pra cada linha da lista e ponha um 1 la na matriz Obrigado e boa semana a todos. -- Grato Augusto C. A. Ribas Site Pessoal: http://augustoribas.heliohost.org Lattes: http://lattes.cnpq.br/7355685961127056

1 1

problema com for loop
by Vitor Aguiar 26 Mar '12

26 Mar '12

Boa noite a todos. Sou novo em programação e estou tentando montar um for loop, então gostaria de pedir ajuda aos que estão familiarizados com isso. A explicação abaixo é longa, mas há poucos códigos de programação. Desde já agradeço aos que puderem ajudar. Isso é muito importante pra mim. ---- # Salve o arquivo "www.datafilehost.com/download-ce4b3fc2.html" no seu diretório de trabalho e leia-o da seguinte forma: > Pop = read.csv(file = "nist.csv", row.names = 1) # Pop = meu conjunto de dados (dataset). Cada 2 colunas é uma unidade de análise, ou seja, preciso analisar os dados de duas em duas colunas. Pode-se observar que as colunas de uma unidade tem o mesmo nome, diferenciadas por .1 e .2 # Nomes = níveis de dados presentes na 1ª unidade (1ª e 2ª colunas). Vou usar isso para nomear as linhas e colunas da matriz que vem adiante. > Nomes = levels(factor(c(Pop[,1], Pop[,2]))) # matrix ObsV = Valores observados para a 1ª unidade. Aqui, estou criando uma matriz de "0" (zeros) e depois usando um loop para preenche-la com os valores do meu dataset. > ObsV = matrix(0, nrow = length(Nomes), ncol = length(Nomes), dimnames = list(Nomes, Nomes)) > for(i in 1:nrow(Pop[,1:2])) { ObsV[paste(Pop[i,1]), paste(Pop[i,2])] = ObsV[paste(Pop[i,1]), paste(Pop[i,2])] + 1 } > ObsV # Muito simples: Leia a 1ª unidade de dados do dataset: quando aparecem na 1ª linha os valores 13 - 14 é somado +1 no spot "linha 13, coluna 14" da matriz ObsV. E assim por diante contando todos os valores observados. Então, se no spot "linha 10, coluna 14" da matrix ObsV eu tenho o valor 11, isso é porque o valor 10 - 14 apareceu 11 vezes no meu dataset. # Isso funciona bem pra mim. Mas eu preciso fazer isso funcionar para todo o dataset de uma vez só, sempre trabalhando de duas em duas colunas. Ou seja, como tenho 30 colunas, preciso produzir 15 matrizes como a ObsV. # Eu pensei em por esse loop dentro de outro loop como for(j in ncol(Pop)) {}, mas não to conseguindo escrever isso e também preciso saber como modificar minhas variáveis, como "N" e "Nomes" para que funcionem para todos os marcadores. #todos os códigos acima podem ser substituídos pelos códigos abaixo, produzindo os mesmo resultados. Se for mais fácil resolver assim. > Pop = read.csv("nist.csv", row.names=1) > Nomes = levels(factor(c(pop[,1], pop[,2]))) > ObsV = table(factor(Pop[,1], levels = Nomes), factor(Pop[,2], levels = Nomes)) > ObsV -------- Obrigado, Vitor Rezende da Costa Aguiar --------------------------------------------------- PhD student in Biotechnology Rede Nordeste de Biotecnologia Universidade Federal do Espírito Santo --------------------------------------------------- Current Adress: Department of Integrative Biology University of California, Berkeley 2033 Valley Life Sciences Building, office 4134 Berkeley, CA - USA 94720 Phone: 1 (510)-643-0060

2 8

Definir efeitos
by andre neves 26 Mar '12

26 Mar '12

Prezados, Bom dia! Estou começando a utilizar o R, já fiz algumas análises e por isso gostaria de parabenizar a todos os coloboradores desta lista pela disponibilidade em responder às perguntas. Pergunta: Como definir os efeitos linear (L), quadrático (Q), cúbico (C) e quártico (Qt) em modelos mistos (melhor ajuste), considerando 5% como nível de significância. Quais comandos utilizar? Um abraço. André =========================== Aviso de confidencialidade: =========================== - Esta mensagem da Empresa Brasileira de Pesquisa Agropecuaria (Embrapa), empresa publica federal regida pelo disposto na Lei Federal n. 5.851, de 7 de dezembro de 1972, eh enviada exclusivamente a seu destinatario e pode conter informacoes confidenciais, protegidas por sigilo profissional. Sua utilizacao desautorizada eh ilegal e sujeita o infrator as penas da lei. Se voce a recebeu indevidamente, queira, por gentileza, reenvia-la ao emitente, esclarecendo o equivoco. ===================== Confidentiality note: ===================== - "This message from Empresa Brasileira de Pesquisa Agropecuaria (Embrapa), a government company established under Brazilian law (5.851/72), is directed exclusively to its addressee and may contain confidential data, protected under professional secrecy rules. Its unauthorized use is illegal and may subject the transgressor to the law's penalties. If you are not the addressee, please send it back, elucidating the failure".

2 3

Análise de covariância
by Gisele Correia 26 Mar '12

26 Mar '12

Bom dia, Alguém teria um exemplo de análise de covariância no R para eu ver como se aplica? Att., -- Gisele Batista Correia Engenheira de Pesca Mestranda em Ciências Pesqueiras nos Trópicos Universidade Federal do Amazonas - UFAM Av. Gen. Rodrigo Octávio Jordão Ramos, 3000, Campus Universitário. Bairro Coroado I. CEP 69077-000. Manaus/AM. Fone:(92) 3305-4063 Celular: (41) 9818-6424 (TIM) (92) 8142-3352 (TIM) (92) 9322-5668 (VIVO) ‎"Agradeço por todas as dificuldades que enfrentei (e ainda enfrento)... se não fosse por elas, não teria saído do lugar. A facilidade nos impede de caminhar. Mesmo as críticas nos auxiliam muito..." Chico Xavier

3 2

Re: [R-br] problema memória
by Benilton Carvalho 26 Mar '12

26 Mar '12

Lucas, por favor, nao mando anexos para a lista... e muito menos para enderecos pessoais sem previa confirmacao (ao ler minhas msgs em dispositivos moveis, eu nem sempre tenho a opcao de nao baixar os anexos, que no seu caso resumem-se a 10MB). as sugestoes que enviei anteriormente continuam validas: 64 bits e mais RAM. boa sorte, b 2012/3/26 Lucas Fernandes de Souza <luc.souza(a)usp.br>: > tenho o windows 7 de 32 bits e 3G de memória ram. > > Segue em anexo as duas planilhas que quero unir. > Já tentei separar por tratamentos, ou seja, primeiro escolho um tratamento e > faço o merge. Com isso, já consegui fazer o merge, porém, os dados aparecem > repetidos. Como coluna em comum (para o merge) já usei o "year", "doy" e > continua aparecendo linhas repetidas... > na verdade, o que me interessa é o "comprimento" da planilha "growth", ou > seja, tenho que adequar os dados da planilha "weatherobs" com o da planilha > "growth". > Não consegui fazer um scrip que resolva isso...se puderem me ajudar, > agradeço. > > > > Citando Benilton Carvalho <beniltoncarvalho(a)gmail.com>: > >> As sugestoes, ja listadas aqui algumas vezes, sao: >> >> 0) De um exemplo reproduzivel >> 1) Instale um SO de 64 bits e o R com 64 bits >> 2) Aumente a RAM disponivel no seu computador >> >> No entremeio, vc pode tentar trabalhar apenas com segmentos do seu >> arquivo... ou ate' jogar tudo num banco de dados SQL e fazer o merge >> por la... (que muito possivelmente se resumira a um INNER JOIN). >> >> O que eh necessario para sugestoes mais uteis eh que vc descreva com >> detalhe o volume de dados e sua estrutura (novamente, encapsulada num >> exemplo reproduzivel) e tambem nos informe informacoes sobre o >> equipamento que voce usa, incluindo RAM e sistema operacional. >> >> b >> _______________________________________________ >> R-br mailing list >> R-br(a)listas.c3sl.ufpr.br >> 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. >> > > > > -- > Lucas F. de Souza > Eng. Agrícola - CONAB > M.Sc. Produção Vegetal - UENF > Doutorando Engenharia de Sistemas Agrícolas - ESALQ/USP > tel: (19) 8160-7068

2 3

problema memória
by Lucas Fernandes de Souza 26 Mar '12

26 Mar '12

Bom dia a todos. Tentei unir duas planilhas, com o comando merge, porém aparece o seguinte erro: Error: cannot allocate vector of size 16.0 Mb In addition: Warning messages: 1: In attr(value, "row.names") <- row.names : Reached total allocation of 1535Mb: see help(memory.size) 2: In attr(value, "row.names") <- row.names : Reached total allocation of 1535Mb: see help(memory.size) 3: In attr(value, "row.names") <- row.names : Reached total allocation of 1535Mb: see help(memory.size) 4: In attr(value, "row.names") <- row.names : Reached total allocation of 1535Mb: see help(memory.size) O que posso fazer para solucionar isso? obrigado -- Lucas F. de Souza Eng. Agrícola - CONAB M.Sc. Produção Vegetal - UENF Doutorando Engenharia de Sistemas Agrícolas - ESALQ/USP tel: (19) 8160-7068

2 2

Re: [R-br] Ajuda com os benditos contrastes.
by Walmes Zeviani 26 Mar '12

26 Mar '12

Olympio, Não é possível fazer testes para níveis de Trat1 dentro de Trat2 com a formula assim ~Trat1/Trat2, isso por causa da codificação das matrizes envolvidas. Você tem que usar Trat2/Trat1, ou Tra3*Trat2 (ou inverso), ou ainda criar um novo fator combinando os níveis desses que pode ser conseguido com a interaction(). No seu caso, você deve ajustar o modelo com Trat2/Tempo/Trat1, o fator de níveis a serem comparados deve ser o último na fórmula do modelo. Se for usar a contrast::contrast() o modelo tem que ter estrutura cruzada e não aninhada. Veja, Id: Unidade Amostral #Estrutura dos dados: set.seed(13032012) d=expand.grid(Trat1=gl(3,18),Tempo=gl(2,1)) d$Trat2=gl(2,9,108) d$Id=as.factor(rep(1:54,2)) d$Y=abs(rnorm(108,c(10,30,8),10)) str(d) #Verificar estrutura para a unidade amostral i=3, por exemplo d[d$Id==3,]: Tratamento1 =1, Tratamento2=1 Tempo=1, 2. #Ajuste que fiz usando lme library(nlme) f=lme(Y~(Trat1/Trat2)*Tempo, random=~1|Id,d) summary(f) m0 <- lme(Y~Trat2/Tempo/Trat1, random=~1|Id,d) fixef(m0) anova(m0) require(aod) # 8 parâmetros = 8 g.l. da interação grep("*.:Tempo*.:Trat1*.", names(fixef(m0)), value=TRUE) idx <- grep("*.:Tempo*.:Trat1*.", names(fixef(m0))) # última linha da anova(m0) wald.test(Sigma=vcov(m0), b=fixef(m0), Terms=idx, df=nrow(d)-length(fixef(m0))) # assim é possível fazer testes de desdobramento de interação usando a wald.test() # é só saber que índices passar para Terms= para ter as hipóteses adequadas require(contrast) levels(d$Trat2) levels(d$Trat1) levels(d$Tempo) c0 <- contrast(m0, list(Trat2="1", Tempo="1", Trat1="1"), list(Trat2="1", Tempo="1", Trat1="2")) # modelo tem que estar na forma cruzada (*) e não aninhada (/) m0 <- lme(Y~Trat2*Tempo*Trat1, random=~1|Id,d) c0 <- contrast(m0, list(Trat2="1", Tempo="1", Trat1="1"), list(Trat2="1", Tempo="1", Trat1="2")) c0 c0$X # vetor de contrastes Depois é só você trocar os niveis na contrast() para demais comparações. Você pode reunir os vetores de contraste e fazer esses contrastes com a multcomp::glht() para ajustar o p-valor dos testes, visto que são muitas hipóteses não ortogonais. Consulte a documentação dessas funções. Aconselho você a sempre enviar perguntas para a R-br pois o sucesso de resposta é maior, sua dúvida pode ajudar os demais além de ficar no arquivo da lista para consulta futura. À 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 twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ========================================================================== 2012/3/25 Olympio Neto <olympiotneto(a)gmail.com> > Walmes, muito obrigado por você estar sempre disponível a ajudar... Olha, > eu deveria te pagar, porque já me ajudou milhares de vezes com suas > intervenções. > > Então, cara, eu tenho um problema. Vou colocar o que fiz aqui, usando > aquele CMR que mandei. > > Id: Unidade Amostral > #Estrutura dos dados: > set.seed(13032012) > d=expand.grid(Trat1=gl(3,18),Tempo=gl(2,1)) > d$Trat2=gl(2,9,108) > d$Id=as.factor(rep(1:54,2)) > d$Y=abs(rnorm(108,c(10,30,8),10)) > > #Verificar estrutura para a unidade amostral i=3, por exemplo > d[d$Id==3,]: Tratamento1 =1, Tratamento2=1 Tempo=1, 2. > > #Ajuste que fiz usando lme > > library(nlme) > f=lme(Y~(Trat1/Trat2)*Tempo, random=~1|Id,d) > > Veja, e "f" é meu modelo final, o melhor e tal. Acontece que quero fazer > comparações múltiplas considerando somente Trat1 e Tempo, desprezando Trat2 > dentro de Trat1. Estou me matando há 3 dias tentando isso e não consigo. > Quando coloco Trat2, ela entra nos contrastes e não quero isso. > Queria fazer isso, mas nem sei como te mostrar, desculpe, mas vai > literalmente. > > Trat1:1 X Tempo:1 - Trat1:2 X Tempo:1 > Trat1:1 X Tempo: 1- Trat1:3 X Tempo:1 > Trat1:2 X Tempo:1 - Trat1:3 X Tempo:1 > > Eu vi um exemplo parecido em > http://www.ats.ucla.edu/stat/R/faq/testing_contrasts.htm, mas lá cria-se > um novo modelo univariado com o comando interaction nos fatores. Poderia > até fazer isso, mas isso iria alterar muito meu modelo, que já é > complicado. Você vislumbra alguma forma de se resolver isso? Não queria > perder a caracterização do aninhamento de Trat2 dentro de Trat1. Já tentei > contrast, multcomp e tudo mais, mas sempre solicitam os níveis de saliva ou > se coloco, eles aumentam as comparações. > > Desculpe te encher com esta dúvida idiota, mas os contrastes no R me > arrebentam, e há anos. Já li um zilhão de livros do R, mas quando preciso > ir às vias de fato, me mato. > > Abs e muito obrigado. > > Olympio >

1 0

chuvapr
by Gilbert Queiroz 25 Mar '12

25 Mar '12

"chuvapr" alguém tem este bd, para pratcar no aRT???? Abs.

3 2

Como ler dados .DAT e CTL
by Victor Eduardo 25 Mar '12

25 Mar '12

Pessoal, estou utilizando os dados do INPE e nele, vem em formato .DAT e um outro com o formato CTL. Como poderia ler no R esses dois formatos? Abraços, Victor Eduardo

6 6

Ajuda com particionamento de soma de quadrados usando nlme
by Olympio Neto 25 Mar '12

25 Mar '12

Pessoal, recorro novamente a vocês para me ajudarem em um problema que já foi objeto de outra mensagem que postei aqui. Estou com um modelo de ANOVA com medidas repetidas em que um dos tratamentos é aninhado a outro. O que gostaria de fazer é particionar o tratamento aninhado como se faz com o comando aov (..., split..) (Vi exemplo no site R-idiculas (excelente, por sinal)). A estrutura dos meus dados é a seguinte: Id: Unidade Amostral #Estrutura dos dados: set.seed(13032012) d=expand.grid(Trat1=gl(3,18),Tempo=gl(2,1)) d$Trat2=gl(2,9,108) d$Id=as.factor(rep(1:54,2)) d$Y=abs(rnorm(108,c(10,30,8),10)) #Verificar estrutura para a unidade amostral i=3, por exemplo d[d$Id==3,]: Tratamento1 =1, Tratamento2=1 Tempo=1, 2. #Ajuste que fiz usando lme library(nlme) f=lme(Y~(Trat1/Trat2)*Tempo, random=~1|Id,d) #Gostaria de realizar algo do tipo, só que com o modelo misto usando o pacote nlme m1=aov(Y~(Trat1/Trat2)*Tempo,d) names(coef(m1))[5:7] grep("Trat11", names(coef(m1))[5:7]) grep("Trat12", names(coef(m1))[5:7]) grep("Trat13", names(coef(m1))[5:7]) ### As partições, mas sem considerar o modelo misto (queria algo semelhante). summary(m1) summary(m1, split=list("Trat1:Trat2"=list( " Trat1-1"=c(1), "Trat2-2"=c(2), "Trat3-3"=c(2) ))) Fico na dúvida, se devo fazer isso na interação tripla, já que nos meus dados reais foi significativa ou se os coeficientes estimados delas também devam entrar no contraste, enfim, eu sempre tenho muitas dúvidas em contrastes no R. Alguém sabe como posso fazer isso? Já tentei os contrastes com multcomp, mas não é isso que quero... gostaria de ter partição da soma dos quadrados, semelhante ao exemplo com o aov. Alguém pode me ajudar? Desde já, muito obrigado. Olympio

2 3

Teste
by Gilbert Queiroz 25 Mar '12

25 Mar '12

Isto aqui é um teste, p/ contato c/ a lista...

1 0

Conectando o R ao TerraView
by Gilbert Queiroz 24 Mar '12

24 Mar '12

Amigos, Estou trabalhando com dados no TerraView. Estava pensando em usar Regessão Espacial. Para isso, iria exportar os dados p/ o R, e trabalhar independentemente, e depois voltar com os resultados prevsitos pelo modelo para o TerraView. Mas encontrei material sobre o aRT. Porém, não estou sabendo fazer a conexão entre os programas. Alguém pode me ajudar nesta tarefa? Abs. Gilbert.

1 0

Teste de comparação de médias e contrastes ortogonais no R
by andre neves 23 Mar '12

23 Mar '12

Prezados, Bom tarde! Gostaria de pedir uma ajuda sobre como rodar (comandos) comparação de médias pelo teste Tukey , quando considero os efeitos fixos e aleatórios, o que constitui um modelo misto. E como faço para rodar (comandos) contrastes ortogonais para definir os efeitos linear (L), quadrático (Q), cúbico (C) e quártico (Qt), considerando 5% como nível de significância. Um abraço a todos. André =========================== Aviso de confidencialidade: =========================== - Esta mensagem da Empresa Brasileira de Pesquisa Agropecuaria (Embrapa), empresa publica federal regida pelo disposto na Lei Federal n. 5.851, de 7 de dezembro de 1972, eh enviada exclusivamente a seu destinatario e pode conter informacoes confidenciais, protegidas por sigilo profissional. Sua utilizacao desautorizada eh ilegal e sujeita o infrator as penas da lei. Se voce a recebeu indevidamente, queira, por gentileza, reenvia-la ao emitente, esclarecendo o equivoco. ===================== Confidentiality note: ===================== - "This message from Empresa Brasileira de Pesquisa Agropecuaria (Embrapa), a government company established under Brazilian law (5.851/72), is directed exclusively to its addressee and may contain confidential data, protected under professional secrecy rules. Its unauthorized use is illegal and may subject the transgressor to the law's penalties. If you are not the addressee, please send it back, elucidating the failure".

2 1

Re: [R-br] Criação de função fortran em pacote (base)
by Junior Beleti 23 Mar '12

23 Mar '12

Olá Benilton, Primeiramente gostaria de agradecer novamente pela paciência e pela atenção, e me desculpar por não ter conhecimento necessário sobre funções e procedimentos estatísticos. O foco da minha pesquisa é o processamento paralelo. Estou utilizando scripts, (tais como de obtenção dos valores preditos utilizando os dados do pivo), pois foi constatado problemas quanto ao tempo de processamento de algumas rotinas desses scripts. Na rotina de obtenção dos valores preditos utilizando os dados do pivo (trecho a seguir), o problema quanto a tempo de processamento está na chamada "preditos.simu <- volta.cokri(md.cov.ck,num.simu=1000,int.conf=0.95)". (lembrando que tais scripts necessitam do pacote geoComp) ------------------------------------------------------------------- require(geoComp) require(MASS) require(statmod) require(geoR) data(pivo) dados <- pivo[,c(6,7,8,1,2)] dados <- as.geoComp(dados) bor <- cbind(c(0,seq(0,200,l=100),0),c(0,sqrt(200^2-seq(0,200,l=100)^2),0)) estima <- mec(dados) gr <- pred_grid(bor, by=4) md.cov.ck <- cokrigagem(estima[[1]]$Estimativas, loc=gr,dados.comp=dados) preditos.gh <- volta.quad(md.cov.ck,n.pontos=7,Variancia=FALSE) preditos.gh <- data.frame(preditos.gh) write.table(preditos.gh,"pred_by4k7ns1000MBM.txt") preditos.simu <- volta.cokri(md.cov.ck,num.simu=1000,int.conf=0.95) #Ponto com problema quanto a processamento preditos.simu <- data.frame(preditos.simu[[1]]) preditos.simu.ic <- data.frame(preditos.simu[[2]]) write.table(preditos.simu,"predsimu_by4k7ns1000MBM.txt") write.table(preditos.simu.ic,"predsimuic_by4k7ns1000MBM.txt") ------------------------------------------------------------------- Visando localizar o ponto em que a execução gasta um maior tempo de processamento, "entrei" na função volta.cokri(código a seguir) e fui investigar a causa da "demora" na execução dessa função. ------------------------------------------------------------------- volta.cokri <- function(mat.cokri, num.simu, retorna.tudo=FALSE, int.conf=0.95){ nlinhas <- dim(mat.cokri[[1]])[1] compos1 <- data.frame(matrix(nrow=nlinhas/2,ncol=3)) compos <- data.frame(matrix(nrow=nlinhas/2,ncol=3)) for(i in 1:num.simu){ g <- mvrnorm(n=1, mat.cokri[[1]],mat.cokri[[2]]) seq1 <- seq(1,nlinhas,by=2) seq2 <- seq(2,nlinhas,by=2) y1 <- g[seq1] y2 <- g[seq2] gerado <- data.frame(y1,y2) compos <- agl(gerado) compos1 <- cbind(compos,compos1) } compos2 <- as.matrix(compos1) dim.vetor <- num.simu*3 sy1 <- seq(1,dim.vetor,by=3) sy2 <- seq(2,dim.vetor,by=3) sy3 <- seq(3,dim.vetor,by=3) amostra1 <- as.matrix(compos2[,sy1],ncol=num.simu) amostra2 <- as.matrix(compos2[,sy2],ncol=num.simu) amostra3 <- as.matrix(compos2[,sy3],ncol=num.simu) if(retorna.tudo == TRUE){ retorna <- list() retorna[[1]] <- amostra1 retorna[[2]] <- amostra2 retorna[[3]] <- amostra3 return(retorna) } if(retorna.tudo == FALSE){ med1 <- apply(amostra1,1,mean) med2 <- apply(amostra2,1,mean) med3 <- apply(amostra3,1,mean) q1 <- t(apply(amostra1,1,quantile,prob=c(1-int.conf,int.conf))) q2 <- t(apply(amostra2,1,quantile,prob=c(1-int.conf,int.conf))) q3 <- t(apply(amostra3,1,quantile,prob=c(1-int.conf,int.conf))) quantis <- cbind(q1,q2,q3) resultado <- list() resultado$preditos <- data.frame(med1,med2,med3) names(resultado$preditos) <- c("Areia","Silte","Argila") resultado$intervalo <- data.frame(quantis) names(resultado$intervalo) <- c("LI Areia", "LS Areia", "LI Silte", "LS Silte", "LI Argila", "LS Argila") return(resultado)} } ------------------------------------------------------------------- Na função volta.cokri verifiquei que o elevado tempo de processamento se concentrava na chamada "g <- mvrnorm(n=1, mat.cokri[[1]],mat.cokri[[2]])". Então desmembrei a função para realizar experimentos com pacotes paralelos afim de criar modelos de execução paralelos. A seguir está o trecho da função desmembrada. ------------------------------------------------------------------- mat.cokri <- md.cov.ck retorna.tudo=FALSE int.conf=0.95 nlinhas <- dim(mat.cokri[[1]])[1] compos1 <- data.frame(matrix(nrow=nlinhas/2,ncol=3)) compos <- data.frame(matrix(nrow=nlinhas/2,ncol=3)) mvrnorm <- function(n = 1, mu, Sigma, tol=1e-6, empirical = FALSE) { p <- length(mu) if(!all(dim(Sigma) == c(p,p))) stop("incompatible arguments") eS <- eigen(Sigma, symmetric = TRUE, EISPACK = TRUE) ev <- eS$values if(!all(ev >= -tol*abs(ev[1L]))) ev <- ev *(-1) X <- matrix(rnorm(p * n), n) if(empirical) { X <- scale(X, TRUE, FALSE) X <- X %*% svd(X, nu = 0)$v X <- scale(X, FALSE, TRUE) } X <- drop(mu) + eS$vectors %*% diag(sqrt(pmax(ev, 0)), p) %*% t(X) nm <- names(mu) if(is.null(nm) && !is.null(dn <- dimnames(Sigma))) nm <- dn[[1L]] dimnames(X) <- list(nm, NULL) if(n == 1) drop(X) else t(X) } m1 <- mat.cokri[[2]][1:2601,1:2601] #divisao da matriz original m2 <- mat.cokri[[2]][2602:5202,1:2601] m3 <- mat.cokri[[2]][1:2601,2602:5202] m4 <- mat.cokri[[2]][2602:5202,2602:5202] #chamada a funcao mclapply (neste exemplo esta com 1:4, mas que na execução real seria 1:1000) mclapply(1:4, function(i) { g1 <- mvrnorm(n=1, mat.cokri[[1]][1:2601],m1) g2 <- mvrnorm(n=1, mat.cokri[[1]][1:2601],m2) g3 <- mvrnorm(n=1, mat.cokri[[1]][1:2601],m3) g4 <- mvrnorm(n=1, mat.cokri[[1]][1:2601],m4) list(g1, g2, g3, g4) },mc.cores=4 ) g <- list(list(1:2601, g1), list(1:2601, g2), list(1:2601, g3),list(1:2601, g4)) seq1 <- seq(1,nlinhas,by=2) seq2 <- seq(2,nlinhas,by=2) y1 <- g[seq1] y2 <- g[seq2] gerado <- data.frame(y1,y2) compos <- agl(gerado) compos1 <- cbind(compos,compos1) compos2 <- as.matrix(compos1) ------------------------------------------------------------------- Esses trechos de código podem reproduzir o problema. Se não tiver o pacote geoComp, posso lhe enviar por email. Meu email é beleti.junior(a)gmail.com Sobre a possivel solução apresentada, consegui aplica-lá. Mais uma vez agradeço a sua disposição e peço desculpas pelo encomodo. Att, Carlos Beleti.

2 3

Posição de um valor em uma variável
by Sérgio Henrique almeida da silva ju 23 Mar '12

23 Mar '12

Olá Estou trabalhando com um banco de mortalidade (SIM) no R e gostaria de uma ajuda. Quero criar uma variável que me diz a posição de um valor em uma variável Exemplo: Banco Original linhaa linhab *T792 *X990 *S069 *X954 *R092*T792 *J189 *R571 *T792 *A301 *T792 *X990*A302 Criando a variável x1 e x2 linhaa linhab x1 x2 *T792 *X990 0 0 *S069 *X954 0 0 *R092*T792 *J189 0 0 *R571 *T792 0 0 *A301 2 0 *T792 *X990*A302 0 7 Ou seja, o valor A30, começa na segunda linha da variável linhaa e na 7 linha da variável linhab Como rodar isso no R, no SPSS o comando é: compute x1 = index(linhaa, "A30"). compute x2 = index(linhab, "A30"). ... execute. Obrigado -- Sérgio Henrique Almeida da Silva Junior Doutorando em Epidemiologia em Saúde Pública Escola Nacional de Saúde Pública Sérgio Arouca - ENSP/FIOCRUZ http://lattes.cnpq.br/1611345552843383 Tel: (21) 94429486/78101651 id: 123*20942

4 12

Re: [R-br] Criação de função fortran em pacote (base)
by Junior Beleti 23 Mar '12

23 Mar '12

Olá Benilton, segue link para download do pacote geoComp: www.din.uem.br/~pg45178/geocomp/geoComp.tar.gz Sobre o problema de desempenho, acredito que as complicações sejam no mvrnorm() pois é lá que o script leva mais tempo na execução. Continuo trabalhando para encontrar solução ou melhoria no código. Obrigado pela atenção. Carlos.

1 0

Msg de Erro em Grafico - In axis(1, x, ...) : "frame.plot" não é um parâmetro gráfico
by Augusto Ribas 22 Mar '12

22 Mar '12

Ola galera. Eu estou com uma duvida que me atormenta aqui. Vejam os seguintes graficos: #CMR demo1 <- read.csv("http://www.ats.ucla.edu/stat/data/demo1.csv") demo1 <- within(demo1, { group <- factor(group) time <- factor(time) id <- factor(id) }) par(mfrow=c(1,2)) par(cex = .6) #grafico with(demo1, interaction.plot(time, group, pulse, ylim = c(5, 20), lty= c(1, 12), lwd = 3, ylab = "mean of pulse", xlab = "time", trace.label = "group",frame.plot=F)) axis(1,at=1:3) title("1") with(demo1, interaction.plot(time, group, pulse, ylim = c(5, 20), lty= c(1, 12), lwd = 3, ylab = "mean of pulse", xlab = "time", trace.label = "group")) title("2") Minha duvida é o seguinte. No grafico 1 eu uso o argumento frame.plot=F, no segundo nao uso nada. Ai no primeiro ele nao desenha o quadrado em volta do grafico, o argumento funciona, mas vem uma msg de erro, "Esse não é um parametro grafico". O grafico é pra ficar assim mesmo, mas eu nao entendo pq ele aceita o argumento e fala que não existe. Alguém saberia dizer o que acontece numa linguagem simples? Que eu ja caçei entender isso e não consegui. Obrigado -- Grato Augusto C. A. Ribas Site Pessoal: http://augustoribas.heliohost.org Lattes: http://lattes.cnpq.br/7355685961127056

2 2

Teste - lista - gentileza ignorar email
by Benilton Carvalho 22 Mar '12

22 Mar '12

benilton

1 0

converter dia do mes em dia juliano
by Lucas Fernandes de Souza 22 Mar '12

22 Mar '12

Bom dia a todos! Tenho um arquivo da seguinte maneira: day year maxt mint rain radn 1 1978 30.7 19.8 0.0 22.6 2 1978 27.2 17.4 6.0 21.4 3 1978 21.4 13.6 0.5 19.1 4 1978 23.2 14.3 0.0 20.4 5 1978 28.8 17.2 0.0 23.2 6 1978 31.5 19.0 16.4 24.1 Onde day é o dia do mes, ou seja 1 a 31 (dependendo do mes) e ano vai de 78 a 2010, incluindo anos bissextos. Preciso converter o day para dia juliano, ou seja, para cada ano, tem q variar de 1 a 365 (anos normais) e 1 a 366 (anos bissextos) Alguem pode me ajudar? Obrigado -- Lucas F. de Souza Eng. Agrícola - CONAB M.Sc. Produção Vegetal - UENF Doutorando Engenharia de Sistemas Agrícolas - ESALQ/USP tel: (19) 8160-7068

4 3

ftp
by Vinicius Brito Rocha 22 Mar '12

22 Mar '12

Caros colegas de lista, estou precisando exportar um resultado .txt para um ftp. Tem como fazer isso via R? abs. -- Vinicius Brito Rocha. Estatístico e Atuário, Mestre em Pesquisa Operacional www.aplicademic.blogspot.com http://twitter.com/viniciusbritor "Não se preocupe muito com as suas dificuldades em Matemática, posso assegurar-lhe que as minhas são ainda maiores." - Albert Einstein.

6 15

E-mails repetidos
by Gisele Correia 22 Mar '12

22 Mar '12

Bom dia, Por algum motivo estou recebendo os e-mails das discussões da lista de forma repetida. Provavelmente devo ter solicitado cadastro tanto na web como no e-mail. Gostaria que o moderador da lista pudesse solucionar isso para mim. Obrigada pela atenção e desculpe o transtorno. -- Gisele

1 0

[Dúvida] Máxima Verossimilhança no R.
by Pedro Rafael 21 Mar '12

21 Mar '12

Pessoal estou com um problema para estimar os parâmetros de dados com distribuição Beta com parâmetros p e q. Gerei aleatoriamente um conjuntos de dados com distribuição densidade de probabilidade beta no intervalo (0,1). Encontrei a distribuição de log verossimilança e como não existe solução fechada para os parâmetros p e q, temos que estimar por métodos numéricos como o de Newton Raphson, BFGS, etc. Código: ######################################## dados = rbeta(1000000,1,3) # onde p = 1 e q = 3. hist(dados) # Cara dos dados. log.verossimilhanca <- function(teta,y){ p = teta[1] q = teta[2] logl <-sum((1-p)*log(y)+(1-q)*log(1-y)+log(gamma(p+q)) - log(gamma(q)) - log(gamma(p))) return(-logl) } optim(theta <- c(1,1),fn = log.verossimilhanca,y=dados,method="BFGS") Como informar o espaço paramétrico de cada um dos parâmetros da distribuição, no caso particular da distribuição beta ambos sendo maiores que zero? -- Saudações, Pedro Rafael Diniz Marinho.

3 5

linhas de batimetria
by Rodrigo Plei 21 Mar '12

21 Mar '12

Pessoal, Alguém sabe como deixar as linhas de batimetria tracejadas quando plotadas a partir de um shapefile usando o maptools? Agradeço desde logo! 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

1 0

Ajuda com Análise de Similaridade - ANOSIM
by Flávia Souza 21 Mar '12

21 Mar '12

Bom dia a todos, parabéns pela iniciativa da lista. Muito bom mesmo. Solicito auxílio quanto a realização de uma análise de Similaridade do tipo ANOSIM, pois preciso realizar e gostaria de usar o R, entretanto, não tenho certeza da entrada dos dados na planilha excel (csv) e os procedimentos quanto a comandos para gerar a análise e respectivo gráfico no software. Minha análise são de três ambientes distintos e três hábitats em cada ambiente. Alguém pode me ajudar, por favor.... Se houver necessidade envio os dados Grata pela atenção, Flávia -- MSc. Flávia Kelly S. Souza Universidade Federal do Amazonas - UFAM Faculdade de Ciências Agrárias - Depto. de Ciências Pesqueiras Manaus - Amazonas flavia.souza(a)piatam.org.br flakel.souza(a)gmail.com flaviasouza(a)ufam.edu.br skype: fla_kel

2 2

[dúvida] Vários plots no mesm frame com margens diferentes
by Marcelo Claro de Souza 21 Mar '12

21 Mar '12

Preciso plotar 8 figuras no mesmo frame em duas colunas e 4 linhas, mas gostaria de juntar as quatro figuras em cada coluna e colocar apenas um eixo na figura inferior. estou fazendo a divisão do frame por layout(matrix(1:8,nrow=4,ncol=2) e a determinação da margem geral do frame por par(mar=c(4.5,4.5,1,1) Como ajusto individualmente as margens de cada figura? Muito obrigado. Marcelo Claro de Souza Biologist, PhD student in Plant Biology Institute of Bioscience - UNESP, Brazil

2 1

Dúvida exportação PNG
by Julianna Trindade 21 Mar '12

21 Mar '12

Pessoal, gostaria que meu plot fosse automaticamente transformado em um arquivo png em paisagem. Estou usando *png(filename="c:/users/Julianna/Documents/SecSaude/produtos_nascvivos/2006.png", height=1360, width=768, bg="white")* e está exportando em retrato. Alguém, por favor, pode me ajudar?? -- Julianna Trindade "Exalta-te, Senhor, na Tua força! Nós cantaremos e louvaremos o Teu poder." Salmos 21; 13

2 5

Ajuda modelo misto ordinal - MCMCglmm.
by Ivan Bezerra Allaman 20 Mar '12

20 Mar '12

#Boa noite senhores! Após vários dias estudando sobre dados categóricos e pesquisando (vários fóruns com respostas do dono da library - todas incipientes) como interpretar o 'output' da função MCMCglmm venho recorrer a ajuda dos senhores, caso alguém já tenha trabalhado com função MCMCglmm no caso de variáveis dependentes ordinais. Já li e reli todos os pdf's do pacote, coursenotes do Jarrod, enfim, estou esgotado. Para esclarecer a base de dados, os tratamentos (denominado FASES) consistem de 3 níveis (1-pre, 2-própolis, 3-vincris), cujo as variável resposta ordinal tem 3 categorias (1-normal, 2-agudo, 3-cronico). Vide tabela! du <- transform(read.table('http://dl.dropbox.com/u/33619290/Dados/Dtest.txt',h=T),FASES=factor(FASES),ALT.RENAIS=ordered(ALT.RENAIS)) summary(du) library('MCMCglmm') du <- subset(du, ALT.RENAIS != 'NA') tabela <- table(du[,c(2,4)]) tabela colnames(tabela) <- c('Normal','Aguda','Crônica') rownames(tabela) <- c('Pre','Propolis','Vincr') tabela #Um modelo misto seria: set.seed(1) mod1 <- MCMCglmm(ALT.RENAIS ~-1+FASES, random= ~ ANIMAIS, family='ordinal',pl=TRUE,data=du) summary(mod1) #Aí começa o sofrimento, pois a documentação é insuficiente neste caso. Segundo o próprio Jarrod (fóruns), as médias a posteriori dos coeficientes das covariáveis estão na escala probit. Segundo o meu estudo, estes coeficientes são os scores da distribuição normal padrão. Mais os escores não deveriam corresponder aos cutpoints? Neste caso, teríamos j(categorias da variável resposta)-1 cutpoints, ou seja, 2 cutpoints. O output só me apresenta um cutpoint. Como é possível então calcular as probabilidades com apenas um cutpoint? Segundo a documentação(Vignettes, pag 22), se P(y=k)=F(yk|l(vlatente),1) - F(yk-1|l,1), esse 'um' provavelmente seria a categoria 'um' da variável dependente? Enfim senhores, como posso extrair as probabilidades para as fases referentes a cada categoria da variável dependente? Desde já agradeço a atenção de todos. #embora um pouco lógico, os resultados são absurdos! latentv <- mean(mod1$Liab) cutpoint <- mean(mod1$CP) pnorm(-(latentv), 0, sqrt(2)) pnorm(cutpoint - (latentv),0, sqrt(2)) - pnorm((latentv),0, sqrt(2)) 1- pnorm(cutpoint - (latentv),0, sqrt(2)) #isso teria algum resultado lógico até certo ponto, mais seria apenas referente a categoria 1 da variável dependente. E as outras? bre <- c(mean(mod1$Liab),mean(mod1$Sol[,1]),mean(mod1$Sol[,2]),mean(mod1$Sol[,3])) pnorm(bre[2])-pnorm(bre[1]) pnorm(bre[3])-pnorm(bre[2]) pnorm(bre[4])-pnorm(bre[3])#probabilidade negativa nãooo! \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

Anova para delineamento hierarquico
by Gisele Correia 20 Mar '12

20 Mar '12

Bom dia a todos, Gostaria de um exemplo de como funciona uma análise em anova num delineamento experimental hierárquico, sendo o primeiro fator fixo e o segundo fator aleatório. Att., -- Gisele Batista Correia Engenheira de Pesca

2 1

dúvida gráfico ponto frequencia relativa
by Lucas Fernandes de Souza 20 Mar '12

20 Mar '12

Boa tarde. Tenho esse intervalo de frequencia relativa de precipitação, conforme abaixo: int freq 0-10 0.8719 10-20 0.06878 20-30 0.03211 30-40 0.01452 40-50 0.00705 50-60 0.00299 60-70 0.00141 70-80 0.00058 80-90 0.00058 90-100 0.00008 quero plotar um gráfico de ponto e linha com esses dados, onde o eixo X fique conforme o intervalo acima. Alguma sugestão? obrigado -- Lucas F. de Souza Eng. Agrícola - CONAB M.Sc. Produção Vegetal - UENF Doutorando Engenharia de Sistemas Agrícolas - ESALQ/USP tel: (19) 8160-7068

2 1

Estouro de Inteiros
by Jodavid Ferreira 20 Mar '12

20 Mar '12

Eu estou realizando um cáculo onde chega uma certa hora, que aparece este mensagem! Warning message: NAs produzidos por estouro de inteiro Alguem sabe com resolver? -- Jodavid Ferreira Responsável pela comissão de Comunicação do MJ da RCCPB - http://jovem.rccpb.com.br Responsável pelo Ministério Jovem e 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

2 1

Ajuda - coeficientes estatisticamente diferentes entre MQO e métodos espaciais.
by Ariene Salgueiro 20 Mar '12

20 Mar '12

Olá pessoal, estou precisando de ajuda para saber se os coeficientes estimados por MQO e por SAR de Kelejian e Prucha são estatisticamente diferentes. E também entre MQO e o SEM e entre o MQO e o SARMA. Alguém pode me ajudar, por favor!? Ariene.

1 0

Erro do shape no linux
by Jodavid Ferreira 20 Mar '12

20 Mar '12

Quando eu chamo o shape no R, acontece este erro da segunda linha: > mapa <- readShapePoly(".../shapes/shape/Municipios.shp") Erro em make.names(onames, unique = TRUE) : string multibyte inválida 6 Alguém sabe dizer o que é, como resolver? -- Jodavid Ferreira Responsável pela comissão de Comunicação do MJ da RCCPB - http://jovem.rccpb.com.br Responsável pelo Ministério Jovem e 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

2 2

Re: [R-br] Dúvida exportação PNG
by Jodavid Ferreira 20 Mar '12

20 Mar '12

Está exportando porque o 'height' está maior do que o 'width', no seu código. O height representa a altura da image, e o width a largura. tente inverter colocar assim *height=**768**, width=**1360* * * *png(filename="c:/users/Julianna/Documents/SecSaude/produtos_nascvivos/2006.png", height=768, width=1360, bg="white")* e está exportando em retrato. Abraço, -- Jodavid Ferreira Responsável pela comissão de Comunicação do MJ da RCCPB - http://jovem.rccpb.com.br Responsável pelo Ministério Jovem e 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 Em 19 de março de 2012 17:21, Julianna Trindade [via R-br] < > ml-node+s2285057n4486243h68(a)n4.nabble.com> escreveu: > > Pessoal, >> >> gostaria que meu plot fosse automaticamente transformado em um arquivo >> png em paisagem. >> Estou usando *png(filename="c:/users/Julianna/Documents/SecSaude/produtos_nascvivos/2006.png", >> height=1360, width=768, bg="white")* e está exportando em retrato. >> >> Alguém, por favor, pode me ajudar?? >> >> >> >> -- >> Julianna Trindade >> >> "Exalta-te, Senhor, na Tua força! Nós cantaremos e louvaremos o Teu >> poder." >> >> Salmos 21; 13 >> >> _______________________________________________ >> R-br mailing list >> [hidden email] <http://user/SendEmail.jtp?type=node&node=4486243&i=0> >> 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-Duvida-exportacao-PNG-tp4486243p4486… >> 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 Responsável pelo Ministério Jovem e 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

1 0

Ajuda parar exibir variação temporal de dados em mapas
by Thiago Veloso 18 Mar '12

18 Mar '12

Pessoal, Possuo três mapas (em formato raster) de desmatamento na Amazônia em diferentes anos. Essas taxas de desmatamento variam de 0 a 1. A minha intenção é mostrar a variação temporal do desmatamento através de histogramas que mostram a contagem de células em diferentes classes (p. ex., de 0 a 0.3, de 0.3 a 0.6 e de 0.6 a 1) em cada ano. Aqui está uma ideia conceitual do que preciso fazer (apenas ilustrativa, criada com dados hipotéticos): http://img24.imageshack.us/img24/6416/examplezxt.png e aqui disponibilizo um arquivo com amostra do arquivo de dados http://db.tt/XXDtRee1. Trata-se um netcdf com três níveis de tempo. Na lista de SIG do R eu já recebi algumas dicas para fazer isso. Entretanto, estou recebendo um erro ao tentar criar fatores nos dados. Será que alguém aqui consegue me dar uma auxílio? Segue o código que estou utilizando e a mensagem de erro: --------------------------------------------> #Load required packages> library (raster)Loading required package: spraster version 1.9-70 (27-February-2012)> library (ggplot2)> library (reshape2)> > #Load data> deforest <- brick("/home/thiago/Dropbox/defor-3-years.nc")Loading required package: ncdf> > #Probe layer names and plot raster data> layerNames(deforest)[1] "X0.1" "X0.2" "X0.3"> plot(deforest)> > #Convert spatial object to data frame (omitting NAs)> meltR <- melt(as.array(deforest),na.rm=TRUE)> > #Define classes used to count grids> meltR$class <- cut(meltR$value, c(0,0.3,0.6,1))> #Define factors to plot individual years> meltR$X3 <- factor(meltR$X3, labels=c('2000','2001','2002'))Error in factor(meltR$X3, labels = c("2000", "2001", "2002")) : invalid labels; length 3 should be 1 or 0-------------------------------------------- O código continuaria assim: --------------------------------------------#Plot datap = ggplot(meltR,aes(x=class))p + geom_histogram(aes(fill=class))+facet_grid(. ~ X3)-------------------------------------------- Caso seja útil, segue o sumário do sessionInfo(): --------------------------------------------> sessionInfo()R version 2.14.1 (2011-12-22)Platform: x86_64-pc-linux-gnu (64-bit) locale: [1] LC_CTYPE=en_US.UTF-8 LC_NUMERIC=C [3] LC_TIME=en_US.UTF-8 LC_COLLATE=en_US.UTF-8 [5] LC_MONETARY=en_US.UTF-8 LC_MESSAGES=en_US.UTF-8 [7] LC_PAPER=C LC_NAME=C [9] LC_ADDRESS=C LC_TELEPHONE=C [11] LC_MEASUREMENT=en_US.UTF-8 LC_IDENTIFICATION=C attached base packages:[1] stats graphics grDevices utils datasets methods base other attached packages:[1] ncdf_1.6.6 reshape2_1.2.1 ggplot2_0.9.0 raster_1.9-70 sp_0.9-97 loaded via a namespace (and not attached): [1] colorspace_1.1-1 dichromat_1.2-4 digest_0.5.2 grid_2.14.1 [5] lattice_0.20-6 MASS_7.3-17 memoise_0.1 munsell_0.3 [9] plyr_1.7.1 proto_0.3-9.2 RColorBrewer_1.0-5 scales_0.2.0 [13] stringr_0.6 -------------------------------------------- Um verdadeiro desafio - estou desde 4a feira quebrando a cabeça com isso. Agradeço desde já qualquer contribuição! Thiago.

2 3

[Dúvida] Colocar figura dentro de gráfico.
by Pedro Rafael 17 Mar '12

17 Mar '12

Pessoal,estou querendo adicionar uma imagem JPG/PNG em um plot. Como fazer isto em R? No meu específico problema estou construindo um mapa e no canto queria colocar uma figura. -- Saudações, Pedro Rafael Diniz Marinho. Estatístico - Secretaria de Estado da Saúde - PB.

2 1

homocedasticidade modelo não linear
by Fernando Antonio de souza 16 Mar '12

16 Mar '12

---------- Mensagem encaminhada ---------- De: Fernando Antonio de souza <nandodesouza(a)gmail.com> Data: 14 de março de 2012 13:50 Assunto: duvida sobre presunção de homocedasticidade Para: r-br(a)listas.c3sl.ufpr.br Caros amigos Estou trabalhando com dados de desenvolcimento de fetos ovinos. Estou modelando o seu crescimento ao longo da gestação. Para tal tenho utilizado o modelo de gompertz da forma log(Y)= a-b*exp(-c*x) sendo x a idade gestacional em dias e Y o peso do feto. Os meus dados atendem o presuposto de normalidade, porém não o faz para o pressuposto de igualdade das variâncias. Tenho lido na literatura a possibilidade de modelar a variância nestes casos. A variância residual observada graficamente se eleva e em seguida rebaixa conforme as médias do modelo se elevam. Na própria literatura comenta que este tipo de comportamento pode ser modelado utilizando o método VarPower() na função glns do pacote nlme. Eu assim o fiz entretanto o novo modelo não obteve melhor ajuste que o anterior quando comparados pela função Anova (p-value e AIC,BIC). Minha dúvida é a seguinte sem a igualdade das variâncias eu não posso validar o modelo. Existe alguma forma de corrigir este problema? abaixo disponho do CMR library(nlme) dput(gestsimp) gompertz1<-nls(log10(PESO.FETO)~a-b*exp(-c*Gest),start=list(a=8.0,b=10.40,c=0.0023),data=gestsimp) gompertz2<-gnls(log10(PESO.FETO)~a-b*exp(-c*Gest),start=list(a=8.0,b=10.40,c=0.0023),weights=varPower(),data=gestsimp) summary(gompertz1) plot(log10(PESO.FETO)~Gest,data=gestsimp,main="Peso do feto em função da idade de cordeiros") curve(5.697742-8.360311*exp(-0.003531*x),add=T,col=2) caso o dput()não funcione segue a saída dele dput(gestsimp)structure(list(Animal = c(3L, 51L, 83L, 89L, 91L, 92L, 93L, 48L, 52L, 55L, 69L, 76L, 78L, 11L, 25L, 41L, 42L, 46L, 50L, 53L, 57L, 10L, 19L, 20L, 30L, 38L, 45L), Manejo = c(1L, 2L, 2L, 2L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 1L, 1L, 2L, 2L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L), Gest = c(90L, 90L, 90L, 90L, 90L, 90L, 90L, 110L, 110L, 110L, 110L, 110L, 110L, 110L, 130L, 130L, 130L, 130L, 130L, 130L, 130L, 140L, 140L, 140L, 140L, 140L, 140L), Fetos = c(1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L ), PV = c(51.1, 46.4, 53.5, 58.2, 30.7, 43.2, 41.5, 42, 43.4, 48.1, 42.2, 47.5, 32, 64.1, 49.9, 40.5, 42, 40.4, 50, 49, 52, 60.1, 55.9, 57, 51.5, 45, 42.8), PCVZ = c(36.93, 36.135, 46.125, 48.02, 23.095, 33.75, 32.575, 34.99, 34.465, 39.578, 30.84, 38.045, 27.455, 54.797, 43.595, 32.12, 33.841, 33.16, 42.06, 40.605, 42.325, 50.436, 46.515, 48.845, 43.075, 36.55, 35.435), PESO.FETO = c(0.215, 0.555, 0.49, 0.53, 0.225, 0.465, 0.61, 1.465, 0.585, 1.755, 0.62, 1.615, 1.045, 0.985, 2.79, 2.42, 2.286, 2.41, 2.3, 3.1, 3.105, 5.38, 4.065, 4.258, 3.878, 3.49, 3.032), FETOGKG = c(60.5065, 50.8621, 67.9386, 32.7586, 41.7391, 46.1815, 47.7506, 51.0421, 66.0087, 63.5711, 50.8569, 60.5784, 49.8817, 62.1197, 60.1357, 0, 0, 62.2886, 0, 55.4768, 53.5301, 68.0826, 58.413, 59.6224, 61.1382, 50.2531, 62.5518 )), .Names = c("Animal", "Manejo", "Gest", "Fetos", "PV", "PCVZ", "PESO.FETO", "FETOGKG"), row.names = c(11L, 13L, 15L, 19L, 20L, 21L, 22L, 30L, 31L, 32L, 35L, 36L, 38L, 42L, 45L, 46L, 47L, 48L, 49L, 51L, 55L, 63L, 66L, 67L, 71L, 76L, 80L), class = "data.frame", na.action = structure(22L, .Names = "62", class = "omit"))

4 5

[dúvida] Como calcular área de um gráfico?
by Marcelo Claro de Souza 16 Mar '12

16 Mar '12

Plotei um gráfico de linhas, e preciso calacular a área da parte de baixo da figura, tem como fazer isso? a figura foi plotada usando plot(x,y) e visualmente se parece com uma gaussiana. Muito obrigado pela ajuda. Marcelo Claro de Souza Biologist, PhD student in Plant Biology Institute of Bioscience - UNESP, Brazil

3 2

Luz em glm() poisson vs quasipoisson
by Marcelo Cardoso mello 16 Mar '12

16 Mar '12

Gostaria de uma luz na análise abaixo. Usando poisson e quasipoisson tenho significância no modelo (ANOVA). Porem, NO SUMMARY NÃO TENHO SIGNIFICÂNCIA PARA OS PARÂMETROS DE QUASIPOISSON. Então é o seguinte, existe ou não diferença entre os tratamentos na quasipoisson? dados<-structure(list(Momento = structure(c(2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L), .Label = c("Alternativa", "Atual"), class = "factor"), CT = c(13000L, 12000L, 31000L, 39000L, 1250000L, 13600000L, 6880000L, 47450000L, 890000L, 33000L, 15400000L, 130000L, 6075L, 850000L, 690000L, 1400000L, 41300000L, 4900000L, 14900000L, 19000L, 31600000L, 1900000L, 600000L, 2000L, 400000L, 840000L, 900L, 900L, 4700000L, 910000L, 1000L, 210L, 90L, 4725L, 200000L, 480000L, 460000L, 107300000L, 8700000L, 3700000L, 4700000L, 15000000L, 20000L, 180000L, 10000L, 30500000L, 59000L, 260000L, 1000L, 100L, 550000L, 300L, 2750L, 46000L, 21500L, 5000L, 1050000L, 9L, 600L, 110L, 310L, 190L, 50L, 18L, 230L, 135L, 100L, 1000L, 108L, 230L, 100L, 3100L, 350L, 30L, 2400L, 19500L, 10000000L, 200L, 10L, 10L, 9L, 9L, 9L, 310000L, 9L, 24000L, 10L, 150L, 9L, 9L, 20L, 80L, 10L, 20L, 9L, 10L, 240L, 40L, 9L, 9L)), .Names = c("Momento", "CT"), row.names = c(NA, -100L), class = "data.frame") m1<-glm(CT~Momento, family=poisson , data=dados) m2<-glm(CT~Momento, family=quasipoisson, data=dados) anova(m1, test= "Chisq") anova(m2, test= "Chisq") summary(m1) summary(m2)

4 4

tcl/tk error : segfault in Rcmdr
by Pedro Paulo Chrispim 16 Mar '12

16 Mar '12

Por algum motivo, o R não conseguia carregar o xwindow. baixei uma versão mais antiga e agora consegui carregar o Rcmdr na versao 2.13.2 as instruções estão em.. http://socserv.socsci.mcmaster.ca/jfox/Misc/Rcmdr/installation-notes.html ab, Pp

1 0

Séries Temporais com frequências distintas
by Taynãna César Simões 16 Mar '12

16 Mar '12

Prezados, Estou analisando uma série de contagens de eventos por semana epidemiológica, entre os anos de 2007 a 2010. No entanto, no ano de 2008 existem 53 semanas epidemiológicas, e não 52 como nos demais anos. Então, o problema é que se crio um série temporal da forma tradicional, com freq=52 semanas, ele considera a semana 53 de 2008 como a semana 1 de 2009. Minha série vai da semana 5 de 2007 à semana 52 de 2010. Usei: serie <- ts(contagens, start(2007,5),freq=52) Minha idéia então era criar uma série para cada ano e juntá-las depois, mas como posso fazer isso? Tentei usar a library zoo para séries irregulares, mas não consegui. Obrigada, Taynãna Fiocruz

2 1