Novembro 2011 - R-br

Tirar o cabeçalho do txt
by Simone D. Sartorio 03 Nov '11

03 Nov '11

Pessoal, boa tarde! Gente, estou com o seguinte problema, apareceu um rapaz aqui com um banco de dados com o formato que eu deixei disponivel no site: http://www.datafilehost.com/download-b2091792.html O programa que ele usa, pelo que ele me disse, solta esse cabeçalho, que não tem como tirar pelo programa: ------------------------------------------------------------------------------------------------ *GROWTH ASPECTS OUTPUT FILE *DSSAT Cropping System Model Ver. 4.5.0.047 APR 11, 2011; 16:16:38 *RUN 1 : nov_01_DAS MZCER045 GOJD1105 1 MODEL : MZCER045 - Maize EXPERIMENT : GOJD1105 SN MILHO_SAFRA_JANDAIA DATA PATH : C:\DSSAT45\Seasonal\ TREATMENT 1 : nov_01_DAS MZCER045 ------------------------------------------------------------------------------------------------ Mas isto ele solta para cada amostra obtida, ou seja, nesse arquivo q deixei no site, que eu salvei em '.txt', tem vários deste tipo, identificando a amostra, mas o cabeçalho é meio q padrão. Alguem sabe se é possível apagar este cabeçalho deste arquivo de uma maneira inteligente, sem ser manualmente?! Pq o cara tem 30 anos de informações, dai pensa no trabalho.... Deletar linha, coluna, ou celula é facil de uma matriz, mas coisa assim eu desconheço se dá p fazer no R. Alguém tem alguma ideia?! Pode ser em outro programa tb! desde já muito obrigada! Simone Simone Daniela Sartorio Licenciada em Matemática - UNESP (Rio Claro); Mestre e atualmente Doutoranda em Estatística e Experimentação Agronômica - ESALQ/USP. Cel: (19)9105-9032 - Claro (19)8182-0586 - Tim ******************************************************************* "Gostaria de te desejar tantas coisas. Mas nada seria suficiente. Então, desejo apenas que você tenha muitos desejos. Desejos grandes. E que eles possam te mover a cada minuto, ao rumo da sua felicidade!" (Carlos Drummond de Andrade) ******************************************************************* Tenha um Bom Dia!!!

4 4

Label teimoso
by Fernando Colugnati 03 Nov '11

03 Nov '11

Pessoal, estou utilizando um pacote que gera alguns gráficos de diagnóstico, conforme o código abaixo. Reparem que no box-plot eu gostaria apenas do label "GBM p-scores", mas ele insiste em colocar o método ks.stat.max junto. Como faço para "apagar" isso? Isso ocorre com vários gráficos, e estou comparando alguns modelos diferentes para estimação dos p-scores (propensity scores). Obrigado, segue o código (retirado do exemplo do pacote, os dados são fornecidos como exemplo): #Pacote twang para estimação de propensity scores library(twang) data(lalonde) ps.lalonde <- ps(treat ~ age + educ + black + hispan + nodegree + married + re74 + re75, data = lalonde, title="Lalonde example", stop.method=stop.methods[c("ks.stat.max")], # generate plots? plots="none", pdf.plots=FALSE, # gbm options n.trees=2000, interaction.depth=3, shrinkage=0.005, perm.test.iters=0, verbose=FALSE) plot(ps.lalonde, label= "GBM p-scores", plots = "ps boxplot") -- Fernando A.B. Colugnati

2 2

Re: [R-br] Programar no R com Emacs ou Vi
by Ivan Bezerra Allaman 03 Nov '11

03 Nov '11

De acordo com o gráfico em anexo o melhor é vim! kkkkkkkkk Só para descontrair! (S,f,P) Allaman \begin{signature} <<>>= Prof. Dr. Ivan Bezerra Allaman Universidade Estadual de Santa Cruz Departamento de Ciências Exatas e Tecnológicas Ilhéus/BA - Brasil Fone: +55 73 3680-5076 E-mail: ivanalaman@yahoo.com.br/ivanalaman@gmail.com @ \end{signature}

1 0

plotando pontos em um mapa
by Rodrigo Plei 03 Nov '11

03 Nov '11

Prezados, Como faço para plotar pontos (coordenadas geográficas) a partir de um arquivo .csv em um mapa .shp usando os pacotes maps and maptools? Gostaria de fazer no R o que se faz utilizando a função "overlay maps" do Surfer! 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

2 2

substituindo valores (data.frame)
by Rogerio R. Silva 03 Nov '11

03 Nov '11

Caros, Tenho dois arquivos. Um deles, com milhares de linhas devido a replicação temporal, inclui três colunas (campanha, transecto, parcela). O segundo, contém informação sobre variáveis medidas em cada parcela (um único valor). Eu preciso substituir os dados das parcelas (primeiro arquivo) por valores das variáveis (segundo arquivo) condicionando a que transecto a parcela pertença. Um exemplo abaixo para uma variável, mas que não trabalha corretamente. coletas<-data.frame (campanha=c(gl(4,8,64,labels=c("C1","C2","C3","C4"))),transecto=c(gl(8,4,64,labels=c("T1","T2","T3","T4","T5","T6","T7","T8"))),parcela=c(gl(4,1,64,labels=c("P1","P2","P3","P4")))) dados<- data.frame (transecto=c(gl(8,4,32,labels=c("T1","T2","T3","T4","T5","T6","T7","T8"))), parcela=c(gl(4,1,32,labels=c("P1","P2","P3","P4"))), alt=rnorm(32)) x.temp<-coletas for (i in 1:length(levels(coletas$transecto))){ x.temp[which(coletas$parcela[i,]==dados$parcela[i,])] <- dados$Alt[i,] } x.temp Error in coletas$parcela[i, ] : incorrect number of dimensions Obrigado.- rogério

