Março 2013 - R-br

sub escrito
by Ze Henrique 13 Mar '13

13 Mar '13

Pessoal como que eu deixo um texto em Italico e/ou negrito e/ou subscrito. Obrigado a todos, desde já. -- *José Henrique Soler Guilhen* Graduando em Agronomia Universidade Federal do Espírito Santo - CCAUFES

4 6

erros no codigo do R
by alanarocha＠sapo.pt 13 Mar '13

13 Mar '13

Boa tarde, alguem me pode ajudar a corrigir estes erros? DADOS: Clínicas X2008 X2009 X2010 X2011 X2012 1 PBarca 82 87 87 87 87 2 Braga 84 85 86 89 87 3 Fafe 85 87 88 86 94 4 Maia 81 86 84 90 95 5 VNGaia 88 93 92 93 95 6 SMFeira 77 83 81 84 92 # Estatística por ano porano<-apply(dados[,-1],2,function(x)c(summary(x),0)[1:7]) desvio padrão<-sd(dados[,-1], na.rm= TRUE) Error: unexpected symbol in "desvio padrão" variância<-diag(var(dados[,-1], na.rm= TRUE)) quantile(dados[,-1], probs = c(0.25, 0.75, NA)/100) ampinterquartil<-Iqr(matrix(as.numeric(dados[,-1], na.rm= TRUE))) Error in quantile(as.numeric(x), c(0.25, 0.75), na.rm = na.rm, names = FALSE, : (list) object cannot be coerced to type 'double' amplitude<-apply(dados[,-1],2,range, na.rm= TRUE) amplitudediff<-c(amplitude[2,] - amplitude[1,]) rownames(porano) <- c('Minimo','Primeiro quartil','Mediana','Média','Terceiro Quartil','Máximo','NAs') poranol<-rbind(porano,desvio padrão,variância,amplitude,amplitudediff,ampinterquartil) > colnames(porano) = 2008:2012 > poranol X2008 X2009 X2010 X2011 X2012 Minimo 66.000000 65.000000 62.000000 59.000000 60.000000 Primeiro quartil 76.000000 84.000000 86.000000 88.250000 91.000000 Mediana 80.000000 87.000000 89.500000 91.500000 94.000000 Média 80.270000 86.290000 88.500000 90.500000 91.840000 Terceiro Quartil 85.000000 89.500000 92.000000 95.000000 96.000000 Máximo 94.000000 98.000000 98.000000 98.000000 99.000000 NAs 4.000000 3.000000 3.000000 3.000000 0.000000 desviopadrão 7.362374 6.254838 6.542402 7.439656 7.879345 variância 54.204545 25.746212 21.530303 25.130682 12.780303 66.000000 65.000000 62.000000 59.000000 60.000000 94.000000 98.000000 98.000000 98.000000 99.000000 amplitudediff 28.000000 33.000000 36.000000 39.000000 39.000000 > porano 2008 2009 2010 2011 2012 Minimo 66.00 65.00 62.0 59.00 60.00 Primeiro quartil 76.00 84.00 86.0 88.25 91.00 Mediana 80.00 87.00 89.5 91.50 94.00 Média 80.27 86.29 88.5 90.50 91.84 Terceiro Quartil 85.00 89.50 92.0 95.00 96.00 Máximo 94.00 98.00 98.0 98.00 99.00 NAs 4.00 3.00 3.0 3.00 0.00 > eu preciso de ter o terceiro quartil menos o primeiro quartil e de ter mínimo menos 1,5 vezes terceiro quartil menos o primeiro quartil e o máximo mais 1,5 vezes o terceiro quartil menos o primeiro quartil ou seja obter a amplitude interquartil e os outliers. alguem tem algum exemplo? e que no grafico me apareça mínimo 1ºquartil mediana media 3ºquartil máximo e as dois tipos de amplitude. boxplot(dados[,-1], axes=FALSE) # O parametro axes=FALSE faz com que os eixos não sejam desenhados title(main= 'indicador eKtv sup 1_2 anual') # Acrescenta o título do gráfico title(xlab= 'Ano') # Acrescenta o titulo do eixo X title(ylab= 'Percentual') # Acrescenta o titulo do eixo Y box() # Desenha a caixa em torno do gráfico axis(2) # Coloca o eixo Y (o 2 ali representa o eixo Y) no gráfico axis(1, 1:5, 2008:2012) # Coloca no eixo X (representado pelo 1), nos pontos 1 a 5, as etiquetas 2008 a 2012. obrigada Ana Rocha

