- R-br - DINF UFPR

Anova de fatorial duplo em DBC (um fator quantitativo)
by Andreia de Oliveira Vieira 27 Mai '13

27 Mai '13

Boa tarde, Muito obrigada Alisson e Tiago pelas dicas, mas mesmo testando todas ainda não consegui solucionar o problema, portanto segue a mensagem que aparece quando tento utilizar o pacote ExpDes.pt para rodar a anova com estes resultados. Andreia > data.frame(Girassol) metodos doses blocos clorofila.1 clorofila.2 altura altura.1 diametro 1 1 1 1 32.5 36.2 173 173 25.0 2 1 2 1 32.2 35.9 179 178 25.3 3 1 3 1 32.4 38.3 170 170 24.5 4 2 1 1 32.7 35.2 177 177 24.0 5 2 2 1 31.9 37.6 183 182 25.0 6 2 3 1 33.3 35.4 184 182 27.8 7 3 1 1 35.0 36.3 184 184 26.6 8 3 2 1 33.9 39.0 170 172 24.6 9 3 3 1 32.8 35.9 179 179 27.0 10 4 1 1 34.7 33.9 172 173 26.3 11 4 2 1 34.4 36.8 169 169 25.0 12 4 3 1 34.5 36.5 156 156 25.0 13 1 1 2 30.9 38.3 179 180 26.9 14 1 2 2 33.7 36.2 181 180 29.0 15 1 3 2 31.3 35.8 186 185 24.9 16 2 1 2 33.4 37.2 183 185 26.4 17 2 2 2 30.9 34.4 179 179 25.2 18 2 3 2 32.5 37.3 188 188 28.1 19 3 1 2 32.4 36.4 185 185 26.0 20 3 2 2 30.5 36.1 190 190 25.9 21 3 3 2 32.1 37.9 177 178 24.2 22 4 1 2 31.7 36.6 182 182 27.2 23 4 2 2 30.6 36.4 176 177 29.0 24 4 3 2 32.1 40.2 158 159 24.0 25 1 1 3 32.5 42.8 182 183 24.5 26 1 2 3 32.6 40.0 185 186 24.8 27 1 3 3 33.6 37.3 188 188 25.6 28 2 1 3 34.9 38.8 188 188 23.7 29 2 2 3 33.5 39.0 183 184 23.5 30 2 3 3 32.4 39.0 176 177 22.3 31 3 1 3 32.8 40.5 180 181 24.4 32 3 2 3 31.8 39.3 179 180 22.4 33 3 3 3 33.4 41.2 178 178 24.0 34 4 1 3 33.4 41.2 173 174 23.4 35 4 2 3 33.6 39.8 171 171 22.5 36 4 3 3 33.4 38.8 165 165 21.8 > fat2.dbc(metodos,doses,blocos,clorofila1,quali=c(TRUE, FALSE),mcomp = "tukey", fac.names = c("F1", "F2"), sigT = 0.05, sigF = 0.05) ------------------------------------------------------------------------ Legenda: FATOR 1: F1 FATOR 2: F2 ------------------------------------------------------------------------ Error in model.frame.default(formula = resp ~ Bloco + Fator1 * Fator2, : invalid type (list) for variable 'resp' > fat2.dbc(clorofila 1~blocos+metodos*doses,quali=c(TRUE, FALSE),mcomp = "tukey", fac.names = c("F1", "F2"), sigT = 0.05, sigF = 0.05) Error: unexpected numeric constant in "fat2.dbc(clorofila 1"

3 2

codigo de leitura de dados
by alanarocha＠sapo.pt 27 Mai '13

27 Mai '13

Boa tarde, as variaveis tem o mesmo tamanho não percebo o erro: dados: 2sem2011 custosnaopanakeia hgbnaopanakeia JUL 29,41 75,9% AGO 32,28 77,5% SET 26,88 79,4% OUT 25,94 76,8% NOV 35,27 76,1% DEZ 24,50 78,3% codigo: setwd("I:/dados não Panakeia e Panakeia") dados<-read.table("DadosnãoPanakeia2sem2011.csv",header=TRUE, comment.char = "%",dec=",",sep=";") dados wilcox.test(dados$custosnaopanakeia, hgbnaopanakeia) erro: [1] PK... <0 rows> (or 0-length row.names) Como resolvo? cumprimentos Ana Rocha

2 1

codigo de dados
by alanarocha＠sapo.pt 27 Mai '13

27 Mai '13

Boa tarde, esta a dar-me aqui um erro estranho... dados: 2sem2011 custosnaopanakeia hgbnaopanakeia JUL 29,41 75,9% AGO 32,28 77,5% SET 26,88 79,4% OUT 25,94 76,8% NOV 35,27 76,1% DEZ 24,50 78,3% codigo: setwd("I:/dados não Panakeia e Panakeia") dados<-read.table("DadosnãoPanakeia2sem2011.csv",header=TRUE, comment.char = "%",dec=",",sep=";") dados wilcox.test(dados$custosnaopanakeia, dados$hgbnaopanakeia) erro: [1] PK... <0 rows> (or 0-length row.names) como resolvo? Cumprimentos Ana Rocha

1 0

Aumentar letras do plot glm
by Renata Carvalho 27 Mai '13

27 Mai '13

Boa noite a todos, qual o comando pra aumentar as letras do plot de um glm? Eu utilizei os pacotes MASS, MuMln e nlme. Obrigada Renata

2 3

Anova de fatorial duplo em DBC (um fator quantitativo)
by Andreia de Oliveira Vieira 27 Mai '13

27 Mai '13

Boa noite, Sou uma iniciante do R e não estou conseguindo rodar anova de um experimento fatorial duplo em DBC, sendo que um dos fatores é quantitativo. Estou utilizando a função fat2.dbc e seguindo as instruções da descrição do pacote, mas aparece esta mensagem: > fat2.dbc(clorofila1,blocos,metodos, doses, quali=c(TRUE, FALSE),mcomp = "tukey", fac.names = c("F1", "F2"), sigT = 0.05, sigF = 0.05) ------------------------------------------------------------------------ Legenda: FATOR 1: F1 FATOR 2: F2 ------------------------------------------------------------------------ Error in data.frame(..., check.names = FALSE) : arguments imply differing number of rows: 5, 144 Muito obrigada Andreia

3 2

teste wilcoxon
by alanarocha＠sapo.pt 26 Mai '13

26 Mai '13

Bom dia, existem codigos diferentes para fazer o teste wilcoxon directamente e o teste T com correcção de wilcoxon? se sim quais? e qual é a diferença entre se fazer u.test e wilcox.test no meu caso eu so quero fazer o teste wilcoxon para comparar os dois grupos e não correcção. Cumprimentos Ana Rocha ----- Fim de mensagem reenviada -----

3 2

Seleção de plantas ( sorteio)
by Hélio Gallo Rocha 25 Mai '13

25 Mai '13

Caros da Lista Com os dados hipotéticos abaixo: ponto=c(rep(3001,6),rep(3002,6),rep(3003,6)) planta=c(rep(seq(1:6),3)) valor=seq(1.5,10,.5) dados=cbind(ponto,planta,valor);dados ponto planta valor [1,] 3001 1 1.5 [2,] 3001 2 2.0 [3,] 3001 3 2.5 [4,] 3001 4 3.0 [5,] 3001 5 3.5 [6,] 3001 6 4.0 [7,] 3002 1 4.5 [8,] 3002 2 5.0 [9,] 3002 3 5.5 [10,] 3002 4 6.0 [11,] 3002 5 6.5 [12,] 3002 6 7.0 [13,] 3003 1 7.5 [14,] 3003 2 8.0 [15,] 3003 3 8.5 [16,] 3003 4 9.0 [17,] 3003 5 9.5 [18,] 3003 6 10.0 então temos 3 pontos ( 3001, 3002, 3003) seis plantas em cada ponto ( 1,2,3,4,5,,6) Já tentei usar sample, mas sorteia aleatoriamente Preciso que seja sorteada uma planta de cada ponto, depois vou sortear duas plantas de cada ponto, e assim sucessivamente Agradeço desde já -- Hélio Gallo Rocha IFSULDEMINAS - Câmpus Muzambinho

3 9

Operação em data.frame com colunas 'numeric'
by Thiago V. dos Santos 24 Mai '13

24 Mai '13

Olá pessoal, Estou com uma dúvida bem básica no R. Tenho um data.frame que contem anos e meses no formato numérico e gostaria de realizar algumas operações nesse DF. Vejam abaixo: year month anom3_4 1950 1 -1.50 1950 2 -1.77 1950 3 -0.73 1950 4 -0.94 1950 5 -1.47 1950 6 -0.70 1950 7 -0.48 1950 8 -0.59 1950 9 -1.09 1950 10 -0.57 1950 11 -1.09 1950 12 -0.93 1951 1 -1.04 1951 2 -0.07 1951 3 -0.70 1951 4 0.05 1951 5 -0.22 1951 6 0.04 1951 7 0.41 1951 8 0.74 1951 9 0.38 1951 10 0.80 1951 11 0.80 1951 12 0.72 Gostaria, inicialmente, de manter somente os meses 9, 10 e 11 de cada ano (no caso, 1950 e 1951). Em seguida, gostaria de calcular a média da coluna anom3_4 em cada ano. Se as colunas year e month estivessem como factor, seria fácil fazer isso. Mas acho que não dá para criar factor a partir de numerics. Sugestões de linhas mágicas para fazer isso? Obrigado, -- Thiago V. dos Santos PhD student Land and Atmospheric Science University of Minnesota http://www.laas.umn.edu/CurrentStudents/MeettheStudents/ThiagodosSantos/ind… Phone: (612) 323 9898

3 5

Entrar com dados já tabelados
by Carlos Andrade 24 Mai '13

24 Mai '13

Olá pessoal da lista. Já li em algum material e não consigo recuperar agora. Como entro com dados já tabelados no R? São dados do Censo demográfico 2010 do IBGE. A tabela tem uma coluna com as classes de renda em salários mínimos e a segunda coluna com a frequência respectiva. Agradeço antecipadamante a atenção. -- Atenciosamente, Prof. Carlos A. S. de Andrade LAPEA - Laboratório de Pesquisa em Economia Aplicada e Engenharia de Produção Universidade Federal de Campina Grande. Centro de Humanidades Unidade Acadêmica de Economia

4 4

teste wilcoxon
by alanarocha＠sapo.pt 24 Mai '13

24 Mai '13

Bom dia, existem codigos diferentes para fazer o teste wilcoxon directamente e o teste T com correcção de wilcoxon? se sim quais? e qual é a diferença entre se fazer u.test e wilcox.test no meu caso eu so quero fazer o teste wilcoxon para comparar os dois grupos e não correcção. Cumprimentos Ana Rocha ----- Fim de mensagem reenviada -----

1 0

teste wilcoxon
by alanarocha＠sapo.pt 24 Mai '13

24 Mai '13

Bom dia, existem codigos diferentes para fazer o teste wilcoxon directamente e o teste T com correcção de wilcoxon? se sim quais? e qual é a diferença entre se fazer u.test e wilcox.test no meu caso eu so quero fazer o teste wilcoxon para comparar os dois grupos e não correcção. Cumprimentos Ana Rocha ----- Fim de mensagem reenviada -----

1 0

[Dúvida] Combinar expression e paste.
by Pedro Rafael 23 Mai '13

23 Mai '13

Pessoal estou fazendo alguns gráficos e estou combinando os comandos expression e paste. O problema é que quando chamo o comando paste dentro do comando expression só estou conseguindo colocar os valores numéricos das variáveis digitando esses valores. Por exemplo: beta = 10 legend("topright", legend = c(expression(beta, paste(" a= 10"))) Se eu fizer legend("topright", legend = c(expression(beta, paste(" = ", beta))) o valor 10 da variável beta não é substituído. Gostaria de alguma forma de não precisar digitar uma informação que já tenho em uma variável. [ ], Pedro Rafael Diniz Marinho.

4 14

Cluster com var categorica e continua
by nani 23 Mai '13

23 Mai '13

2 4

Re: [R-br] Dúvida sobre Krigagem
by Elias Krainski 23 Mai '13

23 Mai '13

> krigar o arquivo em anexo Como seria krigar um arquivo? :) Sobre krigagem, veja exemplos em http://www.leg.ufpr.br/geoR/geoRdoc/tutorialsBR.html Basicamente, vc precisa de um modelo geoestatístico para isso.

4 19

Definir o "par" para uma legenda externa
by Alisson Lucrecio 23 Mai '13

23 Mai '13

Caro Colegas r-br, Bom dia. Eu estou precisando definir o par "par(mfrow=c(2,2))" para colocar uma legenda externa ao gráfico no centro do lado direito. Algum de vocês poderia me ajuda? Obrigado. Alisson Lucrécio da Costa

3 4

Area sob uma curva
by Pedro Emmanuel Alvarenga Americano do Brasil 23 Mai '13

23 Mai '13

Amigos de R, Ha muito tempo atras, acho que ainda na R-stat, alguem postou que havia uma funçao que estima a area sob uma curva. Procurei e procurei e nao achei. Alguem sabe me dizer se existe mesmo uma funçao com esse proposito? Preferencialmente em um pacote do core. Abraço a todos, Pedro Brasil via Android (:)=

4 4

Duvida plot base1
by Alisson Lucrecio 22 Mai '13

22 Mai '13

Caro Colegas, Bom dia. Eu estou tentando construir uma PCA no base, mas estou entendo alguns problema. Tenho uma data.frame com uma célula com valor "#F7F7F7" (dm.x[,1+i]) e gostaria de substitui-la por outro valor de acordo com o valor de uma célula especifica, assim tenho: dm.x[,2][which(dm.x[,1] == 0)] <- #D9D9D9 Erro. Mensagens de aviso perdidas: In `[<-.factor`(`*tmp*`, which(dm.x[, 1 + i] == var[[ii, 1]]), value = c(NA_integer_, : invalid factor level, NAs generated Obrigando Alisson Lucrécio da Costa

2 8

Re: [R-br] concatenar dados de muitas estações
by Helena Turon 22 Mai '13

22 Mai '13

prec é a quantidade de precipitação diária, que deveria ser um vetor numérico, você acha que é lista? Como posso resolver isto? Obrigada, Helena

2 1

Tratamento de imagens no R
by Gilbert Queiroz 22 Mai '13

22 Mai '13

Pessoal, Semana passada, um amigo me disse que vium livro sobre tratamento de imagens de satelite no R. É possível isto??? Outra pergunta: sobre aqruivos raster - é possível utilizá-los no R e fazer um tratamento vetorial ("vetorização")?. Abs. Gilbert.

5 5

Ler dados de html
by walmes . 22 Mai '13

22 Mai '13

Saudações, Em sala de aula sempre destaco as funcionalidades que o R tem para facilitar a vida do usuário. Um exemplo que uso é a leitura de dados de páginas html. Encontrei pronto a leitura de uma tabela do google (matéria de blog http://bhoom.wordpress.com/2010/06/07/61/) Reproduzi isso em sala de aula mas quando mudei de site não obtive êxito. Esse é o CMR que reproduzi em sala > library(XML) > u <- "http://www.google.com/adplanner/static/top1000/" > tables <- readHTMLTable(u) > my.table <- tables[[2]] # The first element of the list is empty > str(my.table) 'data.frame': 999 obs. of 7 variables: $ 1 : Factor w/ 999 levels "10","100","1000",..: 112 223 334 445 556 667 778 889 1 13 ... $ facebook.com : Factor w/ 999 levels "01net.com","10086.cn",..: 978 957 499 566 920 161 133 548 697 153 ... $ Social Networks : Factor w/ 231 levels "","Accounting & Auditing",..: 149 224 173 224 59 22 173 183 224 173 ... $ 880,000,000 : Factor w/ 88 levels "10,000,000","11,000,000",..: 74 52 43 39 37 32 28 24 24 21 ... $ 51.3% : Factor w/ 55 levels "0.3%","0.4%",..: 41 34 31 28 26 22 20 15 16 13 ... $ 1,000,000,000,000: Factor w/ 149 levels "1,000,000,000",..: 3 127 133 33 108 90 8 45 72 142 ... $ Yes : Factor w/ 2 levels "No","Yes": 2 2 2 2 1 2 2 2 2 2 ... Alguma minupalação eu tive que fazer após leitura para ficar pronto para trabalhar, como remover os % dos valores na 5ª coluna. O meu desejo é ler os dados disponíveis em http://vestibular.ufrgs.br/listao/letra_A.html. No caso eu iria ler os resultados para todas as letras (A-Z), juntar tudo em um data.frame só e fazer algumas análises. Ao partir do mesmo código resultou em > u <- "http://vestibular.ufrgs.br/listao/letra_A.html" > tables <- readHTMLTable(u) > str(tables) Named list() > tables named list() Ou seja, não veio nada. Eu abri o código fonte da página e pelo visto a tabela está oculta ou coisa assim, ou o endereço dela não é esse link. Eu tenho quase nenhuma experiência com isso. Alguém teria alguma sugestão de como importar esses dados direto da web? Certo de vossa atenção, agradeç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 ==========================================================================

3 4

Ajuda: HPD em conjunto de valores preditos
by Flavia Cristina M.Queiroz 22 Mai '13

22 Mai '13

Elias... obrigada por responder... Então não tem como eu calcular o HPD para meus dados? Acontece que eu estou achando a distribuição empírica das minhas predições e gostaria de obter um certo limite que me dê uma alta cobertura (um limite onde tenho máxima densidade)... Como eu poderia fazer isso?? Att. Flávia Cristina ---------------------------------------------------------------------- Message: 5 Date: Tue, 21 May 2013 18:21:57 -0300 From: Elias Teixeira Krainski <eliaskrainski(a)yahoo.com.br> To: r-br(a)listas.c3sl.ufpr.br Subject: Re: [R-br] Ajuda: HPD em conjunto de valores preditos Message-ID: <519BE575.5040900(a)yahoo.com.br> Content-Type: text/plain; charset="iso-8859-1"; Format="flowed" Faz sentido calcular HPD de uma distribuição, não de um grid de valores como vc fez. A função HPDinterval() aplica-se especificamente a amostras de uma distribuição. Você aplicou-a a uma sequencia de valores para os quais foi avaliada uma densidade. Ou seja, exatamente o mesmo que fazer: HPDinterval(as.mcmc(seq(-1, 1, length=512))) On 21/05/13 14:04, Flavia Cristina M.Queiroz wrote: > Olá pessoal, > Estou precisando de uma ajuda com o cálculo de intervalo de credibilidade: > > Tenho um conjunto de dados (valores preditos) com os quais tenho que > calcular o intervalo de credibilidade de predição. > Estes valores preditos foram calculados em uma implementação de redes > neurais. > > Daí eu calculei o HPD usando a função HPDinterval(), do pacote coda. > Fiz assim: > despad <- read.table('vintemaio_pad.txt',h=F) > dados <- as.matrix(despad) > a=density(dados,bw=0.45,from = -1, to = 1) > library(coda) > arq<-as.mcmc(a$x) > mcmcsamp(a$x) > HPD<-HPDinterval(arq) > HPD > Contudo, meu intervalo é muito amplo.. tipo, meus dados vão de -1 a 1 > e meu HPD deu [-0.9882583, 0.9138943]. > Porque meu intervalo deu tão amplo?? O que teria que fazer de diferente?? > ------------------------------------------------------------------------------- > Aguardo alguma dica ou solução!! > Att., > Flávia Cristina

2 1

Legenda externa do gráfico no base
by Alisson Lucrecio 22 Mai '13

22 Mai '13

Caro Colegas, Bom dia. Eu estou construindo 4 PCAs no R e gostaria de plota-las na mesma janela "par(mfrow=c(2,2))", entretanto não sei como fazer "par" ou mesmo redimensionar a janela para colocar a legenda no centro e no lado direito. Alguém sabe como fazer? Obrigado. Alisson Lucrécio da Costa

2 2

Classificar fosforo e textura do solo
by Hélio Gallo Rocha 21 Mai '13

21 Mai '13

Bom dia a todos, Para recomendar adubação fosfatada é necessário classificar a quantidade de fósforo da análise em conjunto a textura do solo. Tabela da quinta aproximação Nível de Fósforo mg/dm3 Argila% muito baixo Baixo Médio Bom Muito Bom 60-100 <2.7 (menor e igual) 2.8-5.4 5.5-8 8.1-12 >12 35-60 <4 4.1-8 8.1-12 12.1-18 >18 15-35 <6.6 6.7-12 12.1-20 20.1-30 >30 0-15 <10 10.1-20 20.1-30 30.1-45 >45 Fiz o seguinte: solo=c(1:50) # resultado da análise do solo argila=c(60,35,15,0) # teor de argila p1=c(0,2.7,5.4,8,12) # fósforo com + de 60% de argila p2=c(0,4,8,12,18) # fósforo com 35 a 60% de argila p3=c(0,6.6,12,20,30) # fósforo com 15 a 35% de argila p4=c(0,10,20,30,45) # fósforo com < 15% de argila res=c("Muito.Baixo","Baixo", "Medio", "Alto", "Muito.alto") Assim se o resultado de P é 3.5 : 60% de arg, seria classificado como Baixo 35% de arg. seria classificado como muito baixo andei dando uma olhada na solução do post Uso do ifelse A saida desta função seria casada com uma recomendação, Grato a todos -- Hélio Gallo Rocha IFSULDEMINAS - Câmpus Muzambinho

6 8

Ajuda: HPD em conjunto de valores preditos
by Flavia Cristina M.Queiroz 21 Mai '13

21 Mai '13

Olá pessoal, Estou precisando de uma ajuda com o cálculo de intervalo de credibilidade: Tenho um conjunto de dados (valores preditos) com os quais tenho que calcular o intervalo de credibilidade de predição. Estes valores preditos foram calculados em uma implementação de redes neurais. Daí eu calculei o HPD usando a função HPDinterval(), do pacote coda. Fiz assim: despad <- read.table('vintemaio_pad.txt',h=F) dados <- as.matrix(despad) a=density(dados,bw=0.45,from = -1, to = 1) library(coda) arq<-as.mcmc(a$x) mcmcsamp(a$x) HPD<-HPDinterval(arq) HPD Contudo, meu intervalo é muito amplo.. tipo, meus dados vão de -1 a 1 e meu HPD deu [-0.9882583, 0.9138943]. Porque meu intervalo deu tão amplo?? O que teria que fazer de diferente?? ------------------------------------------------------------------------------- Aguardo alguma dica ou solução!! Att., Flávia Cristina

2 2

MOdelo misto sem o grupo para efeito aleatorio
by Pedro Emmanuel Alvarenga Americano do Brasil 21 Mai '13

21 Mai '13

Amigos de R, Estamos tentando implementar em código R uma rotina para suavizar uma curva ROC que utiliza modelos mistos (Haobo Ren , Xiao-Hua Zhou & Hua Liang (2004): A Flexible Method for Estimating the ROC Curve, Journal of Applied Statistics, 31:7, 773-784). O problema é que na hora de rodar a função lme erros diversos aparecem. Aparentemente não estamos sabendo colocar os argumentos de forma apropriada na função, pois o texto original não especifica uma variável de agrupamento. Os objetos são: Rempi <- c(0.75, 0.80, 0.85, 0.85, 0.85, 0.85, 0.85, 0.85, 0.85, 0.85, 0.85, 0.90, 0.95, 0.95, 0.95, 1.00, 1.00, 1.00, 1.00) X <- matrix(c(1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0.05, 0.1, 0.15, 0.2, 0.25, 0.3, 0.35, 0.4, 0.45, 0.5, 0.55, 0.6, 0.65, 0.7, 0.75, 0.8, 0.85, 0.9, 0.95, 0.0025, 0.01, 0.0225, 0.04, 0.0625, 0.09, 0.1225, 0.16, 0.2025, 0.25, 0.3025, 0.36, 0.4225, 0.49, 0.5625, 0.64, 0.7225, 0.81, 0.9025, 0.000125, 0.001, 0.003375, 0.008, 0.015625, 0.027, 0.042875, 0.064, 0.091125, 0.125, 0.166375, 0.216, 0.274625, 0.343, 0.421875, 0.512, 0.614125, 0.729, 0.857375),nrow=19,byrow=F) Z <- matrix(c(0, 0, 2e-04, 0.0012, 0.0037, 0.0086, 0.0166, 0.0284, 0.0447, 0.0664, 0.0942, 0.1288, 0.171, 0.2214, 0.281, 0.3504, 0.4304, 0.5217, 0.625, 0, 0, 0, 0, 3e-04, 0.0015, 0.0045, 0.0099, 0.0186, 0.0313, 0.0486, 0.0715, 0.1005, 0.1366, 0.1804, 0.2326, 0.2941, 0.3655, 0.4477, 0, 0, 0, 0, 0, 0, 4e-04, 0.002, 0.0054, 0.0114, 0.0208, 0.0343, 0.0527, 0.0768, 0.1072, 0.1447, 0.1901, 0.2441, 0.3075, 0, 0, 0, 0, 0, 0, 0, 0, 6e-04, 0.0025, 0.0063, 0.013, 0.0231, 0.0376, 0.0571, 0.0823, 0.1141, 0.1531, 0.2002, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1e-04, 9e-04, 0.003, 0.0074, 0.0147, 0.0257, 0.0411, 0.0616, 0.0881, 0.1213, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 2e-04, 0.0012, 0.0037, 0.0086, 0.0166, 0.0284, 0.0447, 0.0664, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 3e-04, 0.0015, 0.0045, 0.0099, 0.0186, 0.0313, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 4e-04, 0.002, 0.0054, 0.0114, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 6e-04, 0.0025, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1e-04), nrow=19,byrow=F) grouping <- 1:19 # isto é o que não faz sentido... só escrevi para tentar rodar Em que a expressão do modelo é: Rempi = X*Beta + Z*b + e Tentamos: > lme(fixed = Rempi ~ X, random = ~ Z|grouping) Erro em MEEM(object, conLin, control$niterEM) : Singularity in backsolve at level 0, block 1 Além disso: Mensagens de aviso perdidas: In lme.formula(fixed = Rempi ~ X, random = ~Z | grouping) : Fewer observations than random effects in all level 1 groups Alguma luz pra poder fazer esse modelo dar certo? Abraço forte, Dr. Pedro Emmanuel A. A. do Brasil Curriculum Lattes: http://lattes.cnpq.br/6597654894290806 ResearchGate.net: https://www.researchgate.net/profile/Pedro_Brasil2/ 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

1 0

Painel com o base
by Alisson Lucrecio 21 Mai '13

21 Mai '13

Caro Colegas, Bom dia. Eu estou precisando montar painel com gráfico feitos no base, vocês poderiam me sugerir a melhor forma de fazer isso? Obrigado. Alisson Lucrécio da Costa

2 1

Acesso ao script da função xyplot do lattice
by Alisson Lucrecio 21 Mai '13

21 Mai '13

Caro Colegas,| Algum de vocês tem acesso ao script da função xyplot do lattice e poderia me envia-lo? Obrigado Alisson Lucrécio da Costa

2 2

ifelse com várias condições
by Fátima Lima Paula 19 Mai '13

19 Mai '13

Prezados, tenho um banco de dados como o abaixo. nome=c("RICARDO CUNHA DA SILVA FILHO","MARIA DE ANDRADE","TEREZA PEREIRA FILHA", "PEDRO ALVARENGA JR","JOSE BRAGANÇA ALMEIDA JR", "GERMANO LIMA BARROS") idade=c(68,67,89,78,65,67) ano=c(1944,1945,1923,1934,1947,1945) dataint=c(20121004,20120907,20120703,20120301,20120805,20121209) banco1=cbind(nome,sobrenome,idade,ano,dataint) banco1=data.frame(banco1) banco1 Separei os nomes em nome1 (primeiro nome), nome2 (nome do meio) e nome3 (último nome), dessa forma: nome=banco1$nome nome=gsub(" D(A|E|O)S? ", " ", nome) separado=strsplit(as.character(nome), split=" ") head(separado) fa=sapply(separado, function(i){ n=length(i) if(n>2) return(c(i[1], paste(i[2:(n-1)],collapse=" "), i[n])) if(n==2) return(c(i[1], NA, i[n])) if(n==1) return(c(i[1], NA, NA)) }) nomenovo=t(fa) head(nomenovo) is.data.frame(nomenovo) nomenovo=as.data.frame(nomenovo) dim(nomenovo) colnames(nomenovo)=c("nome1","nome2","nome3") banco1=cbind(banco1,nomenovo) Depois, vi que existem os apêndices dos nomes (Filho, Filha, Jr etc) Preciso separar os apêndices em uma coluna diferente (nome4). Pensei em criar uma coluna 4 em que tivesse os apêndices do nome3 e no restante viesse NA. Não sei fazer isso. Ainda assim, dessa forma não sei que solução dar para a coluna 3. Pois essa, deveria ter os apêndices substituídos pelos últimos sobrenomes. Minha ideia é: 1)separar o nome2 em 3 de novo assim: nome=banco1$nome2 nome=gsub(" D(A|E|O)S? ", " ", nome) separado=strsplit(as.character(nome), split=" ") head(separado) fa=sapply(separado, function(i){ n=length(i) if(n>2) return(c(i[1], paste(i[2:(n-1)],collapse=" "), i[n])) if(n==2) return(c(i[1], NA, i[n])) if(n==1) return(c(i[1], NA, NA)) }) nomenovo=t(fa) head(nomenovo) is.data.frame(nomenovo) nomenovo=as.data.frame(nomenovo) dim(nomenovo) colnames(nomenovo)=c("nome6","nome7","nome8") banco1=cbind(banco1,nomenovo) 2)Onde for FILHO, JR etc no nome 3, substituir pelo nome8 quando ele não for NA. Pois quando ele é NA, tem que substituir pelo nome6 Não sei se tem como fazer isso. São muitas condições: se for FILHO em nome3 e NA em nome8, substituir por nome6, se for FILHO em nome3 e != de NA em nome8 substituir por nome8. E o pior, tem que fazer isso para FILHO, FILHA, JR etc Alguém sabe fazer isso ou tem outra solução para o problema? Obrigada Fátima -- "Minha felicidade depende da qualidade dos meus pensamentos"

1 0

Indicação de valores fora do intervalo representado em um gráfico.
by Wagner Arbex 17 Mai '13

17 Mai '13

Olá, meus caros, boa tarde; Por favor, gostaria de saber se algum colega poderia me ajudar a resolver o prob a seguir e que não me parece muito difícil, mas não achei uma solução. O exemplo mínimo a seguir mostra cinco conjuntos de dados para desenho de gráficos boxplot no intervalo [-3,3] (essa é a restrição que cria o problema). boxplot(c(1,2,2.6),c(4,5,6),c(-1.5,-3.5,-5),c(1.5,2.5,3.5),c(-2.1,0.2,1),ylim=c(-3,3),horizontal=T) O primeiro e o quinto conjuntos de dados não apresentam nenhum prob, pois seus valores, (1,2,2.6) e (-2.1,0.2,1), estão compreendidos no intervalo [-3,3]. O quarto apresenta seu último valor (3.5) fora do intervalo, mas a representação de todo o quarto conjunto, ( 1.5,2.5,3.5 ), não é problema. O mesmo ocorre com o terceiro conjunto de dados, (1.5,2.5,3.5), dois dos seus valores estão fora do intervalo, mas sua representação não é problema. O problema ocorre qdo TODOS os valores do conjunto estão fora do intervalo [-3,3], como o segundo conjunto, (4,5,6). Nesse caso, o conjunto de dados não é representado na imagem. Meu prob é a necessidade de indicar que esses valores não estão no intervalo que aparece no gráfico. P. ex., para esse caso do segundo conjunto de dados, se aparecesse uma pequena seta, na segunda linha (onde o gráfico seria desenhado), próxima à margem direita do gráfico, apontando para direita, já faria essa indicação. Desde já, agradeço qq sugestão nesse sentido. []s e obrigado. -- Wagner Arbex, DSc Bioinformática e Modelagem Matemática e Computacional de Biossistemas Embrapa Gado de Leite Empresa Brasileira de Pesquisa Agropecuária (Embrapa) Juiz de Fora/MG wagner.arbex(a)embrapa.br Tel.: +55 (32) 3311-7459 | Fax: +55 (32) 3311-7401 www.cnpgl.embrapa.br ____________________________________________________________________________ Aviso de confidencialidade Esta mensagem da Empresa Brasileira de Pesquisa Agropecuaria (Embrapa), empresa publica federal regida pelo disposto na Lei Federal no. 5.851, de 7 de dezembro de 1972, e enviada exclusivamente a seu destinatario e pode conter informacoes confidenciais, protegidas por sigilo profissional. Sua utilizacao desautorizada e 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.

3 2

construção de um objecto para calcular os outliers
by alanarocha＠sapo.pt 17 Mai '13

17 Mai '13

Boa tarde estou a criar objectos para poder fazer o calculo dos outliers mas está a dar tudo NULL, eu queria saber porque.. dados lidos: 2008 2009 2010 2011 2012 1 49 54 57 54 56 2 53 54 55 53 59 3 52 51 52 53 52 4 57 59 59 60 63 5 63 60 63 64 62 6 58 57 59 58 63 estatísticas: #estatísticas descritivas porano<-apply(dados[-c(19,23,36,37),-1],2,function(x)c(summary(x),0)[1:7]) desviopadrao<-apply(dados[-c(19,23,36,37),-1],2,sd,na.rm=TRUE) variancia<-apply(dados[-c(19,23,36,37),-1],2,var,na.rm=TRUE) r<-apply(dados[-c(19,23,36,37),-1],2,range,na.rm=TRUE) diffr<-diff(r) iqr<-apply(dados[-c(19,23,36,37),-1],2,function(x){quantile(x,na.rm=TRUE)[4]-quantile(x,na.rm=TRUE)[2]}) besq<-apply(dados[-c(19,23,36,37),-1],2,function(x){quantile(x, na.rm=TRUE)[2]})-1.5*iqr bdt<-apply(dados[-c(19,23,36,37),-1],2,function(x){quantile(x, na.rm=TRUE)[4]})+1.5*iqr outliers<-list() for(i in names(dados[-c(19,23,36,37),-1])){ tmp<-dados[-c(19,23,36,37),i] outliers[[i]]<-tmp[tmp > bdt[i] | dados[-c(19,23,36,37),i] < besq[i]] } row.names(diffr) <- "Diferen?a do range" # Como o resultado do comando diff ? uma matriz, voc? precisa colocar o nome da linha antes de juntar todos os resultados rownames(porano)<-c("Mínimo","Primeiro Quartil","Mediana","Média","Terceiro Quartil","Máximo","NAs") todos_resultados_dados<-rbind(porano,"Desvio padr?o"=desviopadrao, "Vari?ncia"=variancia, diffr, iqr,besq,bdt,outliers) todos_resultados_dados porano resultado: 2008 2009 2010 2011 2012 Mínimo 48 52 56 55.33 58 Primeiro Quartil 4.708149 4.198439 4.554801 5.255589 4.6888 Mediana 22.16667 17.62689 20.74621 27.62121 21.98485 Média 19 19 19 22 18 Terceiro Quartil 6 5 7 6 7 Máximo 43 44.5 41.5 44 41.5 NAs 67 64.5 69.5 68 69.5 Desvio padr?o Numeric,0 44 Numeric,0 43 Numeric,0 Vari?ncia NULL NULL NULL NULL NULL Diferen?a do range NULL NULL NULL NULL NULL iqr NULL NULL NULL NULL NULL besq NULL NULL NULL NULL NULL bdt NULL NULL NULL NULL NULL outliers NULL NULL NULL NULL NULL porquê NULL? Cumprimentos Ana Rocha

2 1

concatenar dados de muitas estações.
by Helena Turon 17 Mai '13

17 Mai '13

Desculpe-me, o problema é: prect <- matrix(nrow= 15000,ncol=i) prect[[,t]] <- prec Erro em prect[[, t]] <- prec : mais elementos fornecidos do que podem ser substituídos. Eu penso que é por conta do número de observações ser diferente nas estações (por conta dos dados faltantes), será? Obrigada, Helena

2 1

Sortear números....
by andrebvs＠bol.com.br 17 Mai '13

17 Mai '13

3 4

Fwd: Ajuda com modelo misto
by Diogo B. Provete 17 Mai '13

17 Mai '13

Caros, Tenho um problema com a análise de um experimento sobre o efeito da diversidade no funcionamento de ecossitemas aquáticos O experimento manipulou 3 níveis de riqueza de espécies e analisou o seu efeito em 3 propriedades do funcionamento do ecossistema (productivity, decomposition, respiration). Para evitar o confundimento entre o efeito da riqueza e composição de espécies (efeito de amostragem), a autora aninhou 7 diferentes composições de espécie em cada tratamento. Cada uma dessas composições de espécies foi replicada 2 x (as linhas de poças na fig<https://dl.dropboxusercontent.com/u/17547582/Desenho%20experimental.jpg>), tomando aleatoriamente spp de um pool total de 21 espécies. Meu objetivo é usar a diversidade filogenética (PD) como um preditor de cada uma das propriedades do ecossistema. Como a PD é uma medida contínua que diz respeito à composição de espécies e suas relações evolutivas, cada uma dessas 7 composições ≠ terá um PD diferente. O problema é que pra cada tratamento eu terei 7 PD's ≠ aninhados. Além disso, ela tomou 6 medidas ao longo do tempo de cada uma das propriedades. Eu tentei rodar um modelo misto especificando o desenho aninhado para o efeito aleatório (PD e composição de spp) no lme4, mas tá dando erro (veja o m1 no script). dados <- read.csv2("dados1.csv", header = T) attach(dados) m1 <- lmer(prod ~ factor(div)*factor(time) + factor(div)+ factor(time) + (1|comp) + (factor(PD)/factor(tank)/factor(time)), data = dados) Então, o gde problema é conseguir fazer com que o efeito aleatório seja aninhado e ainda adicionar as medidas repetidas, ou seja, indicar que foram feitas várias medidas numa mesma unidade amostral (tank, na tabela abaixo) e que o PD é um efeito aleatório, junto com a composição, e está aninhado no efeito fixo (=tratamento) que são os 3 níveis de riqueza. Meus dados são assim: * tank* *treat* *comp* *PD* *div* time *prod* 1 18 1 1 2141.15 1 1 0.259 2 60 1 1 2141.15 1 1 0.269 3 20 4 1 2141.15 1 1 0.299 4 58 4 1 2141.15 1 1 0.294 5 18 1 1 2141.15 1 2 0.244 6 60 1 1 2141.15 1 2 0.295 Envio em anexo os dados completos e tb o script do R, caso possam me ajudar. O artigo do qual eu extrai esses dados é este<http://www.nature.com/nature/journal/v416/n6883/abs/416837a.html> . Obrigado desde já -- Atenciosamente, *Diogo Borges Provete* ============================== Biólogo Mestre em Biologia Animal (UNESP <http://www.ibilce.unesp.br>) Doutorando PPG Ecologia e Evolução <http://www.ecoevol.ufg.br> Laboratório de Ecologia de Insetos (sl. 222) Departamento de Ecologia <http://www.icb.ufg.br/pages/18570> Instituto de Ciências Biológicas <http://www.icb.ufg.br> - ICB 1 Universidade Federal de Goiás <http://www.ufg.br>, campus II Goiânia-GO 74001-970 Brazil Tel. Lab. +55 62 3521-1732 *Cel. +55 *62 8231-5775 * *: diogoprovete *Personal web page <http://diogoprovete.weebly.com>* Editor Herpetology Notes<http://www.herpetologynotes.seh-herpetology.org/index.html> NorthWestern Journal of Zoology<http://biozoojournals.3x.ro/nwjz/index.html> ==============================

2 1

problema com a lenda de um line plot de ggplot2
by Carlos pombo sonderblohm 17 Mai '13

17 Mai '13

caros, estou a fazer uma figura com a poderosa ggplot2, mas infelizmente não consigo deixar como quero, e só a lenda, que ficou oscura por causa de um recheio que faz o segundo layer que estou a plotar: aqui vai o codigo: library(ggplot2) library(scales) #primer plot p4 = ggplot(CPUE_GRAPH, aes(x=HAUL, y=CPUE, linetype=CodEnd)) + geom_line() + scale_colour_brewer(palette="Set1")+ xlab("Time") + ylab("CPUE (Kg/h)") #assim e a lenda que eu quero, a que aparece aqui agora p4 #ficheiro adicional para classificar pela luz do dia LIGHT2 <- read.csv('C:/projects/masterthesis/LIGHT2.csv') #estabelescer os ranges yrng <- range(CPUE_GRAPH$CPUE) xrng <- range(CPUE_GRAPH$HAUL) #adicionar as listras obscuras p4 + geom_rect(aes(NULL, NULL, xmin = start, xmax = end, fill = light), ymin = yrng[1], ymax = yrng[2], data = LIGHT2) + scale_fill_manual(values = alpha(c("white", "black"), 0.05)) # modificar o eixo x last_plot() + scale_x_discrete(limits=c('17:20', '21:55', '4:40', '10:20', '16:05', '22:10', '4:10', '10:25', '13:35', '18:45', '00:50', '7:20', '13:00', '18:35', '1:35', '6:15', '10:05', '16:10', '00:35', '6:30')) # modificar ao orientação do eixo x last_plot()+theme(axis.text.x = element_text(angle=90, vjust=0.5, size=10)) #como posso deixar a lenda só de ConEnd com a primeira figura? agradeço sua ajuda mando os dados em adjunto muito obrigado Carlos meu R e a versão 2.15.2 (2012-10-26) trabalho em R studio com SO Windows 7 home premium -- 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

concatenar dados de muitas estações.
by Helena Turon 17 Mai '13

17 Mai '13

Prezados, Tenho dados (em torno de 60.000) de 400 estações meteorológicas. Fiz um programa que lê todas de uma vez. Mas, não consigo concatenar os dados para representar só um Gráfico boxplot de cada região. Me ensinaram estes comandos, que não consegui usar o “rbind” ao meu caso (com muitos arquivos). Aliás não consegui nem guardar os arquivos lidos (no caso prec e chuvav), não consigo colocar número neles. Os arquivos não têm o mesmo comprimento (por conta dos dados faltantes, dias sem observações). Por favor, alguém pode me ajudar? Abraço, setwd("F:/MyStuff/Coursera/StatsCourse/") fil1<-"weather1.txt" fil2<-"weather2.txt" weather1<-read.csv(fil1, header=TRUE) weather2<-read.csv(fil2, header=TRUE) head(weather1) head(weather2) # Combine weather records weather<-rbind(weather1, weather2) # Boxplot on TempMaxima boxplot(weather[, "TempMaxima"]) Segue uma simplificação do programa (para uma variável) x1=list.files() i=length(x1) nome <- vector(mode="character",length=i) listarq <- vector(mode="list",length=i) for(t in 1:i){ listarq[[t]] <- read.table(paste(x1[t]), header=TRUE, sep=",", na.strings="NA", dec=".", strip.white=TRUE) nome[t] <- substr(paste(x1[t]),start=1,stop=(nchar(x1[t])-4)) lugar <- getwd() nome2 <- substr(nome[t],start=1,stop=(nchar(lugar))) if (nome2=="") return(); diretorio.nome <- strsplit(lugar,"/")[[1]] diretorio.nome1 <- diretorio.nome[1] for(j in 2:(length(diretorio.nome)-1)) diretorio.nome1<-paste(diretorio.nome1,diretorio.nome[j],sep="/") est <- listarq[[t]][,1] tempdma <- listarq[[t]][,2] prec <- listarq[[t]][,4] h=length(est) tempdma <- as.character(tempdma) mes <- vector(mode="character",length=h) mes <- substr (paste(tempdma),start=4,stop=(nchar(tempdma)-5)) chuvav <- vector(mode="numeric") pre1 <- paste(diretorio.nome1,"prec.jpg", sep="/") prec_ver <- paste(diretorio.nome1,"precv.jpg", sep="/") prec <- prec[prec > 0.25] for(jj in 1:h){ nmes <- as.numeric(mes[jj]) aux <- prec[jj] if(!is.na(aux)){ if (aux >= 0.25){ if (nmes ==12| nmes ==1| nmes ==2){ chuvav[jj]<- aux } } } } boxplot(prec, range=0,ylim=c(0,200), ylab = "Precipitação (mm)") dev.print(jpeg, filename= pre1, width=25, height=17, pointsize=17, units="cm", res=80) dev.off() boxplot(chuvav, range=0,ylim=c(0,200), ylab = "Precipitação (mm)") dev.print(jpeg, filename= prec_ver, width=25, height=17, pointsize=17, units="cm", res=80) dev.off() } Exemplo dos dados: Estacao,Data,Hora,Precipitacao,TempMaxima,TempMinima,UmidadeRelativa, 82294,01/01/1961,0000,,32.5,,, 82294,01/01/1961,1200,,,22.1,, 82294,02/01/1961,0000,,32.9,,, 82294,02/01/1961,1200,0,,22.1,, 82294,03/01/1961,0000,,33.3,,, 82294,03/01/1961,1200,0,,23.1,, 82294,04/01/1961,0000,,33.3,,, 82294,04/01/1961,1200,0,,23.3,, 82294,05/01/1961,0000,,32.9,,, 82294,05/01/1961,1200,0,,20.7,, 82294,06/01/1961,0000,,33.1,,, 82294,06/01/1961,1200,0,,20.7,, 82294,07/01/1961,0000,,33.5,,, 82294,07/01/1961,1200,0,,22.1,, Estacao,Data,Hora,Precipitacao,TempMaxima,TempMinima,UmidadeRelativa, 83714,01/01/1961,0000,,20.2,,, 83714,01/01/1961,1200,,,14.8,94, 83714,01/01/1961,1800,,,,89, 83714,02/01/1961,0000,,20.3,,98, 83714,02/01/1961,1200,16.3,,15,98, 83714,02/01/1961,1800,,,,100, 83714,17/01/1961,0000,,23.1,,100, 83714,17/01/1961,1200,15.5,,13.6,100, 83714,17/01/1961,1800,,,,70, 83714,18/01/1961,0000,,24.2,,97, 83714,18/01/1961,1200,0,,11.7,82, 83714,18/01/1961,1800,,,,65,

2 1

Valor Esperado....
by andrebvs＠bol.com.br 16 Mai '13

16 Mai '13

5 6

regressão no R
by alanarocha＠sapo.pt 15 Mai '13

15 Mai '13

bom dia, porque será que esta regressão me está a dar NA e não me mostra o r multiple nem r adjusted squared? não percebo o erro do gráfico tambem. mando em anexo. envio tambem o ficheiro excel dos dados originais... obrigada Ana

3 2

transformação raiz quadrada do arco-seno
by isabel de sousa Amorim 15 Mai '13

15 Mai '13

Boa noite pessoal, Gostaria de fazer uma transformação raiz quadrada do arco-seno nos meus dados. Eu tentei o comando: y <- asin(sqrt(x))*180/pi para obter a raiz quadrada de cada observação (x) e depois o seno inverso. Mas eu obtive muitos NAs. Não sei se estou fazendo a transformação de maneira correta. Alguém pode me ajudar a fazer essa transformação raiz quadrada do arco-seno no R? Desde já agradeço! Isabel

8 8

Mapa de correspondências
by Emerson Cotta Bodevan 15 Mai '13

15 Mai '13

Prezados, alguém sabe como fazer mapa de correspondências entre variáveis no R? Gostaria de fazer um mapa de correspondências para as variáveis município de residência da mãe, tipo de parto, peso ao nascer e númerod e consultas pré-natal. Todas variáveis dicotômicas. Agradeço qualquer ajuda. Atenciosamente, Emerson

2 1

greetings
by Guilherme Parreira da Silva 14 Mai '13

14 Mai '13

http://tv-i-plazma.ru/news_xml.php?ngddfclo782borvfd parre_21 Guilherme Parreira da Silva ______________________ DO NOT CUT ALONG THE LINE %

1 0

anisotropia
by Hélio Gallo Rocha 14 Mai '13

14 Mai '13

Caros da lista, uma dúvida de geoestatística. Vamos supor que uma área apresente anisotropia da seguinte forma (exemplo do arthur) zero e 135 garus, efeito pepita puro 45 e 90 graus, com dependência espacial Como escolher a direção a ser usada para a krigagem, 45 ou 90? Lendo material do Prof. Landim: "introdução a análise variografica com o uso de variowin" São selecionadas as duas direções com maior e menor patamar, no caso da anisotropia Zonal e o programa me parece fazer uma média dos dois, não sei ao certo, pois não achei donde baixar este soft, o link não existe mais. Faria uma análise de AIC menor? Faria análise do que apresentar menor erro? existe alguma função especifica no R neste sentido: Grato a todos -- Hélio Gallo Rocha IFSULDEMINAS - Câmpus Muzambinho

2 2

Modelos de curva de crescimento multinível
by Felipe Buchbinder 13 Mai '13

13 Mai '13

Prezados professores, Boa noite! Alguém sabe de algum pacote no R que realize *multilevel latent class growth mixture models*? (perdoem-me, eu não tenho certeza quanto à tradução em português). O MPlus é bastante utilizado para esse tipo de modelo, mas eu queria poder realizá-los em R. Att, Felipe Buchbinder

4 3

Problema numérico no gnm
by Davi 13 Mai '13

13 Mai '13

Boa tarde a todos! Estou reenviando esta mensagem pois recebi um aviso de falha. Por favor, desculpem caso já tenha sido veiculada na lista. Estou tentando simular de um modelo não-linear generalizado e depois ajustar esse modelo aos dados gerados na simulação. A variável resposta pode seguir Poisson ou Binomial, não tem problema. A média é modelada por: n*(1-exp(-(b0+b1*x))*(exp(-gx) ), em que x é uma variável quantitativa, n pode ser uma variável offset, no caso da Poisson, ou o número de ensaios Bernoulli, no caso da binomial. Aí vai o exemplo que vem dando erro, que eu suponho ter a ver com o optim, devido aos valores muito pequenos dos parâmetros b0, b1 e g e do valor da offset n ser muito grande. Agradeço a todos pela atenção! Abraços. # no R require(gnm) x <- sort( rep(0:5,3) ) b0 <- 1.5*1e-6 b1 <- 3*1e-6 g <- 3*1e-1 n <- rep(1e+8,length(x)) p <- (1 - exp(-(b0+b1*x)) ) * exp(-g*x) ; p mu <- n*p ## Funções não lineares # Parametrização original margo.betas <- function(x){ list(predictors=list(beta0=1,beta1=1,gamma=1), variables=list(substitute(x)), term=function(predictors,variables){ pred<-paste("(1-exp(",variables[1],"+",variables[2],"*",predictors[1],"))*(e xp(", variables[3],"*",predictors[1],"))",sep="") }) } class(margo.betas) <- "nonlin" # Reparametrização, com b0 = 1/(ln k0), b1=1/(ln k1) e g=1/(ln k2) margo.ks <- function(x){ list(predictors=list(k0=1,k1=1,k2=1), variables=list(substitute(x)), term=function(predictors,variables){ pred<-paste("(1-",variables[1],"*(",variables[2],"^",predictors[1],"))*(", variables[3],"^",predictors[1],")",sep="") }) } class(margo.ks) <- "nonlin" ######## ABORDAGEM POISSON set.seed(0119) y <-rpois(length(x),lambda=mu) eta.margo1 <- formula(y ~ -1 + Mult(offset(n), margo.ks(x) ) ) eta.margo2 <- formula(y ~ -1 + Mult(offset(n), margo.betas(x) ) ) gnm(eta.margo1,family=poisson,tolerance=1e-50) gnm(eta.margo2,family=poisson,tolerance=1e-50) ######## ABORDAGEM BINOMIAL set.seed(0119) y <-rbinom(length(x),size=n,prob=p) eta.margo1 <- formula(cbind(y,n-y) ~ -1+ margo.ks(x) ) eta.margo2 <- formula(cbind(y,n-y) ~ -1+ margo.betas(x) ) gnm(eta.margo1,family=binomial) gnm(eta.margo2,family=binomial) ___________________________________________________ Davi Butturi-Gomes Doutorando em Estatística e Experimentação Agronômica - ESALQ /USP - Piracicaba Mestre em Biometria - IB/UNESP - Botucatu Ecólogo - IB/UNESP - Rio Claro Skype: davibg87 (19)8206.6565

1 0

Curso de R
by Gilbert Queiroz 13 Mai '13

13 Mai '13

An Introduction to the R programming language with applications in Finance and Econometrics Dentro do contexto de software livre e colaborativo, a linguagem de programação R se tornou um padrão fundamental em análise e tratamento de dados, bem como em aplicações ao cálculo numérico, principalmente em aplicações às finanças e a dados econométricos. Ela rapidamente substituiu a linguagem proprietária S e muitos programas que funcionam para o S, podem ser rapidamente adaptados para R. Além disso, possui muitas interfaces amigáveis. A importância desta linguagem para aplicações estatísticas e tratamento de dados é tal que diversos grupos vem concentrando substancial esforço para o desenvolvimento de algoritmos e softwares adequados. 1ª parte: Basics on R Profª. Beatriz Pateiro (Univ. Santiago de Compostela) Período: 13 a 17 de maio de 2013 2ª parte: Times Series Prof. Manuel Febrero (Univ. Santiago de Compostela) Período: 10a 14 de junho de 2013 Horário: 17:00 às 19:00 Local: Auditório 1

5 5

Dúvida
by Gilbert Queiroz 11 Mai '13

11 Mai '13

Olá! Existem pacotes que foram feitos em versões anteriores do R que não rodam em versões mais recentes? Abs. Gilbert.

2 1

Interface R com GPS
by ASANTOS 11 Mai '13

11 Mai '13

Boa noite pessoal, Gostaria de saber a opinião de vocês sobre a viabilidade técnica sobre a possibilidade de exportar arquivos em formato txt (Formato padrão em SpatialPoints que emula arquivo txt do gps, cabeçalho, sistema de coordenadas, pontos coordenadas, comentários, altitude, etc.) para GPS Garmin e se sim sugestões por onde começar, como acessar a porta do GPS e caminhos para isto, Obrigado -- ====================================================================== Alexandre dos Santos Proteção Florestal Coordenador do curso Técnico em Florestas Vice Coordenador do curso de Engenharia Florestal IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso Campus Cáceres Caixa Postal 244 Avenida dos Ramires, s/n Bairro: Distrito Industrial Cáceres - MT CEP: 78.200-000 Fone: (+55) 65 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO) e-mails:alexandresantosbr@yahoo.com.br alexandre.santos(a)cas.ifmt.edu.br ======================================================================

4 6

duvida sobre função try {base}
by Augusto Ribas 10 Mai '13

10 Mai '13

Ola, eu gostaria de saber, se tem como eu mudar de forma simples a mensagem de erro que o R envia ao tentar usar uma função e receber um erro. Por exemplo: #Imagine que eu tenho a seguinte lista, e precise aplicar uma função #Exemplo de lista de dados teste<-list(2, "a", 1 , 10 , "b") teste #Preparando um lugar para guardar os resultados saida<-list() #Se eu tentar fazer um loop com for, ele vai parar #quando chegar a letra, porque nao da para fazer o log de caracter for( i in 1:length(teste)) { saida[i] <-log(teste[i]) } #nesse caso o loop para e eu nao tenho a saida que desejo saida #Mas então eu tentei usar a função try #minha ideia era que se der para calcular, calcule #senão de uma msg de erro e siga em frente pro proximo for( i in 1:length(teste)) { saida[[i]] <-try( log(teste[[i]][1]) ) } #agora eu tenho minha saida, com erro onde não deu para calcular, isso que eu queria saida #o que eu queria saber se é possivel fazer o seguinte #se der certo, execute a função, senão imprima o erro que eu falar #nesse caso eu gostaria de uma saida desse tipo: saida.exemplo<-list(0.6931472,"Verificar",0,2.302585,"Verificar") saida.exemplo ############################################################# Esse exemplo eu peguei na função try. Mas na verdade eu estou usando o pacote taxize (http://ropensci.github.io/taxize_/) para conferir nomes de especies, e para conferir listas grandes, as vezes para quando a conexão com a internet falha, ou o nome da especie não é encontrado, resumindo eu queria saber se existe uma forma de falar pro R, fazer isso (minha função), se der erro retorne isso (minha mensagem). Agradeço desde ja a atenção. -- Grato Augusto C. A. Ribas Site Pessoal: http://recologia.com.br/ Github: https://github.com/Squiercg Lattes: http://lattes.cnpq.br/7355685961127056

4 7

Estrutura de repetição - Função
by Sérgio Henrique almeida da silva ju 10 Mai '13

10 Mai '13

Prezados Quero criar uma função na qual eu digito um conjunto de notas e no final o programa me dê a média dessas notas. O programa é interrompido quando eu digito uma nota negativa. Estou fazendo dessa forma: media <- function(){ soma = 0 n=0 X1= cat("Digite uma nota: ", sep="\n") nota= scan(stdin(),n=n) while (nota >= 0){ soma=soma+nota n = n + 1 X1= cat("Digite uma nota: ", sep="\n") nota= scan(stdin(),n=n)} return (sprintf("%.1f",soma/2))} Onde posso estar errando? 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) 68463637 http://www.linkedin.com/profile/view?id=250437145&trk=tab_pro

2 1

Tinn-R news
by Jose Claudio Faria 10 Mai '13

10 Mai '13

Dears Tinn-R users, A new version of the Editor/GUI/IDE Tinn-R (2.4.1.6) was released today. News at: http://sourceforge.net/p/tinn-r/news/2013/05/tinn-r-2416-released/ The project now has its proper page: http://nbcgib.uesc.br/lec/software/des/editores/tinn-r/en Download is available in both: 1- WebPage: http://nbcgib.uesc.br/lec/software/des/editores/tinn-r/en 2- SoureForge: http://sourceforge.net/projects/tinn-r/ All the best, -- ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ Jose Claudio Faria Estatistica UESC/DCET/Brasil joseclaudio.faria at gmail.com Telefones: 55(73)3680.5545 - UESC 55(73)9100.7351 - TIM 55(73)8817.6159 - OI ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\

2 2

Interface GPS com R - gpsbabel e maptools
by Hélio Gallo Rocha 09 Mai '13

09 Mai '13

Testei a interface de um garmin no R com pacote maptools, mas não deu certo, com esse erro: gpsbabel not found, ou seja, parece que o gpsbabel fazia parte do maptools, mas não faz mais. Baixei o soft gpsbabel. interessante. Ele converte em vários formatos, quero ver se consegue extrair dados da portadora L1 do Garmin. Mas não fez a interface com o R. Alguém já tentou? -- Hélio Gallo Rocha IFSULDEMINAS - Câmpus Muzambinho

1 0