2 2

Carregar muitos arquivos txt
by Eder David Borges da Silva 03 Nov '11

03 Nov '11

Bom dia pessoal, venho mais uma vez pedir ajuda de vocês, O problema é o seguinte, tenho 7 pastas(Cenários), dentro dos mesmo tenho 949 arquivos TXT (1.5 Mb) cada, preciso carregar estes 7*949 arquivos e guardar em um array, tenho que carregar tudo pois tudo isso será plotado. meu código é o seguinte: cenarios <- c(1:7) con <- 1 for (i in cenarios){ setwd(paste('C:\\Users\\l01199\\Desktop\\Modelo\\Resultados\\Modelo_Cenario',i,'_Q1300',sep='')) files <- sort(list.files(pattern = "m000.TXT"))[1:50] tempo <- as.numeric(gsub("m000.TXT","",gsub("CON", "", files))) dim.a <- c(dim(read.table(files[1],head=T)),length(tempo)) if(con==1){ dados <- array(NA,c(length(cenarios),dim.a)) nomes.col <- names(read.table(files[1],head=T)) } con1 <- 1 for (j in files){ dados[con,,,con1] <- as.matrix(read.table(j,head=TRUE)) con1 <- con1+1 } print(con) con <- con+1 } tenho 4 gb de memoria na maquina, porem quando crio o objeto dados que é um array de tamanho c(length(cenarios),dim.a) [1] 7 11069 11 949 Não tenho memoria para tudo isso, alguém tem uma sugestão para guardar isso, dei uma olhada no pacote SOAR, mas como ele guarda e depois tenho que acessar novamente não sei se é a melhor saída(Posso estar equivocado em minha interpretação). Alguém possui um dica para contornar este problema. Att

2 5

Script para atualização de versão
by Marcio Nicolau 02 Nov '11

02 Nov '11

Olá Pessoas! Eu sei que já passaram várias formas para se atualizar o R para uma nova versão, mas gostaria de deixar aqui os meus $0,002 O processo é o seguinte: 1.) Crie uma pasta chamada UpdateR no home do usuário 2.) Copie os arquivos oldVersion.R e newVersion.R 3.) Na versão antiga (antes da atualização) execute o arquivo oldVersion.R no Passo (3) o script copia os arquivos (Rprofile.site e Rconsole) com as "configurações do usuário" - caso tenha outros arquivos, altere o script - para a pasta atual "UpdateR" e salva uma lista com as bibliotecas instaladas. 4.) Desinstale o R 5.) Instale a nova versão 6.) Execute o arquivo newVersion.R no Passo (6) o script copia os arquivos (Rprofile.site e Rconsole) e instala as bibliotecas da lista. Espero que seja util! Abraços e bom feriado! MN (@marcionicolau)

4 6

teste de médias para interação
by Daniele Portela de oliveira 02 Nov '11

02 Nov '11

Caros colegas gostaria de fazer um teste de médias para interação pelo pacote multcomp. Onde eu poderia encontrar um exemplo para distribuição de poisson ou binomial? Daniele Portela de Oliveira Mestranda Produção Animal - UEM Cel: 44 9949 3525 67 9617 0965

4 6

Duvida com métodos de plot pra função stripchart
by Augusto Ribas 02 Nov '11

02 Nov '11

Ola pessoal. Estava usando o comando stripchart para fazer alguns gráficos e gostaria de saber se alguém tem alguma forma simples de fazer um grafico similar a esse da figura do systat. http://www.intesoft.com/produits/tech/systat/images/density/pic7.jpg aqui tem outra imagem mas com um boxplot do lado, mas eu quero so as bolinhas alinhadas centralizadas mesmo. http://www.solutions4u-asia.com/PDT/SYSTAT/sigmaplot/SP1110/DotdensityGraph… Acredito que seria alguma coisa tipo o method="stack" mas que ele não alinha-se os pontos a esquerda mas centralizado. Não encontrei mais métodos fora o stack e o jitter. Espero ter conseguido me explicar. Bem é isso, segue um cmr de exemplo: #dados Densidade<-round(c(rnorm(30,2),rnorm(30,5))) Local<-rep(c("A","B"),each=30) #grafico stripchart(Densidade~Local,vertical=T,method="stack",offset=0.6,at=c(1.2,2.2),frame.plot=F,pch=16,cex=1.4) -- Grato Augusto C. A. Ribas Site Pessoal: http://augustoribas.heliohost.org Lattes: http://lattes.cnpq.br/7355685961127056

4 5

Fwd: Windows binary of 2.14.0
by Leonard de Assis 01 Nov '11

01 Nov '11

-------- Mensagem original -------- Assunto: Windows binary of 2.14.0 Data: Tue, 01 Nov 2011 16:03:48 -0400 De: Duncan Murdoch <murdoch.duncan(a)gmail.com> Para: r-announce(a)lists.r-project.org <r-announce(a)r-project.org> To all Windows users: The binary build of 2.14.0 that was uploaded yesterday was missing Cairo support. I have rebuilt it, and uploaded a new copy. You can tell which one you have by running "svg()", which works on the new one, but not the old one. You can tell which one is on your CRAN mirror by looking at the last changed date: it is today (November 1) for the corrected build. Sorry for the inconvenience.... Duncan Murdoch _______________________________________________ R-announce(a)r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-announce

1 0