2 1

código do R
by alanarocha＠sapo.pt 13 Mar '13

13 Mar '13

Eu encontrei a forma de calcular os quantis que precisava mas não conheço nem percebo os erros que me deram e os resultados de "IRQ" não me aparecem. código: # Estatística por ano porano<-apply(dados[,-1],2,function(x)c(summary(x),0)[1:7]) desvio padrão<-sd(dados[,-1], na.rm= TRUE) Error: unexpected symbol in "desvio padrão" variância<-diag(var(dados[,-1], na.rm= TRUE)) IQR) <- quantile(dados[,-1], 3/4) - Error: unexpected ')' in "IQR)" quantile(dados[,-1], 1/4) Error in quantile.default(dados[, -1], 1/4) : missing values and NaN's not allowed if 'na.rm' is FALSE amplitude<-apply(dados[,-1],2,range, na.rm= TRUE) amplitudediff<-c(amplitude[2,] - amplitude[1,]) rownames(porano) <- c('Minimo','Primeiro quartil','Mediana','Média','Terceiro Quartil','Máximo','NAs') poranol<-rbind(porano,desvio padrão,variância,amplitude,amplitudediff,IQR) Error: unexpected symbol in "poranol<-rbind(porano,desvio padrão" > colnames(porano) = 2008:2012 > poranol X2008 X2009 X2010 X2011 X2012 Minimo 66.000000 65.000000 62.000000 59.000000 60.000000 Primeiro quartil 76.000000 84.000000 86.000000 88.250000 91.000000 Mediana 80.000000 87.000000 89.500000 91.500000 94.000000 Média 80.270000 86.290000 88.500000 90.500000 91.840000 Terceiro Quartil 85.000000 89.500000 92.000000 95.000000 96.000000 Máximo 94.000000 98.000000 98.000000 98.000000 99.000000 NAs 4.000000 3.000000 3.000000 3.000000 0.000000 desviopadrão 7.362374 6.254838 6.542402 7.439656 7.879345 variância 54.204545 25.746212 21.530303 25.130682 12.780303 66.000000 65.000000 62.000000 59.000000 60.000000 94.000000 98.000000 98.000000 98.000000 99.000000 amplitudediff 28.000000 33.000000 36.000000 39.000000 39.000000 > porano 2008 2009 2010 2011 2012 Minimo 66.00 65.00 62.0 59.00 60.00 Primeiro quartil 76.00 84.00 86.0 88.25 91.00 Mediana 80.00 87.00 89.5 91.50 94.00 Média 80.27 86.29 88.5 90.50 91.84 Terceiro Quartil 85.00 89.50 92.0 95.00 96.00 Máximo 94.00 98.00 98.0 98.00 99.00 NAs 4.00 3.00 3.0 3.00 0.00

2 1

cor no titulo
by Ze Henrique 13 Mar '13

13 Mar '13

Pessoal gostaria de saber se alguem sabe como mudar a cor do título do gráfico e dos eixos x e y. Obrigado. -- *José Henrique Soler Guilhen* Graduando em Agronomia Universidade Federal do Espírito Santo - CCAUFES

3 3

gráficos no R
by alanarocha＠sapo.pt 13 Mar '13

13 Mar '13

Olá, eu queria fazer um gráfico boxplot do sumário das estatísticas descritivas, tenho um exemplo mas não aparece a média no gráfico, o que terei de acrescentar ao seguinte código: boxplot(dados[,-1], axes=FALSE) # O parametro axes=FALSE faz com que os eixos não sejam desenhados title(main= 'indicador mortalidade anual') # Acrescenta o título do gráfico title(xlab= 'Ano') # Acrescenta o titulo do eixo X title(ylab= 'Percentual') # Acrescenta o titulo do eixo Y box() # Desenha a caixa em torno do gráfico axis(2) # Coloca o eixo Y (o 2 ali representa o eixo Y) no gráfico axis(1, 1:5, 2008:2012) # Coloca no eixo X (representado pelo 1), nos pontos 1 a 5, as etiquetas 2008 a 2012. para que a média apareça no gráfico juntamente com as outras estatísticas. beijinhos Ana

1 0

Formatação de números
by Edmar Caldas 13 Mar '13

13 Mar '13

Boa Tarde, Pessoal preciso de uma ajuda, tenho um arquivo em txt dessa forma ,100 preciso de algum comando que faz com que ele fique assim 0,100 com o ZERO na frente. Edmar

3 2

Criar objeto com operações distintas
by ASANTOS 12 Mar '13

12 Mar '13

Bom dia pessoal, Estou tentando obter um novo objeto com alguns dados que sera ó resultado de várias operações, sendo: # dados<-expand.grid(especie=c("sp1","sp2","sp3","sp4","sp5", "sp6","sp7","sp8"),##8 espécies de inseto dossel=rnorm(5,75,15), ## Cobertura de dossel elevacao=rnorm(5,750,50),##Elevação do terreno argila=rnorm(5,45,15),## Teor de argila presenca=rbinom(5,1,0.6)) ## Presença ou ausência do inseto # Gostaria de explicitar que a cada variação em dados$especie, fosse realizado a soma de dados$presenca e a média de dados$dossel, dados$elevacao e dados$argila, na criação do novo objeto, mas com o comando tapply(), não estou tendo sucesso em uma única linha de comando, alguém teria alguma sugestão para fazer esta operação de maneira direta, 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 ======================================================================

3 8

Re: [R-br] Atualização R 2.15.3i
by Sérgio Henrique almeida da silva ju 12 Mar '13

12 Mar '13

Oi Fernando Obrigado pela ajuda. Não entendi muito bem o que você quis dizer... Abraços Em 11 de março de 2013 21:05, FHRB Toledo <fernandohtoledo(a)gmail.com>escreveu: > Use outro repositório, nesse caso o repositório da versão mais atual e não > o padrão da sua distro! > > 2013/3/11 Sérgio Henrique almeida da silva ju <sergio.edfisica(a)gmail.com> > >> Prezados >> >> Estou tentando atualizar meu R para a versão 2.15.3, porém não estou >> conseguindo. Já utilizei os tutorias anteriores aqui da lista mais não deu >> certo. >> >> Utilizo o Ubuntu 12.10 (quantal) 64 bit e meu R está na versão 2.15.1. >> >> Alguém tem alguma ideia de como posso atualizar? >> >> Att >> >> -- >> 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 / 94429486 >> >> _______________________________________________ >> 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. >> > > > _______________________________________________ > 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. > -- 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 / 94429486

2 1

Outra dica de compartilhamento de arquivo
by Benilton Carvalho 12 Mar '13

12 Mar '13

Pessoal, aqui vai uma outra dica (além das que já discutimos antes) de como compartilhar arquivos com outros membros da lista. A) File Dropper: http://www.filedropper.com vcs não precisam assinar... basta escolher o arquivo e fazer o upload... um link será gerado e vocês compartilham apenas o link. Nao consegui gerar um link direto para o arquivo. B) Dropbox: http://www.dropbox.com precisa de conta e é possível obter um link direto para leitura a partir do R (desde q o arquivo esteja em Public) C) Ge.tt: http://www.ge.tt nao precisa de conta e vc tbm pode obter um link direto para o arquivo (clique o botao direito sobre "Download" e use o link dado)... Por exemplo: cat(paste(readLines('http://ge.tt/api/1/files/9EuKLua/0/blob?download'), collapse='\n')) Poderia usar mais tempo descobrindo outros serviços, mas acho que por aí basta... Sintam-se à vontade para adicionar outras sugestões. b

1 0

aRT, show de bola
by Gilbert Queiroz 12 Mar '13

12 Mar '13

Pessoal, estou utilizando o pacote aRT para trabalhar com dados de um db criado no TerraView. Consegui chegar até os dados das tabelas do layer. Porém, surgiu uma dúvida: é possivel modelar esses dados no R e criar uma nova coluna, ainda no R, com os resultados dessa modelagem, para ser introduzida no db original, via R ? Abs. Gilbert.

1 0